How much does the amphioxus genome represent the ancestor of chordates?

Fast cycling culture of the annelid model Platynereis dumerilii

Platynereis dumerilii, a marine annelid, is a model animal that has gained popularity in various fields such as developmental biology, biological rhythms, nervous system organization and physiology, behaviour, reproductive biology, and epigenetic regulation. The transparency of P. dumerilii tissues at all developmental stages makes it easy to perform live microscopic imaging of all cell types. In addition, the slow-evolving genome of P. dumerilii and its phylogenetic position as a representative of the vast branch of Lophotrochozoans add to its evolutionary significance. Although P. dumerilii is amenable to transgenesis and CRISPR-Cas9 knockouts, its relatively long and indefinite life cycle, as well as its semelparous reproduction have been hindrances to its adoption as a reverse genetics model. To overcome this limitation, an adapted culturing method has been developed allowing much faster life cycling, with median reproductive age at 13-14 weeks instead of 25-35 weeks using the traditional protocol. A low worm density in boxes and a strictly controlled feeding regime are important factors for the rapid growth and health of the worms. This culture method has several advantages, such as being much more compact, not requiring air bubbling or an artificial moonlight regime for synchronized sexual maturation and necessitating only limited water change. A full protocol for worm care and handling is provided.

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.

Tunicates Illuminate the Enigmatic Evolution of Chordate Metallothioneins by Gene Gains and Losses, Independent Modular Expansions, and Functional Convergences

To investigate novel patterns and processes of protein evolution, we have focused in the metallothioneins (MTs), a singular group of metal-binding, cysteine-rich proteins that, due to their high degree of sequence diversity, still represents a "black hole" in Evolutionary Biology. We have identified and analyzed more than 160 new MTs in nonvertebrate chordates (especially in 37 species of ascidians, 4 thaliaceans, and 3 appendicularians) showing that prototypic tunicate MTs are mono-modular proteins with a pervasive preference for cadmium ions, whereas vertebrate and cephalochordate MTs are bimodular proteins with diverse metal preferences. These structural and functional differences imply a complex evolutionary history of chordate MTs-including de novo emergence of genes and domains, processes of convergent evolution, events of gene gains and losses, and recurrent amplifications of functional domains-that would stand for an unprecedented case in the field of protein evolution.

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.

Reconstruction of proto-vertebrate, proto-cyclostome and proto-gnathostome genomes provides new insights into early vertebrate evolution

Ancient polyploidization events have had a lasting impact on vertebrate genome structure, organization and function. Some key questions regarding the number of ancient polyploidization events and their timing in relation to the cyclostome-gnathostome divergence have remained contentious. Here we generate de novo long-read-based chromosome-scale genome assemblies for the Japanese lamprey and elephant shark. Using these and other representative genomes and developing algorithms for the probabilistic macrosynteny model, we reconstruct high-resolution proto-vertebrate, proto-cyclostome and proto-gnathostome genomes. Our reconstructions resolve key questions regarding the early evolutionary history of vertebrates. First, cyclostomes diverged from the lineage leading to gnathostomes after a shared tetraploidization (1R) but before a gnathostome-specific tetraploidization (2R). Second, the cyclostome lineage experienced an additional hexaploidization. Third, 2R in the gnathostome lineage was an allotetraploidization event, and biased gene loss from one of the subgenomes shaped the gnathostome genome by giving rise to remarkably conserved microchromosomes. Thus, our reconstructions reveal the major evolutionary events and offer new insights into the origin and evolution of vertebrate genomes.

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.

Evolutionary Analysis of the Zinc Finger and Homeoboxes Family of Proteins Identifies Multiple Conserved Domains and a Common Early Chordate Ancestor

