A Phylogenomic Framework and Divergence History of Cephalochordata Amphioxus

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.

Improved Understanding of the Role of Gene and Genome Duplications in Chordate Evolution With New Genome and Transcriptome Sequences

Comparative approaches to understanding chordate genomes have uncovered a significant role for gene duplications, including whole genome duplications (WGDs), giving rise to and expanding gene families. In developmental biology, gene families created and expanded by both tandem and WGDs are paramount. These genes, often involved in transcription and signalling, are candidates for underpinning major evolutionary transitions because they are particularly prone to retention and subfunctionalisation, neofunctionalisation, or specialisation following duplication. Under the subfunctionalisation model, duplication lays the foundation for the diversification of paralogues, especially in the context of gene regulation. Tandemly duplicated paralogues reside in the same regulatory environment, which may constrain them and result in a gene cluster with closely linked but subtly different expression patterns and functions. Ohnologues (WGD paralogues) often diversify by partitioning their expression domains between retained paralogues, amidst the many changes in the genome during rediploidisation, including chromosomal rearrangements and extensive gene losses. The patterns of these retentions and losses are still not fully understood, nor is the full extent of the impact of gene duplication on chordate evolution. The growing number of sequencing projects, genomic resources, transcriptomics, and improvements to genome assemblies for diverse chordates from non-model and under-sampled lineages like the coelacanth, as well as key lineages, such as amphioxus and lamprey, has allowed more informative comparisons within developmental gene families as well as revealing the extent of conserved synteny across whole genomes. This influx of data provides the tools necessary for phylogenetically informed comparative genomics, which will bring us closer to understanding the evolution of chordate body plan diversity and the changes underpinning the origin and diversification of vertebrates.

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.

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.

Differential expression pattern of two Brachyury genes in amphioxus embryos

Cephalochordate amphioxus contain two Brachyury genes (AmphiBra1 and AmphiBra2). Using probes from the highly conserved coding regions, a summation of their expression profiles in amphioxus embryos have been investigated by several previous studies. However, their respective expression patterns have not been determined up to date. We here address this issue using both qRT-PCR and in situ hybridization methods (with probes from the divergent untranslated regions). qRT-PCR detected a very low maternal expression for AmphiBra2, but not for AmphiBra1. Zygotic expression of both genes are activated around early gastrula stage and change in a similar pattern at subsequent stages. However, compared to AmphiBra1, the expression level of AmphiBra2 is much higher in all examined stages of embryos; in some extreme cases an over fifty-times difference is observed. In situ hybridization and embryonic sections reveal that while AmphiBra2 is highly expressed in the blastopore, the tail bud and the notochord, AmphiBra1 is weakly transcribed only in the notochord. Our results show that the two Brachyury genes, resulted from a lineage-specific duplication in amphioxus, have evolved different embryonic expression profiles.

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.

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

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

The Complete Mitogenome of Pyrrhocoris tibialis (Hemiptera: Pyrrhocoridae) and Phylogenetic Implications

We determined the complete mitogenome of Pyrrhocoris tibialis (Hemiptera: Heteroptera: Pyrrhocoridae) to better understand the diversity and phylogeny within Pentatomomorpha, which is the second largest infra-order of Heteroptera. Gene content, gene arrangement, nucleotide composition, codon usage, ribosomal RNA (rRNA) structures, and sequences of the mitochondrial transcription termination factor were well conserved in Pyrrhocoroidea. Different protein-coding genes have been subject to different evolutionary rates correlated with the G + C content. The size of control regions (CRs) was highly variable among mitogenomes of three sequenced Pyrrhocoroidea species, with the P. tibialis CR being the largest. All the transfer RNA genes found in Pyrrhocoroidea had the typical clover leaf secondary structure, except for trnS1 (AGN), which lacked the dihydrouridine arm and possessed an unusual anticodon stem (9 bp vs. the normal 5 bp). A total of three different phylogenetic relationships among the five super-families of Pentatomomorpha were obtained using three analytical methods (MrBayes and RAxML under site-homogeneous models and PhyloBayes under a site-heterogeneous CAT + GTR model) and two mitogenomic datasets (nucleotides and amino acids). The tree topology test using seven methods statistically supported a phylogeny of (Aradoidea + (Pentatomoidea + (Lygaeoidea + (Pyrrhocoroidea + Coreoidea)))) as the best topology, as recognized by both RAxML and MrBayes based on the two datasets.