Molecular and mineral responses of corals grown under artificial Calcite Sea conditions

The formation of skeletal structures composed of different calcium carbonate polymorphs (e.g. aragonite and calcite) appears to be both biologically and environmentally regulated. Among environmental factors influencing aragonite and calcite precipitation, changes in seawater conditions-primarily in the molar ratio of magnesium and calcium during so-called 'Calcite' (mMg:mCa below 2) or 'Aragonite' seas (mMg:mCa above 2)-have had profound impacts on the distribution and performance of marine calcifiers throughout Earth's history. Nonetheless, the fossil record shows that some species appear to have counteracted such changes and kept their skeleton polymorph unaltered. Here, the aragonitic octocoral Heliopora coerulea and the aragonitic scleractinian Montipora digitata were exposed to Calcite Sea-like mMg:mCa with various levels of magnesium and calcium concentration, and changes in both the mineralogy (i.e. CaCO3 polymorph) and gene expression were monitored. Both species maintained aragonite deposition at lower mMg:mCa ratios, while concurrent calcite presence was only detected in M. digitata. Despite a strong variability between independent experimental replicates for both species, the expression for a set of putative calcification-related genes, including known components of the M. digitata skeleton organic matrix (SkOM), was found to consistently change at lower mMg:mCa. These results support the previously proposed involvements of the SkOM in counteracting decreases in seawater mMg:mCa. Although no consistent expression changes in calcium and magnesium transporters were observed, down-regulation calcium channels in H. coerulea in one experimental replicate and at an mMg:mCa of 2.5, pointing to a possible active calcium uptake regulation by the corals under altered mMg:mCa.

Transcriptomic response of Cinachyrella cf. cavernosa sponges to spatial competition

Spatial competition in the intertidal zones drives the community structure in marine benthic habitats. Organisms inhabiting these areas not only need to withstand fluctuations in temperature, water level, pH, and salinity but also need to compete for the best available space. Sponges are key members of the intertidal zones, and their life history processes (e.g. growth, reproduction, and regeneration) are affected by competition. Here, we used transcriptomics to investigate the effects of interspecific competition between the tetillid sponge Cinachyrella cf. cavernosa, the zoantharid Zoanthus sansibaricus and the macroalgae Dictyota ciliolata in the field. The analysis of differentially expressed genes showed that Z. sansibaricus was the more stressful competitor to C. cf. cavernosa, which showed an upregulation of cellular respiration under stress of competition. Similarly, an upregulation of energy metabolism, lipid metabolism and the heat-shock protein (HSP) 70 was also observed along with an increase in viral load and decreased ability to synthesize protein. A downregulation of purine and pyrimidine metabolism indicated a reduction in the physiological activities of the competing sponges. Moreover, a putative case of possible kleptocnidism, not previously reported in C. cf. cavernosa, was also observed. This study offers a glimpse into the inner workings of marine organisms competing for spatial resources using transcriptome data.

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.

Description of four Millepora spp. transcriptomes and their potential to delimit the Caribbean fire coral species

Millepora is a relatively species-rich genus of hydrocorals, with 16 species distributed around the globe. It is considered an important reef building cnidarian. The current diversity of Caribbean Millepora species consists of Millepora complanata, M. alcicornis, M. squarrosa and M. striata. Here, we report the de novo transcriptome assembly and phylotranscriptomic analysis of M. alcicornis, M. complanata, M. squarrosa and a undescribed morphotype (Millepora sp.) found in exposed Thalassia beds and mangrove areas in southwest Puerto Rico. Over 345 million sequence reads were obtained for the analysis of the Millepora transcriptomes (Illumina HiSeq4000; 2x150bp). The analysis pipeline consisted of assembly with Trinity, BUSCO, RSEM and ORFs calling for each transcriptome, followed by ontology (Blast2GO) and phylogenetic analysis. The phylogenetic analysis was performed after selecting homologous genes among the transcriptomes, resulting in 10,596 sequences. Concatenation analysis (Maximum Likelihood and Bayesian inference) and a coalescence-based analysis were performed to the dataset too. Concatenation analysis yielded a topology supporting a clade of M. complanata and M. alcicornis, with Millepora sp. outside this clade and M. squarrosa as an outgroup. The coalescence-based tree estimation analysis (ASTRAL-II), presented a different topology placing M. alcicornis and Millepora sp. as sister taxa, rather than grouping with M. alcicornis with M. complanata. Our coalescence analysis indicated that there is a high degree of incomplete lineage sorting, suggesting a very recent time of species emergence among three out of the four Caribbean Millepora species. Calculations of ABBA-BABA statistics derived from transcriptome-wide SNP data indicate the possible presence of introgression between Millepora complanata and M. alcicornis.

Tracing animal genomic evolution with the chromosomal-level assembly of the freshwater sponge Ephydatia muelleri

The genomes of non-bilaterian metazoans are key to understanding the molecular basis of early animal evolution. However, a full comprehension of how animal-specific traits, such as nervous systems, arose is hindered by the scarcity and fragmented nature of genomes from key taxa, such as Porifera. Ephydatia muelleri is a freshwater sponge found across the northern hemisphere. Here, we present its 326 Mb genome, assembled to high contiguity (N50: 9.88 Mb) with 23 chromosomes on 24 scaffolds. Our analyses reveal a metazoan-typical genome architecture, with highly shared synteny across Metazoa, and suggest that adaptation to the extreme temperatures and conditions found in freshwater often involves gene duplication. The pancontinental distribution and ready laboratory culture of E. muelleri make this a highly practical model system which, with RNAseq, DNA methylation and bacterial amplicon data spanning its development and range, allows exploration of genomic changes both within sponges and in early animal evolution.

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.

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

The plating coral, Agaricia lamarcki is a widely distributed species inhabiting reefs across the Caribbean basin and Florida. This species is of interest since it is considered a depth-generalist, found from 10 to 70 m. Given the scope of contemporary studies on this coral's population dynamics and physiology, as well as, the potential of mesophotic reefs to be refuge habitats for deteriorated shallow water reefs, we present the first de novo transcriptome assembly of an important mesophotic coral. Using next-generation paired-end sequencing (Illumina Hiseq4000; 2 × 150 bp), we obtained a total of 82,506,058 raw reads. The novel transcriptome assembly strategy included the recently developed National Center for Genome Analysis Support de novo transcriptome assembly pipeline. Assembly produced a total of 101,322 biologically true, non-redundant transcripts with an average contig length of 959 and N50 of 1830. EvidentialGene and TransDecoder were used to identify open reading frames (ORFs) with homology insight provided by the UniProtKb and PFAM databases. ORF prediction resulted in 38,517 putative ORFs of which 12,107 ORFs were annotated as genes dealing with molecular function, 1266 with biological processes and 416 with cellular components.