Next Article in Journal
Higher Mutation Burden in High Proliferation Compartments of Heterogeneous Melanoma Tumors
Next Article in Special Issue
Bioinformatics Applications to Reveal Molecular Mechanisms of Gene Expression Regulation in Model Organisms
Previous Article in Journal
Genetic Screen in Adult Drosophila Reveals That dCBP Depletion in Glial Cells Mitigates Huntington Disease Pathology through a Foxo-Dependent Pathway
Previous Article in Special Issue
Somatic Copy Number Alterations and Associated Genes in Clear-Cell Renal-Cell Carcinoma in Brazilian Patients
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 22
Issue 8
10.3390/ijms22083882
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparative Transcriptomes of the Body Wall of Wild and Farmed Sea Cucumber Isostichopus badionotus

by
Roberto Martín-Hernández
1,†,
Rossanna Rodríguez-Canul
2,
Nuvia Kantún-Moreno
2,
Miguel A. Olvera-Novoa
3,
Oscar Medina-Contreras
4,
Cristobal Garikoitz-Legarda
5,
Juan Carlos Triviño
5,
Jesús Alejandro Zamora-Briseño
2,
Víctor May-Solis
3,
Alicia Poot-Salazar
6,
Juan Antonio Pérez-Vega
2,
Judit Gil-Zamorano
7,
George Grant
8,*,
Alberto Dávalos
7,* and
Leticia Olivera-Castillo
3,*,†
1
Bioinformatics and Biostatistics Unit, IMDEA Food Institute, CEI UAM+CSIC, Carretera de Cantoblanco 8, 28049 Madrid, Spain
2
Laboratorio de Inmunología y Biología Molecular, Centro de Investigación y de Estudios Avanzados del IPN-Unidad Mérida, Antigua Carretera a Progreso Km. 6, Mérida 97310, Yucatán, Mexico
3
Laboratorio de Nutrición Acuícola, Centro de Investigación y de Estudios Avanzados del IPN-Unidad Mérida, Antigua Carretera a Progreso Km. 6, Mérida 97310, Yucatán, Mexico
4
Unidad de Investigación Epidemiológica en Endocrinología y Nutrición, Hospital Infantil de México “Federico Gómez”, Mexico City 06720, Mexico
5
Bioinformatics Department, Sistemas Genómicos S.L., Ronda de Guglielmo Marconi 6, 46980 Paterna, Spain
6
Centro Regional de Investigaciones Acuícola y Pesqueras en Yucalpetén, Instituto Nacional de Pesca y Acuacultura, Boulevard del Pescador S/N, Puerto de Abrigo, Progreso 97320, Yucatán, Mexico
7
Laboratory of Epigenetics of Lipid Metabolism, IMDEA Food Institute, CEI UAM+CSIC, Carretera de Cantoblanco 8, 28049 Madrid, Spain
8
School of Medicine, Medical Sciences and Nutrition, University of Aberdeen, Aberdeen AB25 2ZD, UK
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this study.
Int. J. Mol. Sci. 2021, 22(8), 3882; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22083882
Submission received: 18 February 2021 / Revised: 27 March 2021 / Accepted: 6 April 2021 / Published: 9 April 2021
(This article belongs to the Special Issue Molecular Mechanisms of Gene Expression: “Bioinformatics of Gene Regulations and Structure”)
Browse Figures
Versions Notes

Abstract

:
Overfishing of sea cucumber Isostichopus badionotus from Yucatan has led to a major population decline. They are being captured as an alternative to traditional species despite a paucity of information about their health-promoting properties. The transcriptome of the body wall of wild and farmed I. badionotus has now been studied for the first time by an RNA-Seq approach. The functional profile of wild I. badionotus was comparable with data in the literature for other regularly captured species. In contrast, the metabolism of first generation farmed I. badionotus was impaired. This had multiple possible causes including a sub-optimal growth environment and impaired nutrient utilization. Several key metabolic pathways that are important in effective handling and accretion of nutrients and energy, or clearance of harmful cellular metabolites, were disrupted or dysregulated. For instance, collagen mRNAs were greatly reduced and deposition of collagen proteins impaired. Wild I. badionotus is, therefore, a suitable alternative to other widely used species but, at present, the potential of farmed I. badionotus is unclear. The environmental or nutritional factors responsible for their impaired function in culture remain unknown, but the present data gives useful pointers to the underlying problems associated with their aquaculture.

1. Introduction

Sea cucumbers (Phylum Echinodermata, Class Holothuroidea) are a diverse group of benthic marine invertebrates that inhabit intertidal to abyssal zones in temperate and tropical oceans worldwide [1]. There is a high demand for them in Asian markets, where sea cucumber are widely used as a luxury seafood and in traditional medicines. They are also a rich source of potential health-promoting factors such as fucosylated chondroitin sulphates, collagens, peptides, and sphingolipids, among others [2,3]. Overall, this has led to over-fishing of wild stocks and depletion of some species [4], including Isostichopus badionotus from the Caribbean Sea and Gulf of Mexico. In fact, capture of I. badionotus from off the coast of the Yucatan State has been so intense in recent years that federal regulatory authorities have implemented protective measures to prevent the regional extinction of the species [5,6]. This was vital not only to preserve stocks but also to protect the local marine environment since sea cucumbers recycle nutrients and break down detritus and other organic matter on the sea bottom and have a critical role in maintaining the marine ecosystems [4,7].
Farming of I. badionotus could in part alleviate this problem by offering a readily available commercial source of mature sea cucumbers or alternatively by providing seed animals for restocking and regeneration of wild populations. Cultivation studies are ongoing [8,9,10]. In these, animals are kept in indoor fiberglass tanks in a closed recirculating system with a temperature maintained at a range between 24–26 °C and a 13:11 (light:dark) photoperiod [8]. The main feed offered is macroalgae-based, with most diets including a mixture of Ulva sp., Sargassum sp. and Macrocystis sp., with fatty acid-rich supplements [8]. Under these conditions most I. badionotus survive, go through larval development, and slowly grow and mature as juveniles and adults [8,9,10].
Much of the therapeutic and commercial potential of sea cucumber lies in the range of health-modulating factors that are expressed by it [2,3]. However, the profile, yield, and reactivity of constituent bioactive factors are dependent on the exact metabolic state of the sea cucumber and individual compounds may be reduced, lost, or produced in inactive forms if host metabolism is disturbed. While I. badionotus has been shown to grow in culture, little is known about its metabolic state during cultivation and hence its capacity to produce functional health-promoting components.
Omics approaches, such as high-throughput RNA sequencing and subsequent bioinformatics analyses, have provided useful insights into the biology of sea cucumber species. For example, the de novo transcriptomic assembly of A. japonicus has identified key genes associated with important traits such as disease response, color variation, growth, and development [11,12,13,14].
In the present study, the transcriptome of body wall tissue from I. badionotus from Yucatan has been described for the first time. In addition, a comparison between the transcriptomes of wild (captured) and farmed age-equivalent I. badionotus has been done to assess the similarity/dissimilarity in their metabolic profiles.

2. Results

2.1. Samples

Adult Isostichopus badionotus (Holothuria) were collected off the coast of the state of Yucatan, Mexico. The organisms were placed in individual plastic bags containing sea water from the collection site and kept at 22 to 24 °C, the average water temperature at the collection site during transport to the laboratory. The individual lengths and body weights were measured (Supplementary Table S1) and four representative animals (14–25 cm, 160–420 g), whose ages were estimated at approximately two years by the procedures of Poot-Salazar et al., 2014 [15], were selected for transcriptome analysis.
I. badionotus brood-stock had been collected previously from the same location and transported to the Sea Cucumber Aquaculture Laboratory, Telchac Puerto Marine Station, CINVESTAV. The animals were kept in indoor fiberglass tanks at 24–26 °C under a photoperiod of 13:11 (light:dark), to simulate summer conditions and fed algae-based meals in accordance with recommended procedures [8].
Brood-stock maintained under these culture conditions spawned spontaneously during the normal reproductive season of July to November without induction [9]. When spawning occurred, gametes were collected for artificial fertilization and incubated in 250 L fiberglass tanks. Larval feeding was done with live microalgae produced in situ, and commercial microalgae concentrates [8,9].
Early juveniles were transferred to fresh fiberglass tanks for subsequent growth. They were fed a mix of algae supplemented with periodic additions of live and concentrated diatoms, following published protocols [10]. Four representative first generation individuals of approximately two years of age were taken for transcriptome analysis (5–7 cm, 17–43 g) (Supplementary Table S1).

2.2. Transcriptome Sequencing

