Single symbiotic cell transcriptome sequencing of coral

Abstract

Zooxanthellae and coral can form an intracellular symbiotic system. Yet, little is known about the molecular mechanism underlying this symbiosis. In this study, we characterized the symbiosis based on analyses of gene expression at the single-cell level. Among 9110 single coral cells, we identified 4871 symbiotic cells based on the detection of both coral and zooxanthellae gene transcripts within a single cell. Using the bioinformatics tool Seurat, symbiotic cells were further clustered into five groups, 52 genes exhibited differential expression between groups. We proposed an index called the "symbiosis index", to indicate the degree of gene expression of both species in a single symbiotic cell. Interestingly, the index differed distinctly among the five groups. The symbiosis index was highly correlated with the expression of the coral gene gfas1.m1.6761 (ANKRD40), which encodes ankyrin repeat domain-containing protein 40 and is involved in DNA replication (r = 0.76). Two metabolism-related genes, DAGLA and betaGlu, were highly expressed in cells with a high symbiosis index. Four zooxanthellae genes, PRPF19, ATRN, aAA-ATPases and AK812-SmicGene44833, exhibited substantial changes in expression levels when zooxanthellae lived within coral. A trajectory analysis suggested that cells with a higher symbiosis index may be derived from those with a lower index during coral colony development. Taken together, our results provide evidence for zooxanthellae residing within coral, forming a symbiotic system. The symbiosis index is an effective indicator of different cell groups, with lineage relationships among groups. Additionally, we identified specific genes that exhibit expression changes in the symbiotic system.

Keywords: Coral; Single-cell RNA sequencing; Symbiosis; Zooxanthellae.