The Zinc Fingers and Homeoboxes (Zhx) proteins, Zhx1, Zhx2, and Zhx3, comprise a small family of proteins containing two amino-terminal C₂–H₂ zinc fingers and four or five carboxy-terminal homeodomains. These multiple homeodomains make Zhx proteins unusual because the majority of homeodomain-containing proteins contain a single homeodomain. Studies in cultured cells and mice suggest that Zhx proteins can function as positive or negative transcriptional regulators. Zhx2 regulates numerous hepatic genes, and all three Zhx proteins have been implicated in different cancers. Because Zhx proteins contain multiple predicted homeodomains, are associated with interesting physiological traits, and seem to be only present in the vertebrate lineage, we investigated the evolutionary history of this small family by comparing Zhx homologs from a wide range of chordates. This analysis indicates that the zinc finger motifs and homeodomains are highly similar among all Zhx proteins and also identifies additional Zhx-specific conserved regions, including a 13 amino acid amino-terminal motif that is nearly identical among all gnathostome Zhx proteins. We found single Zhx proteins in the sea lamprey (Petromyzon marinus) and in the nonvertebrate chordates sea squirt (Ciona intestinalis) and lancelet (Branchiostoma floridae); these Zhx proteins are most similar to gnathostome Zhx3. Based on our analyses, we propose that a duplication of the primordial Zhx gene gave rise to Zhx3 and the precursor to Zhx1 and Zhx2. A subsequent tandem duplication of this precursor generated Zhx1 and Zhx2 found in gnathostomes.

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.

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.

Molecular evolution of CRH and CRHR subfamily before the evolutionary origin of vertebrate

Corticotropin-releasing hormone (CRH) is well-cited for its important role in governing the stress responses via neuroendocrine system in vertebrates. After the identification of homologs of CRH receptor (CRHR) in both deuterostome and arthropod lineages, it was suggested that the ancestral homolog of CRH-CRHR molecular system is present in the bilaterian. However, homolog sequences from arthropods differ considerably from vertebrate CRH-like peptide sequences. Due to the significant difference between the biological system, as well as the gene regulatory network, of protostome and that of vertebrate, physiological studies on the protostomes may not provide important insight into the evolutionary history of vertebrate CRH system, while tunicate and amphioxus, two close relatives to vertebrate, which have diverged before two rounds of whole genome duplication (2WGDs) do. Given the identification of amphioxus CRH-like peptide by our group, this review aims to reexamine the current hypotheses on the evolution of CRH subfamily. It is generally accepted that paralogs of CRH and CRHR have been produced through 2WGDs, which occurred during the early vertebrate evolution. The identification of a single crh-like gene in amphioxi and tunicates by in silico search and the presence of two paralogons with a total of 5 crh-like genes in gnathostomes has shown that an additional duplication event might have happened to the ancestral crh-like gene before 2WGDs. On the other hand, the evolution of crhr gene subfamily appears to be mainly influenced by 2WGDs and only two receptor genes have been retained in the genomes of jawed vertebrates.

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.

Functional lability of RNA-dependent RNA polymerases in animals

RNA interference (RNAi) requires RNA-dependent RNA polymerases (RdRPs) in many eukaryotes, and RNAi amplification constitutes the only known function for eukaryotic RdRPs. Yet in animals, classical model organisms can elicit RNAi without possessing RdRPs, and only nematode RNAi was shown to require RdRPs. Here we show that RdRP genes are much more common in animals than previously thought, even in insects, where they had been assumed not to exist. RdRP genes were present in the ancestors of numerous clades, and they were subsequently lost at a high frequency. In order to probe the function of RdRPs in a deuterostome (the cephalochordate Branchiostoma lanceolatum), we performed high-throughput analyses of small RNAs from various Branchiostoma developmental stages. Our results show that Branchiostoma RdRPs do not appear to participate in RNAi: we did not detect any candidate small RNA population exhibiting classical siRNA length or sequence features. Our results show that RdRPs have been independently lost in dozens of animal clades, and even in a clade where they have been conserved (cephalochordates) their function in RNAi amplification is not preserved. Such a dramatic functional variability reveals an unexpected plasticity in RNA silencing pathways.

RNA interference (RNAi) is a conserved gene regulation system in eukaryotes. In non-animal eukaryotes, it necessitates RNA-dependent RNA polymerases (“RdRPs”). Among animals, only nematodes appear to require RdRPs for RNAi. Yet additional animal clades have RdRPs and it is assumed that they participate in RNAi. Here, we find that RdRPs are much more common in animals than previously thought, but their genes were independently lost in many lineages. Focusing on a species with RdRP genes (a cephalochordate), we found that it does not use them for RNAi. While RNAi is the only known function for eukaryotic RdRPs, our results suggest additional roles. Eukaryotic RdRPs thus have a complex evolutionary history in animals, with frequent independent losses and apparent functional diversification.

My Favorite Animal, Amphioxus: Unparalleled for Studying Early Vertebrate Evolution

