- Research article
- Open Access
Deep whole-genome resequencing sheds light on the distribution and effect of amphioxus SNPs
BMC Genomic Data volume 23, Article number: 26 (2022)
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.
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.
Our work provides explanation for “amphioxus polymorphism paradox” in a micro view, and generates an enhanced genomic dataset for amphioxus research.
The amphioxus, also known as lancelet, is the modern representative of the subphylum Cephalochordata, providing evolutionary insight into the origin of vertebrates [1, 2]. It is recognized that cephalochordates, urochordates, and vertebrates belonging to the phylum Chordata, evolved from a common ancestor that lived about 550 million years ago. Amphioxus has a vertebrate-like but simpler body plan, and different to most chordates, their genomes remained intact without any WGD events . Hence, they are considered to be intermediate between invertebrates and vertebrates, and widely utilized as model organism for exploring vertebrate origin [1,2,3,4,5,6].
Previous studies have revealed a contradictory phenomenon that extreme phenotypic similarity and high-level genetic diversity [2, 3] [7,8,9,10,11] co-exist in amphioxus populations. Early publication claimed that the polymorphism rate of amphioxus might be as high as 5.37% , while they are commonly observed to share similar phenotypic characteristics , including body length, asymmetric body shape, number of oral cirri, etc. . Taking Branchiostoma belcheri inhabiting Xiamen waters as an example, most adult individuals among these amphioxi present no apparent morphological differences, even females and males only differ in the reproductive organs , however genomic analysis revealed there might be up to 1 mutation site per 30 bases on their genome . The “polymorphism paradox” has raised interests in studying the actual impact of amphioxus genomic mutation, which may provide not only evolutionary insights but also guidance for several potential applications, such as breeding .
Variant calling based on resequencing of main amphioxus species is an approach to seek the reasons behind “amphioxus polymorphism paradox”, for it helps researchers to get whole-genome variant distribution, where affected genes and regulatory elements can be identified for functional analysis. High-quality reference genome is the prerequisite of variant calling, while it is certain that representative amphioxus genomes have been successfully assembled [2, 3] [15, 16], some of which were assembled to chromosome level. In 2021, haploid genomes of three amphioxus species, B. belcheri (20 chromosomes), B. japonicum (18 chromosomes), and B. floridae (19 chromosomes), got completely resolved , further contributing to the amphioxus genomic toolbox. These works offer more options for subsequent resequencing experiments as researchers have higher chance to access the reference genome best representing their samples. On the basis above, several attempts have been made to build amphioxus variant datasets [10, 11], yet the low sequencing depth causing loss of polymorphisms  might make them difficult to provide reliable evidence.
In addition, the booming protein structure prediction algorithms have paved the way for exploring detailed variant effect on the protein level. As evidenced by the results of the biennial Critical Assessment of protein Structure Prediction (CASP), structure prediction has seen substantial progress in recent years . In CASP14 (2020), the program AlphaFold2  achieved a record score of 92.4, vastly more accurate than competing methods. Up to 2021, AlphaFold2 has predicted over 98.5% of human protein structures, and its results are recognized as reliable structure sources by some public protein databases such as UniProt. Compared to other prediction algorithms including RoseTTAFold , AlphaFold2 presents to be single in function and huge in performance overhead, while its unmatched accuracy maintains its “gold medal” in this field. As the rapid structure modelling solution challenging x-ray crystallography and cryo-electron microscopy, AlphaFold2 can directly transform sequences to structures, enabling researchers to observe SNP-brought structural difference intuitively.
Here we resequenced 21 amphioxus (B. belcheri) genomes to over 100X depth using the Illumina HiSeq X Ten platform, selected the new haploid genome  as the reference for joint SNP calling, and utilized AlphaFold2 to detect structural variant brought by SNPs. 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. Our work provides explanation for “amphioxus polymorphism paradox” in a micro view, and generates an enhanced genomic dataset for amphioxus research.
Reads passing quality control were mapped to the haploid reference genome assembled into 20 chromosomes using BWA , then GATK  was employed for SNP calling on each sample. The genomic mapping results are listed in Table 1. The offline read bases are 1,268.5G totally, and average sequencing depth is 124.66X. The polymorphism ratio is about 1 SNP per 30 bases, which does not differ much between individuals.
Joint calling by GATK and PLINK  was launched after all individual SNP callings, where common SNPs (minor allele frequency > 0.05) were selected and got annotation from SnpEff . There are totally 59,895,836 SNPs among 21 individuals, of which 11,541,148 are identified as common SNPs, approximately one per 34 bases in B. belcheri genome. Statistics of SNP distribution and effect are illustrated in Fig. 1. Figure 1a presents the SNP distribution centres on chromosomes in the form of heatmap. It’s obvious that for each chromosome there are one or two relatively wide distribution peaks. There is little difference among mutation rates of each chromosome, except for Chr1 and Chr3, whose polymorphism rates are relatively lower (Table S1).
Figure 1b-c reveal that most common SNPs are located in non-coding regions. 33.15% of them hit intron regions and 14.85% hit intergenic regions, and these two region types were found with most variants in previous studies . Meanwhile it’s conspicuous that up to 34.2% get identified in up/downstream regions to genes, where cis-regulatory elements are located. SNPs hitting exon regions account for 8.07% (Fig. 1b), among which synonymous variants turn out to be slightly more than missense variants (Fig. 1c). Nonsynonymous variants having potential high impact on gene function are counted as a limited proportion of the total (Table S2). It should be mentioned that intron variants, synonymous variants as well as regulatory element variants mainly affect gene expression  rather than structure, thus numerous as they are, their impacts are considered not as high as nonsynonymous variants . In addition, differences of SNP distributions among female, hermaphrodite and male amphioxus individuals are observed to be far from significant (Fig S1, Table S3), indicating that SNP distribution is not closely related to their sex.
GO enrichment analysis was separately performed on top-1000 genes with largest number of intron variants, synonymous variants and missense variants, results of which are summarized as Fig. 2. The top enriched terms of biological process and molecular function are related to ion transport, signal transduction and homophilic cell adhesion via plasma membrane adhesion molecules, where highly missense-mutated genes outnumber those synonymous-mutated. It indicates that variant selectivity might exist in related amphioxus gene families. Besides, enriched cellular component terms show that mutation frequently appears in amphioxus cytoskeleton along with its related complex. Enrichment results of these genes are generally consistent with their function distribution according to COG annotations (Fig S2).
The pN/pS ratios of genes with a minimum of five synonymous SNPs were calculated, distribution of which is displayed in Fig. 3a. Genes with pN/pS > 1.5 account for about 2.14%, presenting to be a relatively large positive tail. KEGG enrichment analysis was performed on them, and results (Fig. 3b) reveal that these genes be involved in hormone regulation, mainly including growth-related pathways (MAPK, PI3K-Akt, etc.) and reproduction-related pathways (GnRH, Estrogen, Oxytocin, etc.). Meanwhile they might participate in response to hypoxic stress (HIF-1 and FoxO).
AlphaFold2 was employed to help evaluate protein structural variation brought by SNPs. Based on gene annotation, SNP distribution statistics and enrichment results, 100 domains were picked out as samples for structure analysis. Comparison between each raw-mutated structure pair reveals that only few SNPs could cause significant change in protein spatial structure, examples of which are displayed in Fig. 4. The largest variation is found in the SEA domain of a gene from Chr2, where SNPs at the beginning of the first α-helix interfere with the formation of a crucial hydrogen bond, triggering to conformational change. Several SNPs can also affect secondary structure such as α-helix and β-sheet, as observed in genes from Chr10 and Chr14. Except for them, other SNPs seem to modify intermolecular interaction at most, failing to exert great effects on protein structure, let alone the function (Data S2).
The polymorphism rate of B. belcheri in our results (about 3%) turns out to be larger than that of previous researches using low-depth sequencing data [10, 11], though still smaller than several reported high numbers (5.37% for B. belcheri and 4% for B. floridae of the same genus). The improvement of reference genome and sequencing method can exert positive influence on SNP calling, for it helps researchers to capture more mutation sites, thus reducing false negative rate and laying foundation for deeper studies on a larger sample size. In fact, genomic polymorphism rates of various marine lives are identified as high, for example, 1–1.5% in sea squirts (urochordate) , 4–5% in sea urchin (echinoderm) , 2–2.5% in oyster (lophotrochozoan) , etc. Generally, high mutation rate might offer them better natural physiological tolerance, enabling them to extend their species lifespan . The actual impact of such numerous variants on these fantastic animals is a complex but valuable scientific question.
As for amphioxus, it indeed owns widely-distributed SNPs, while two factors limit their functional effect. For one thing, the majority of SNPs (more than 90%) hit non-coding regions, and those located in protein-coding regions mainly trigger to synonymous variants, in that pN/pS ratios of most genes (more than 75%) are below 1 (Fig S3), the symbol of neutral or purifying selection. This is consistent with other species for coding regions commonly maintain highly conserved to perform relevant functions and avoid potential damaging. For another, structure modelling reveals that most SNP-affected sequences are translated to proteins without significant structural changes. It seems that inside these living fossils, the role of translation is not only conversion, but also protection, in that effects of DNA-level mutation get weakened on the protein structure level to a certain extent. These above together shape amphioxus into phenotypically convergent species. Whether there is a similar situation in other marine organisms with high mutation rate needs to be further explored.
While it’s still discovered that in contrast to the overall variant distribution, genes identified as highly missense-mutated turn out to outnumber those synonymous-mutated in several pathways, mainly involved in signal transduction, ion transport and cell adhesion (See Table S4 for details). As the filter feeder, amphioxus utilize cilia to filter seawater and gather a variety of algae as food. They engulf food particles via phagocytic intracellular digestion mechanism [29, 30], which means they need to efficiently degrade algal toxin while intracellularly transforming algae into small molecules. Ion channels are the targets of many algal toxins , hence related genes inevitably suffer from high selective pressure, meanwhile their immune system and digestive system have to maintain active almost all the time. These may all contribute to the functional variation for survival. The SEA domain shown in Fig. 4, one of few samples detected with significant structural change in this work, has been recorded to participate in immune regulation and digestive enzyme production [32, 33], serving as evidence for the assumption above. Furthermore, pN/pS analysis sheds light on the positive selection acting in amphioxus endocrine-related genes, as illustrated in Fig. 3b. Previous studies indicate that amphioxus own several primitive regulatory axes, for example, their GH-IGF signalling system consisting of Hatschek’s pit and hepatic cecum regulates their growth and osmotic pressure [34, 35], analogous to pituitary-hepatic axis in vertebrates; brain vesicle, Hatschek’s pit and gonad constitute their reproductive endocrine axis, sharing functional similarity with HPG (hypothalamic-pituitary–gonadal) axis [36,37,38]; The coordination of their endostyle and Hatschek’s pit is potentially the early version of TSH-TH signalling system [39,40,41]; etc. Positive selection upon related genes might help to complicate and refine their signal transduction networks, so as to promote evolutionary processes in regulatory mechanisms.
In view of their extremely long history of existence and evolution, there might be few spaces for amphioxus to gain spontaneous genetic optimization now. From another perspective, it indicates the great threats from booming human activities, even the construction of a seawall could lead to the rapid disappearance of a hundreds-year-old B. belcheri fishing ground . Therefore, the aim of studying amphioxus SNPs shouldn’t be confined to supplementing contents of biological treatise. Several variants hidden in our large dataset are probably the key to save amphioxus, of which scientists can take advantages to breed individuals and enlarge their population. Nevertheless, it’s bound to be an effort-consuming work to mine useful information from such seas of data, challenging the knowledge base and experimental infrastructure of any single group, besides an expanded sample size is required for subsequent confirmatory and translational researches.
Our work not only sheds light on the distribution and effect of amphioxus SNPs but also makes progress in the explanation for “amphioxus polymorphism paradox”. It is proposed that influence of high-level genetic polymorphism be sharply weakened on the protein level. Numerous and selective as they are, amphioxus SNPs lack the macromolecular basis of impact on phenotypic characteristics.
It is expected that the upgraded sequencing technology with much deeper coverage, the reference genome with higher assembly level, and the advanced protein structure prediction algorithm make our new genomic dataset a valuable resource for exploring amphioxus biology as well as the origin of vertebrates. Due to limitation of computing power and storage, there remain lots of unmined information such as rare variants in it, thus joint efforts from further research are always welcome.
Additional annotation for reference genome
The haploid genome of Branchiostoma belcheri is provided by Fujian Key Laboratory of Special Marine Bio-resources Sustainable Utilization , including full sequences, protein-coding sequences, and annotation in GFF3 format (GCA_019207075.1, Data S1).
HMMER 3.1 package was employed for Pfam  annotation of reference genome. GO  annotations were converted from Pfam annotation via the external2go service at http://current.geneontology.org/ontology/external2go/pfam2go. KEGG and COG annotations were acquired from online services KAAS (https://www.genome.jp/kegg/kaas/) and eggNOG-mapper (http://eggnog-mapper.embl.de) [45,46,47].
Sample preparation, DNA extraction, and sequencing
Twenty-one B. belcheri individuals, 10 male, 10 female and 1 hermaphrodite, were obtained from shoals (within 1km2) in Zhanjiang, Guangdong province during the summer breeding season, and kept at 28–32 °C at the Beihai Marine Station of Nanjing University, Guangxi, China.
Twenty-one samples were anesthetized and fixed by ethanol (100%, AR) before muscles were dissected. Genomic DNA was extracted separately using a QIAamp® DNA mini kit (Qiagen, Germany) following the standard manufacturer's protocol . The purity and concentration of total DNA were determined with a NanoDrop spectrophotometer (NanoDrop, Wilmington, DE). DNA integrity was assessed by agarose gel electrophoresis. Briefly, the DNA sample was fragmented using a Covaris ultrasonic processor (Covaris, USA) to a size of ~ 350 bp, then the fragmented DNA was end repaired, “A”-tailed, and ligated with the full-length adaptor for Illumina sequencing with further PCR amplification. The concentrations of the constructed libraries were initially measured and diluted to 1 ng/μl by Qubit®2.0 (Life technologies, USA). Then, an Agilent Bioanalyzer 2100 system (Agilent, USA) was used to check the insert size of the libraries. To ensure the quality of these constructed libraries, the SYBR green qRT-PCR protocol was used with a Kapa Probe Fast qPCR kit (Kapa Biosystems, USA) to accurately dose the effective concentrations of the libraries. Finally, these libraries were sequenced on the Illumina HiSeq X Ten platform (Illumina, USA) by the Novogene Bioinformatics Institute, Beijing, China.
Quality control criteria  were applied to remove low-quality reads:
Removal of reads with more than 10% unidentified nucleotides (N);
Removal of reads containing more than 50% of bases with a Phred score ≤ 5;
Removal of putative PCR duplicate reads generated by PCR amplification.
Then sequences got further filtering via fastp  (-q 20 -c).
Individual SNP calling
After removing low-quality reads, the clean paired-end reads were mapped to the haploid reference genome using BWA  (-M -k 19), then secondary or supplementary alignments were filtered by sambamba  (- “ "not (secondary_alignment or supplementary”)" -p -l 9).
The remaining mapped reads were sorted and converted into BAM format files using SAMtools . PCR duplicates were marked using the GATK MarkDuplicates module.
Two rounds of individual SNP calling were performed using GATK, including following steps:
RealignerTargetCreator and IndelRealigner modules applied run to reduce the false-positive variants where alignment error occurred across overlapping reads;
HaplotypeCaller and VariantFiltration modules were applied to detect SNPs (QD < 10.0, MQ < 50.0, FS > 10.0, MQRankSum < − 5.0, ReadPosRankSum < − 8.0);
BaseRecalibrator and ApplyBQSR modules were applied to generate recalibrated bam files for each individual;
HaplotypeCaller module was applied again to detect variants (–emit-ref-confidence GVCF).
Joint calling and SNP annotation
All GVCF files were processed by CombineGVCFs and GenotypeGVCFs modules to generate the population genotype file, then SNPs were selected using SelectVariants and VariantFiltration modules (QD < 2.0, MQ < 40.0, FS > 60.0, MQRankSum < − 12.5, ReadPosRankSum < − 8.0, QUAL < 30). Common SNPs in the population were selected by PLINK  (–maf 0.05, –geno 0.05, –hwe 1e-4), and got annotation by SnpEff . In addition to annotation in VCF format, SnpEff generated the statistics of SNPs’ genomic locations and coding effect defined in SO  terms, which were extracted to plot SNP distribution heatmap and other statistic charts (Fig. 1, Fig S1).
Gene functional enrichment analysis
GO enrichment analysis was performed on top-1000 (P value < 1e-5) genes with largest number of intron variant (SO:0,001,627), synonymous variant (SO:0,001,819), and missense variant (SO:0,001,583) using clusterProfiler  (pvalueCutoff = 0.05, pAdjustMethod = 'BH', qvalueCutoff = 0.2).
The pN/pS ratios of genes were calculated referring to existing approach , and KEGG enrichment analysis using the same tool and parameters as GO analysis were performed on those with pN/pS > 1.5 and a minimum of five synonymous SNPs.
Protein structural variation analysis
Sequences of missense-variant-enriched genes were extracted based on enrichment analysis and SNP statistics above, then 100 crucial protein domains according to Pfam annotation were chosen for structure analysis. For each raw-mutated sequence pair, AlphaFold2  (default parameters) was run to get the highest-score structures (ranked_0.pdb) to make comparison.
Availability of data and materials
Raw data from our whole-genome sequencing are available at NCBI (PRJNA742127).
Total SNP set is available at figshare (https://doi.org/10.6084/m9.figshare.18833234.v1).
Single Nucleotide Polymorphism
Holland LZ, Albalat R, Azumi K, et al. The amphioxus genome illuminates vertebrate origins and cephalochordate biology. Genome Res. 2008;18(7):1100–11. https://doi.org/10.1101/gr.073676.107.
Putnam NH, Butts T, Ferrier DE, et al. The amphioxus genome and the evolution of the chordate karyotype. Nature. 2008;453(7198):1064–71. https://doi.org/10.1038/nature06967.
Huang S, Chen Z, Yan X, et al. Decelerated genome evolution in modern vertebrates revealed by analysis of multiple lancelet genomes. Nat Commun. 2014;5(1):5896. https://doi.org/10.1038/ncomms6896.
Marlétaz F, Firbas PN, Maeso I, et al. Amphioxus functional genomics and the origins of vertebrate gene regulation. Nature. 2018;564(7734):64–70. https://doi.org/10.1038/s41586-018-0734-6.
Yuan S, Ruan J, Huang S, Chen S, Xu A. Amphioxus as a model for investigating evolution of the vertebrate immune system. Dev Comp Immunol. 2015;48(2):297–305. https://doi.org/10.1016/j.dci.2014.05.004.
Zhang QL, Zhang GL, Yuan ML, et al. A Phylogenomic Framework and Divergence History of Cephalochordata Amphioxus. Front Physiol. 2018;9:1833. https://doi.org/10.3389/fphys.2018.01833.
Li WY, Zhong J, Xu W, Wang Y-Q. Microsatellite DNA marker development and genetic diversity of Branchiostoma belcheri in Xiamen waters. Mar Biol Res. 2011;7(8):826–31. https://doi.org/10.1080/17451000.2011.569553.
Yue JX, Yu J-K, Putnam NH, Holland LZ. The transcriptome of an amphioxus, asymmetron lucayanum, from the Bahamas: A window into chordate evolution. Genome Biol Evol. 2014;6(10):2681–96. https://doi.org/10.1093/gbe/evu212.
Yue JX, Kozmikova I, Ono H, et al. Conserved noncoding elements in the most distant genera of cephalochordates: the goldilocks principle. Genome Biol Evol. 2016;8(8):2387–405. https://doi.org/10.1093/gbe/evw158.
Bi C, Lu N, Han T, et al. Whole-genome resequencing of twenty Branchiostoma belcheri individuals provides a brand-new variant dataset for Branchiostoma. Biomed Res Int. 2020;2020:3697342. https://doi.org/10.1155/2020/3697342.
Bi C, Lu N, Huang Z, Chen J, He C, Lu Z. Whole-genome resequencing reveals the pleistocene temporal dynamics of Branchiostoma belcheri and Branchiostoma floridae populations. Ecol Evol. 2020;10(15):8210–24. https://doi.org/10.1002/ece3.6527.
Zhang S. Evolutionary Biology of Amphioxus: Tracing Origin of Vertebrate. 1st ed. Beijing: Science Press; 2020. p. 3–16.
Jin D. Amphioxus. Fujian: Fujian People’s Publishing House; 1957. p. 1–2.
Zhang X, Yuan J, Sun Y, et al. Penaeid shrimp genome provides insights into benthic adaptation and frequent molting. Nat Commun. 2019;10(1):356. https://doi.org/10.1038/s41467-018-08197-4.
Simakov O, Marlétaz F, Yue JX, et al. Deeply conserved synteny resolves early events in vertebrate evolution. Nat Ecol Evol. 2020;4(6):820–30. https://doi.org/10.1038/s41559-020-1156-z.
Huang Z, Xu L, Cai C, Zhou Y, Liu J, Zhu Z, Kang W, Chen D, Pei S, Xue T, Cen W, Shi C, Wu X, Huang Y, Xu C, Yan Y, Yang Y, He W, Hu X, Zhang Y, Chen Y, Bi C, He C, Xue L, Xiao S, Yue Z, Jiang Y, Yu J-K, Jarvis ED, Li G, Lin G, Zhang Q, Zhou Q. Three amphioxus reference genomes reveal gene and chromosome evolution of chordates [Internet]. bioRxiv. 2022. Available from: https://doi.org/10.1101/2022.01.04.475009.
Tunyasuvunakool K, Adler J, Wu Z, et al. Highly accurate protein structure prediction for the human proteome. Nature. 2021;596(7873):590–6. https://doi.org/10.1038/s41586-021-03828-1.
Jumper J, Evans R, Pritzel A, et al. Highly accurate protein structure prediction with AlphaFold. Nature. 2021;596(7873):583–9. https://doi.org/10.1038/s41586-021-03819-2.
Baek M, DiMaio F, Anishchenko I, et al. Accurate prediction of protein structures and interactions using a three-track neural network. Science. 2021;373(6557):871–6. https://doi.org/10.1126/science.abj8754.
Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25(14):1754–60. https://doi.org/10.1093/bioinformatics/btp324.
McKenna A, Hanna M, Banks E, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297–303. https://doi.org/10.1101/gr.107524.110.
Chang CC, Chow CC, Tellier LC, Vattikuti S, Purcell SM, Lee JJ. Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience. 2015;4:7. https://doi.org/10.1186/s13742-015-0047-8.
Cingolani P, Platts A, Wang le L, Coon M, Nguyen T, Wang L, Land SJ, Lu X, Ruden DM. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly (Austin). 2012;6(2):80-92. https://doi.org/10.4161/fly.19695.
Rasal KD, Chakrapani V, Pandey AK, et al. Status and future perspectives of single nucleotide polymorphisms (SNPs) markers in farmed fishes: Way ahead using next generation sequencing. Gene Reports. 6:81–86 https://doi.org/10.1016/j.genrep.2016.12.004
Dehal P, Satou Y, Campbell RK, et al. The draft genome of Ciona intestinalis: insights into chordate and vertebrate origins. Science. 2002;298(5601):2157–67. https://doi.org/10.1126/science.1080049.
Sea Urchin Genome Sequencing Consortium, Sodergren E, Weinstock GM, et al. The genome of the sea urchin Strongylocentrotus purpuratus. Science. 2006;314(5801):941-52. https://doi.org/10.1126/science.1133609. Erratum in: Science. 2007 Feb 9;315(5813):766.
Sauvage C, Bierne N, Lapègue S, Boudry P. Single Nucleotide polymorphisms and their relationship to codon usage bias in the Pacific oyster Crassostrea gigas. Gene. 2007;406(1–2):13–22. https://doi.org/10.1016/j.gene.2007.05.011.
Raup DM. The role of extinction in evolution. Proc Natl Acad Sci U S A. 1994;91(15):6758–63. https://doi.org/10.1073/pnas.91.15.6758.
He C, Han T, Liao X, et al. Phagocytic intracellular digestion in amphioxus (Branchiostoma) [published correction appears in Proc Biol Sci. 2018 Jun 27;285(1881):]. Proc Biol Sci. 2018;285(1880):20180438 https://doi.org/10.1098/rspb.2018.0438
The sea creature that swallows its food whole - twice. Nature. 2018;558(7709):165 https://doi.org/10.1038/d41586-018-05396-3
Van Dolah FM. Marine algal toxins: origins, health effects, and their increased occurrence. Environ Health Perspect. 2000;108 Suppl 1(Suppl 1):133–141 doi:https://doi.org/10.1289/ehp.00108s1133
Kitamoto Y, Yuan X, Wu Q, McCourt DW, Sadler JE. Enterokinase, the initiator of intestinal digestion, is a mosaic protease composed of a distinctive assortment of domains. Proc Natl Acad Sci U S A. 1994;91(16):7588–92. https://doi.org/10.1073/pnas.91.16.7588.
Palmai-Pallag T, Khodabukus N, Kinarsky L, Leir SH, Sherman S, Hollingsworth MA, Harris A. The role of the SEA (sea urchin sperm protein, enterokinase and agrin) module in cleavage of membrane-tethered mucins. FEBS J. 2005;272(11):2901-11. https://doi.org/10.1111/j.1742-4658.2005.04711.x.
Li M, Jiang C, Zhang Y, Zhang S. Activities of Amphioxus GH-Like Protein in Osmoregulation: Insight into Origin of Vertebrate GH Family. Int J Endocrinol. 2017;2017:9538685. https://doi.org/10.1155/2017/9538685.
Liu M, Zhang S. Amphioxus IGF-like peptide induces mouse muscle cell development via binding to IGF receptors and activating MAPK and PI3K/Akt signaling pathways. Mol Cell Endocrinol. 2011;343(1–2):45–54. https://doi.org/10.1016/j.mce.2011.06.004.
Wang P, Wang M, Ji G, Yang S, Zhang S, Liu Z. Demonstration of a Functional Kisspeptin/Kisspeptin Receptor System in Amphioxus With Implications for Origin of Neuroendocrine Regulation. Endocrinology. 2017;158(5):1461–73. https://doi.org/10.1210/en.2016-1848.
Tello JA, Sherwood NM. Amphioxus: beginning of vertebrate and end of invertebrate type GnRH receptor lineage. Endocrinology. 2009;150(6):2847–56. https://doi.org/10.1210/en.2009-0028.
Tjoa LT, Welsch U. Electron microscopical observations on Kölliker’s and Hatschek’s pit and on the wheel organ in the head region of Amphioxus (Branchiostoma lanceolatum). Cell Tissue Res. 1974;153(2):175–87. https://doi.org/10.1007/BF00226606.
Wang P, Liu S, Yang Q, Liu Z, Zhang S. Functional Characterization of Thyrostimulin in Amphioxus Suggests an Ancestral Origin of the TH Signaling Pathway. Endocrinology. 2018;159(10):3536–48. https://doi.org/10.1210/en.2018-00550.
Wang S, Zhang S, Zhao B, Lun L. Up-regulation of C/EBP by thyroid hormones: a case demonstrating the vertebrate-like thyroid hormone signaling pathway in amphioxus. Mol Cell Endocrinol. 2009;313(1–2):57–63. https://doi.org/10.1016/j.mce.2009.08.024.
Paris M, Hillenweck A, Bertrand S, Delous G, Escriva H, Zalko D, Cravedi JP, Laudet V. Active metabolism of thyroid hormone during metamorphosis of amphioxus. Integr Comp Biol. 2010;50(1):63-74. https://doi.org/10.1093/icb/icq052. Epub 24 May 2010.
Jin D, Cheng Z, Deng Y. Branchiostoma belcheri is on the verge of extinction in Liuwudian. Fujian Aquatic Products. 1987;1:32–3. https://doi.org/10.14012/j.cnki.fjsc.1987.01.008.
El-Gebali S, Mistry J, Bateman A, Eddy SR, Luciani A, Potter SC, Qureshi M, Richardson LJ, Salazar GA, Smart A, Sonnhammer ELL, Hirsh L, Paladin L, Piovesan D, Tosatto SCE, Finn RD. The Pfam protein families database in 2019. Nucleic Acids Res. 2019;47(D1):D427-32. https://doi.org/10.1093/nar/gky995.
The Gene Ontology Consortium. The Gene Ontology Resource: 20 years and still GOing strong. Nucleic Acids Res. 2019;47(D1):D330–8. https://doi.org/10.1093/nar/gky1055.
Kanehisa M, Furumichi M, Tanabe M, Sato Y, Morishima K. KEGG: new perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. 2017;45(D1):D353–61. https://doi.org/10.1093/nar/gkw1092.
Cantalapiedra CP, Hernández-Plaza A, Letunic I, Bork P, Huerta-Cepas J. eggNOG-mapper v2: functional annotation, orthology assignments, and domain prediction at the metagenomic scale. Mol iolE. 2021;38(12):5825–9. https://doi.org/10.1093/molbev/msab293.
Huerta-Cepas J, Szklarczyk D, Heller D, et al. eggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses. Nucleic Acids Res. 2019;47(D1):D309-D314 https://doi.org/10.1093/nar/gky1085
Yan L, Yang M, Guo H, et al. Single-cell RNA-Seq profiling of human preimplantation embryos and embryonic stem cells. Nat Struct Mol Biol. 2013;20(9):1131–9. https://doi.org/10.1038/nsmb.2660.
Chen S, Zhou Y, Chen Y, Gu J. fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 2018;34(17):i884–90. https://doi.org/10.1093/bioinformatics/bty560.
Tarasov A, Vilella AJ, Cuppen E, Nijman IJ, Prins P. Sambamba: fast processing of NGS alignment formats. Bioinformatics. 2015;31(12):2032–4. https://doi.org/10.1093/bioinformatics/btv098.
Li H, Handsaker B, Wysoker A, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25(16):2078–9. https://doi.org/10.1093/bioinformatics/btp352.
Eilbeck K, Lewis SE, Mungall CJ, et al. The sequence ontology: a tool for the unification of genome annotations. Genome Biol. 2005;6(5):R44. https://doi.org/10.1186/gb-2005-6-5-r44.
Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7. https://doi.org/10.1089/omi.2011.0118.
Hao Y, Washburn JD, Rosenthal J, et al. Patterns of population variation in two paleopolyploid eudicot lineages suggest that dosage-based selection on homeologs is long-lived. Genome Biol Evol. 2018;10(3):999–1011. https://doi.org/10.1093/gbe/evy061.
We are grateful for Prof. Ping Zhou and Dr. Ke Xiao’s technical support in the computing environment.
This work was supported by the open research fund of State Key Laboratory of Bioelectronics, Southeast University [Sklb2021-k02].
Ethics approval and consent to participate
All amphioxus samples were collected and processed in accordance with local laws for aquatic animal protection and approved by the Ethics Committee of Institutional Animal Care and Use Committee of Nanjing Medical University (protocol code IACUC-1910003 and date of approval is 10 October 2019).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Protein-coding sequences and annotation of the haploid reference genome. Data S2. Samples for AlphaFold2 analysis. Fig S1. Pie charts of SNPs detected in each type of genomic region in amphioxus of different sexes. Fig S2. Statistics of COG functional classes of top-1000 genes with largest number of synonymous variants (a), missense variants (b), and intron variants (c). A: RNA processing and modification; B: Chromatin structure and dynamics; C: Energy production and conversion; D: Cell cycle control, cell division, chromosome partitioning; E: Amino acid transport and metabolism; F: Nucleotide transport and metabolism; G: Carbohydrate transport and metabolism; H: Coenzyme transport and metabolism; I: Lipid transport and metabolism; J: Translation, ribosomal structure and biogenesis; K: Transcription; L: Replication, recombination and repair; M: Cell wall/membrane/envelope biogenesis; N: Cell motility; O: Posttranslational modification, protein turnover, chaperones; P: Inorganic ion transport and metabolism; Q: Secondary metabolites biosynthesis, transport and catabolism; R: General function prediction only; S: Function unknown; T: Signal transduction mechanisms; U: Intracellular trafficking, secretion, and vesicular transport; V: Defense mechanisms; W: Extracellular structures; X: Unnamed protein; Y: Nuclear structure; Z: Cytoskeleton. Fig S3. Box plot of total pN/pS ratio distribution. Table S1. Number of common SNPs on each chromosome. Table S2. Number of variants affecting each gene calculated via SnpEff. Table S3. Statistics of SNPs detected in each type of genomic region in each amphioxus sample. Table S4. Full tables of GO enrichment results.
About this article
Cite this article
Zhu, Y., Lu, N., Chen, JY. et al. Deep whole-genome resequencing sheds light on the distribution and effect of amphioxus SNPs. BMC Genom Data 23, 26 (2022). https://doi.org/10.1186/s12863-022-01038-w
- Whole-genome resequencing