Eight cDNA libraries were generated from the I. badionotus body wall mRNA for comparative transcriptomic analysis between the two sea cucumber groups: wild adult (~two-year-old) specimens (4 cDNA libraries) and two-year-old farmed animals born and reared in captivity (n = 4).
Total RNA was extracted from individual body wall samples using TRIzol Reagent. Preparations with an RNA integrity number (RIN) >7.0 were taken forward for cDNA library preparation. mRNA isolation, cDNA synthesis, and library preparation were done using the TruSeq Stranded mRNA Sample Preparation Kit. The quality and size distribution of the cDNA libraries was validated using a TapeStation Genomic DNA system.
The libraries were processed with Illumina deep-sequencing. After paired-end sequencing, a total of 320,540,514 raw reads with a 76 bp average length were acquired, with four libraries corresponding to wild adult specimens: 74,199,152 (Wild1), 83,116,920 (Wild2), 77,041,560 (Wild3) and 86,182,882 (Wild4). These were used for de novo transcriptome assembly. Raw data quality is described in Supplementary Table S2. The raw reads were deposited in the Sequence Read Archive (SRA) in the GenBank database, linked to BioProject PRJNA639785 (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/bioproject/PRJNA639785 accessed on 16 June 2020), which also includes the sequences resulting from transcriptome assembly. A total of 293,055,064 raw reads with a 150 bp average length were obtained from farmed individuals: 39,511,252. The raw reads from the farmed individuals were used to evaluate differential gene expression and are available in the GEO database with ID GSE157183 (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/geo/query/acc.cgi?acc=GSE157183 accessed on 1 September 2020).

2.3. De Novo Transcriptome Assembly

Processed paired-end reads obtained from the wild sea cucumber specimens were used for de novo assembly. As low k-mer values are prone to introduce false contigs, and it is recommended to use a maximum k-mer size corresponding to 2/3 of the read length to allow enough overlapping, a range of k-mers between 23 and 51 was tested in order to capture short and longer transcripts with sufficient coverage (Figure 1). Following selection of the best assemblies based on N50 values (Supplementary Table S3), and merging and filtering contigs shorter than 100 bp, the merged assembly produced with Velvet included 148,431 contigs. Contig size ranged from 101 to 6957 bp, with a 384.9 bp average length (Supplementary Figure S1A). The relationship between contig size and contig coverage exhibited a positive trend, with a 5.52 mean coverage value (Supplementary Figure S1B).
Oases was used for the scaffolding step. Therefore, after filtering redundant sequences exploiting pairing information of reads to further connect contigs, and producing transcript isoforms, the final assembly was reduced to 132,257 consensus sequences, representing a total of 81,784 unique loci or unigenes (Supplementary Table S4). Average unigene length was 679 bp, with a maximum of 10,560 bp and a N50 value of 1219.
Quantitative measurement of transcriptome assembly process completeness was done with the BUSCO software, which is based on sets of Benchmarking Universal Single-Copy Orthologs from OrthoDB [16]. The utilized lineage-specific library corresponds to Eukaryota. We were able to completely recover a total of 278 BUSCO genes, thus achieving a 65% of completeness. In addition, this tool was used to compare our de novo transcriptome assembly with two other sea cucumber species belonging to Apostichopus japonicus transcriptome shotgun assemblies (Acc. IDs: GFKU00000000.1 and GHCH00000000.1) available online on Genbank, part of the NCBI Nucleotide database.
Two of the GenBank assemblies slightly outperformed our assembly in terms of complete identified BUSCOs (Figure 2), but neither was able to recover the full set of 429 eukaryote orthologs. Our assembly also exhibits a significant reduction in the number of duplicated BUSCOs, probably due to removal of duplicated reads during the de novo assembly process. Further differences may be due to the tool chosen for de novo assembly, since the two GenBank assemblies were produced with the Trinity de novo assembly tool [17].

2.4. Transcriptome Annotation

BLAST analyses of the 132,257 non-redundant assembled transcripts were run against different databases. A total of 5275 assembled transcripts had hits with 3599 unique sequences from the NT database, while 21,783 had hits with 12,609 unique protein sequences from the Uniprot database. GO functional annotation resulted in annotation of 11,425 unique transcripts with 7390 highly specific GO terms and KEGG enzyme IDs (Figure 3). Comparing these results to those reported by a previous transcriptome assembly of a closely related species [14], the amount of functionally annotated genes against GO database is very similar. However, the above-mentioned study reported higher annotation rates against NT and Uniprot databases. This fact might be explained by the low level of gene duplication and higher gene missingness of the transcriptome assembly generated in this research (Figure 2), as well as by the use of a different assembler.
RNAcentral analysis produced 218 potential RNAs, of which 84 corresponded to ribosomal RNAs, 80 to long non-coding RNAs, 34 to miscellaneous RNAs, 8 to transfer RNAs and 4 to antisense RNAs. A small fraction was identified as small nucleolar RNAs (n = 3). One RNA was a precursor and another a bacterial transfer-messenger RNA, while three sequences remained as unclassified RNAs (data not shown)

2.5. Gene Ontology (GO) Classification

Assigned GO terms for the wild sea cucumber assembled transcriptome were distributed among highly specific biological processes (3071), molecular functions (2077), and cellular components (736) (Supplementary Table S5). Among the biological processes, the top three most frequent terms were GO:1901671 (“positive regulation of superoxide dismutase activity”) (149), GO:1904247 (“positive regulation of polynucleotide adenylyltransferase activity”) (94) and GO:0043547 (“positive regulation of GTPase activity”) (59). For molecular functions, GO:0015020 (“glucuronosyltransferase activity”) (57), GO:0003887 (“DNA-directed DNA polymerase activity”) (51) and GO:0004722 (“protein serine/threonine phosphatase activity”) (41) were the most frequent terms. For cellular components, GO:0030964 (“NADH dehydrogenase complex”) (12), GO:1990316 (“ATG1/ULK1 kinase complex”) (10) and GO:0045275 (“respiratory chain complex III”) (9) were the most frequent (Figure 4).

2.6. Differentially Expressed Genes between Wild and Farmed Sea Cucumbers

Sequenced reads from each individual were mapped to the reference transcriptome to quantify gene expression of body wall samples in both groups (Table 1). In addition, reads obtained from wild specimens were mapped against the two previously tested transcriptome assemblies (GFKU00000000.1 and GHCH00000000.1), and against a publicly available Apostichopus japonicus genome assembly (ASM275485v1), thus obtaining spurious overall alignment rates (Supplemental Table S6). Results from this comparative analysis highlight a significant genomic difference between sea cucumbers from the Apostichopus and Isostichopus genus. Furthermore, a sample (W2) showing the higher alignment rate on ASM275485v1 genome assembly was selected to evaluate the most highly expressed genes which were detected (Supplemental Table S7).
In transcript counting per analyzed sample, only reads showing a unique match to the transcriptome reference were counted as valid alignments (MAPQ = 60). For differential expression assessment, only transcripts detected in both groups were considered for differential expression.
Of the 132,257 newly assembled consensus transcript sequences from wild I. badionotus, consistent expression of 18,449 transcripts was detected across all the analyzed samples from both groups. From these, a total of 3760 transcripts (20.4%) were found to be differentially expressed between the wild and farmed groups (Figure 5a); the False Discovery Rate threshold was significant (FDR < 0.05) and a log2 fold change threshold was used. Of these transcripts, only 303 were known sequences showing a hit in the NT database (Supplementary Table S8). Gene names for each annotated transcript were retrieved when available. Principal components analysis (PCA) clearly separated the farmed (reared in culture or maintained in culture after capture) from the wild animals (Figure 5b). The contribution of each variable and the samples to the PCA are shown in Supplementary Figure S1.
The identified differentially expressed genes (DEGs) were mapped against the KEGG KO (KEGG Orthology) database, which contains molecular functions represented in terms of functional orthologs. The identified genes were mapped to 83 pathways (table KEGG_ko mapping, Supplementary Table S9), which are implicated in organism development and growth.
Several annotated genes involved in host responses to metabolic stress, including COXI, EIF2AK3, ATPA, MAPK14-1, HSP40, hsc70, and TRPM2, were upregulated in farmed I. badionotus (Table 2). The EIF2AK3 gene mapped to the Mitophagy pathway (ko04137).
Of the annotated downregulated transcripts, Techylectin 5A, CREB1 and collagen mRNAs, were of particular relevance (Table 2). Sea cucumber lectins are involved in host immune responses to infection, inflammation, cell-cell or cell-extracellular matrix interactions and self- recognition [18,19]. If lectin expression was low in the farmed animals, they were likely to be prone to immunological compromise and disease. CREB1 is a gluconeogenesis regulatory factor and linked to the Longevity regulating pathway. Collagen is a major component of the body wall of sea cucumber and any reduction in its synthesis, as suggested by the under-expression of collagen mRNAs in farmed I badionotus, is likely to greatly impair their health and well-being.
Preliminary collagen extraction studies done on body wall from wild and farmed I. badionotus do indeed indicate that the constituent collagens differ significantly (Olivera-Castillo, unpublished data). Collagen from farmed animals was more fibrous and difficult to extract by standard procedures [20] than that from wild counterparts and was of reduced proportion of body wall mass.

2.7. Functional Enrichment between Wild and Farmed Sea Cucumbers