Amphioxus represents the most basally divergent group in chordates and probably the best extant proxy to the ancestor of all chordates including vertebrates. The amphioxus, or lancelets, are benthic filter feeding marine animals and their interest as a model in research is due to their phylogenetic position and their anatomical and genetic stasis throughout their evolutionary history. From the first works in the 19th century to the present day, enormous progress is made mainly favored by technical development at different levels, from spawning induction and husbandry techniques, through techniques for studies of gene function or of the role of different signalling pathways through embryonic development, to functional genomics techniques. Together, these advances foretell a plethora of interesting developments in the world of research with the amphioxus model. Here, the discovery and development of amphioxus as a superb model organism in evolutionary and evolutionary-developmental biology are reviewed.

Chromosome evolution at the origin of the ancestral vertebrate genome

BACKGROUND

It has been proposed that more than 450 million years ago, two successive whole genome duplications took place in a marine chordate lineage before leading to the common ancestor of vertebrates. A precise reconstruction of these founding events would provide a framework to better understand the impact of these early whole genome duplications on extant vertebrates.

RESULTS

We reconstruct the evolution of chromosomes at the beginning of vertebrate evolution. We first compare 61 extant animal genomes to reconstruct the highly contiguous order of genes in a 326-million-year-old ancestral Amniota genome. In this genome, we establish a well-supported list of duplicated genes originating from the two whole genome duplications to identify tetrads of duplicated chromosomes. From this, we reconstruct a chronology in which a pre-vertebrate genome composed of 17 chromosomes duplicated to 34 chromosomes and was subject to seven chromosome fusions before duplicating again into 54 chromosomes. After the separation of the lineage of Gnathostomata (jawed vertebrates) from Cyclostomata (extant jawless fish), four more fusions took place to form the ancestral Euteleostomi (bony vertebrates) genome of 50 chromosomes.

CONCLUSIONS

These results firmly establish the occurrence of two whole genome duplications in the lineage that precedes the ancestor of vertebrates, resolving in particular the ambiguity raised by the analysis of the lamprey genome. This work provides a foundation for studying the evolution of vertebrate chromosomes from the standpoint of a common ancestor and particularly the pattern of duplicate gene retention and loss that resulted in the gene composition of extant vertebrate genomes.

Review: Structure, function and evolution of GnIH

Neuropeptides that possess the Arg-Phe-NH₂ motif at their C-termini (i.e., RFamide peptides) have been characterized in the nervous system of both invertebrates and vertebrates. In vertebrates, RFamide peptides make a family and consist of the groups of gonadotropin-inhibitory hormone (GnIH), neuropeptide FF (NPFF), prolactin-releasing peptide (PrRP), kisspeptin (kiss1 and kiss2), and pyroglutamylated RFamide peptide/26RFamide peptide (QRFP/26RFa). It now appears that these vertebrate RFamide peptides exert important neuroendocrine, behavioral, sensory, and autonomic functions. In 2000, GnIH was discovered as a novel hypothalamic RFamide peptide inhibiting gonadotropin release in quail. Subsequent studies have demonstrated that GnIH acts on the brain and pituitary to modulate reproductive physiology and behavior across vertebrates. To clarify the origin and evolution of GnIH, the existence of GnIH was investigated in agnathans, the most ancient lineage of vertebrates, and basal chordates, such as tunicates and cephalochordates (represented by amphioxus). This review first summarizes the structure and function of GnIH and other RFamide peptides, in particular NPFF having a similar C-terminal structure of GnIH, in vertebrates. Then, this review describes the evolutionary origin of GnIH based on the studies in agnathans and basal chordates.

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.

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.

Gene family innovation, conservation and loss on the animal stem lineage

Choanoflagellates, the closest living relatives of animals, can provide unique insights into the changes in gene content that preceded the origin of animals. However, only two choanoflagellate genomes are currently available, providing poor coverage of their diversity. We sequenced transcriptomes of 19 additional choanoflagellate species to produce a comprehensive reconstruction of the gains and losses that shaped the ancestral animal gene repertoire. We identified ~1944 gene families that originated on the animal stem lineage, of which only 39 are conserved across all animals in our study. In addition, ~372 gene families previously thought to be animal-specific, including Notch, Delta, and homologs of the animal Toll-like receptor genes, instead evolved prior to the animal-choanoflagellate divergence. Our findings contribute to an increasingly detailed portrait of the gene families that defined the biology of the Urmetazoan and that may underpin core features of extant animals.

