Evolutionary history of the extant amphioxus lineage with shallow-branching diversification

Asymmetric Segregation of Maternal mRNAs and Germline-related Determinants in Cephalochordate Embryos: Implications for the Evolution of Early Patterning Events in Chordates

How animal embryos determine their early cell fates is an important question in developmental biology. In various model animals, asymmetrically localized maternal transcripts play important roles in axial patterning and cell fate specification. Cephalochordates (amphioxus), which have three living genera (Asymmetron, Epigonichthys, and Branchiostoma), are an early branching chordate lineage and thus occupy a key phylogenetic position for understanding the evolution of chordate developmental mechanisms. It has been shown that in the zygote of Branchiostoma amphioxus, which possesses bilateral gonads flanking both sides of their trunk region, maternal transcripts of germline determinants form a compact granule. During early embryogenesis, this granule is inherited by a single blastomere, which subsequently gives rise to a cluster of cells displaying typical characteristics of primordial germ cells (PGC). These PGCs then come to lie in the tailbud region and proliferate during posterior elongation of the larvae to join in the gonad anlagen at the ventral tip of the developing myomeres in amphioxus larvae. However, in Asymmetron and Epigonichthys amphioxus, whose gonads are present only on the right side of their bodies, nothing is known about their PGC development or the cellular/morphogenetic processes resulting in the asymmetric distribution of gonads. Using conserved germline determinants as markers, we show that similarly to Branchiostoma amphioxus, Asymmetron also employs a preformation mechanism to specify their PGCs, suggesting that this mechanism represents an ancient trait dating back to the common ancestor of Cephalochordates. Surprisingly, we found that Asymmetron PGCs are initially deposited on both sides of the body during early larval development; however, the left-side PGCs cease to exist in young juveniles, suggesting that PGCs are eliminated from the left body side during larval development or following metamorphosis. This is reminiscent of the PGC development in the sea urchin embryo, and we discuss the implications of this observation for the evolution of developmental mechanisms.

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.

Nuclear import of IRF11 via the importin α/β pathway is essential for its antiviral activity

Interferon regulatory factor 11 (IRF11), an intriguing IRF member found only in fish species, has recently been shown to have antiviral properties that are dependent on its nuclear entry and DNA binding affinity. However, the mechanisms by which IRF11 enters the nucleus are unknown. In the present study, we found orthologs of IRF11 in lamprey and lancelet species by combining positional, phylogenetic and structural comparison data, showing that this gene has an ancient origin. The IRF11 gene (AjIRF11) from the Japanese eel, Anguilla japonica, was subsequently characterized, and it was found that AjIRF11 has antiviral activities against spring viremia of carp virus (SVCV), which are accomplished by regulating the production of type I IFN and IFN-stimulated genes. In addition to its known DNA binding residues in the α3 helix, two residues in Loop 1, His40 and Trp46, are also involved in DNA binding and activation of the IFN promoter. Using immunofluorescence microscopy and site-directed mutagenesis analysis, we confirmed that full nuclear localization of AjIRF11 requires the bipartite nuclear localization sequence (NLS) spanning residues 75 to 101, as well as the monopartite NLS situated between residues 119 and 122. Coimmunoprecipitation assays confirmed that AjIRF11 interacts with importin α via its NLSs and can also bind to importin β directly, implying that IRF11 can be imported to the nucleus by one or more transport receptors.

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.

An Updated Staging System for Cephalochordate Development: One Table Suits Them All

Chordates are divided into three subphyla: Vertebrata, Tunicata, and Cephalochordata. Phylogenetically, the Cephalochordata, more commonly known as lancelets or amphioxus, constitute the sister group of Vertebrata and Tunicata. Lancelets are small, benthic, marine filter feeders, and their roughly three dozen described species are divided into three genera: Branchiostoma, Epigonichthys, and Asymmetron. Due to their phylogenetic position and their stereotypical chordate morphology and genome architecture, lancelets are key models for understanding the evolutionary history of chordates. Lancelets have thus been studied by generations of scientists, with the first descriptions of adult anatomy and developmental morphology dating back to the 19th century. Today, several different lancelet species are used as laboratory models, predominantly for developmental, molecular and genomic studies. Surprisingly, however, a universal staging system and an unambiguous nomenclature for developing lancelets have not yet been adopted by the scientific community. In this work, we characterized the development of the European lancelet (Branchiostoma lanceolatum) using confocal microscopy and compiled a streamlined developmental staging system, from fertilization through larval life, including an unambiguous stage nomenclature. By tracing growth curves of the European lancelet reared at different temperatures, we were able to show that our staging system permitted an easy conversion of any developmental time into a specific stage name. Furthermore, comparisons of embryos and larvae from the European lancelet (B. lanceolatum), the Florida lancelet (Branchiostoma floridae), two Asian lancelets (Branchiostoma belcheri and Branchiostoma japonicum), and the Bahamas lancelet (Asymmetron lucayanum) demonstrated that our staging system could readily be applied to other lancelet species. Although the detailed staging description was carried out on developing B. lanceolatum, the comparisons with other lancelet species thus strongly suggested that both staging and nomenclature are applicable to all extant lancelets. We conclude that this description of embryonic and larval development will be of great use for the scientific community and that it should be adopted as the new standard for defining and naming developing lancelets. More generally, we anticipate that this work will facilitate future studies comparing representatives from different chordate lineages.