Enrichment analysis allows identification of overrepresented functional categories among differentially expressed transcripts. Some of the most important biological processes and molecular functions induced by captivity (Figure 6) were ATP synthesis (GO:0015986), proton-transporting ATP (GO:0046933), proton transport (GO:0015992) and the oxidoreductase process (GO:0016705). These functions are associated with energy synthesis for organism growth and proliferation. Other interesting significant functions were gluconeogenesis (GO:0006094), cytoskeleton (GO:0005200), and the vesicle process (GO:0030127). Other biological processes related to protein synthesis included protein folding (GO:0006457), protein glycosylation (GO:0006486), protein N-linked glycosylation via asparagine (GO:0018279), and unfolded protein binding (GO:0051082).

2.8. RT-qPCR Validation of Differentially Expressed Genes Induced by Captivity

We validated the DEG results produced in the wild vs. farmed (reared in culture) sea cucumber comparisons by sampling five sea cucumbers from the coast of the Yucatan state. We randomly selected eleven genes (oligo sequences in Supplementary Table S10) from the 303 transcripts (NT database hits) for validation analysis by RT-qPCR. In most cases, the RT-qPCR results exhibit behavior similar to the RNA-seq analysis (Figure 7), suggesting that our de novo assembly of the I. badionotus transcriptome is reliable.

3. Discussion

Sea cucumber (Isostichopus badionotus) is vital for maintaining the marine microenvironment and water quality due to its ability to scavenge calcareous algae, foraminifera, molluscs, and clear sedimentary detritus from the seabed [8,20]. However, a rapid decline in I. badionotus stocks from Yucatan’s coast due to overfishing has seriously compromised this ecological protection [4,5]. This exploitation of I. badionotus is due to its use as a substitute luxury seafood and as a potential source of health-promoting factors to replace traditional sea cucumber species whose wild stocks are already depleted [2,3,5,6]. Despite the high usage of I. badionotus, little is known about its general metabolism or indeed, its comparability to other sea cucumber species. To address this in part, the functional transcriptome of the body wall of wild-type I. badionotus from Yucatan has herein been described for the first time and compared with data available in the literature for other sea cucumber species.
The transcriptome of wild-type I. badionotus yielded 132,257 consensus sequences, representing a total of 81,784 unique loci or unigenes. From these, 5275 assembled transcripts had hits with 3599 unique sequences from the NT database, while 21,783 had hits with 12,609 unique protein sequences from the Uniprot database. GO functional annotation identified 11,425 unique transcripts with 7390 highly specific GO terms and KEGG enzyme IDs. Amongst the biological processes, molecular functions, and cellular components identified those linked to cellular respiratory chain complex activity, positive regulation of cell growth and protection from generated reactive oxygen species were prominent. The transcriptome of wild I. badionotus had a high degree of similarity with that for Apostichopus japonicus [12,13,14,21,22,23,24], and Stichopus horrens [25] and the proteome reported for Stichopus japonicus [26]. The immune molecule progenitors identified in wild I. badionotus also closely matched those reported for other sea cucumber species [27,28].
Overall, the functional metabolic profile of wild I. badionotus from Yucatan is comparable with that of A. japonicus and other captured species. Recent studies have also shown that I. badionotus contains bioactive factors, such as glycosaminoglycans, peptides, saponins, and lectins, which have health-modulating properties comparable to those products from other sea cucumber species [2,3,29,30]. Therefore, I. badionotus from Yucatan is a proper alternative to traditionally used sea cucumber species. However, its ongoing capture from the wild is non- sustainable. To supply the commercial demands for this species a viable alternative to fishing, such as aquaculture, is needed.
Studies on farming of I. badionotus are ongoing [9,10]. However, as with many sea cucumber species [8,20,31,32], rearing I. badionotus in culture has proved problematic with development and growth being variable and slow. Here we now show that the active transcriptome of farmed I. badionotus was much altered compared to its wild counterpart. Of 18,449 transcripts common to both sets of animals, twenty-one percent were differentially expressed in farmed I. badionotus (4% downregulated by 2.5 to 8.0-fold and 8% by 1.5 to 2.5-fold; 7% upregulated by 1.5 to 2.5-fold and 2% by 2.5 to 9.0-fold).
The farmed I. badionotus studied here were first generation offspring derived from I. badionotus brood-stock collected from the same location off the coast of the state of Yucatan as the wild animals that were tested. The differential changes in their transcriptome were, therefore, due mainly to adaptational changes to the laboratory environment in which they were raised rather than longer-term effects of domestication or selection. Indeed, the findings of a recent study [33] indicate that these growth-limiting changes may be reversed if farmed juveniles are transferred to sea shelters close to the area of parental capture. In this case, growth of the juveniles was greatly increased compared to that laboratory environment [33]. This would be consistent with a previous report of little molecular variance between sea-based wild and cultured A. japonicus populations [34].
The present study’s data demonstrates that multiple aspects of normal host metabolism were disrupted in farmed I. badionotus. In the wild, the local environment, including season, pressure, temperature, light, salinity, and nutrient and micronutrient availability, greatly influence sea cucumber physiology and development [4,12,13,14,35,36,37]. While farmed I. badionotus were reared under conditions that approximated those off the Yucatan coast and received appropriate nutrients, our present findings indicate that some of these environmental or nutritional parameters were sub-optimal. Indeed, some of the differential changes observed could be associated with low pressure, intense lighting, inappropriate habitat, stress, microbiota, as well as nutrient or micronutrient deficiencies in the laboratory setting [12,38,39,40,41].
Major deregulated processes included the Mitophagy pathway, the Longevity regulating pathway and ATP synthesis coupled proton transport. Of the proportion of the transcripts that were annotated, the most significant changes were in collagen mRNAs (−4–5 to −5.8-fold) and cytochrome c oxidase subunit I (COI) gene (+5.5 to +9.0-fold).
The EIF2AK3 gene of the Mitophagy pathway was overexpressed in farmed I. badionotus. Mitophagy is the process by which damaged or stressed mitochondria are degraded and cleared by autophagy. EIF2AK3 is linked with the inactivation of the eukaryotic translation-initiation factor’s alpha subunit, reduced translational initiation, and repression of global protein synthesis [42]. In contrast to EIF2AK3, the gluconeogenesis regulatory gene CREB1 [43] of the Longevity regulating pathway was downregulated in the farmed animals.
Mitochondrial function and viability may be significantly impaired in farmed I. badionotus. In addition to EIF2AK3, cytochrome c oxidase subunit 1 (COX1) and ATP synthase subunit alpha (ATPA), which are involved in the mitochondrial inner membrane electron transport chain were also upregulated. This change may be associated with the accumulation of reactive oxygen species, cytoplasmic acidification, caspase activation, and cell death [44].
Collagen mRNAs were considerably downregulated in farmed I. badionotus. Since collagen is the major component (~70%) of the body wall of sea cucumber [45], this finding indicates that collagen production was significantly impaired or suppressed in these cultured animals. Preliminary extraction studies tend to confirm this in that the collagen of farmed I. badionotus was more fibrous, less readily extractable and a lower proportion body wall mass than that from wild counterparts.
The diet offered to farmed I. badionotus was formulated per present knowledge of their requirements [8,9,10]. However, it cannot be excluded that important, yet unidentified nutrients or micronutrients were missing or deficient in the food. Alternatively, the digestion, absorption, and systemic processing/utilization of nutrients and micronutrients, including important organic and inorganic compounds, may be impaired due to modifications to digestive enzyme activity in the intestine and the constituent bacterial and fungal populations resulting from the culturing conditions [46,47,48,49,50,51].
Alternatively, the capacity of farmed I. badionotus to deal with and clear toxins and other deleterious cellular factors and metabolites may be compromised. The increased activity of ATP-dependent glucuronosyltransferases, which are responsible for the formation of glucuronides with a large variety of cytotoxic and genotoxic compounds and aid in their clearance [52], and increased expression of EIF2AK3 gene, which encodes a protein with critical roles in detection and initiation of cellular responses to endoplasmic reticulum stress [26], may be indicators of such disruption to systemic metabolism. The farmed animal may thus need to restrict metabolic activity to survive and remain healthy [31,53].
The transcriptome data presented here show the functional metabolism in wild I. badionotus from Yucatan was like that of A. japonicus and other widely captured sea cucumber species. However, when I. badionotus was raised and maintained in culture its metabolism was significantly disrupted and production of at least one major body component impaired. The exact factors responsible for these adverse changes in farmed I. badionotus and their modes of action need further detailed study particularly considering the urgent need to find a sustainable alternative to capture from the wild to supply the commercial market.

4. Materials and Methods

4.1. Collection of Wild Sea Cucumbers

Adult Isostichopus badionotus (Holothuria) were collected by coastal scuba diving off the state of Yucatan, México during June–July 2019. Collection was authorized via a permit (DGOPA/1009/210809/08761) issued by the National Aquaculture and Fishing Commission (Comisión Nacional de Acuacultura y Pesca).
Immediately after removal from the sea floor, the organisms were individually placed in plastic bags containing sea water from the collection site. At the surface, the bags were submerged in cold water inside coolers and kept at 22 to 24 °C, the water temperature range at the collection sites. This temperature control prevented proteolysis or autolysis. They were transported live to the laboratory. Before processing, individual length was measured with a plastic ruler in centimeter units and body weight measured with a spring scale in grams to 0.1 accuracy. This group was called “wild” sea cucumbers (14–25 cm, 160–420 g; Supplementary Table S1).