All animals, from sea sponges and reef-building corals to elephants and humans, share a single common ancestor that lived over half a billion years ago. This single-celled predecessor evolved the ability to develop into a creature made up of many cells with specialized jobs. Reconstructing the steps in this evolutionary process has been difficult because the earliest animals were soft-bodied and microscopic and did not leave behind fossils that scientists can study.

Though their bodies have since disintegrated, many of the instructions for building the first animals live on in genes that were passed on to life forms that still exist. Scientists are trying to retrace those genes back to the first animal by comparing the genomes of living animals with their closest relatives, the choanoflagellates. Choanoflagellates are single-celled, colony-forming organisms that live in waters around the world. Comparisons with choanoflagellates may help scientists identify which genes were necessary to help animals evolve and diversify into so many different species. So far, 1,000 animal and two choanoflagellate genomes have been sequenced. But the gene repertoires of most species of choanoflagellates have yet to be analyzed.

Now, Richter et al. have cataloged the genes of 19 more species of choanoflagellates. This added information allowed them to recreate the likely gene set of the first animal and to identify genetic changes that occurred during animal evolution. The analyses showed that modern animals lost about a quarter of the genes present in their last common ancestor with choanoflagellates and gained an equal number of new genes. Richter et al. identified several dozen core animal genes that were gained and subsequently preserved throughout animal evolution. Many of these are necessary so that an embryo can develop properly, but the precise roles of some core genes remain a mystery. Most other genes that emerged in the first animals have been lost in at least one living animal.

The study of Richter et al. also showed that some very important genes in animals, including genes essential for early development and genes that help the immune system detect pathogens, predate animals. These key genes trace back to animals’ last common ancestor with choanoflagellates and may have evolved new roles in animals.

The conserved ancient role of chordate PIAS as a multilevel repressor of the NF-κB pathway

In vertebrates, PIAS genes encode versatile cellular regulators, with functions extremely complex and redundant. Here we try to understand their functions from an evolutionary perspective. we evaluate the sequences, expression and molecular functions of amphioxus PIAS genes and compare them with their vertebrate counterparts. Phylogenetic analysis suggests a single PIAS gene in ancestral chordates, which has been duplicated into four families (PIAS1-4) in vertebrates by 2R-WGD but remained single in a basal chordate (amphioxus). Amphioxus PIAS encodes two variants with and without a Serine/Threonine-rich tail, which are retained in human PIAS1-3 but lost in PIAS4. We show that amphioxus PIAS binds C-terminus of NF-κB Rel and blocks the DNA binding activity. In humans, such function is retained in PIAS1, altered in PIAS4, and lost in PIAS2-3. Instead, PIAS3 has evolved new ability to inhibit Rel by binding RHD and promoting SUMOylation. We show that amphioxus PIAS also inhibits NF-κB by binding with upstream signalling adaptor TICAM-like and MyD88. Finally, we verify that human PIAS1, 3 and 4, but not 2, were capable of these newly-discovered functions. Our study offers insight into the sub- and neo-functionalization of PIAS genes and suggests a conserved ancient role for chordate PIAS in NF-κB signalling.

Exploring gene expression changes in the amphioxus gill after poly(I:C) challenge using digital expression profiling

Amphioxus, a cephalochordate, is a key model animal for studying the evolution of vertebrate immunity. Recently, studies have revealed that microRNA (miRNA) expression profiles change significantly in the amphioxus gill after immune stimulation, but it remains largely unknown how gene expression responds to immune stress. Elucidating gene expression changes in the amphioxus gill will provide a deeper understanding of the evolution of gill immunity in vertebrates. Here, we used high-throughput RNA sequencing technology (RNA-seq) to conduct tag-based digital gene expression profiling (DGE) analyses of the gills of control Branchiostoma belcheri and of those exposed to the viral mimic, poly(I:C) (pIC). Six libraries were created for the control and treatment groups including three biological replicates per group. A total of 1999 differently expressed genes (DEGs) were obtained, with 571 and 1428 DEGs showing up- or down-regulation, respectively, in the treatment group. Enrichment analysis of gene ontology (GO) terms and pathways revealed that the DEGs were primarily related to immune and defense response, apoptosis, human disease, cancer, protein metabolism, enzyme activity, and regulatory processes. In addition, eight DEGs were randomly selected to validate the RNA-seq data using real-time quantitative PCR (qRT-PCR), and the results confirmed the accuracy of the RNA-seq approach. Next, we screened eight key responding genes to examine the dynamic changes in expression levels at different time points in more detail. The results indicated that expressions of TRADD, MARCH, RNF31, NF-κb, CYP450, TNFRSF6B, IFI and LECT1 were induced to participate in the antiviral response against pIC. This study provides a valuable resource for understanding the role of the amphioxus gill in antiviral immunity and the evolution of gill immunity in vertebrates.