Molecular taxonomy confirms that the northeastern Atlantic and Mediterranean Sea harbor a single lancelet, Branchiostoma lanceolatum (Pallas, 1774) (Cephalochordata: Leptocardii: Branchiostomatidae)

Branchiostomatidae (lancelets or amphioxus) comprises about 30 species, several of which are well-established models in evolutionary development. Our zoological and ecological knowledge of the family is nonetheless limited. Despite evident differences can be found among known populations, the taxonomy of Branchiostoma lanceolatum (type species of the genus Branchiostoma) has never been investigated with modern methods through its range in the northeastern Atlantic and Mediterranean Sea. We address this via a multilocus molecular approach and comparing specimens collected from different European populations. Results obtained here confirm the presence of a single species inhabiting the range between the topotypical localities of B. lanceolatum (Atlantic Ocean) and of its junior synonym B. lubricum (Mediterranean Sea), without evincing geographical structure between populations. This suggests that environment most likely drives the characteristics observed in different geographic areas. The long larval phase and the slow mutation rate in lancelets may have played a key role in the evolutionary history of this iconic species.

Dynamic Evolution of the Cthrc1 Genes, a Newly Defined Collagen-Like Family

Collagen triple helix repeat containing protein 1 (Cthrc1) is a secreted glycoprotein reported to regulate collagen deposition and to be linked to the Transforming growth factor β/Bone morphogenetic protein and the Wnt/planar cell polarity pathways. It was first identified as being induced upon injury to rat arteries and was found to be highly expressed in multiple human cancer types. Here, we explore the phylogenetic and evolutionary trends of this metazoan gene family, previously studied only in vertebrates. We identify Cthrc1 orthologs in two distant cnidarian species, the sea anemone Nematostella vectensis and the hydrozoan Clytia hemisphaerica, both of which harbor multiple copies of this gene. We find that Cthrc1 clade-specific diversification occurred multiple times in cnidarians as well as in most metazoan clades where we detected this gene. Many other groups, such as arthropods and nematodes, have entirely lost this gene family. Most vertebrates display a single highly conserved gene, and we show that the sequence evolutionary rate of Cthrc1 drastically decreased within the gnathostome lineage. Interestingly, this reduction coincided with the origin of its conserved upstream neighboring gene, Frizzled 6 (FZD6), which in mice has been shown to functionally interact with Cthrc1. Structural modeling methods further reveal that the yet uncharacterized C-terminal domain of Cthrc1 is similar in structure to the globular C1q superfamily domain, also found in the C-termini of collagens VIII and X. Thus, our studies show that the Cthrc1 genes are a collagen-like family with a variable short collagen triple helix domain and a highly conserved C-terminal domain structure resembling the C1q family.

Filter feeding, deviations from bilateral symmetry, developmental noise, and heterochrony of hemichordate and cephalochordate gills

We measured gill slit fluctuating asymmetry (FA), a measure of developmental noise, in adults of three invertebrate deuterostomes with different feeding modes: the cephalochordate Branchiostoma floridae (an obligate filter feeder), the enteropneusts Protoglossus graveolens (a facultative filter feeder/deposit feeder) and Saccoglossus bromophenolosus (a deposit feeder). FA was substantially and significantly low in B. floridae and P. graveolens and high in S. bromophenolosus. Our results suggest that the gills of species that have experienced a relaxation of the filter feeding trait exhibit elevated FA. We found that the timing of development of the secondary collagenous gill bars, compared to the primary gill bars, was highly variable in P. graveolens but not the other two species, demonstrating an independence of gill FA from gill bar heterochrony. We also discovered the occasional ectopic expression of a second set of paired gills posterior to the first set of gills in the enteropneusts and that these were more common in S. bromophenolosus. Moreover, our finding that gill slits in enteropneusts exhibit bilateral symmetry suggests that the left-sidedness of larval cephalochordate gills, and the directional asymmetry of Cambrian stylophoran echinoderm fossil gills, evolved independently from a bilaterally symmetrical ancestor.

Ferritins in Chordata: Potential evolutionary trajectory marked by discrete selective pressures: History and reclassification of ferritins in chordates and geological events' influence on their evolution and radiation