4.2. Cultivation of Sea Cucumbers

Isostichopus badionotus brood-stock had been collected from the same location and transported to the Sea Cucumber Aquaculture Laboratory, Telchac Puerto Marine Station, Center for Research and Advanced Studies (Centro de Investigación y de Estudios Avanzados—CINVESTAV). In the laboratory, the animals were kept in indoor fiberglass tanks forming part of a closed recirculating system including mechanical and biological filters, skimmer, and UV water-sterilization. Water temperature was maintained within a 24–26 °C range using chillers and air-conditioning. Laboratory photoperiod was a constant 13:11 (light:dark) to simulate summer conditions. The organisms were fed daily using a formulated mixture of algae meals (1:1:1; Sargassum sp.:Ulva sp.:Macrocystis piryfera) [8] mixed with washed and chlorinated beach sand. Before feeding, feces were removed by siphoning.
Brood-stock maintained under these culture conditions spawned spontaneously during the normal reproductive season of July to November without induction [9]. When spawning occurred, gametes were collected for artificial fertilization and incubated in 250 L fiberglass tanks. Larval feeding was done following Zacarías-Soto et al. [9], using live microalgae produced in situ and commercial microalgae concentrates (Instant Algae®, Reed Mariculture, Campbell, CA, USA).
Early juveniles were transferred to fresh fiberglass tanks (see above) for subsequent growth. They were fed a mix of algae enriched with Algamac® (Aquafauna Bio-Marine, Hawthorne, CA, USA), and supplemented with periodic additions of live and concentrated diatoms, following published protocols [8,9,10]. First generation individuals of approximately two years of age were collected for the transcriptome analysis; these were called “farmed” sea cucumbers (5–7 cm, 17–43 g) (Supplementary Table S1).

4.3. RNA Isolation, cDNA Library, and Illumina Sequencing

Both the wild and farmed sea cucumbers were processed by first washing with cold sterile distilled water to remove excess sea salt followed by immediate and rapid evisceration. Each body wall was cut into several pieces (approx. 1 cm2), and each piece immersed in 1 mL RNAlater™ stabilization solution (Ambion©, (Thermo Fisher Scientific), Mexico City, MX)©, Cat. AM7021), snap frozen in liquid N2 and stored at −80 °C until use. Total RNA was extracted using TRIzol Reagent® (Ambion©, Cat. 15596-026) essentially in accordance with Puch-Hau et al., 2019 [54]. Quantification of total RNA was done with a Qubit 2.0 (Thermo Fisher Scientific, Mexico City, MX, Mexico), and its quality [RNA integrity number (RIN)] assessed using an Agilent 2100 Bioanalyzer® with the RNA 6000 Nano Lab Chip Kit® (Agilent Technologies©, Santa Clara, CA, USA). Total RNA with a RIN score >7.0 was accepted for cDNA library preparation. Isolation of mRNA, cDNA synthesis, and library preparation were done using the TruSeq Stranded mRNA Sample Preparation Kit® (Illumina©, San Diego, CA, USA). Briefly, the poly(A)+ mRNA fraction was isolated from total RNA, using poly(T) oligo attached to magnetic beads. Following purification, the RNA was fragmented and primed with random hexamers, and reverse transcribed using SuperScript II Reverse Transcriptase (Invitrogen©, Thermo Fisher Scientific, Mexico City, MX) with addition of actinomycin D. Second-strand synthesis was done using polymerase I and RNaseH, with replacement of dTTP for dUTP to generate the double-stranded cDNA (ds-cDNA). The ds-cDNA fragments were 3′-end adenylated and ligated to Illumina paired-end sequencing adapters. Finally, the products were purified and enriched with PCR to create the final indexed double-stranded cDNA library. The quality and size distribution of the cDNA libraries was validated using a TapeStation Genomic DNA system® (Agilent©). Library construction and sequencing were done at the National Institute of Genomic Medicine (Instituto Nacional de Medicina Genómica—INMEGEN, México City, MX). Paired-end sequencing was conducted on an Illumina NextSeq 500 platform to generate 2 × 76 or 2 × 150 paired-end (PE) reads.

4.4. Short Reads De Novo Assembly

Quality of the raw reads was checked using FASTQC tools [55] to identify any deviations or problems in the sequencing process. The resulting raw reads were pre-processed by consistently trimming Illumina sequencing adaptors to discard low quality read pairs using a Q24 threshold Trim Galore v.0.5.0 tool. Duplicated reads were removed to reduce redundancy and improve assembly process computational efficiency, using the same software.
The Velvet v.1.2.1 sequence assembler tool [56] was used for de novo transcriptome assembly with a multi-k assembly pipeline. Length of k-mers ranged from 51 to 23, and the corresponding estimated insert length for each sequenced library was used as a parameter. A total of nine assemblies, obtained with higher and lower k-mer lengths, were selected based on the higher N50 values obtained, and then merged using Oases v.0.2.08 [57] to obtain the final assembly. Finally, the CD-HIT v.4.7 tool [58] was used to remove redundant transcripts, with a 90% similarity threshold.
Quantitative measures of the de novo assembled transcriptome were generated by comparing the generated assembly to previous assemblies from closely related species using the BUSCO software v.4.0.2 [59].

4.5. Functional Annotation

The de novo assembled transcriptome was annotated using the Basic Local Alignment Search Tool (BLAST) [60] locally, with parameters such as max target seqs set to 1, best hit overhang set to 0.1 and best hit score edge set to 0.1. Queried sequence databases were the NCBI nucleotide (NT), Uniprot and RNACentral databases [61,62]. BLAST analyses results were filtered against nucleotide and protein databases, and only hits with an e-value score equal to or lower than 0.001 retained. For sequence annotation with GO terms the annotated sequences against the Uniprot database were introduced into the Blast2Go v.5.1 tool [63]. Default parameters were used which allow a compromise between coverage and correctness of the assigned functional GO term. The predicted proteins were functionally annotated according to three different functional categories available in the GO database: biological process, molecular function, and cellular component.

4.6. Functional Enrichment

A functional enrichment analysis was run using the GO functional information, based on hypergeometric distribution and using the R ver. 3.2.3 statistical software. Following a Fisher’s exact test, the false discovery rate (FDR) adjusted p-value was set at a 0.05 cut-off.

4.7. Differential Expression between Sea Cucumber Body Wall Transcriptomes, and Pathway Analysis

The newly assembled transcriptome from wild I. badionotus was used as reference for sequence alignment. The sequenced paired-end reads from each sample of both groups (wild and farmed) were mapped against the reference transcriptome using the Hisat2 software [64] setting rna-strandness parameter to FR. Only reads showing a unique match were counted as valid alignments. The resulting count table was analyzed using the edgeR (v 3.20.9) [65] library in the R environment. After calculating normalization factors to scale the raw library sizes, exact tests were applied to evaluate differences in the mean counts for each transcript between both groups, and setting the False Discovery rate (FDR) statistical significance level below 0.05. Transcripts exhibiting a minimum expression level value of 3 CPM in at least three samples in both groups were considered as consistently expressed and were used in the differential expression analysis. Differentially expressed genes (DEGs) were mapped to their corresponding pathways using the KO (KEGG Orthology) database, a database of molecular functions represented in terms of functional orthologs [66].

4.8. Gene Validation

Eleven DEGs from the comparative transcriptome analysis were used for validation. Eight wild and eight farmed animals were chosen for gene validation (see Supplementary Table S10 for the complete set of primers used herein). After weighing and measuring, the animals were washed, quickly eviscerated and the epidermis removed from the body wall. The body wall was snap frozen in liquid nitrogen and quickly macerated into powder. A ≈ 100 mg sample of this powder was placed in 1 mL Trizol® (Invitrogen©). Total RNA was isolated with RNA Direct-zol RNA (Thermo Scientific©). The RNA was analyzed by 1.2% agarose gel electrophoresis, quantified using a NanoDrop 2000c® (Thermo Scientific) and RNA quality analyzed with the Qubit system. Genomic DNA (gDNA) was removed using the TURBO DNA-free™ Kit (Ambion©, Cat. AM1907) per manufacturer instructions and confirmed by PCR analysis of the 16S rRNA gene. Synthesis of cDNA was done with the RevertAid H Minus First Strand cDNA kit® (Thermo Scientific©). The RT-qPCR was done with a Rotor Gene Q (2-plex) real-time PCR detector (Qiagen©, Hilden, Germany) using Quantinova SYBR Green PCR master® (Qiagen©). Four RT-qPCR reactions were run, including non-template reactions as negative controls. Thermo-cycling conditions were 5 min at 95 °C, followed by 40 cycles of 15 s at 94 °C, and 45 s at an annealing temperature (Ta), depending on each primer pair. Dissociation curves were generated to verify amplified product specificity. The expression levels of selected DEGs were normalized by comparison to MTB2 internal reference genes. Values for Ct were exported, the fold-change was estimated as relative gene expression levels using the ddCt method, transformed to log2 scale, and compared with transcriptome data [67,68].