Evidence for Conservation of the Calcitonin Superfamily and Activity-regulating Mechanisms in the Basal Chordate Branchiostoma floridae: INSIGHTS INTO THE MOLECULAR AND FUNCTIONAL EVOLUTION IN CHORDATES

The calcitonin (CT)/CT gene-related peptide (CGRP) family is conserved in vertebrates. The activities of this peptide family are regulated by a combination of two receptors, namely the calcitonin receptor (CTR) and the CTR-like receptor (CLR), and three receptor activity-modifying proteins (RAMPs). Furthermore, RAMPs act as escort proteins by translocating CLR to the cell membrane. Recently, CT/CGRP family peptides have been identified or inferred in several invertebrates. However, the molecular characteristics and relevant functions of the CTR/CLR and RAMPs in invertebrates remain unclear. In this study, we identified three CT/CGRP family peptides (Bf-CTFPs), one CTR/CLR-like receptor (Bf-CTFP-R), and three RAMP-like proteins (Bf-RAMP-LPs) in the basal chordate amphioxus (Branchiostoma floridae). The Bf-CTFPs were shown to possess an N-terminal circular region typical of the CT/CGRP family and a C-terminal Pro-NH2. The Bf-CTFP genes were expressed in the central nervous system and in endocrine cells of the midgut, indicating that Bf-CTFPs serve as brain and/or gut peptides. Cell surface expression of the Bf-CTFP-R was enhanced by co-expression with each Bf-RAMP-LP. Furthermore, Bf-CTFPs activated Bf-CTFP-R·Bf-RAMP-LP complexes, resulting in cAMP accumulation. These results confirmed that Bf-RAMP-LPs, like vertebrate RAMPs, are prerequisites for the function and translocation of the Bf-CTFP-R. The relative potencies of the three peptides at each receptor were similar. Bf-CTFP2 was a potent ligand at all receptors in cAMP assays. Bf-RAMP-LP effects on ligand potency order were distinct to vertebrate CGRP/adrenomedullin/amylin receptors. To the best of our knowledge, this is the first molecular and functional characterization of an authentic invertebrate CT/CGRP family receptor and RAMPs.

Evolution of the CNS myelin gene regulatory program

Myelin is a specialized subcellular structure that evolved uniquely in vertebrates. A myelinated axon conducts action potentials many times faster than an unmyelinated axon of the same diameter; for the same conduction speed, the unmyelinated axon would need a much larger diameter and volume than its myelinated counterpart. Hence myelin speeds information transfer and saves space, allowing the evolution of a powerful yet portable brain. Myelination in the central nervous system (CNS) is controlled by a gene regulatory program that features a number of master transcriptional regulators including Olig1, Olig2 and Myrf. Olig family genes evolved from a single ancestral gene in non-chordates. Olig2, which executes multiple functions with regard to oligodendrocyte identity and development in vertebrates, might have evolved functional versatility through post-translational modification, especially phosphorylation, as illustrated by its evolutionarily conserved serine/threonine phospho-acceptor sites and its accumulation of serine residues during more recent stages of vertebrate evolution. Olig1, derived from a duplicated copy of Olig2 in early bony fish, is involved in oligodendrocyte development and is critical to remyelination in bony vertebrates, but is lost in birds. The origin of Myrf orthologs might be the result of DNA integration between an invading phage or bacterium and an early protist, producing a fusion protein capable of self-cleavage and DNA binding. Myrf seems to have adopted new functions in early vertebrates - initiation of the CNS myelination program as well as the maintenance of mature oligodendrocyte identity and myelin structure - by developing new ways to interact with DNA motifs specific to myelin genes. This article is part of a Special Issue entitled SI: Myelin Evolution.