Ferritins (FTs) are iron storage proteins that are involved in managing iron-oxygen balance. In our work, we present a hypothesis on the putative effect of geological changes that have affected the evolution and radiation of ferritin proteins. Based on sequence analysis and phylogeny reconstruction, we hypothesize that two significant factors have been involved in the evolution of ferritin proteins: fluctuations of atmospheric oxygen concentrations, altering redox potential, and changing availability of water rich in bioavailable ferric ions. Fish, ancient amphibians, reptiles, and placental mammals developed the broadest repertoire of singular FTs, attributable to embryonic growth in aquatic environments containing low oxygen levels and abundant forms of soluble iron. In contrast, oviparous land vertebrates, like reptiles and birds, that have developed in high oxygen levels and limited levels of environmental Fe²⁺ exhibit a lower diversity of singular FTs, but display a broad repertoire of subfamilies, particularly notable in early reptiles.

Cephalochordates: A window into vertebrate origins

How vertebrates evolved from their invertebrate ancestors has long been a central topic of discussion in biology. Evolutionary developmental biology (evodevo) has provided a new tool-using gene expression patterns as phenotypic characters to infer homologies between body parts in distantly related organisms-to address this question. Combined with micro-anatomy and genomics, evodevo has provided convincing evidence that vertebrates evolved from an ancestral invertebrate chordate, in many respects resembling a modern amphioxus. The present review focuses on the role of evodevo in addressing two major questions of chordate evolution: (1) how the vertebrate brain evolved from the much simpler central nervous system (CNS) in of this ancestral chordate and (2) whether or not the head mesoderm of this ancestor was segmented.

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.

The evolution of zebrafish RAG2 protein is required for adapting to the elevated body temperature of the higher endothermic vertebrates

The recombination activating gene (RAG or RAG1/RAG2 complex)-mediated adaptive immune system is a hallmark of jawed vertebrates. It has been reported that RAG originated in invertebrates. However, whether RAG further evolved once it arose in jawed vertebrates remains largely unknown. Here, we found that zebrafish RAG (zRAG) had a lower activity than mouse RAG (mRAG). Intriguingly, the attenuated stability of zebrafish RAG2 (zRAG2), but not zebrafish RAG1, caused the reduced V(D)J recombination efficiency compared to mRAG at 37 °C which are the body temperature of most endotherms except birds. Importantly, the lower temperature 28 °C, which is the best temperature for zebrafish growth, made the recombination efficiency of zRAG similar to that of mRAG by improving the stability of zRAG2. Consistent with the prementioned observation, the V(D)J recombination of Rag2^KI/KI mice, which zRAG2 was substituted for mRAG2, was also severely impaired. Unexpectedly, Rag2^KI/KI mice developed cachexia syndromes accompanied by premature death. Taken together, our findings illustrate that the evolution of zebrafish RAG2 protein is required for adapting to the elevated body temperature of the higher endothermic vertebrates.

Asymmetron lucayanum: How many species are valid?

The cephalochordates amphioxus or lancelets are benthic marine animals representing the earliest divergent evolutionary lineage within chordates. Although amphioxus are present in most of the world's tropical and temperate oceans, only about thirty different species grouped into three different genera, Branchiostoma, Epigonichthys and Asymmetron have been described. In the genus Asymmetron, only two species have been characterized, although for one of them, A. lucayanum, several cryptic lineages exist. In this work we have sequenced and analyzed the mitogenome of an A. lucayanum population previously described in the Red Sea. The phylogenetic study using this complete mitogenome as well as the analysis of COI gene sequences of several individuals of this Red Sea population show that the Red Sea population is a new cryptic species. We propose to call this new species Asymmetron rubrum.

Using high-resolution annotation of insect mitochondrial DNA to decipher tandem repeats in the control region

In this study, we used a small RNA sequencing (sRNA-seq) based method to annotate the mitochondrial genome of the insect Erthesina fullo Thunberg at 1 bp resolution. The high-resolution annotations cover both entire strands of the mitochondrial genome without any gaps or overlaps. Most of the new annotations were consistent with the previous annotations which had been obtained using PacBio full-length transcripts. Two important findings were that animals transcribe both entire strands of mitochondrial genomes and the tandem repeats in the control region of the E. fullo mitochondrial genome contains the repeated Transcription Initiation Sites (TISs) of the heavy strand. In addition, we found that the copy numbers of tandem repeats showed a great diversity within an individual, suggesting that mitochondrial DNA recombination occurs in an individual. In conclusion, the sRNA-seq based method uses 5' and 3' end small RNAs to annotate nuclear non-coding and mitochondrial genes at 1 bp resolution, and can be used to identify new steady RNAs, particularly long non-coding RNAs (lncRNAs). The high-resolution annotations of mitochondrial genomes can also be used to study the molecular phylogenetics and evolution of animals or to investigate mitochondrial gene transcription, RNA processing, RNA maturation and several other related topics. The complete mitochondrial genome sequence of E. fullo with the new annotations using the sRNA-seq based method is available at the NCBI GenBank database under the accession number MK374364. We publish our theories, methods, the high quality sRNA-seq and RNA-seq data (SRA: SRP174926) for extensive use.