4.9. Ethics Statement

The capture of wild organisms was done with the authorization, participation, and supervision of scientific personnel of the National Fisheries Institute, Yucalpeten, Yucatan, Mexico (INAPESCA, CRIP-Yucalpeten). Sea cucumber handling and euthanasia comply with the Mexican Official Standard for the Care and Use of Laboratory Animals (NOM-062-ZOO-1999) and with the rules of the Internal Committee for the Care and Use of Laboratory Animals of CINVESTAV.

5. Conclusions

De novo transcriptomes for the body wall from adult I. badionotus from Yucatan captured in the wild or grown in culture have been generated for the first time. The functional metabolic profile of wild I. badionotus was comparable with that in the literature for A. japonicus and other regularly captured species. In conjunction with recent findings that I. badionotus contains bioactive factors, such as glycosaminoglycans, peptides, etc., which have health-modulating properties, the present data indicates that I. badionotus from the wild can be used as a substitute for traditionally captured sea cucumber species.
The data for first generation farmed I. badionotus revealed that their functional transcriptome was greatly altered compared to that of their wild counterparts. This disruption to host metabolism had multiple possible causes. Despite the animals being reared under conditions approximating those off the Yucatan coast and having received appropriate nutrition, the data indicates that aspects of their growth environment and availability or utilization of nutrients and micronutrients were sub-optimal. An important indicator of the combined adverse effects was that several key metabolic pathways that are important in effective handling and accretion of nutrients and energy, or clearance of harmful cellular metabolites were disrupted or dysregulated. For example, collagen mRNAs were greatly reduced, and collagen protein deposition was impaired.
The nature of the environmental or nutritional factors responsible for the disruption of metabolism in farmed I. badionotus remains unknown, but the present data gives useful insights into the difficulties associated with aquaculture of this important species.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/ijms22083882/s1. Supplementary Figure S1: Contig and PCA analysis of differentially expressed genes between wild and farmed sea cucumber Isostichopus badionotus. (A) Contig size distribution generated by Velvet. (B) Log-log plots showing the relation of contig lengths to the number of reads assembled into each. (C) Correlation among samples. (D) Multidimensional scaling plot of the analyzed RNA-Seq samples. (E) Hierarchical cluster analysis of differentially expressed genes. Supplementary Table S1: Summary of weight, length, age and life status of each Isostichopus badionotus collected for RNAseq analysis, Supplementary Table S2: Summary of data quality of sequencing data used for de-novo assembly before and after pre-processing, Supplementary Table S3: Summary of obtained assemblies using different k-mer length. Based on the best N50 values, selected assemblies for merging are highlighted, Supplementary Table S4: Summary statistics of de novo assembly for Isostichopus badionotus transcriptome and its functional annotation, Supplementary Table S5: Gene Ontology (GO) specific term frequency, Supplementary Table S6: Overall alignment rates obtained after mapping wild I. badionotus samples against 3 different genomic references, Supplementary Table S7: Top 10 Genes detected after mapping Wild2 sample (W2), Supplementary Table S8: Differentially expressed transcripts showing a hit with NT database. Data includes gene symbols without determined orthologs (LOC) and locus tags obtained from sequencing projects. Supplementary Table S9: KEGG_ko mapping, Supplementary Table S10: Primers used in the expression analysis by RT-qPCR.

Author Contributions

Conceptualization, L.O.-C. and A.D.; methodology, L.O.-C. and N.K.-M.; validation, L.O.-C., O.M.-C. and V.M.-S.; formal analysis, N.K.-M., R.M.-H., J.C.T., C.G.-L., L.O.-C., V.M.-S., J.A.Z.-B. and J.A.P.-V.; resources, R.R.-C., L.O.-C., A.P.-S., M.A.O.-N. and G.G.; data curation, R.M.-H., J.G.-Z., J.C.T. and C.G.-L.; writing—original draft preparation, R.M.-H., A.D., G.G., L.O.-C., R.R.-C., A.P.-S., M.A.O.-N.; writing—review and editing, G.G. and L.O.-C.; supervision, L.O.-C. and A.D.; project administration, R.R.-C.; funding acquisition, L.O.-C. and R.R.-C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Consejo Nacional de Ciencia y Tecnología (CONACYT) through four grants: CONACyT Fondos Mixtos Yucatán (M0023-A1) “El pepino de mar como un alimento funcional: obtención de sus principios bióactivos, caracterización biológica, efectos sobre el metabolismo y sistema inmune utilizando un modelo murino”; CONACyT Ciencia Básica (2015, 013/221734CB-2: 221734) “Actividad antiinflamatoria y cicatrizante del pepino de mar (Isostichopus badionotus) en un modelo murino: caracterización de la actividad farmacológica y los mecanismos moleculares involucrados”; CONACyT Fondos Mixtos Yucatán (FOMIX YUC-2011-C09-169961) “Aportaciones al conocimiento para el desarrollo tecnológico del cultivo de pepino de mar en Yucatán”; CONACyT Problemas Nacionales 2015 (2015-01-1551) “Bases biológicas y de manejo para la integración multitrófica del pepino de mar Isostichopus badionotus y macroalgas en sistemas de maricultivo comercial de corvina roja (Sciaenops ocellatus)”.

Institutional Review Board Statement

The capture of wild I. badionotus was done with the authorization, participation, and supervision of the National Fisheries Institute, Yucalpeten, Yucatan, Mexico (INAPESCA, CRIP-Yucalpeten). Sea cucumber handling and euthanasia was in accordance the Mexican Official Standard for the Care and Use of Laboratory Animals (NOM-062-ZOO-1999) and the rules of the Internal Committee for the Care and Use of Laboratory Animals of CINVESTAV.

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw reads for wild animals were deposited in the Sequence Read Archive in the GenBank database. Data from corresponding samples are available in the following links: Wild1 (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/sra/SRX8584713), Wild2 (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/sra/SRX8584714, accessed on 9 April 2021) Wild3 (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/sra/SRX8584711, accessed on 9 April 2021), Wild4 (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/sra/SRX8584712, accessed on 9 April 2021). BioProject PRJNA639785 (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/bioproject/PRJNA639785, accessed on 17 February 2021). The raw reads from the farmed individuals were used to evaluate differential gene expression and are available in the GEO database with ID GSE157183 (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/geo/query/acc.cgi?acc=GSE157183, accessed on 9 April 2021).

Acknowledgments

The authors thank John Lindsay-Edwards for assistance with manuscript editing, proofreading, and formatting.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ATPAATP synthase subunit alpha
BLASTBasic Local Alignment Search Tool
COX1cytochrome c oxidase subunit 1
CREB1CAMP responsive element binding protein 1
DEGDifferentially expressed gene
EIF2AK3Eukaryotic Translation Initiation Factor 2 Alpha Kinase 3
FDRFalse discovery rate
GOGene Ontology
KOKEGG Orthology
RT-qPCRReverse transcription qualitative real-time polymerase chain rection
SRASequence Read Archive