Differential activation of kiss receptors by Kiss1 and Kiss2 peptides in the sea bass

Two forms of kiss gene (kiss1 and kiss2) have been described in the teleost sea bass. This study assesses the cloning and characterization of two Kiss receptor genes, namely kissr2 and kissr3 (known as gpr54-1b and gpr54-2b, respectively), and their signal transduction pathways in response to Kiss1 and Kiss2 peptides. Phylogenetic and synteny analyses indicate that these paralogs originated by duplication of an ancestral gene before teleost specific duplication. The kissr2 and kissr3 mRNAs encode proteins of 368 and 378 amino acids, respectively, and share 53.1% similarity in amino acid sequences. In silico analysis of the putative promoter regions of the sea bass Kiss receptor genes revealed conserved flanking regulatory sequences among teleosts. Both kissr2 and kissr3 are predominantly expressed in brain and gonads of sea bass, medaka and zebrafish. In the testis, the expression levels of sea bass kisspeptins and Kiss receptors point to a significant variation during the reproductive cycle. In vitro functional analyses revealed that sea bass Kiss receptor signals are transduced both via the protein kinase C and protein kinase A pathway. Synthetic sea bass Kiss1-15 and Kiss2-12 peptides activated Kiss receptors with different potencies, indicating a differential ligand selectivity. Our data suggest that Kissr2 and Kissr3 have a preference for Kiss1 and Kiss2 peptides, respectively, thus providing the basis for future studies aimed at establishing their physiologic roles in sea bass.

Characterization of peptide QRFP (26RFa) and its receptor from amphioxus, Branchiostoma floridae

A peptide ending with RFamide (Arg-Phe-amide) was discovered independently by three different laboratories in 2003 and named 26RFa or QRFP. In mammals, a longer version of the peptide, 43 amino acids, was identified and found to bind to the orphan G protein-coupled receptor GPR103. We searched the genome database of Branchiostoma floridae (Bfl) for receptor sequences related to those that bind peptides ending with RFa or RYa (including receptors for NPFF, PRLH, GnIH, and NPY). One receptor clustered in phylogenetic analyses with mammalian QRFP receptors. The gene has 3 introns in Bfl and 5 in human, but all intron positions differ, implying that the introns were inserted independently. A QRFP-like peptide consisting of 25 amino acids and ending with RFa was identified in the amphioxus genome. Eight of the ten last amino acids are identical between Bfl and human. The prepro-QRFP gene in Bfl has one intron in the propeptide whereas the human gene lacks introns. The Bfl QRFP peptide was synthesized and the receptor was functionally expressed in human cells. The response was measured as inositol phosphate (IP) turnover. The Bfl QRFP peptide was found to potently stimulate the receptor's ability to induce IP turnover with an EC50 of 0.28nM. Also the human QRFP peptides with 26 and 43 amino acids were found to stimulate the receptor (1.9 and 5.1nM, respectively). Human QRFP with 26 amino acids without the carboxyterminal amide had dramatically lower potency at 1.3μM. Thus, we have identified an amphioxus QRFP-related peptide and a corresponding receptor and shown that they interact to give a functional response.

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

Vertebrates diverged from other chordates ~500 Myr ago and experienced successful innovations and adaptations, but the genomic basis underlying vertebrate origins are not fully understood. Here we suggest, through comparison with multiple lancelet (amphioxus) genomes, that ancient vertebrates experienced high rates of protein evolution, genome rearrangement and domain shuffling and that these rates greatly slowed down after the divergence of jawed and jawless vertebrates. Compared with lancelets, modern vertebrates retain, at least relatively, less protein diversity, fewer nucleotide polymorphisms, domain combinations and conserved non-coding elements (CNE). Modern vertebrates also lost substantial transposable element (TE) diversity, whereas lancelets preserve high TE diversity that includes even the long-sought RAG transposon. Lancelets also exhibit rapid gene turnover, pervasive transcription, fastest exon shuffling in metazoans and substantial TE methylation not observed in other invertebrates. These new lancelet genome sequences provide new insights into the chordate ancestral state and the vertebrate evolution.

The lancelet, or amphioxus, is an extant basal chordate that diverged from other chordate lineages about 550 million years ago. Here the authors sequence and assemble the diploid genome of a male adult of the Chinese lancelet, B. belcheri, and highlight genomic features that may have played an important role in the origin and evolution of vertebrates.

