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

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.

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

Background

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

Results

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

Conclusions

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

Supplementary Information

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

Functional Conservation and Genetic Divergence of Chordate Glycinergic Neurotransmission: Insights from Amphioxus Glycine Transporters

Glycine is an important neurotransmitter in vertebrates, performing both excitatory and inhibitory actions. Synaptic levels of glycine are tightly controlled by the action of two glycine transporters, GlyT1 and GlyT2, located on the surface of glial cells and neurons, respectively. Only limited information is available on glycinergic neurotransmission in invertebrates, and the evolution of glycinergic neurotransmission is poorly understood. Here, by combining phylogenetic and gene expression analyses, we characterized the glycine transporter complement of amphioxus, an important invertebrate model for studying the evolution of chordates. We show that amphioxus possess three glycine transporter genes. Two of these (GlyT2.1 and GlyT2.2) are closely related to GlyT2 of vertebrates, whereas the third (GlyT) is a member of an ancestral clade of deuterostome glycine transporters. GlyT2.2 expression is predominantly non-neural, whereas GlyT and GlyT2.1 are widely expressed in the amphioxus nervous system and are differentially expressed, respectively, in neurons and glia. Vertebrate glycinergic neurons express GlyT2 and glia GlyT1, suggesting that the evolution of the chordate glycinergic system was accompanied by a paralog-specific inversion of gene expression. Despite this genetic divergence between amphioxus and vertebrates, we found strong evidence for conservation in the role glycinergic neurotransmission plays during larval swimming, the implication being that the neural networks controlling the rhythmic movement of chordate bodies may be homologous.

Detection of the First Epoxyalcohol Synthase/Allene Oxide Synthase (CYP74 Clan) in the Lancelet (Branchiostoma belcheri, Chordata)

The CYP74 clan cytochromes (P450) are key enzymes of oxidative metabolism of polyunsaturated fatty acids in plants, some Proteobacteria, brown and green algae, and Metazoa. The CYP74 enzymes, including the allene oxide synthases (AOSs), hydroperoxide lyases, divinyl ether synthases, and epoxyalcohol synthases (EASs) transform the fatty acid hydroperoxides to bioactive oxylipins. A novel CYP74 clan enzyme CYP440A18 of the Asian (Belcher’s) lancelet (Branchiostoma belcheri, Chordata) was biochemically characterized in the present work. The recombinant CYP440A18 enzyme was active towards all substrates used: linoleate and α-linolenate 9- and 13-hydroperoxides, as well as with eicosatetraenoate and eicosapentaenoate 15-hydroperoxides. The enzyme specifically converted α-linolenate 13-hydroperoxide (13-HPOT) to the oxiranyl carbinol (9Z,11R,12R,13S,15Z)-11-hydroxy-12,13-epoxy-9,15-octadecadienoic acid (EAS product), α-ketol, 12-oxo-13-hydroxy-9,15-octadecadienoic acid (AOS product), and cis-12-oxo-10,15-phytodienoic acid (AOS product) at a ratio of around 35:5:1. Other hydroperoxides were converted by this enzyme to the analogous products. In contrast to other substrates, the 13-HPOT and 15-HPEPE yielded higher proportions of α-ketols, as well as the small amounts of cyclopentenones, cis-12-oxo-10,15-phytodienoic acid and its higher homologue, dihomo-cis-12-oxo-3,6,10,15-phytotetraenoic acid, respectively. Thus, the CYP440A18 enzyme exhibited dual EAS/AOS activity. The obtained results allowed us to ascribe a name “B. belcheri EAS/AOS” (BbEAS/AOS) to this enzyme. BbEAS/AOS is a first CYP74 clan enzyme of Chordata species possessing AOS activity.

The population genomics of adaptive loss of function

Discoveries of adaptive gene knockouts and widespread losses of complete genes have in recent years led to a major rethink of the early view that loss-of-function alleles are almost always deleterious. Today, surveys of population genomic diversity are revealing extensive loss-of-function and gene content variation, yet the adaptive significance of much of this variation remains unknown. Here we examine the evolutionary dynamics of adaptive loss of function through the lens of population genomics and consider the challenges and opportunities of studying adaptive loss-of-function alleles using population genetics models. We discuss how the theoretically expected existence of allelic heterogeneity, defined as multiple functionally analogous mutations at the same locus, has proven consistent with empirical evidence and why this impedes both the detection of selection and causal relationships with phenotypes. We then review technical progress towards new functionally explicit population genomic tools and genotype-phenotype methods to overcome these limitations. More broadly, we discuss how the challenges of studying adaptive loss of function highlight the value of classifying genomic variation in a way consistent with the functional concept of an allele from classical population genetics.