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.

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.

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.

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.

Impacts of microplastics and the associated plastisphere on physiological, biochemical, genetic expression and gut microbiota of the filter-feeder amphioxus

Oceanic plastic pollution is of major concern to marine organisms, especially filter feeders. However, limited is known about the toxic effects of the weathered microplastics instead of the pristine ones. This study evaluates the effects of weathered polystyrene microplastic on a filter-feeder amphioxus under starvation conditions via its exposure to the microplastics previously deployed in the natural seawater allowing for the development of a mature biofilm (so-called plastisphere). The study focused on the integration of physiological, histological, biochemical, molecular, and microbiota impacts on amphioxus. Overall, specific alterations in gene expression of marker genes were observed to be associated with oxidative stresses and immune systems. Negligible impacts were observed on antioxidant biochemical activities and gut microbiota of amphioxus, while we highlighted the potential transfer of 12 bacterial taxa from the plastisphere to the amphioxus gut microbiota. Moreover, the classical perturbation of body shape detected in control animals under starvation conditions (a slim and curved body) but not for amphioxus exposed to microplastic, indicates that the microorganisms colonizing plastics could serve as a nutrient source for this filter-feeder, commitment with the elevated proportions of goblet cell-like structures after the microplastic exposure. The multidisciplinary approach developed in this study underlined the trait of microplastics that acted as vectors for transporting microorganisms from the plastisphere to amphioxus.

miR-92b-5p negatively regulates IKK through targeting its ORF region in the innate immune responses of amphioxus (Branchiostoma belcheri)

Cephalochordate (Amphioxus), situated at a key phylogenetic position in the phylum Chordata, serves as a model organism for studying the origin and evolution of the vertebrate innate immune. In this study, five members of precursor miR-92 family (miR-92a-1, miR-92a-2, miR-92b, miR-92c and miR-92d) are identified in Branchiostoma belcheri, and their evolutionary conservation and potential molecular functions in innate immunity are analyzed. Among them, miR-92b-5p was validated in HEK293 cells to target the coding region but not classic 3'UTR of IKK (inhibitor of nuclear factor kappa-B kinase) mRNA, one integral component of MAPK and TLR4 immune signaling. Furthermore, the spatiotemporal expression patterns of miR-92b-5p and IKK were examined in different tissues or different time points (2 h, 4 h, 8 h, 12 h, 24 h and 48 h) post LPS stimulation at RNA and protein level in vivo. The seemingly inverse expression pattern between miR-92b-5p and IKK supports the involvement of miR-92b-5p in Branchiostoma belcheri innate immune response. In conclusion, our work not only illustrates the evolutionary pattern of Branchiostoma belcheri miR-92 family across chordates, but also reveals that miR-92b-5p could target IKK expression to regulate innate immune response.

miR-2013 negatively regulates phylogenetically conserved PIP5K involved in TLR4 mediated immune responses of amphioxus (Branchiostoma belcheri Tsingtaunese)

Phosphatidylinositol 4-phosphate 5-kinase (PIP5K) is a catalytic kinase that performs multiple functions in organisms. Recent studies have shown that PIP5Kα in mammals can directly participate in the TLR-mediated innate immune regulation by controlling the production of PIP2. However, the PIP5K homologous gene has not been identified in Cephalochordata to date. In this study, we firstly identify and characterize a new PIP5K family member (designed as AmphiPIP5K) from Cephalochordata amphioxus (Branchiostoma belcheri tsingtaunese), particularly AmphiPIP5K is orthologous with vertebrate PIP5Kα and paralogous with PIP5Kβ and PIP5Kγ. Secondly, we find that the AmphiPIP5K is involved in amphioxus innate immune response to LPS stimulation. Thirdly, our results demonstrate that miR-2013 can inhibit AmphiPIP5K expression by binding to the CDS and 3' UTR regions of AmphiPIP5K. Collectively, our work not only demonstrates the evolutionary pattern of amphioxus PIP5K, but also reveals that miR-2013 negatively regulate PIP5K expression to involve in amphioxus innate immune response.

Antifungal Activity of NP20 Derived from Amphioxus Midkine/Pleiotrophin Homolog Against Aspergillus niger and Aspergillus fumigatus