Two apextrin-like proteins mediate extracellular and intracellular bacterial recognition in amphioxus

Animals exploit different germ-line-encoded proteins with various domain structures to detect the signature molecules of pathogenic microbes. These molecules are known as pathogen-associated molecular patterns (PAMPs), and the host proteins that react with PAMPs are called pattern recognition proteins (PRPs). Here, we present a novel type of protein domain structure capable of binding to bacterial peptidoglycan (PGN) and the minimal PGN motif muramyl dipeptide (MDP). This domain is designated as apextrin C-terminal domain (ApeC), and its presence was confirmed in several invertebrate phyla and subphyla. Two apextrin-like proteins (ALP1 and ALP2) were identified in a basal chordate, the Japanese amphioxus Branchiostoma japonicum (bj). bjALP1 is a mucosal effector secreted into the gut lumen to agglutinate the Gram-positive bacterium Staphylococcus aureus via PGN binding. Neutralization of secreted bjALP1 by anti-bjALP1 monoclonal antibodies caused serious damage to the gut epithelium and rapid death of the animals after bacterial infection. bjALP2 is an intracellular PGN sensor that binds to TNF receptor-associated factor 6 (TRAF6) and prevents TRAF6 from self-ubiquitination and hence from NF-κB activation. MDP was found to compete with TRAF6 for bjALP2, which released TRAF6 to activate the NF-κB pathway. BjALP1 and bjALP2 therefore play distinct and complementary functions in amphioxus gut mucosal immunity. In conclusion, discovery of the ApeC domain and the functional analyses of amphioxus ALP1 and ALP2 allowed us to define a previously undocumented type of PRP that is represented across different animal phyla.

Genome-wide survey and expression analysis of the bHLH-PAS genes in the amphioxus Branchiostoma floridae reveal both conserved and diverged expression patterns between cephalochordates and vertebrates

Background

The bHLH-PAS transcription factors are found in both protostomes and deuterostomes. They are involved in many developmental and physiological processes, including regional differentiation of the central nervous system, tube-formation, hypoxia signaling, aromatic hydrocarbon sensing, and circadian rhythm regulation. To understand the evolution of these genes in chordates, we analyzed the bHLH-PAS genes of the basal chordate amphioxus (Branchiostoma floridae).

Results

From the amphioxus draft genome database, we identified ten bHLH-PAS genes, nine of which could be assigned to known orthologous families. The tenth bHLH-PAS gene could not be assigned confidently to any known bHLH family; however, phylogenetic analysis clustered this gene with arthropod Met family genes and two spiralian bHLH-PAS-containing sequences, suggesting that they may share the same ancestry. We examined temporal and spatial expression patterns of these bHLH-PAS genes in developing amphioxus embryos. We found that BfArnt, BfNcoa, BfSim, and BfHifα were expressed in the central nervous system in patterns similar to those of their vertebrate homologs, suggesting that their functions may be conserved. By contrast, the amphioxus BfAhr and BfNpas4 had expression patterns distinct from those in vertebrates. These results imply that there were changes in gene regulation after the divergence of cephalochordates and vertebrates.

Conclusions

We have identified ten bHLH-PAS genes from the amphioxus genome and determined the embryonic expression profiles for these genes. In addition to the nine currently recognized bHLH-PAS families, our survey suggests that the BfbHLHPAS-orphan gene along with arthropod Met genes and the newly identified spiralian bHLH-PAS-containing sequences represent an ancient group of genes that were lost in the vertebrate lineage. In a comparison with the expression patterns of the vertebrate bHLH-PAS paralogs, which are the result of whole-genome duplication, we found that although several members seem to retain conserved expression patterns during chordate evolution, many duplicated paralogs may have undergone subfunctionalization and neofunctionalization in the vertebrate lineage. In addition, our survey of amphioxus bHLH-PAS gene models from genome browser with experimentally verified cDNA sequences calls into question the accuracy of the current in silico gene annotation of the B. floridae genome.

Punctuated emergences of genetic and phenotypic innovations in eumetazoan, bilaterian, euteleostome, and hominidae ancestors