Wnt evolution and function shuffling in liberal and conservative chordate genomes

BACKGROUND

What impact gene loss has on the evolution of developmental processes, and how function shuffling has affected retained genes driving essential biological processes, remain open questions in the fields of genome evolution and EvoDevo. To investigate these problems, we have analyzed the evolution of the Wnt ligand repertoire in the chordate phylum as a case study.

RESULTS

We conduct an exhaustive survey of Wnt genes in genomic databases, identifying 156 Wnt genes in 13 non-vertebrate chordates. This represents the most complete Wnt gene catalog of the chordate subphyla and has allowed us to resolve previous ambiguities about the orthology of many Wnt genes, including the identification of WntA for the first time in chordates. Moreover, we create the first complete expression atlas for the Wnt family during amphioxus development, providing a useful resource to investigate the evolution of Wnt expression throughout the radiation of chordates.

CONCLUSIONS

Our data underscore extraordinary genomic stasis in cephalochordates, which contrasts with the liberal and dynamic evolutionary patterns of gene loss and duplication in urochordate genomes. Our analysis has allowed us to infer ancestral Wnt functions shared among all chordates, several cases of function shuffling among Wnt paralogs, as well as unique expression domains for Wnt genes that likely reflect functional innovations in each chordate lineage. Finally, we propose a potential relationship between the evolution of WntA and the evolution of the mouth in chordates.

Morphological Stasis and Proteome Innovation in Cephalochordates

Lancelets, extant representatives of basal chordates, are prototypic examples of evolutionary stasis; they preserved a morphology and body-plan most similar to the fossil chordates from the early Cambrian. Such a low level of morphological evolution is in harmony with a low rate of amino acid substitution; cephalochordate proteins were shown to evolve slower than those of the slowest evolving vertebrate, the elephant shark. Surprisingly, a study comparing the predicted proteomes of Chinese amphioxus, Branchiostoma belcheri and the Florida amphioxus, Branchiostoma floridae has led to the conclusion that the rate of creation of novel domain combinations is orders of magnitude greater in lancelets than in any other Metazoa, a finding that contradicts the notion that high rates of protein innovation are usually associated with major evolutionary innovations. Our earlier studies on a representative sample of proteins have provided evidence suggesting that the differences in the domain architectures of predicted proteins of these two lancelet species reflect annotation errors, rather than true innovations. In the present work, we have extended these studies to include a larger sample of genes and two additional lancelet species, Asymmetron lucayanum and Branchiostoma lanceolatum. These analyses have confirmed that the domain architecture differences of orthologous proteins of the four lancelet species are because of errors of gene prediction, the error rate in the given species being inversely related to the quality of the transcriptome dataset that was used to aid gene prediction.

Pax3/7 duplicated and diverged independently in amphioxus, the basal chordate lineage

The Pax3/7 transcription factor family is integral to developmental gene networks contributing to important innovations in vertebrate evolution, including the neural crest. The basal chordate lineage of amphioxus is ideally placed to understand the dynamics of the gene regulatory network evolution that produced these novelties. We report here the discovery that the cephalochordate lineage possesses two Pax3/7 genes, Pax3/7a and Pax3/7b. The tandem duplication is ancestral to all extant amphioxus, occurring in both Asymmetron and Branchiostoma, but originated after the split from the lineage leading to vertebrates. The two paralogues are differentially expressed during embryonic development, particularly in neural and somitic tissues, suggesting distinct regulation. Our results have implications for the study of amphioxus regeneration, neural plate and crest evolution, and differential tandem paralogue evolution.

The role of transposable elements in functional evolution of amphioxus genome: the case of opsin gene family

Transposable elements (TEs) are able to jump to new locations (transposition) in the genome, usually after replication. They constitute the so-called selfish or junk DNA and take over large proportions of some genomes. Due to their ability to move around they can change the DNA landscape of genomes and are therefore a rich source of innovation in genes and gene regulation. Surge of sequence data in the past years has significantly facilitated large scale comparative studies. Cephalochordates have been regarded as a useful proxy to ancestral chordate condition partially due to the comparatively slow evolutionary rate at morphological and genomic level. In this study, we used opsin gene family from three Branchiostoma species as a window into cephalochordate genome evolution. We compared opsin complements in terms of family size, gene structure and sequence allowing us to identify gene duplication and gene loss events. Furthermore, analysis of the opsin containing genomic loci showed that they are populated by TEs. In summary, we provide evidence of the way transposable elements may have contributed to the evolution of opsin gene family and to the shaping of cephalochordate genomes in general.