With the emergence of antifungal resistance, systematic infections with Aspergillus are becoming the major cause of the clinical morbidity. The development of novel antifungal agents with high efficacy, low drug tolerance, and few side effects is urgent. In response to that need, we have identified NP20. Here we demonstrate clearly that NP20 has antifungal activity, capable of killing the spores of Aspergillus niger and Aspergillus fumigatus as well as causing direct damage to the surface, membrane, cytoplasm, organelle, and nucleus of the fungal spores. Interestingly, NP20 is active under temperature stress and a wide range of pH. Subsequently, MTT assay, assay for binding of NP20 to fungal cell wall components, membrane depolarization assay, confocal microscopy, ROS assay, DNA replication, and protein synthesis assay are performed to clarify the mechanisms underlying NP20 against Aspergillus. The results show that NP20 can bind with and pass through the fungal cell wall, and then interfere with the lipid membrane. Moreover, NP20 can induce intracellular ROS production, DNA fragmentation, and protein synthesis inhibition of the fungal cells. These together indicate that NP20 is a novel antifungal peptide, which has considerable potential for future development as novel peptide antibiotics against Aspergillus.

Autophagy participates in innate immune defense in lamprey

Autophagy is a homeostatic process which degrades cytoplasmic constituents to maintain the balance of organs when they were challenged with nutrient stress. It also participates in cancer, neurodegenerative disorders, aging and innate immune defense. In order to reveal how autophagy participates in innate immune response when invertebrates evolved into vertebrates. Firstly, we performed a systematic analysis of Atg genes and found that they are highly conserved among lancelet, lamprey and zebrafish. Then, we observed autophagosomes upon starvation by TEM in lancelet, lamprey and zebrafish and found that the morphology of autophagosome is similar to that was observed in yeast and mammals. In addition, rapamycin can induce autophagy in lamprey leukocytes and the deficiency of human Beclin1 protein can be rescued by lancelet and lamprey Beclin1 proteins. When lamprey leukocytes were treated with polyI:C and LPS, autophagy was induced. Moreover, when lamprey leukocytes were challenged with live E. coli, phagocytosis along with autophagy was triggered to degrade pathogenic bacteria. In all, our study here indicated that autophagy is highly conserved during evolution and plays a key role in innate defense when invertebrates evolved into vertebrates.

The promiscuous estrogen receptor: Evolution of physiological estrogens and response to phytochemicals and endocrine disruptors

Many actions of estradiol (E2), the principal physiological estrogen in vertebrates, are mediated by estrogen receptor-α (ERα) and ERβ. An important physiological feature of vertebrate ERs is their promiscuous response to several physiological steroids, including estradiol (E2), Δ⁵-androstenediol, 5α-androstanediol, and 27-hydroxycholesterol. A novel structural characteristic of Δ⁵-androstenediol, 5α-androstanediol, and 27-hydroxycholesterol is the presence of a C19 methyl group, which precludes the presence of an aromatic A ring with a C3 phenolic group that is a defining property of E2. The structural diversity of these estrogens can explain the response of the ER to synthetic chemicals such as bisphenol A and DDT, which disrupt estrogen physiology in vertebrates, and the estrogenic activity of a variety of plant-derived chemicals such as genistein, coumestrol, and resveratrol. Diversity in the A ring of physiological estrogens also expands potential structures of industrial chemicals that can act as endocrine disruptors. Compared to E2, synthesis of 27-hydroxycholesterol and Δ⁵-androstenediol is simpler, leading us, based on parsimony, to propose that one or both of these steroids or a related metabolite was a physiological estrogen early in the evolution of the ER, with E2 assuming this role later as the canonical estrogen. In addition to the well-studied role of the ER in reproductive physiology, the ER also is an important transcription factor in non-reproductive tissues such as the cardiovascular system, kidney, bone, and brain. Some of these ER actions in non-reproductive tissues appeared early in vertebrate evolution, long before the emergence of mammals.

Is the microplastic selective according to the habitat? Records in amphioxus sands, Mäerl bed habitats and Cymodocea nodosa habitats

This study estimated for the first time the total loads of plastic litter (macro- meso- and micro-plastics) in sediments of different habitat types from the Northern Adriatic Sea. Samples were collected in March 2016. The sampling sites were settled in shoreline, on the C. nodosa bottoms, Amphioxus sands, and Mäerl bed habitats. Microplastics items were present in all sampling site and ranging within 137-703 items/kg d.w. from Mäerl bed habitat to the shoreline. In C. nodosa bottoms 170 items/kg d.w. were found, while in Amphioxus sands were recorded on average 194 items/kg d.w. Due to the absence of statistical associations among litter levels and abundance of B. lanceolatum in the study area, this research present the needs to develop a new method and more research to for the evaluation of how much the interrelation between sensible habitats and microplastic exist.

Structural similarities and functional differences clarify evolutionary relationships between tRNA healing enzymes and the myelin enzyme CNPase