Phenotypic traits derive from the selective recruitment of genetic materials over macroevolutionary times, and protein-coding genes constitute an essential component of these materials. We took advantage of the recent production of genomic scale data from sponges and cnidarians, sister groups from eumetazoans and bilaterians, respectively, to date the emergence of human proteins and to infer the timing of acquisition of novel traits through metazoan evolution. Comparing the proteomes of 23 eukaryotes, we find that 33% human proteins have an ortholog in nonmetazoan species. This premetazoan proteome associates with 43% of all annotated human biological processes. Subsequently, four major waves of innovations can be inferred in the last common ancestors of eumetazoans, bilaterians, euteleostomi (bony vertebrates), and hominidae, largely specific to each epoch, whereas early branching deuterostome and chordate phyla show very few innovations. Interestingly, groups of proteins that act together in their modern human functions often originated concomitantly, although the corresponding human phenotypes frequently emerged later. For example, the three cnidarians Acropora, Nematostella, and Hydra express a highly similar protein inventory, and their protein innovations can be affiliated either to traits shared by all eumetazoans (gut differentiation, neurogenesis); or to bilaterian traits present in only some cnidarians (eyes, striated muscle); or to traits not identified yet in this phylum (mesodermal layer, endocrine glands). The variable correspondence between phenotypes predicted from protein enrichments and observed phenotypes suggests that a parallel mechanism repeatedly produce similar phenotypes, thanks to novel regulatory events that independently tie preexisting conserved genetic modules.

At the transition from invertebrates to vertebrates, a novel GnRH-like peptide emerges in amphioxus

Gonadotropin-releasing hormone (GnRH) is a critical reproductive regulator in vertebrates. Homologous peptides are also found in invertebrates, with a variety of characterized functions. In the amphioxus, an invertebrate that provides the best model for the transition to vertebrates, four GnRH receptors (GnRHRs) were previously described, but their native ligands were not identified. Using a more sensitive search methodology with hidden Markov models, we identified the first GnRH-like peptide confirmed in the amphioxus Branchiostoma floridae. This peptide specifically activated one of the four GnRHRs. Although the primary structure of this peptide was divergent from any previously isolated GnRH peptide, the minimal conserved residues found in all other GnRH superfamily members were retained. The peptide was immunolocalized in proximity of the central canal of the anterior nerve cord, a region where other neuropeptides and receptors have been found. Additionally, the amphioxus GnRH-like gene was positioned in a locus surrounded by syntenic homologs of the human GnRH paralogon. The amphioxus GnRH-like peptide, with its distinct primary structure, activated a receptor with equal potency to multiple ligands that span the GnRH superfamily.

The evolutionary portrait of metazoan NAD salvage

Nicotinamide Adenine Dinucleotide (NAD) levels are essential for cellular homeostasis and survival. Main sources of intracellular NAD are the salvage pathways from nicotinamide, where Nicotinamide phosphoribosyltransferases (NAMPTs) and Nicotinamidases (PNCs) have a key role. NAMPTs and PNCs are important in aging, infection and disease conditions such as diabetes and cancer. These enzymes have been considered redundant since either one or the other exists in each individual genome. The co-occurrence of NAMPT and PNC was only recently detected in invertebrates though no structural or functional characterization exists for them. Here, using expression and evolutionary analysis combined with homology modeling and protein-ligand docking, we show that both genes are expressed simultaneously in key species of major invertebrate branches and emphasize sequence and structural conservation patterns in metazoan NAMPT and PNC homologues. The results anticipate that NAMPTs and PNCs are simultaneously active, raising the possibility that NAD salvage pathways are not redundant as both are maintained to fulfill the requirement for NAD production in some species.

Evolution of new characters after whole genome duplications: insights from amphioxus

Additional copies of genes resulting from two whole genome duplications at the base of the vertebrates have been suggested as enabling the evolution of vertebrate-specific structures such as neural crest, a midbrain/hindbrain organizer and neurogenic placodes. These structures, however, did not evolve entirely de novo, but arose from tissues already present in an ancestral chordate. This review discusses the evolutionary history of co-option of old genes for new roles in vertebrate development as well as the relative contributions of changes in cis-regulation and in protein structure. Particular examples are the FoxD, FGF8/17/18 and Pax2/5/8 genes. Comparisons with invertebrate chordates (amphioxus and tunicates) paint a complex picture with co-option of genes into new structures occurring both after and before the whole genome duplications. In addition, while cis-regulatory changes are likely of primary importance in evolution of vertebrate-specific structures, changes in protein structure including alternative splicing are non-trivial.