ECB-ART-45806
PLoS Biol
2017 Oct 12;1510:e2003790. doi: 10.1371/journal.pbio.2003790.
Show Gene links
Show Anatomy links
The sea cucumber genome provides insights into morphological evolution and visceral regeneration.
Zhang X
,
Sun L
,
Yuan J
,
Sun Y
,
Gao Y
,
Zhang L
,
Li S
,
Dai H
,
Hamel JF
,
Liu C
,
Yu Y
,
Liu S
,
Lin W
,
Guo K
,
Jin S
,
Xu P
,
Storey KB
,
Huan P
,
Zhang T
,
Zhou Y
,
Zhang J
,
Lin C
,
Li X
,
Xing L
,
Huo D
,
Sun M
,
Wang L
,
Mercier A
,
Li F
,
Yang H
,
Xiang J
.
???displayArticle.abstract???
Apart from sharing common ancestry with chordates, sea cucumbers exhibit a unique morphology and exceptional regenerative capacity. Here we present the complete genome sequence of an economically important sea cucumber, A. japonicus, generated using Illumina and PacBio platforms, to achieve an assembly of approximately 805 Mb (contig N50 of 190 Kb and scaffold N50 of 486 Kb), with 30,350 protein-coding genes and high continuity. We used this resource to explore key genetic mechanisms behind the unique biological characters of sea cucumbers. Phylogenetic and comparative genomic analyses revealed the presence of marker genes associated with notochord and gill slits, suggesting that these chordate features were present in ancestral echinoderms. The unique shape and weak mineralization of the sea cucumber adult body were also preliminarily explained by the contraction of biomineralization genes. Genome, transcriptome, and proteome analyses of organ regrowth after induced evisceration provided insight into the molecular underpinnings of visceral regeneration, including a specific tandem-duplicated prostatic secretory protein of 94 amino acids (PSP94)-like gene family and a significantly expanded fibrinogen-related protein (FREP) gene family. This high-quality genome resource will provide a useful framework for future research into biological processes and evolution in deuterostomes, including remarkable regenerative abilities that could have medical applications. Moreover, the multiomics data will be of prime value for commercial sea cucumber breeding programs.
???displayArticle.pubmedLink??? 29023486
???displayArticle.pmcLink??? PMC5638244
???displayArticle.link??? PLoS Biol
Species referenced: Echinodermata
Genes referenced: egfIII fgf fras1 hoxa7l lars1 LOC100887844 LOC100892375 LOC115920572 LOC115923268 LOC579909 LOC586799 LOC590297 tubgcp2
???attribute.lit??? ???displayArticles.show???
![]() |
Fig 1. A schematic representation of the genomic characteristics of A. japonicus.Track 1: 22 linkage groups of the genome. Track 2: anchored scaffolds to each linkage group. Long bars represent scaffolds with a length > 500 Kb; short bars represent scaffolds with a length ⤠500 Kb. Track 3: protein-coding genes located on scaffolds. Red stands for genes on the forward strand, and green stands for genes on the reverse strand. Track 4: distributions of 5 significantly expanded gene families in the genome. The 5 gene families are fibrinogen-related proteins (FREPs) (F), retrovirus-related Pol polyprotein from transposon (R), nucleotide-binding oligomerization domain-like receptors (NLR) family caspase recruitment domain (CARD) domain-containing protein (N), tyrosine-protein kinase receptor (T), and zinc finger CysCysHisCys (CCHC) domain-containing protein (Z). These expanded gene families show a clustered distribution in the genome. FREPs mainly clustered at LG02, and NLR chiefly accumulated at LG04. Track 5: distribution of microRNA (miRNA) in the genome. Most miRNAs are distributed in clusters across the genome (LG01, LG10, LG11, LG12, and LG22). The data underlying Fig 1 can be found in S3 Data. |
![]() |
Fig 2. Comparative genomic analysis between A. japonicus and other metazoans.(A) Phylogenetic placement of A. japonicus within the metazoan tree. The numbers on the branches indicate the number of gene gains (+) and the number of gene losses (â), which are also displayed as bar plots: gene gain (in green), gene loss (in red), and the remaining gene families (in blue). The divergence times were estimated and displayed below the phylogenetic tree. Image credits: Robert Michniewicz; Kobie Mercury-Clarke; Vector Open Stock, Patrick Narbonne, David E. Simpson, John B. Gurdon; Lars Simonsen; Freshwater and Marine Image Bank; Encyclopædia Britannica; public domain; Jerry Kirkhart; authors' own; Johny Ha; Virginia Gewin; public domain; public domain; Cnidaria; Michael Eitel, Hans-Jürgen Osigus, Rob DeSalle, Bernd Schierwater; Maja Adamska. (B) The shared and unique gene families in 4 species of Ambulacraria are shown in the Venn diagram. There are 763 gene families shared by 3 echinoderms. (C) Comparison of the gene repertoire of 10 metazoan genomes. Here, "1:1" indicates single-copy genes; "X:X" indicates orthologous genes present in multiple copies in all the 9 species, where X means 1 or more orthologs per species; and "patchy" indicates the existence of other orthologs that are present in at least 1 genome. A. japonicus and S. purpuratus show a similar distribution of gene repertoire. The data underlying Fig 2 can be found in S3 Data. |
![]() |
Fig 3. Gene families and gene clusters related to notochord and pharyngeal gill slit formation.(A) Distribution pattern of genes related to notochord and stomochord formation in metazoans. These genes are all present in deuterostomes, but some of them are present in nondeuterostomes. (B) Phylogenetic analysis of the fibroblast growth factor (FGF) gene family in metazoans. Different genes of the FGF gene family are distinguished with different colors. The FGF gene family is important for transcriptional activation of bruchury in notochord formation, but only 1 such gene exists in the genome of A. japonicus (red arrow) and the other 2 echinoderms (yellow background), S. purpuratus (Spur_15908_340) and A. planci (Apla_2356_405). The data underlying Fig 3A and 3B can be found in S3 Data. (C) A gene cluster related to the pharyngeal gill slit. The genes connected on a line indicate that they are clustered in order in the genome. This gene cluster is conserved among deuterostomes, whereas it is incomplete and shows poor synteny within nondeuterostomes. Image credits: Robert Michniewicz; Kobie Mercury-Clarke; Patrick Narbonne, David E. Simpson, John B. Gurdon; Vector Open Stock; Lars Simonsen; Freshwater and Marine Image Bank; Encyclopædia Britannica; public domain; Jerry Kirkhart; authors' own; Martin Cooper; public domain; public domain; Cnidaria; Freshwater and Marine Image Bank; Maja Adamska; Johny Ha; Virginia Gewin; Michael Eitel, Hans-Jürgen Osigus, Rob DeSalle, Bernd Schierwater. |
![]() |
Fig 4. The conserved homeobox gene clusters in the genome of A. japonicus.A. planci, Amphiura filiformis, and Florometra serratissima belong to Echinodermata class Asteroidea, Ophiuroidea, and Crinoidea, respectively. The general phylogenetic tree shows the Hox and ParaHox cluster order and genes. Hox4 and Hox6 are not found in the genome of A. japonicus. Different colors indicate different Hox gene groups, including anterior (blue), central (green), and posterior (red). The arrow indicates the direction of transcription. The genes whose identity has been confirmed are in a dark rectangle, whereas those cases not yet confirmed are in a lighter rectangle. The Hox cluster images of Branchiostoma floridae, S. purpuratus, A. filiformis, and F. serratissima were modified from Byrne et al. 2016 [34]. Image credits: "Introduction to zoology; a guide to the study of animals, for the use of secondary schools" (1900) Macmillan; Freshwater and Marine Image Bank; Encyclopædia Britannica; Jerry Kirkhart; authors' own; public domain; Freshwater and Marine Image Bank; National Oceanic and Atmospheric Administration; Martin Cooper. |
![]() |
Fig 5. Gene families related to skeletogenesis and biomineralization in the genome of A. japonicus.(A) Comparison of the skeleton from A. japonicus and S. purpuratus. (B) The comparison of the copy numbers of biomineralization-related genes in 3 echinoderms, A. japonicus, S. purpuratus, and A. planci, and the hemichordate S. kowalevskii. Biomineralization-related genes are apparently contracted in A. japonicus. (C) The expression levels of biomineralization-related genes at early developmental stages of A. japonicus (red) and S. purpuratus (black). The expression level of these genes at different development stages is calculated relative to a housekeeping gene (tubulin). Compared with sea urchin, there is a significant gene loss in the A. japonicus genome, and the expression levels for those genes at different development stages are relatively low. The data underlying Fig 5B and 5C can be found in S3 Data. |
![]() |
Fig 6. Prostatic secretory protein of 94 amino acids (PSP94)-like gene family in the genome of A. japonicus.(A) A heatmap showing the expression profile of tandem-duplicated gene clusters during the intestinal regeneration of A. japonicus. It depicts the cluster ID, the repeat number, and the expression profile. The color represents the relative expression level (fragments per kilobase of transcript per million mapped reads [FPKM] value) of the gene clusters during the intestinal regeneration. The Scaffold889 with 11 tandem-duplicated genes (blue asterisk) shows significant up-regulation post evisceration. (B) The detailed expression of 11 PSP94-like genes located on Scaffold889 at different time points post evisceration. The scale covers log expression values. Genes are clustered by Euclidean distance of the log10(FPKM+1) value and grouped with average-linkage clustering. (C) The distribution and alignment of the PSP94-like gene family. Eleven PSP94-like genes are tandemly arranged on Scaffold889. Arrows denote the transcriptional orientation. The amino acid sequence alignment and domain prediction of proteins coded by these PSP94-like genes show highly conserved cysteine residues in the PSP94 domain. The data underlying Fig 6B and 6C can be found in S3 Data. |
![]() |
Fig 7. Genes involved in the intestine regeneration of A. japonicus.(A) Regeneration diagram showing A. japonicus (a) undergoing evisceration, (b) immediately post evisceration, and (c) after complete recovery. (B) A heatmap showing the expression profile of molecular events applicable to the intestinal regeneration of A. japonicus. It depicts the molecular events identified in this study in relation to the morphological diversification of regeneration in a diagrammatic sketch. The height and the color of trapezium represent the relative expression level (fragments per kilobase of transcript per million mapped reads [FPKM] value per gene) during the regeneration. High-expressed and low-expressed genes are labeled in red and green, respectively. Extracellular matrix (ECM)-related genes (collagen and fibropellin), signaling pathways (Wnt, bone morphogenetic protein-related genes, and epidermal growth factor-related genes), and myogenesis-related genes (tubulin) are activated at the early stage of regeneration (0â3 days post evisceration [dpe]). During the middle stage of regeneration (3â7 dpe), factors related to hormonal regulation are up-regulated. In the late stages of regeneration, ECM-related genes (tenascin, FRAS1, and collagen) and myogenesis-related genes (actin and myosin) show up-regulated expression. The physiological events happening during regeneration are shown below the x-axis. The data underlying Fig 7B can be found in S3 Data. FREP, fibrinogen-related protein. |
References [+] :
Baughman,
Genomic organization of Hox and ParaHox clusters in the echinoderm, Acanthaster planci.
2014, Pubmed,
Echinobase
Baughman, Genomic organization of Hox and ParaHox clusters in the echinoderm, Acanthaster planci. 2014, Pubmed , Echinobase
Benson, Tandem repeats finder: a program to analyze DNA sequences. 1999, Pubmed
Birney, GeneWise and Genomewise. 2004, Pubmed
Birney, Using GeneWise in the Drosophila annotation experiment. 2000, Pubmed
Blair, Molecular phylogeny and divergence times of deuterostome animals. 2005, Pubmed
Burke, A genomic view of the sea urchin nervous system. 2006, Pubmed , Echinobase
Byrne, Evolution of a pentameral body plan was not linked to translocation of anterior Hox genes: the echinoderm HOX cluster revisited. 2016, Pubmed , Echinobase
Cameron, Unusual gene order and organization of the sea urchin hox cluster. 2006, Pubmed , Echinobase
Cao, Hepassocin regulates cell proliferation of the human hepatic cells L02 and hepatocarcinoma cells through different mechanisms. 2011, Pubmed
Chapman, The dynamic genome of Hydra. 2010, Pubmed
Cheers, P16 is an essential regulator of skeletogenesis in the sea urchin embryo. 2005, Pubmed , Echinobase
Cisternas, Expression of Hox4 during development of the pentamerous juvenile sea star, Parvulastra exigua. 2009, Pubmed , Echinobase
Conesa, Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. 2005, Pubmed
Darriba, jModelTest 2: more models, new heuristics and parallel computing. 2012, Pubmed
Davey, Genome-wide genetic marker discovery and genotyping using next-generation sequencing. 2011, Pubmed
De Bie, CAFE: a computational tool for the study of gene family evolution. 2006, Pubmed
Dehal, The draft genome of Ciona intestinalis: insights into chordate and vertebrate origins. 2002, Pubmed
Dominguez, Paired gill slits in a fossil with a calcite skeleton. 2002, Pubmed , Echinobase
Douzery, The timing of eukaryotic evolution: does a relaxed molecular clock reconcile proteins and fossils? 2004, Pubmed
Ebert, Longevity and lack of senescence in the red sea urchin Strongylocentrotus franciscanus. 2008, Pubmed , Echinobase
Edgar, MUSCLE: multiple sequence alignment with high accuracy and high throughput. 2004, Pubmed
Ettensohn, Lessons from a gene regulatory network: echinoderm skeletogenesis provides insights into evolution, plasticity and morphogenesis. 2009, Pubmed , Echinobase
Friedländer, miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. 2012, Pubmed
Fritzenwanker, The Fox/Forkhead transcription factor family of the hemichordate Saccoglossus kowalevskii. 2014, Pubmed , Echinobase
Godwin, Regeneration, tissue injury and the immune response. 2006, Pubmed
Grabherr, Full-length transcriptome assembly from RNA-Seq data without a reference genome. 2011, Pubmed
Gross, The role of Brachyury (T) during gastrulation movements in the sea urchin Lytechinus variegatus. 2001, Pubmed , Echinobase
Guindon, A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. 2003, Pubmed
Haas, Automated eukaryotic gene structure annotation using EVidenceModeler and the Program to Assemble Spliced Alignments. 2008, Pubmed
Hall, The crown-of-thorns starfish genome as a guide for biocontrol of this coral reef pest. 2017, Pubmed , Echinobase
Hedges, Tree of life reveals clock-like speciation and diversification. 2015, Pubmed
Huang, CAP3: A DNA sequence assembly program. 1999, Pubmed
Imokawa, A critical role for thrombin in vertebrate lens regeneration. 2004, Pubmed
Jo, Draft genome of the sea cucumber Apostichopus japonicus and genetic polymorphism among color variants. 2017, Pubmed , Echinobase
Karandikar, beta-Microseminoprotein-related molecules may participate in formation of the mesoderm in the chick embryo. 2003, Pubmed
Koran, Impact of family structure and common environment on heritability estimation for neuroimaging genetics studies using Sequential Oligogenic Linkage Analysis Routines. 2014, Pubmed
Kotov, Mesozoic fossils (>145 Mya) suggest the antiquity of the subgenera of Daphnia and their coevolution with chaoborid predators. 2011, Pubmed
Krzywinski, Circos: an information aesthetic for comparative genomics. 2009, Pubmed
Langmead, Fast gapped-read alignment with Bowtie 2. 2012, Pubmed
Li, OrthoMCL: identification of ortholog groups for eukaryotic genomes. 2003, Pubmed
Li, Fast and accurate short read alignment with Burrows-Wheeler transform. 2009, Pubmed
Li, The Sequence Alignment/Map format and SAMtools. 2009, Pubmed
Li, De novo assembly of human genomes with massively parallel short read sequencing. 2010, Pubmed
Liu, The role of the Pax1/9 gene in the early development of amphioxus pharyngeal gill slits. 2015, Pubmed
Livesey, Netrins and netrin receptors. 1999, Pubmed
Livingston, A genome-wide analysis of biomineralization-related proteins in the sea urchin Strongylocentrotus purpuratus. 2006, Pubmed , Echinobase
Love, Fossil steroids record the appearance of Demospongiae during the Cryogenian period. 2009, Pubmed
Lynch, The evolutionary fate and consequences of duplicate genes. 2000, Pubmed
Lyons, Specification to biomineralization: following a single cell type as it constructs a skeleton. 2014, Pubmed , Echinobase
Marçais, A fast, lock-free approach for efficient parallel counting of occurrences of k-mers. 2011, Pubmed
Mashanov, Radial glial cells play a key role in echinoderm neural regeneration. 2013, Pubmed , Echinobase
Mashanov, Postembryonic organogenesis of the digestive tube: why does it occur in worms and sea cucumbers but fail in humans? 2014, Pubmed , Echinobase
Mashanov, Gut regeneration in holothurians: a snapshot of recent developments. 2011, Pubmed , Echinobase
Mashanov, Transcriptomic changes during regeneration of the central nervous system in an echinoderm. 2014, Pubmed , Echinobase
Mathers, Multiple global radiations in tadpole shrimps challenge the concept of 'living fossils'. 2013, Pubmed
McCauley, Development of an embryonic skeletogenic mesenchyme lineage in a sea cucumber reveals the trajectory of change for the evolution of novel structures in echinoderms. 2012, Pubmed , Echinobase
McCauley, A conserved gene regulatory network subcircuit drives different developmental fates in the vegetal pole of highly divergent echinoderm embryos. 2010, Pubmed , Echinobase
Mofarrahi, Angiopoietin-1 enhances skeletal muscle regeneration in mice. 2015, Pubmed
Moriya, KAAS: an automatic genome annotation and pathway reconstruction server. 2007, Pubmed
Nagase, Antithrombin III-independent effect of depolymerized holothurian glycosaminoglycan (DHG) on acute thromboembolism in mice. 1997, Pubmed , Echinobase
Nakamura-Ishizu, Extracellular matrix protein tenascin-C is required in the bone marrow microenvironment primed for hematopoietic regeneration. 2012, Pubmed
Nawrocki, Infernal 1.1: 100-fold faster RNA homology searches. 2013, Pubmed
Ortiz-Pineda, Gene expression profiling of intestinal regeneration in the sea cucumber. 2009, Pubmed , Echinobase
Parra, CEGMA: a pipeline to accurately annotate core genes in eukaryotic genomes. 2007, Pubmed
Pascual-Anaya, Evolution of Hox gene clusters in deuterostomes. 2013, Pubmed
Patel, NGS QC Toolkit: a toolkit for quality control of next generation sequencing data. 2012, Pubmed
Purcell, The cost of being valuable: predictors of extinction risk in marine invertebrates exploited as luxury seafood. 2014, Pubmed , Echinobase
Röttinger, FGF signals guide migration of mesenchymal cells, control skeletal morphogenesis [corrected] and regulate gastrulation during sea urchin development. 2008, Pubmed , Echinobase
Ruta, Brief review of the stylophoran debate. 1999, Pubmed
Sanderson, r8s: inferring absolute rates of molecular evolution and divergence times in the absence of a molecular clock. 2003, Pubmed
Satoh, How was the notochord born? 2012, Pubmed
Satoh, On a possible evolutionary link of the stomochord of hemichordates to pharyngeal organs of chordates. 2014, Pubmed
Shukeir, Prostate secretory protein PSP-94 decreases tumor growth and hypercalcemia of malignancy in a syngenic in vivo model of prostate cancer. 2003, Pubmed
Simakov, Hemichordate genomes and deuterostome origins. 2015, Pubmed , Echinobase
Slater, Automated generation of heuristics for biological sequence comparison. 2005, Pubmed
Sodergren, The genome of the sea urchin Strongylocentrotus purpuratus. 2006, Pubmed , Echinobase
Stamatakis, RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. 2014, Pubmed
Sun, Signal-dependent regulation of the sea urchin skeletogenic gene regulatory network. 2014, Pubmed , Echinobase
Sun, Large scale gene expression profiling during intestine and body wall regeneration in the sea cucumber Apostichopus japonicus. 2011, Pubmed , Echinobase
Sun, RNA-Seq reveals dynamic changes of gene expression in key stages of intestine regeneration in the sea cucumber Apostichopus japonicus. [corrected]. 2013, Pubmed , Echinobase
Sutcliffe, MSMB variation and prostate cancer risk: clues towards a possible fungal etiology. 2014, Pubmed
Suzek, UniRef clusters: a comprehensive and scalable alternative for improving sequence similarity searches. 2015, Pubmed
Suzuki, Antithrombotic and anticoagulant activity of depolymerized fragment of the glycosaminoglycan extracted from Stichopus japonicus Selenka. 1991, Pubmed , Echinobase
Suzuki, Genomic approaches reveal unexpected genetic divergence within Ciona intestinalis. 2005, Pubmed
Swalla, Deciphering deuterostome phylogeny: molecular, morphological and palaeontological perspectives. 2008, Pubmed , Echinobase
Takahashi, Induction of pluripotent stem cells from mouse embryonic and adult fibroblast cultures by defined factors. 2006, Pubmed
Tamura, MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. 2007, Pubmed
Tanaka, A target of thrombin activation promotes cell cycle re-entry by urodele muscle cells. 1998, Pubmed
Tarailo-Graovac, Using RepeatMasker to identify repetitive elements in genomic sequences. 2009, Pubmed
Telford, Phylogenomic analysis of echinoderm class relationships supports Asterozoa. 2014, Pubmed , Echinobase
Trapnell, TopHat: discovering splice junctions with RNA-Seq. 2009, Pubmed
Trapnell, Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. 2012, Pubmed
Wang, The draft genome of the grass carp (Ctenopharyngodon idellus) provides insights into its evolution and vegetarian adaptation. 2015, Pubmed
Yan, Cloning and characterization of a mouse liver-specific gene mfrep-1, up-regulated in liver regeneration. 2002, Pubmed
Yu, PTGBase: an integrated database to study tandem duplicated genes in plants. 2015, Pubmed
Zhang, The oyster genome reveals stress adaptation and complexity of shell formation. 2012, Pubmed