De novo transcriptome assembly of the coral Agaricia lamarcki (Lamarck's sheet coral) from mesophotic depth in southwest Puerto Rico

Analysis of Crassostrea gasar transcriptome reveals candidate genes involved in metal metabolism

Oysters have been extensively employed for monitoring of metal pollution in dynamic aquatic ecosystems. Therefore, the use of specific biomarkers can assist in discriminating the ecotoxicological implications of different elements in such complex environments. In this study, we revisited the sequencing data of gills and digestive glands transcripts in the mangrove oyster Crassostrea gasar and generated a reference transcriptome assembly from multiple assemblers, seven in total. Overall, we were able to identify a total of 11,917 transcripts, with 86.6% of them being functionally annotated and 1.4 times more than the first annotation. We screened the annotated transcripts to identify genes potentially involved in metals' transport, storage, and detoxification. Our findings included genes related to Zn distribution in cells (Zn transporters - ZIP, ZnT), metallothionein (MT-I and MT-IV), GSH biosynthesis, Ca⁺ transporter (NCX and ATP2B), and Cu distribution in cells (ATP7, ATOX1, CCS, and laccase-like). These results provided a reference transcriptome for additional insights into the transcriptional profile of C. gasar and other bivalves to better understand the molecular pathways underpinning metal tolerance and susceptibility. The study also provided an auxiliary tool for biomonitoring metal contamination in dynamic environments as estuaries.

TransPi - a comprehensive TRanscriptome ANalysiS PIpeline for de novo transcriptome assembly

The use of RNA-Seq data and the generation of de novo transcriptome assemblies have been pivotal for studies in ecology and evolution. This is distinctly true for non-model organisms, where no genome information is available. In such organisms, studies of differential gene expression, DNA enrichment baits design, and phylogenetics can all be accomplished with de novo transcriptome assemblies. Multiple tools are available for transcriptome assembly, however, no single tool can provide the best assembly for all datasets. Therefore, a multi assembler approach, followed by a reduction step, is often sought to generate an improved representation of the assembly. To reduce errors in these complex analyses while at the same time attaining reproducibility and scalability, automated workflows have been essential in the analysis of RNA-Seq data. However, most of these tools are designed for species where genome data is used as reference for the assembly process, limiting their use in non-model organisms. We present TransPi, a comprehensive pipeline for de novo transcriptome assembly, with minimum user input but without losing the ability of a thorough analysis. A combination of different model organisms, k-mer sets, read lengths, and read quantities were used for assessing the tool. Furthermore, a total of 49 non-model organisms, spanning different phyla, were also analysed. Compared to approaches using single assemblers only, TransPi produces higher BUSCO completeness percentages, and a concurrent significant reduction in duplication rates. TransPi is easy to configure and can be deployed seamlessly using Conda, Docker and Singularity.

A comparative genomics study of neuropeptide genes in the cnidarian subclasses Hexacorallia and Ceriantharia

Background

Nervous systems originated before the split of Proto- and Deuterostomia, more than 600 million years ago. Four animal phyla (Cnidaria, Placozoa, Ctenophora, Porifera) diverged before this split and studying these phyla could give us important information on the evolution of the nervous system. Here, we have annotated the neuropeptide preprohormone genes of twenty species belonging to the subclass Hexacorallia or Ceriantharia (Anthozoa: Cnidaria), using thirty-seven publicly accessible genome or transcriptome databases. Studying hexacorals is important, because they are versatile laboratory models for development (e.g., Nematostella vectensis) and symbiosis (e.g., Exaiptasia diaphana) and also are prominent reef-builders.

Results

We found that each hexacoral or ceriantharian species contains five to ten neuropeptide preprohormone genes. Many of these preprohormones contain multiple copies of immature neuropeptides, which can be up to 50 copies of identical or similar neuropeptide sequences. We also discovered preprohormones that only contained one neuropeptide sequence positioned directly after the signal sequence. Examples of them are neuropeptides that terminate with the sequence RWamide (the Antho-RWamides). Most neuropeptide sequences are N-terminally protected by pyroglutamyl (pQ) or one or more prolyl residues, while they are C-terminally protected by an amide group. Previously, we isolated and sequenced small neuropeptides from hexacorals that were N-terminally protected by an unusual L-3-phenyllactyl group. In our current analysis, we found that these N-phenyllactyl-peptides are derived from N-phenylalanyl-peptides located directly after the signal sequence of the preprohormone. The N-phenyllactyl- peptides appear to be confined to the hexacorallian order Actiniaria and do not occur in other cnidarians. On the other hand, (1) the neuropeptide Antho-RFamide (pQGRFamide); (2) peptides with the C-terminal sequence GLWamide; and (3) tetrapeptides with the X₁PRX₂amide consensus sequence (most frequently GPRGamide) are ubiquitous in Hexacorallia.

Conclusions

We found GRFamide, GLWamide, and X₁PRX₂amide peptides in all tested Hexacorallia. Previously, we discovered these three neuropeptide classes also in Cubozoa, Scyphozoa, and Staurozoa, indicating that these neuropeptides originated in the common cnidarian ancestor and are evolutionarily ancient. In addition to these ubiquitous neuropeptides, other neuropeptides appear to be confined to specific cnidarian orders or subclasses.

De novo transcriptome assembly of the digitate morphotype of Briareum asbestinum (Octocorallia: Alcyonacea) from the southwest shelf of Puerto Rico

Octocorals have now become the most visually dominant metazoan benthic taxa of most Caribbean reefs, following the precipitous decline of scleractinian corals. Yet taxonomic issues because of their extensive phenotypic plasticity are still abound. Briareum asbestinum one of the iconic octocorals of the shallow Caribbean coral reefs exhibits a biform morphology, the digitate and the encrusting one. The taxonomic status of each form has not been clarified, yet. Until recently, there were few genetic resources for non-model metazoans, however, affordable high-throughput DNA sequencing has removed this hindrance. We present the first transcriptome of the digitate form of Briareum asbestinum from southwest Puerto Rico. We used paired-end sequencing (Illumina NextSeq 500), with a total yield of 159,754,702 raw reads. De novo assembly was performed utilizing a multi-assembler approach generating 371,554 biologically true, non-redundant transcripts. Open reading frame analysis identified 102,839 putative ORFs of which 78,607 were with annotations. BUSCO analysis indicated a total of 96.4% complete orthologous genes from the metazoan dataset. The assembly presented here serves as an important new genomic reference for the Briareum genus that will facilitate future population and phylogenetic studies aiming to better understand the molecular basis of phenotypic plasticity exhibited throughout the genus.

Global Neuropeptide Annotations From the Genomes and Transcriptomes of Cubozoa, Scyphozoa, Staurozoa (Cnidaria: Medusozoa), and Octocorallia (Cnidaria: Anthozoa)

During animal evolution, ancestral Cnidaria and Bilateria diverged more than 600 million years ago. The nervous systems of extant cnidarians are strongly peptidergic. Neuropeptides have been isolated and sequenced from a few model cnidarians, but a global investigation of the presence of neuropeptides in all cnidarian classes has been lacking. Here, we have used a recently developed software program to annotate neuropeptides in the publicly available genomes and transcriptomes from members of the classes Cubozoa, Scyphozoa, and Staurozoa (which all belong to the subphylum Medusozoa) and contrasted these results with neuropeptides present in the subclass Octocorallia (belonging to the class Anthozoa). We found three to six neuropeptide preprohormone genes in members of the above-mentioned cnidarian classes or subclasses, each coding for several (up to thirty-two) similar or identical neuropeptide copies. Two of these neuropeptide preprohormone genes are present in all cnidarian classes/subclasses investigated, so they are good candidates for being among the first neuropeptide genes evolved in cnidarians. One of these primordial neuropeptide genes codes for neuropeptides having the C-terminal sequence GRFamide (pQGRFamide in Octocorallia; pQWLRGRFamide in Cubozoa and Scyphozoa; pQFLRGRFamide in Staurozoa). The other primordial neuropeptide gene codes for peptides having RPRSamide or closely resembling amino acid sequences. In addition to these two primordial neuropeptide sequences, cnidarians have their own class- or subclass-specific neuropeptides, which probably evolved to serve class/subclass-specific needs. When we carried out phylogenetic tree analyses of the GRFamide or RPRSamide preprohormones from cubozoans, scyphozoans, staurozoans, and octocorallia, we found that their phylogenetic relationships perfectly agreed with current models of the phylogeny of the studied cnidarian classes and subclasses. These results support the early origins of the GRFamide and RPRSamide preprohormone genes.