Since the last century, episodes of coral reef bleaching and mortality have occurred almost annually in tropical or subtropical seas. When the temperature exceeds the tolerant limit of a coral–zooxanthellae holobiont, it induces physiological stress and disrupts the vulnerable fine-tuned balance between the two partners, leading to bleaching. The gene expression profiles of a scleractinian coral and its symbiotic zooxanthellae can offer important information with which to decipher this balanced relationship at the functional level of genes. Here, we sequence a full-length transcriptome of a well-known, common and frequently dominant reef-building coral, Pocillopora damicornis, to acquire gene expression information for the coral–zooxanthellae holobiont. To this end, we identify 21,926 and 465 unique genes in the coral and algal symbiont, respectively, and examine the functional enrichment among these genes based on GO (gene ontology) terms and KEGG (the Kyoto Encyclopedia of Genes and Genomes) pathways. The results show that the zooxanthellae provide for their coral host through energy and nutrition metabolism by photosynthesis, and that both the coral host and zooxanthellae have an anti-stress molecular mechanism, though the two parties have independent abilities to survive in the short term. This work sheds light on the valuable gene expression profile of a coral–zooxanthellae holobiont and provides grounds for further molecular biological research to support ecological protection work.

Coral–zooxanthellae holobionts are one of the most productive ecosystems in the ocean. With global warming and ocean acidification, coral ecosystems are facing unprecedented challenges. To save the coral ecosystems, we need to understand the symbiosis of coral–zooxanthellae. Although some Scleractinia (stony corals) transcriptomes have been sequenced, the reliable full-length transcriptome is still lacking due to the short-read length of second-generation sequencing and the uncertainty of the assembly results. Herein, PacBio Sequel II sequencing technology polished with the Illumina RNA-seq platform was used to obtain relatively complete scleractinian coral M. foliosa transcriptome data and to quantify M. foliosa gene expression. A total of 38,365 consensus sequences and 20,751 unique genes were identified. Seven databases were used for the gene function annotation, and 19,972 genes were annotated in at least one database. We found 131 zooxanthellae transcripts and 18,829 M. foliosa transcripts. A total of 6328 lncRNAs, 847 M. foliosa transcription factors (TFs), and 2 zooxanthellae TF were identified. In zooxanthellae we found pathways related to symbiosis, such as photosynthesis and nitrogen metabolism. Pathways related to symbiosis in M. foliosa include oxidative phosphorylation and nitrogen metabolism, etc. We summarized the isoforms and expression level of the symbiont recognition genes. Among the membrane proteins, we found three pathways of glycan biosynthesis, which may be involved in the organic matter storage and monosaccharide stabilization in M. foliosa. Our results provide better material for studying coral symbiosis.