BACKGROUND

Eukaryotic tRNA splicing is an essential process in the transformation of a primary tRNA transcript into a mature functional tRNA molecule. 5'-phosphate ligation involves two steps: a healing reaction catalyzed by polynucleotide kinase (PNK) in association with cyclic phosphodiesterase (CPDase), and a sealing reaction catalyzed by an RNA ligase. The enzymes that catalyze tRNA healing in yeast and higher eukaryotes are homologous to the members of the 2H phosphoesterase superfamily, in particular to the vertebrate myelin enzyme 2',3'-cyclic nucleotide 3'-phosphodiesterase (CNPase).

RESULTS

We employed different biophysical and biochemical methods to elucidate the overall structural and functional features of the tRNA healing enzymes yeast Trl1 PNK/CPDase and lancelet PNK/CPDase and compared them with vertebrate CNPase. The yeast and the lancelet enzymes have cyclic phosphodiesterase and polynucleotide kinase activity, while vertebrate CNPase lacks PNK activity. In addition, we also show that the healing enzymes are structurally similar to the vertebrate CNPase by applying synchrotron radiation circular dichroism spectroscopy and small-angle X-ray scattering.

CONCLUSIONS

We provide a structural analysis of the tRNA healing enzyme PNK and CPDase domains together. Our results support evolution of vertebrate CNPase from tRNA healing enzymes with a loss of function at its N-terminal PNK-like domain.

Identification and functional characterization of viperin of amphioxus Branchiostoma japonicum: Implications for ancient origin of viperin-mediated antiviral response

Viperin, an antiviral protein, has been shown to be active against a wide range of DNA and RNA viruses, but no information is available regarding functional characterization of viperin in invertebrate species. In this study, we clearly demonstrate that amphioxus (Branchiostoma japonicum) viperin, BjVip, has features in common with those of vertebrate viperin, including the presence of the SAM superfamily domain with the characteristic CNYKCGFC motif, syntenic conservation, and predicted 3D structure. Bjvip exhibits a tissue-specific expression with abundant levels in the hepatic cecum, hind-gut, gill and muscle, and following challenge with the viral mimic poly I:C, its expression is significantly up-regulated, suggesting an involvement of BjVip in immune response of amphioxus against viral infection. Importantly, we show that the cells transfected with Bjvip is able to kill LCDV or inhibiting its propagation, and co-incubation of rBjVip with WSSV markedly attenuates its infectivity. Thus, we provide the first evidences that amphioxus viperin, like that of vertebrates, is capable of promoting resistance against viral infection in vitro and in vivo, indicating that viperin-mediated antiviral response already emerged in the primitive chordate. We also prove that amphioxus viperin has evolved under positive selection.

Activity Augmentation of Amphioxus Peptidoglycan Recognition Protein BbtPGRP3 via Fusion with a Chitin Binding Domain

Peptidoglycan recognition proteins (PGRPs), which have been identified in most animals, are pattern recognition molecules that involve antimicrobial defense. Resulting from extraordinary expansion of innate immune genes, the amphioxus encodes many PGRPs of diverse functions. For instance, three isoforms of PGRP encoded by Branchiostoma belcheri tsingtauense, termed BbtPGRP1~3, are fused with a chitin binding domain (CBD) at the N-terminus. Here we report the 2.7 Å crystal structure of BbtPGRP3, revealing an overall structure of an N-terminal hevein-like CBD followed by a catalytic PGRP domain. Activity assays combined with site-directed mutagenesis indicated that the individual PGRP domain exhibits amidase activity towards both DAP-type and Lys-type peptidoglycans (PGNs), the former of which is favored. The N-terminal CBD not only has the chitin-binding activity, but also enables BbtPGRP3 to gain a five-fold increase of amidase activity towards the Lys-type PGNs, leading to a significantly broadened substrate spectrum. Together, we propose that modular evolution via domain shuffling combined with gene horizontal transfer makes BbtPGRP1~3 novel PGRPs of augmented catalytic activity and broad recognition spectrum.

The cult of amphioxus in German Darwinism; or, our gelatinous ancestors in Naples' blue and balmy bay

Biologists having rediscovered amphioxus, also known as the lancelet or Branchiostoma, it is time to reassess its place in early Darwinist debates over vertebrate origins. While the advent of the ascidian-amphioxus theory and challenges from various competitors have been, documented, this article offers a richer account of the public appeal of amphioxus as a primitive ancestor. The focus is on how the 'German Darwin' Ernst Haeckel persuaded general magazine and newspaper readers to revere this "flesh of our flesh and blood of our blood", and especially on Das neue Laienbrevier des Haeckelismus (The new lay breviary of Haeckelism) by Moritz Reymond with cartoons by Fritz Steub. From the late 1870s these successful little books of verse introduced the Neapolitan discoveries that made the animal's name and satirized Haeckel's rise as high priest of its cult. One song is reproduced and translated here, with a contemporary "imitation" by the Canadian palaeontologist Edward John Chapman, and extracts from others. Predating the American "It's a long way from amphioxus" by decades, these rhymes dramatize neglected 'species politics' of Darwinism and highlight the roles of humour in negotiating evolution.

Crystallization and preliminary X-ray diffraction analysis of a single variable domain of the immunoglobulin superfamily in amphioxus, Amphi-IgSF-V

Amphioxus is regarded as an essential animal model for the study of immune evolution. Discovery of new molecules with the immunoglobulin superfamily (IgSF) variable (V) domain in amphioxus would help in studying the evolution of IgSF V molecules in the immune system. A protein was found which just contains only one IgSF V domain in amphioxus, termed Amphi-IgSF-V; it has over 30% sequence identity to the V domains of human immunoglobulins and mammalian T-cell receptors. In order to clarify the three-dimensional structure of this new molecule in amphioxus, Amphi-IgSF-V was expressed, purified and crystallized, and diffraction data were collected to a resolution of 1.95 Å. The crystal belonged to space group P3221, with unit-cell parameters a = b = 53.9, c = 135.5 Å. The Matthews coefficient and solvent content were calculated to be 2.58 Å(3) Da(-1) and 52.38%, respectively. The results will provide structural information to study the evolution of IgSF V molecules in the immune system.

Structure of the red fluorescent protein from a lancelet (Branchiostoma lanceolatum): a novel GYG chromophore covalently bound to a nearby tyrosine

A key property of proteins of the green fluorescent protein (GFP) family is their ability to form a chromophore group by post-translational modifications of internal amino acids, e.g. Ser65-Tyr66-Gly67 in GFP from the jellyfish Aequorea victoria (Cnidaria). Numerous structural studies have demonstrated that the green GFP-like chromophore represents the `core' structure, which can be extended in red-shifted proteins owing to modifications of the protein backbone at the first chromophore-forming position. Here, the three-dimensional structures of green laGFP (λex/λem = 502/511 nm) and red laRFP (λex/λem ≃ 521/592 nm), which are fluorescent proteins (FPs) from the lancelet Branchiostoma lanceolatum (Chordata), were determined together with the structure of a red variant laRFP-ΔS83 (deletion of Ser83) with improved folding. Lancelet FPs are evolutionarily distant and share only ∼20% sequence identity with cnidarian FPs, which have been extensively characterized and widely used as genetically encoded probes. The structure of red-emitting laRFP revealed three exceptional features that have not been observed in wild-type fluorescent proteins from Cnidaria reported to date: (i) an unusual chromophore-forming sequence Gly58-Tyr59-Gly60, (ii) the presence of Gln211 at the position of the conserved catalytic Glu (Glu222 in Aequorea GFP), which proved to be crucial for chromophore formation, and (iii) the absence of modifications typical of known red chromophores and the presence of an extremely unusual covalent bond between the Tyr59 C(β) atom and the hydroxyl of the proximal Tyr62. The impact of this covalent bond on the red emission and the large Stokes shift (∼70 nm) of laRFP was verified by extensive structure-based site-directed mutagenesis.

Structure of the longitudinal body muscles of amphioxus

The structure of the longitudinal body muscles of Branchiostoma caribaeum has been studied by light and electron microscopy. These muscles are shown to be composed of fibers in the form of flat lamellae about 0.8µ in thickness, more than 100 µ wide, and reaching in length from one intermuscular septum to the next, a distance of about 0.6 mm. Each flat fiber is covered by a plasma membrane and contains a single myofibril consisting of myofilaments packed in the interdigitating hexagonal array characteristic of vertebrate striated muscle. Little or no sarcoplasmic reticulum is present. Mitochondria are found infrequently and have a tubular internal structure. These morphological observations are discussed in relation to a proposed hypothesis of excitation-contraction coupling. It is pointed out that the maximum distance from surface to myofilament in these muscles is about 0.5 µ and that diffusion of an "activating" substance over this distance would essentially be complete in less than 0.5 msec. after its release from the plasma membrane. It is concluded that the flat form of amphioxus muscle substitutes for the specialized mechanisms of excitation-contraction coupling thought possibly to involve the sarcoplasmic reticulum in higher vertebrate muscles.