References

  1. Pawson, D.L. Phylum Echinodermata. In Linnaeus Tercentenary: Progress in Invertebrate Taxonomy, Zootaxa 1668; Zhang, Z.-Q., Shear, W.A., Eds.; Magnolia Press: Auckland, New Zealand, 2007; pp. 649–764. [Google Scholar] [CrossRef]
  2. Khotimchenko, Y. Pharmacological potential of sea cucumbers. Int. J. Mol. Sci. 2018, 19, 1342. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Pangestuti, R.; Arifin, Z. Medicinal and health benefit effects of functional sea cucumbers. J. Tradit. Complement. Med. 2018, 8, 341–351. [Google Scholar] [CrossRef] [PubMed]
  4. Baker-Médard, M.; Ohl, K.N. Sea cucumber management strategies: Challenges and opportunities in a developing country context. Environ. Conserv. 2019, 46, 267–277. [Google Scholar] [CrossRef]
  5. Gamboa-Álvarez, M.Á.; López-Rocha, J.A.; Poot-López, G.R.; Aguilar-Perera, A.; Villegas-Hernández, H. Rise and decline of the sea cucumber fishery in Campeche Bank, Mexico. Ocean Coast. Manag. 2020, 184, 105011. [Google Scholar] [CrossRef]
  6. López-Rocha, J.A.; Velázquez-Abunader, I. Fast decline of the sea cucumber Isostichopus badionotus as a consequence of high exploitation in Yucatan, Mexico. Reg. Stud. Mar. Sci. 2019, 27, 100547. [Google Scholar] [CrossRef]
  7. Anderson, S.C.; Flemming, J.M.; Watson, R.; Lotze, H.K. Serial exploitation of global sea cucumber fisheries. Fish Fish. 2010, 12, 317–339. [Google Scholar] [CrossRef]
  8. Zacarías-Soto, M.; Tec-Tec, P.; Olvera-Novoa, M. Effect of diet on growth and body biochemical composition of juvenile four-sided sea cucumber Isostichopus badionotus (Selenka, 1867). Aquac. Res. 2018, 49, 939–946. [Google Scholar] [CrossRef]
  9. Zacarías-Soto, M.; Olvera-Novoa, M.A.; Pensamiento-Villarauz, S.; Sánchez-Tapia, I. Spawning and larval development of the four-sided sea cucumber, Isostichopus badionotus (Selenka 1867), under controlled conditions. J. World Aquac. Soc. 2013, 44, 694–705. [Google Scholar] [CrossRef]
  10. Sánchez-Tapia, I.A.; Slater, M.; Olvera-Novoa, M.A. Evaluation of the growth and survival rate of the Caribbean Sea cucumber, Isostichopus badionotus (Selenka, 1867), early juveniles produced in captivity. J. World Aquac. Soc. 2019, 50, 763–773. [Google Scholar] [CrossRef]
  11. Boyko, A.V.; Girich, A.S.; Tkacheva, E.S.; Dolmatov, I.Y. The Eupentacta fraudatrix transcriptome provides insights into regulation of cell transdifferentiation. Sci. Rep. 2020, 10, 1–11. [Google Scholar] [CrossRef] [Green Version]
  12. Liang, L.; Chen, J.; Li, Y.; Zhang, H. Insights into high-pressure acclimation: Comparative transcriptome analysis of sea cucumber Apostichopus japonicus at different hydrostatic pressure exposures. BMC Genom. 2020, 21, 68. [Google Scholar] [CrossRef]
  13. Gao, K.; Wang, Z.; Qiu, X.; Song, J.; Wang, H.; Zhao, C.; Wang, X.; Chang, Y. Transcriptome analysis of body wall reveals growth difference between the largest and smallest individuals in the pure and hybrid populations of Apostichopus japonicus. Comp. Biochem. Physiol. Part D Genom. Proteom. 2019, 31, 100591. [Google Scholar] [CrossRef]
  14. Xing, L.; Sun, L.; Liu, S.; Li, X.; Zhang, L.; Yang, H. De Novo assembly and comparative transcriptome analyses of purple and green morphs of Apostichopus japonicus during body wall pigmentation process. Comp. Biochem. Physiol. Part D Genom. Proteom. 2018, 28, 151–161. [Google Scholar] [CrossRef]
  15. Poot-Salazar, A.; Hernández-Flores, Á.; Ardisson, P.-L. Use of the SLW index to calculate growth function in the sea cucumber Isostichopus badionotus. Sci. Rep. 2014, 4, 5151. [Google Scholar] [CrossRef] [Green Version]
  16. Seppey, M.; Manni, M.; Zdobnov, E.M. BUSCO: Assessing genome assembly and annotation completeness. In Gene Prediction. Methods in Molecular Biology; Kollmar, M., Ed.; Humana: New York, NY, USA, 2019; Volume 1962, pp. 227–245. [Google Scholar] [CrossRef]
  17. Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.Z.; Thompson, D.A.; Amit, I.; Adiconis, X.; Fan, L.; Raychowdhury, R.; Zeng, Q.; et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 2011, 29, 644–652. [Google Scholar] [CrossRef] [Green Version]
  18. Wei, X.; Liu, X.; Yang, J.; Wang, S.; Sun, G.; Yang, J. Critical roles of sea cucumber C-type lectin in non-self recognition and bacterial clearance. Fish Shellfish. Immunol. 2015, 45, 791–799. [Google Scholar] [CrossRef]
  19. Lv, Z.; Guo, M.; Li, C.; Shao, Y.; Zhao, X.; Zhang, W. Divergent proteomics response of Apostichopus japonicus suffering from skin ulceration syndrome and pathogen infection. Comp. Biochem. Physiol. Part D Genom. Proteom. 2019, 30, 196–205. [Google Scholar] [CrossRef]
  20. Zhong, M.; Chen, T.; Hu, C.; Ren, C. Isolation and Characterization of Collagen from the Body Wall of Sea Cucumber Stichopus monotuberculatus. J. Food Sci. 2015, 80, C671–C679. [Google Scholar] [CrossRef]
  21. Du, H.; Bao, Z.; Hou, R.; Wang, S.; Su, H.; Yan, J.; Tian, M.; Li, Y.; Wei, W.; Lu, W.; et al. Transcriptome sequencing and characterization for the sea cucumber Apostichopus japonicus (Selenka, 1867). PLoS ONE 2012, 7, e33311. [Google Scholar] [CrossRef]
  22. Ding, K.; Zhang, L.; Sun, L.; Lin, C.; Feng, Q.; Zhang, S.; Yang, H.; Brinkman, R.; Lin, G.; Huang, Z. Transcriptome analysis provides insights into the molecular mechanisms responsible for evisceration behavior in the sea cucumber Apostichopus japonicus. Comp. Biochem. Physiol. Part D Genom. Proteom. 2019, 30, 143–157. [Google Scholar] [CrossRef]
  23. Li, C.; Fang, H.; Xu, D. Effect of seasonal high temperature on the immune response in Apostichopus japonicus by transcriptome analysis. Fish Shellfish. Immunol. 2019, 92, 765–771. [Google Scholar] [CrossRef]
  24. Gao, L.; He, C.; Bao, X.; Tian, M.; Ma, Z. Transcriptome analysis of the sea cucumber (Apostichopus japonicus) with variation in individual growth. PLoS ONE 2017, 12, e0181471. [Google Scholar] [CrossRef] [Green Version]
  25. Liu, H.; Kong, X.; Chen, J.; Zhang, H. De novo sequencing and transcriptome analysis of Stichopus horrens to reveal genes related to biosynthesis of triterpenoids. Aquaculture 2018, 491, 358–367. [Google Scholar] [CrossRef]
  26. Feng, J.; Zhang, L.; Xia, X.; Hu, W.; Zhou, P. Effect of geographic variation on the proteome of sea cucumber (Stichopus japonicus. Food Res. Int. 2020, 136, 109498. [Google Scholar] [CrossRef]
  27. Liu, R.; Liu, J.; Zhang, H. Positive selection analysis reveals the deep-sea adaptation of a hadal sea cucumber (Paelopatides sp.) to the Mariana Trench. J. Oceanol. Limnol. 2021, 39, 266–281. [Google Scholar] [CrossRef]
  28. Xue, Z.; Li, H.; Wang, X.; Li, X.; Liu, Y.; Sun, J.; Liu, C. A review of the immune molecules in the sea cucumber. Fish Shellfish. Immunol. 2015, 44, 1–11. [Google Scholar] [CrossRef]
  29. Xu, D.; Zhou, S.; Sun, L. RNA-seq based transcriptional analysis reveals dynamic genes expression profiles and immune-associated regulation under heat stress in Apostichopus japonicus. Fish Shellfish. Immunol. 2018, 78, 169–176. [Google Scholar] [CrossRef]
  30. Olivera-Castillo, L.; Grant, G.; Kantún-Moreno, N.; Barrera-Pérez, H.A.; Montero, J.; Olvera-Novoa, M.A.; Carrillo-Cocom, L.M.; Acevedo, J.J.; Puerto-Castillo, C.; Solís, V.M.; et al. A glycosaminoglycan-rich fraction from sea cucumber Isostichopus badionotus has potent anti-inflammatory properties in vitro and in vivo. Nutrients 2020, 12, 1698. [Google Scholar] [CrossRef]
  31. Olivera-Castillo, L.; Barrientos, R.G.; Legarreta, I.G.; Sámano, A.H.; Chi, Y.C. Sea Cucumber as a Source of Bioactive Compounds: Current Research on Isostichopus badionotus and Isostichopus fuscus from Mexico. Bioact. Compd. Mar. Foods Plant Anim. Sources 2013, 329–342. [Google Scholar] [CrossRef]
  32. Laguerre, H.; Raymond, G.; Plan, P.; Améziane, N.; Bailly, X.; Le Chevalier, P. First description of embryonic and larval development, juvenile growth of the black sea-cucumber Holothuria forskali (Echinodermata: Holothuroidea), a new species for aquaculture in the north-eastern Atlantic. Aquaculture 2020, 521, 734961. [Google Scholar] [CrossRef]
  33. Purcell, S.W.; Hair, C.A.; Mills, D.J. Sea cucumber culture, farming and sea ranching in the tropics: Progress, problems and opportunities. Aquaculture 2012, 368-369, 68–81. [Google Scholar] [CrossRef]
  34. Gamboa-Álvarez, M.Á.; Poot-López, G.R.; Olvera-Novoa, M.A.; Aguilar-Perera, A.; Ponce-Márquez, M.A.; López-Rocha, J.A. Survival and growth of wild-translocated individuals and released-cultured juveniles of sea cucumber Isostichopus badionotus off the northern Yucatan Peninsula, Mexico. Estuar. Coast. Shelf Sci. 2021, 252, 107273. [Google Scholar] [CrossRef]
  35. Dong, Y.; Li, Q.; Zhong, X.; Kong, L. Development of Gene-derived SNP Markers and Their Application for the Assessment of Genetic Diversity in Wild and Cultured Populations in Sea Cucumber, Apostichopus japonicus. J. World Aquac. Soc. 2016, 47, 873–888. [Google Scholar] [CrossRef]
  36. Dong, G.; Dong, S.; Wang, F.; Tian, X. Effects of light intensity on daily activity rhythm of juvenile sea cucumber, Apostichopus japonicus (Selenka). Aquac. Res. 2010, 41, 1640–1647. [Google Scholar] [CrossRef]
  37. Liu, J.; Li, J.; Cao, X.; Song, S.; Li, J.; Wang, F.; Wang, L.; Liu, S. Effect of temperature on growth, digestion and immunity of green, white and purple sea cucumbers, Apostichopus japonicus (Selenka). Aquac. Res. 2019, 50, 219–225. [Google Scholar] [CrossRef] [Green Version]
  38. Palomar-Abesamis, N.; Juinio-Meñez, M.A.; Slater, M.J. Effects of light and microhabitat on activity pattern and behaviour of wild and hatchery-reared juveniles of Stichopus cf. horrens. J. Mar. Biol. Assoc. UK 2018, 98, 1703–1713. [Google Scholar] [CrossRef]
  39. Gao, L.; Yuan, Z.; Yu, S.; Yang, Y.; Li, Y.; He, C. Genome-wide identification of HSP70/110 genes in sea cucumber Apostichopus japonicus and comparative analysis of their involvement in aestivation. Comp. Biochem. Physiol. Part D Genom. Proteom. 2018, 28, 162–171. [Google Scholar] [CrossRef]
  40. Zhang, L.; Feng, Q.; Sun, L.; Fang, Y.; Xu, D.; Zhang, T.; Yang, H. Differential gene expression in the body wall of the sea cucumber (Apostichopus japonicus) under strong lighting and dark conditions. Acta Oceanol. Sin. 2018, 37, 54–66. [Google Scholar] [CrossRef]
  41. Wang, Q.-L.; Yu, S.-S.; Qin, C.-X.; Dong, S.-L.; Dong, Y.-W. Combined effects of acute thermal and hypo-osmotic stresses on osmolality and hsp70, hsp90 and sod expression in the sea cucumber Apostichopus japonicus Selenka. Aquac. Int. 2014, 22, 1149–1161. [Google Scholar] [CrossRef]
  42. Dong, Y.-W.; Yu, S.-S.; Wang, Q.-L.; Dong, S.-L. Physiological responses in a variable environment: Relationships between metabolism, hsp and thermotolerance in an intertidal-subtidal species. PLoS ONE 2011, 6, e26446. [Google Scholar] [CrossRef] [Green Version]
  43. Leprivier, G.; Remke, M.; Rotblat, B.; Dubuc, A.; Mateo, A.-R.F.; Kool, M.; Agnihotri, S.; El-Naggar, A.; Yu, B.; Somasekharan, S.P.; et al. The eEF2 kinase confers resistance to nutrient deprivation by blocking translation elongation. Cell 2013, 153, 1064–1079. [Google Scholar] [CrossRef] [Green Version]
  44. Wang, Q.; Zhang, X.; Chen, M.; Li, W.; Zhang, P. Comparison of intestinal microbiota and activities of digestive and immune-related enzymes of sea cucumber Apostichopus japonicus in two habitats. J. Oceanol. Limnol. 2017, 36, 990–1001. [Google Scholar] [CrossRef]
  45. Xia, S.; Wang, X. Nutritional and medicinal value. Dev.Aquac. Fish. Sci. 2015, 39, 353–365. [Google Scholar] [CrossRef]
  46. Senadheera, T.R.; Dave, D.; Shahidi, F. Sea cucumber derived type i collagen: A comprehensive review. Mar. Drugs 2020, 18, 471. [Google Scholar] [CrossRef]
  47. Zhao, Z.; Jiang, J.; Pan, Y.; Dong, Y.; Chen, Z.; Zhang, G.; Gao, S.; Sun, H.; Guan, X.; Wang, B.; et al. Temporal dynamics of bacterial communities in the water and sediments of sea cucumber (Apostichopus japonicus) culture ponds. Aquaculture 2020, 528, 735498. [Google Scholar] [CrossRef]
  48. Günay, D.; Emiroğlu, D.; Suzer, C. Seasonal variations of digestive enzymes in sea cucumbers (Holothuria tubulosa, G. 1788) under culture conditions. J. Exp. Zool. A Ecol. Integr. Physiol. 2020, 333, 144–150. [Google Scholar] [CrossRef]
  49. Li, L.; Zhang, Y.; Dong, S.; Wang, F.; Gao, Q. Effects of light intensity on larval development and juvenile growth of sea cucumber Apostichopus japonicus. Aquac. Res. 2019, 50, 2333–2340. [Google Scholar] [CrossRef]
  50. Sun, J.; Zhang, L.; Pan, Y.; Lin, C.; Wang, F.; Kan, R.; Yang, H. Feeding behavior and digestive physiology in sea cucumber Apostichopus japonicus. Physiol. Behav. 2015, 139, 336–343. [Google Scholar] [CrossRef]
  51. Yingst, J.Y. The utilization of organic matter in shallow marine sediments by an epibenthic deposit-feeding holothurian. J. Exp. Mar. Biol. Ecol. 1976, 23, 55–69. [Google Scholar] [CrossRef]
  52. Fujiwara, R.; Yokoi, T.; Nakajima, M. Structure and protein–protein interactions of human UDP-glucuronosyltransferases. Front. Pharmacol. 2016, 7, 388. [Google Scholar] [CrossRef] [Green Version]
  53. Gullian Klanian, M.; Terrats Preciat, M. Effect of pH on temperature controlled degradation of reactive oxygen species, heat shock protein expression, and mucosal immunity in the sea cucumber Isostichopus badionotus. PLoS ONE 2017, 12, e0175812. [Google Scholar] [CrossRef] [Green Version]
  54. Puch-Hau, C.; Sánchez-Tapia, I.A.; Patiño-Suárez, V.; Olvera-Novoa, M.A.; Klanian, M.G.; Quintanilla-Mena, M.; Zapata-Pérez, O. Evaluation of two independent protocols for the extraction of DNA and RNA from different tissues of sea cucumber Isostichopus badionotus. MethodsX 2019, 6, 1627–1634. [Google Scholar] [CrossRef]
  55. Andrews, S. FastQC: A Quality Control Tool for High throughput Sequence Data. 2010. Available online: http://www.bioinformatics.babraham.ac.uk/projects/fastqc (accessed on 9 April 2019).
  56. Zerbino, D.R.; Birney, E. Velvet: Algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008, 18, 821–829. [Google Scholar] [CrossRef] [Green Version]
  57. Schulz, M.H.; Zerbino, D.R.; Vingron, M.; Birney, E. Oases: Robust de novo RNA-seq assembly across the dynamic range of expression levels. Bioinform 2012, 28, 1086–1092. [Google Scholar] [CrossRef]
  58. Fu, L.; Niu, B.; Zhu, Z.; Wu, S.; Li, W. CD-HIT: Accelerated for clustering the next-generation sequencing data. Bioinform 2012, 28, 3150–3152. [Google Scholar] [CrossRef]
  59. Waterhouse, R.M.; Seppey, M.; Simão, F.A.; Manni, M.; Ioannidis, P.; Klioutchnikov, G.; Kriventseva, E.V.; Zdobnov, E.M. BUSCO applications from quality assessments to gene prediction and phylogenomics. Mol. Biol. Evol. 2018, 35, 543–548. [Google Scholar] [CrossRef] [Green Version]
  60. Camacho, C.; Coulouris, G.; Avagyan, V.; Ma, N.; Papadopoulos, J.; Bealer, K.; Madden, T.L. BLAST+: Architecture and applications. BMC Bioinform. 2009, 10, 421. [Google Scholar] [CrossRef] [Green Version]
  61. Sayers, E.W.; Agarwala, R.; Bolton, E.; Brister, J.R.; Canese, K.; Clark, K.; Connor, R.; Fiorini, N.; Funk, K.; Hefferon, T.; et al. Database resources of the National Center for Biotechnology Information. Nucleic Acids Res. 2018, 47, D23–D28. [Google Scholar] [CrossRef] [Green Version]
  62. Petrov, A.I.; Kay, S.; Kalvari, I.; Howe, K.L.; Gray, K.A.; Bruford, E.A.; Kersey, P.J.; Cochrane, G.; Finn, R.D.; Bateman, A.; et al. RNAcentral: A comprehensive database of non-coding RNA sequences. Nucleic Acids Res. 2017, 45, D128–D134. [Google Scholar] [CrossRef]
  63. Götz, S.; García-Gómez, J.M.; Terol, J.; Williams, T.D.; Nagaraj, S.H.; Nueda, M.J.; Robles, M.; Talón, M.; Dopazo, J.; Conesa, A. High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008, 36, 3420–3435. [Google Scholar] [CrossRef] [PubMed]
  64. Kim, D.; Paggi, J.M.; Park, C.; Bennett, C.; Salzberg, S.L. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat. Biotechnol. 2019, 37, 907–915. [Google Scholar] [CrossRef] [PubMed]
  65. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinform 2010, 26, 139–140. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Kanehisa, M.; Sato, Y.; Kawashima, M.; Furumichi, M.; Tanabe, M. KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res. 2016, 44, D457–D462. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Pfaffl, M.W. A new mathematical model for relative quantification in real-time RT–PCR. Nucleic Acids Res. 2001, 29, e45. [Google Scholar] [CrossRef]
  68. Zamora-Briseño, J.A.; Ruiz-May, E.; Elizalde-Contreras, J.M.; Hernández-Velázquez, I.M.; Hernández-Pérez, A.; Fuentes-García, A.G.; Herrera-Salvatierra, N.; Briones-Fourzán, P.; Pascual-Jiménez, C.; Lozano-Álvarez, E.; et al. iTRAQ-Based Proteomic Profile Analysis of the Hepatopancreas of Caribbean Spiny Lobsters Infected With Panulirus argus Virus 1: Metabolic and Physiological Implications. Front. Microbiol. 2020, 11, 1084. [Google Scholar] [CrossRef]
Ijms 22 03882 g001 550
Figure 1. Bioinformatics pipeline used for transcriptome de novo assembly.
Figure 1. Bioinformatics pipeline used for transcriptome de novo assembly.
Ijms 22 03882 g001
Ijms 22 03882 g002 550
Figure 2. Comparison of BUSCO results among three assemblies. Quantitative assessment of transcriptome completeness for the newly assembled wild cucumber transcriptome, and comparison with related publicly available transcriptome assemblies.
Figure 2. Comparison of BUSCO results among three assemblies. Quantitative assessment of transcriptome completeness for the newly assembled wild cucumber transcriptome, and comparison with related publicly available transcriptome assemblies.
Ijms 22 03882 g002
Ijms 22 03882 g003 550
Figure 3. Distribution of annotated unigenes across databases. Venn diagram for common annotated unigenes in NT, UniProt and Gene Ontology databases.
Figure 3. Distribution of annotated unigenes across databases. Venn diagram for common annotated unigenes in NT, UniProt and Gene Ontology databases.
Ijms 22 03882 g003
Ijms 22 03882 g004 550
Figure 4. Gene Ontology (GO) classification of wild sea cucumber body wall assembled transcripts.
Figure 4. Gene Ontology (GO) classification of wild sea cucumber body wall assembled transcripts.
Ijms 22 03882 g004
Ijms 22 03882 g005 550
Figure 5. Transcriptomic comparisons between wild and farmed sea cucumbers. (a) MA plots of differentially expressed genes (DEGs) between wild and farmed sea cucumbers. Red dots show transcripts with a false discovery rate (FDR) < 0.05. (b) PCA and variable contribution to the transcriptomic data. (c) Differentially altered transcripts. (d) Percentage of altered transcripts.
Figure 5. Transcriptomic comparisons between wild and farmed sea cucumbers. (a) MA plots of differentially expressed genes (DEGs) between wild and farmed sea cucumbers. Red dots show transcripts with a false discovery rate (FDR) < 0.05. (b) PCA and variable contribution to the transcriptomic data. (c) Differentially altered transcripts. (d) Percentage of altered transcripts.
Ijms 22 03882 g005
Ijms 22 03882 g006 550
Figure 6. Gene Ontology (GO) analysis of sea cucumber differentially expressed genes modulated in farmed animals.
Figure 6. Gene Ontology (GO) analysis of sea cucumber differentially expressed genes modulated in farmed animals.
Ijms 22 03882 g006
Ijms 22 03882 g007 550
Figure 7. Reverse-transcription real-time PCR validation of eleven genes differentially modulated by captivity. COX1, cytochrome c oxidase subunit 1; TBA2, tubulin α chain 2; TBB, Tubulin β chain; ZN709, zinc finger protein 709; ATPA, ATP synthase subunit α; OST48, dolichyl-diphosphooligosaccharide-protein glycosyltransferase non-catalytic subunit; TBA, tubulin alpha chain; DMBT1, deleted in malignant brain tumors 1; RL39, ribosomal protein L39; RS5, ribosomal protein S5; RL18A, ribosomal protein L18a. Data is expressed as relative expression compared to that of wild sea cucumbers (n = 5).
Figure 7. Reverse-transcription real-time PCR validation of eleven genes differentially modulated by captivity. COX1, cytochrome c oxidase subunit 1; TBA2, tubulin α chain 2; TBB, Tubulin β chain; ZN709, zinc finger protein 709; ATPA, ATP synthase subunit α; OST48, dolichyl-diphosphooligosaccharide-protein glycosyltransferase non-catalytic subunit; TBA, tubulin alpha chain; DMBT1, deleted in malignant brain tumors 1; RL39, ribosomal protein L39; RS5, ribosomal protein S5; RL18A, ribosomal protein L18a. Data is expressed as relative expression compared to that of wild sea cucumbers (n = 5).
Ijms 22 03882 g007
Table
Table 1. Summary statistics of sequence alignment process.
Table 1. Summary statistics of sequence alignment process.
Sample *Number of Paired ReadsAligned Concordantly
0 Times		Aligned Concordantly Exactly 1 TimeAligned Concordantly
>1 Time		Overall Alignment Rate
W143,091,44114,423,783 (33.5%)21,501,037 (49.9%)7,166,621 (16.6%)83.90%
W238,520,78015,694,808 (40.7%)16,972,681 (44.1%)5,853,290 (15.2%)75.50%
W341,558,46016,415,528 (39.5%)18,538,978 (44.6%)6,603,954 (15.9%)76.60%
W437,099,57615,847,351 (42.7%)16,619,310 (44.8%)4,632,915 (12.5%)84.30%
F122,251,49011,211,107 (50.4%)8,994,180 (40.4%)2,046,203 (9.2%)63.60%
F222,603,02510,767,121 (47.6%)9,890,391 (43.7%)1,945,513 (8.6%)67.40%
F329,004,58712,871,831 (44.4%)12,775,028 (44.0%)3,357,728 (11.6%)67.40%
F419,755,6269,614,938 (48.7%)8,230,706 (41.7%)1,909,982 (9/7%)65.90%
* W = Wild; F = Farmed.
Table
Table 2. Selection of transcripts differentially expressed in farmed I. badionotus.
Table 2. Selection of transcripts differentially expressed in farmed I. badionotus.
GeneFold Change
Cytochrome c oxidase subunit I (COXI)+6-fold
Eukaryotic translation initiation factor 2-alpha kinase 3 (EIF2AK3)+2-fold
ATP synthase subunit alpha (APTA)+2.2-fold
Mitogen-activated protein kinase 14 (MAPK14-1)+2.9-fold
Heat Shock Protein-40 (HSP40)+2.3-fold
Heat Shock protein cognate70+1.7-fold
* TRPM2+2.9-fold
Cyclic AMP-responsive element-binding protein 1 (CREB1)−1.9-fold
Techylectin 5A−2.9-fold
Suppressors of cytokine signaling 2 (SOCS2)−2.6-fold
I. badionotus isolate FAO44 16S ribosomal RNA (mitochondrial)−3.9-fold
Collagen mRNAs~−5-fold
* Transient receptor potential cation channel subfamily M member 2-like.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Martín-Hernández, R.; Rodríguez-Canul, R.; Kantún-Moreno, N.; Olvera-Novoa, M.A.; Medina-Contreras, O.; Garikoitz-Legarda, C.; Triviño, J.C.; Zamora-Briseño, J.A.; May-Solis, V.; Poot-Salazar, A.; et al. Comparative Transcriptomes of the Body Wall of Wild and Farmed Sea Cucumber Isostichopus badionotus. Int. J. Mol. Sci. 2021, 22, 3882. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22083882

AMA Style

Martín-Hernández R, Rodríguez-Canul R, Kantún-Moreno N, Olvera-Novoa MA, Medina-Contreras O, Garikoitz-Legarda C, Triviño JC, Zamora-Briseño JA, May-Solis V, Poot-Salazar A, et al. Comparative Transcriptomes of the Body Wall of Wild and Farmed Sea Cucumber Isostichopus badionotus. International Journal of Molecular Sciences. 2021; 22(8):3882. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22083882

Chicago/Turabian Style

Martín-Hernández, Roberto, Rossanna Rodríguez-Canul, Nuvia Kantún-Moreno, Miguel A. Olvera-Novoa, Oscar Medina-Contreras, Cristobal Garikoitz-Legarda, Juan Carlos Triviño, Jesús Alejandro Zamora-Briseño, Víctor May-Solis, Alicia Poot-Salazar, and et al. 2021. "Comparative Transcriptomes of the Body Wall of Wild and Farmed Sea Cucumber Isostichopus badionotus" International Journal of Molecular Sciences 22, no. 8: 3882. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22083882

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop