Transcriptome of the Caribbean stony coral Porites astreoides from three developmental stages
GigaScience volume 5, Article number: 33 (2016)
Porites astreoides is a ubiquitous species of coral on modern Caribbean reefs that is resistant to increasing temperatures, overfishing, and other anthropogenic impacts that have threatened most other coral species. We assembled and annotated a transcriptome from this coral using Illumina sequences from three different developmental stages collected over several years: free-swimming larvae, newly settled larvae, and adults (>10 cm in diameter). This resource will aid understanding of coral calcification, larval settlement, and host–symbiont interactions.
A de novo transcriptome for the P. astreoides holobiont (coral plus algal symbiont) was assembled using 594 Mbp of raw Illumina sequencing data generated from five age-specific cDNA libraries. The new transcriptome consists of 867 255 transcript elements with an average length of 685 bases. The isolated P. astreoides assembly consists of 129 718 transcript elements with an average length of 811 bases, and the isolated Symbiodinium sp. assembly had 186 177 transcript elements with an average length of 1105 bases.
This contribution to coral transcriptome data provides a valuable resource for researchers studying the ontogeny of gene expression patterns within both the coral and its dinoflagellate symbiont.
With an increasing focus on threats such as climate change in recent years, there is a growing body of research on the mechanisms underlying coral calcification, and coral response to environmental change [1–6]. Evidence suggests that corals regulate the movement of ions such as bicarbonate, calcium, and hydrogen to facilitate calcification [1, 2, 6, 7], and that some species are more tolerant of changes in their environment [8, 9], yet the mechanisms behind these important processes and their molecular components are unknown. In particular, details of how the symbiotic dinoflagellates (zooxanthellae) enhance calcification and their role in skeleton formation have not been identified to date. A library of gene transcripts from key developmental stages such as settlement can provide valuable information about which genes are important for processes that are turned on at a particular stage, in this case at the onset of calcification. Previous studies have focused individually on the adult stage [3, 5, 10] or the larval stage ; few have examined the holobiont, particularly during the early developmental stages [5, 12–14].
As a model, we used the species Porites astreoides, a rapidly growing stony coral that is ubiquitous in the Caribbean and relatively tolerant to anthropogenic stresses [15, 16]. Thus, this species may provide insight into changes at the molecular level that have allowed it to resist or acclimate to environmental change. Additionally, the availability of brooded, rapidly settling symbiotic larvae that lack a skeleton [9, 17–19] make this an ideal species for studying early events in the development of the machinery used in calcification, and the role of symbionts in skeleton formation. Although 454 sequence data are already available for Acropora millepora larval stages , here we present the first transcriptomes for a coral holobiont at multiple life stages with much higher depth and breadth of coverage than previously published Porites transcriptomes [3, 11].
Here we used deep sequencing to build a combined pre-/post-settling larval and adult coral metatranscriptome, which includes the zooxanthellae symbionts. The addition of early life history stages and the algal symbiont to currently available assemblies will allow events important to settlement and calcification to be studied, as well as how they relate to simultaneous changes in the holobiont’s gene expression [12, 20]. To aid in the capture of rare holobiont transcripts, including those from organisms besides coral and zooxanthellae (e.g., bacteria ), we used relaxed assembly parameters. Though this can potentially include some erroneous assemblies, it is optimized for potentially low abundance organisms. We have included the isolated coral and zooxanthellae transcriptomes through homology filtering to aid in the identification of genes expressed exclusively in the animal or plant components of the holobiont. The annotated transcriptomes presented here should therefore provide a valuable resource for researchers interested in coral calcification, larval settlement, host–symbiont interactions, and development of the complex holobiont.
All corals were collected from El Mario Cay, La Parguera, Puerto Rico (N 17.95258, W 67.0563) in 2012, 2013 and 2014 under permits issued by the Puerto Rico Departamento de Recursos Naturales y Ambientales (Department of Natural and Environmental Resources; DRNA), numbers 2012-IC-051 and #2014-IC-075. Cone traps, developed by Wade Cooper of the Fish and Wildlife Research Institute, Florida (Fig. 1a), were placed and retrieved by SCUBA divers to collect larvae from adult P. astreoides colonies in 2012 and 2014. During the time of maximum larval release (April–July), ten healthy adult colonies between 10 and 25 cm in diameter were removed from the reef and placed on racks in an adjacent sandy area to facilitate anchoring the traps, and to minimize damage to the reef (Fig. 1b). Actively swimming planula larvae were collected from each trap the morning after trap placement (Fig. 1c), transported to the Magueyes Island Marine Laboratory in a large cooler, and immediately processed for sampling. Larvae from each trap were kept separate and divided into groups for RNA extraction, settlement, or further studies. Samples for RNA extraction were placed in a 2 ml Eppendorf tube and completely immersed in RNAlater (Thermofisher Scientific, USA) for extraction at a later date. Samples for settlement were added to 500 ml glass bowls containing settlement tiles (glass, ceramic tile, dead coral skeleton, and shell fragments) that were seasoned for 1–12 months in the field to induce rapid settlement . Once skeleton formation could be seen by eye (3–5 days after settlement; Fig. 1d), individuals were carefully scraped from the surface using a scalpel and placed in RNAlater. A total of 10–25 swimming larvae or newly settled larvae were pooled for extraction to provide enough RNA for sequencing. All samples were from different individuals, except in 2014 where we collected both swimming and settled larvae from the same adult colony. A sample of adult tissue, taken from a single colony in April 2013, contained 4–5 polyps (less than 5 mm diameter plug), and a minimal amount of the underlying skeleton. The adult was different from those that were used to collect larvae. All adult colonies were returned to the reef after sampling and monitored for survivorship. Colony mortality was less than 5 %.
Sample preparation and sequencing
Total RNA was isolated from swimming and recently settled larvae using the RNAqueous kit (Thermo Fisher Scientific, USA). The adult sample was homogenized directly into RNAqueous kit homogenization buffer and processed as per the manufacturer’s instructions. RNA was quantified, checked for quality with an agarose gel, and then sent to the Genomics Core Facility of the Research Technology Support Facilities of Michigan State University. Samples were quantified and quality checked using a BioAnalyzer (RNA integrity number > 8), and then converted into libraries using poly-A selection. Several methods were used to prepare and sequence samples (Table 1). For larval samples from 2012 (SRX1045048 and SRX1045047), the Illumina TruSeq mRNA version 2 library preparation kit was used. For adult samples from 2013 (SRX1045045 and SRX1045046), and larval samples from 2014 (SRX1045052, SRX1045051, SRX1045050, and SRX1045049; two runs per sample of swimming and settled larvae, respectively), the Illumina TruSeq Stranded mRNA library preparation kit LT was used. Libraries were pooled prior to sequencing. Samples from 2012 were sequenced using an Illumina HiSeq 2000 flow cell using TruSeq SBS version 3 reagents. For samples from 2013 and 2014, an Illumina HiSeq 2500 Rapid Run flow cell (v1) was used with Rapid SBS reagents. In 2012, larvae samples were multiplexed to run in one lane. The adult sample was sequenced in different runs on two lanes. Sequencing was done using 2 × 100 bp paired-end sequencing cycles producing 133 268 440, 109 623 674 and 104 806 800 read pairs from swimming larvae, settled larvae, and adult tissues, respectively. In 2014, another two samples from swimming and newly settled larvae were treated in the same way as the samples from 2013 but were sequenced with 2 × 150 bp paired-end sequencing cycles producing 25 183 766 and 20 628 911 read pairs, respectively. Base calling was done by Illumina Real Time Analysis (RTA) v22.214.171.124, and RTA output was demultiplexed and converted to FastQ format with Illumina Bcl2fastq v1.8.4.
Transcriptome preprocessing, assembly and annotation
Sequencing data from all five samples (393 511 591 paired-end reads) were pooled for assembly. Data went through two stages of read trimming: first, we performed quality-based trimming using Trimmomatic v0.33 . A sliding window of 4 bp and trimming threshold of phred score equal to 2 were chosen to maximize sensitivity , followed by K-mer spectral analysis to remove low abundance k-mers using ‘filter-abund.py –V’ from the Khmer 2.0 package . FastQC v0.11.3 was used to check data quality before and after trimming . After filtering, the remaining 391 297 779 high-quality read pairs were used for de novo transcriptome assembly using Trinity v6.0.2 producing 881 402 transcripts . SeqClean was used to trim poly-A tails and remove low complexity sequences . To enhance the quality of the assembly, 1077 short transcripts (<200 bp) were excluded. To exclude 11 380 uncovered isoforms, we also back-mapped input sequencing reads to the assembly using Salmon software, which allows an unambiguous alignment . For functional annotation, assembled transcripts were blasted against the Swiss-Prot database and best hits with p-values less than 1 × e10−3 were selected. Assembly and annotation statistics are listed in Table 2.
To identify zooxanthellae transcripts, the assembled sequences were compared against a homemade library of publically available zooxanthellae transcriptomes using BLASTN for sequences with e-values less than 10−5. The library included transcriptomes for Symbiodinium clades A and B , S. minutum , S. kawagutii, Symbiodinium sp. (Clades C1, C15, CCMP2430, and Mp), Alexandrium fundyense, A. monilatum, A. temarense, Peridinium aciculiferum, Karenia brevis (CCMP2229, SP1, SP3, and Wilson), and Prorocentrum minimum (CCMP1329 and CCMP2233) . To identify a high confidence P. astreoides transcriptome, the zooxanthellae-free transcripts went through a second BLASTN screen using two homemade libraries of publically available Cnidarian transcriptomes and genomic sequences for sequences with e-values less than 10−5. The Cnidarian transcriptome library included Acropora digitifera , A. millepora , A. hyacinthus, A. tenuis , P. astreoides , Nematostella vectensis , and Hydra vulgaris . The genomic library included the genomes of A. digitifera  and N. vectensis , in addition to genomic sequences from the Trace archive  for A. millepora, A. palmata, P. lobata, and M. faveolata.
To assess the differences between our P. astreoides transcriptome and the previously published one , Transrate (v1.0.1)  was used to examine the contigs of both assemblies. We used Transrate to run a conditional reciprocal best BLAST analysis using the older, smaller assembly as a query against our new assembly as a reference . A p-value of 1x10−5 was used as a threshold. Then, we compared the change of transcript length between reciprocal hits and calculated the total gain or loss in the transcript lengths. The conditional reciprocal BLAST using the original assembly (30 740 isoforms) as a query against our assembly as a reference was more sensitive than the classical reciprocal BLAST. This identified 21 232 reciprocal hits (17 382 of which were unique), and resulted in a total gain of approximately 30 Mbp and a loss of 2.5 Mbp in the transcript lengths.
Availability of supporting data
The raw data used for assembly are deposited into the National Center for Biotechnology Information (NCBI) Sequence Reads Archive (SRA) under accession number SRX1045045-52, which is associated with BioProject number PRJNA283441. This Transcriptome Shotgun Assembly project has been deposited at DDBJ/ENA/GenBank under the accession GEHP0000000. Further supporting data can be found in the GigaScience repository, GigaDB .
DRNA, Departamento de Recursos Naturales y Ambientales (Department of Natural and Environmental Resources); NCBI, National Center for Biotechnology Information; RTA, real time analysis; SRA, Sequence Read Archive
Gattuso J-P, Allemand D, Frankignoulle M. Photosynthesis and calcification at cellular, organismal and community levels in coral reefs: a review on interactions and control by carbonate chemistry. Amer Zool. 1999;39:160–83.
Goreau T. The physiology of skeleton formation in corals. I. A method for measuring the rate of calcium deposition by corals under different conditions. Biol Bull. 1959;116:59–75.
Kenkel CD, Meyer E, Matz MV. Gene expression under chronic heat stress in populations of the mustard hill coral (Porites astreoides) from different thermal environments. Mol Eco. 2013;22(16):4322–34.
Moya A. Study of calcification during a daily cycle of the coral Stylophora pistillata: implications for ‘light-enhanced calcification’. J Exp Biol. 2006;209:3413–9.
Shinzato C, Inoue M, Kusakabe M. A snapshot of a coral “holobiont”: a transcriptome assembly of the scleractinian coral, Porites, captures a wide variety of genes from both the host and symbiotic zooxanthellae. PLoS One. 2014;9:e85182.
Zoccola D, Ganot P, Bertucci A, Caminiti-Segonds N, Techer N, Voolstra CR, Aranda M, Tambutté E, Allemand D, Casey JR, et al. Bicarbonate transporters in corals point towards a key step in the evolution of Cnidarian calcification. Sci Rep. 2015;5:9983.
Cai W-J, Ma Y, Hopkinson BM, Grottoli AG, Warner ME, Ding Q, Hu X, Yuan X, Schoepf V, Xu H, et al. Microelectrode characterization of coral daytime interior pH and carbonate chemistry. Nature Comm. 2016;7:11144.
Hughes TP, Baird AH, Bellwood DR, Card M, Connolly SR, Folke C, Grosberg R, Hoegh-Guldberg O, Jackson JBC, Kleypas J, et al. Climate change, human impacts, and the resilience of coral reefs. Science. 2003;301:929–33.
Edmunds PJ. Population biology of Porites astreoides and Diploria strigosa on a shallow Caribbean reef. Mar Ecol Prog Series. 2010;418:87–104.
Mayfield AB, Wang Y-B, Chen C-S, Lin C-Y, Chen S-H. Compartment-specific transcriptomics in a reef-building coral exposed to elevated temperatures. Mol Ecol. 2014;23:5816–30.
Meyer E, Aglyamova GV, Wang S, Buchanan-Carter J, Abrego D, Colbourne JK, Willis BL, Matz MV. Sequencing and de novo analysis of a coral larval transcriptome using 454 GSFlx. BMC Gen. 2009;10:219.
Kaniewska P, Chan C-KK, Kline D, Ling EYS, Rosic N, Edwards D, Hoegh-Guldberg O, Dove S. Transcriptomic changes in coral holobionts provide insights into physiological challenges of future climate and ocean change. PLoS One. 2015;10:e0139223.
Littman R, Willis BL, Bourne DG. Metagenomic analysis of the coral holobiont during a natural bleaching event on the Great Barrier Reef. Envir Micro Rep. 2011;3(6):651-60.
Thurber R, Willner-Hall D, Rodriguez-Mueller B, Desnues C, Edwards RA, Angly F, Dinsdale E, Kelly L, Rohwer F. Metagenomic analysis of stressed coral holobionts. Enviro Micro. 2009;11:2148–63.
Edmunds PJ. Patterns in the distribution of juvenile corals and coral reef community structure in St. John, US Virgin Islands. Mar Ecol Prog Series. 2000;202:113–24.
Edmunds PJ. Long-term dynamics of coral reefs in St. John, US Virgin Islands. Coral Reefs. 2002;21:357–67.
Chornesky EA, Peters EC. Sexual reproduction and colony growth and colony growth in the scleractinian coral Porites astreoides. Biol Bull. 1987;172:161–77.
McGuire MP. Timing of larval release by Porites astreoides in the northern Florida Keys. Coral Reefs. 1998;17:369–75.
Edmunds PJ, Gates R, Gleason D. The biology of larvae from the reef coral Porites astreoides and their response to temperature disturbances. Mar Bio. 2001;139:981–9.
Gust KA, Najar FZ, Habib T, Lotufo GR, Piggot AM, Fouke BW, Laird JG, Wilbanks MS, Rawat A, Indest KJ, Roe BA. Coral-zooxanthellae meta-transcriptomics reveals integrated response to pollutant stress. BMC Genomics. 2014;15(1):1.
Ainsworth TD, Krause L, Bridge T, Torda G, Raina JB, Zakrzewski M, Gates RD, Padilla-Gamino JL, Spalding HL, Smith C, et al. The coral core microbiome identifies rare bacterial taxa as ubiquitous endosymbionts. ISME J. 2015;9(10):2261–74.
Tomascik T. Settlement patterns of Caribbean scleractinian corals on artificial substrata along a eutrophication gradient, Barbados, West Indies. Mar Ecol Prog Series. 1991;77(2):261–9.
MacManes MD. On the optimal trimming of high-throughput mRNA sequence data. Front Genet. 2014;5:13.
Crusoe MR, Alameldin HF, Awad S, Boucher E, Caldwell A, Cartwright R, Charbonneau A, Constantinides B, Edvenson G, Fay S, et al. The khmer software package: enabling efficient nucleotide sequence analysis. F1000Res. 2015;4:900.
Boguski MS, Lowe TMJ, Tolstoshev CM. dbEST — database for “expressed sequence tags”. Nat Genet. 1993;4(4):332–3.
NCBI Trace Archive. Available from: http://www.ncbi.nlm.nih.gov/Traces. Accessed 18 Jun 2016.
Smith-Unna RD, Boursnell C, Patro R, Hibberd JM, Kelly S. TransRate: reference free quality assessment of de-novo transcriptome assemblies. Biorxiv. 2015. http://dx.doi.org/10.1101/021626.
Patro R, Duggal G, Kingsford C. Accurate, fast, and model-aware transcript expression quantification with Salmon. Biorxiv. 2015. http://dx.doi.org/10.1101/021592.
Bayer T, Aranda M, Sunagawa S, Yum LK, Desalvo MK, Lindquist E, Coffroth MA, Voolstra CR, Medina M. Symbiodinium transcriptomes: genome insights into the dinoflagellate symbionts of reef-building corals. PLoS One. 2012;7:e35269.
Shoguchi E, Shinzato C, Kawashima T, Gyoja F, Mungpakdee S, Koyanagi R, Takeuchi T, Hisata K, Tanaka M, Fujiwara M, et al. Draft assembly of the Symbiodinium minutum nuclear genome reveals dinoflagellate gene structure. Curr Biol. 2013;23(15):1399–408.
Marine Microbial Eukaryote Transcriptome Sequencing Project. Available from: http://marinemicroeukaryotes.org/. Accessed 18 Jun 2016.
Shinzato C, Shoguchi E, Kawashima T, Hamada M, Hisata K, Tanaka M, Fujie M, Fujiwara M, Koyanagi R, Ikuta T, et al. Using the Acropora digitifera genome to understand coral responses to environmental change. Nature. 2011;476:320–3.
Moya A, Huisman L, Ball EE, Hayward DC, Grasso LC, Chua CM, Woo HN, Gattuso J-P, Forêt S, Miller DJ. Whole transcriptome analysis of the coral Acropora millepora reveals complex responses to CO2-driven acidification during the initiation of calcification. Mol Eco. 2012;21:2440–54.
matzLAB. Available from: http://www.bio.utexas.edu/research/matz_lab/matzlab/Data.html. Accessed 18 Jun 2016.
Putnam NH, Srivastava M, Hellsten U, Dirks B, Chapman J, Salamov A, Terry A, Shapiro H, Lindquist E, Kapitonov VV, et al. Sea anemone genome reveals ancestral eumetazoan gene repertoire and genomic organization. Science. 2007;317(5834):86–94.
Aubry S, Kelly S, Kumpers BM, Smith-Unna RD, Hibberd JM. Deep evolutionary comparison of gene expression identifies parallel recruitment of trans-factors in two independent origins of C4 photosynthesis. PLoS Genet. 2014;10(6):e1004365.
Mansour TA, Rosenthal JJC, Brown CT, Roberson LM. Supporting data for “Transcriptome of the Caribbean stony coral Porites astreoides from three developmental stages”. GigaScience Database. 2016. http://dx.doi.org/10.5524/100207
The authors thank Magueyes Island Marine Laboratory, University of Puerto Rico Mayagüez, and the Caribbean Coastal Ocean Observing System (CariCOOS), with special thanks to Jorge Corredor and Julio Morell for marine operations and laboratory support in La Parguera. Mayra Sanchez, Melissa Melendez and Team Larvae provided essential assistance in the field collections.
Bioinformatic analysis was performed in part on computing resources at the University of Puerto Rico (UPR) Puerto Rico Center for Environmental Neuroscience (PRCEN)’s High Performance Computing Facility, which is supported by: Institutional Development Award Networks of Biomedical Research Excellent (INBRE) grant P20GM103475 from the National Institute of General Medical Sciences, National Institutes of Health; the Institute for Functional Nanomaterials (IFN) award from the Experimental Program to Stimulate Competitive Research (EPSCoR) Track 1 program of the National Science Foundation (NSF); and EPSCoR Track 2 awards for computational nanoscience (EPS 1002410, EPS 1010094). Funding and support of the research was provided by PRCEN thanks to an NSF Centers of Research Excellent in Science and Technology (CREST) award, number HRD-1137725.
LR, JR, and TM conceived and designed the experiments; TM, LR, and JR performed the experiments; TM, LR, JR, and CTB analyzed the data; LR and JR evaluated the conclusions; and LR, JR and TM wrote the paper. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
About this article
Cite this article
Mansour, T.A., Rosenthal, J.J.C., Brown, C.T. et al. Transcriptome of the Caribbean stony coral Porites astreoides from three developmental stages. GigaSci 5, 33 (2016). https://doi.org/10.1186/s13742-016-0138-1
- Porites astreoides
- Swimming larvae
- Larval settlement