ORIGINAL RESEARCH article

Front. Mar. Sci., 28 February 2025

Sec. Marine Biology

Volume 12 - 2025 | https://doi.org/10.3389/fmars.2025.1551763

Single-cell transcriptome sequencing reveals the cell populations and functions associated with stony coral Pocillopora damicornis

  • 1School of Computer and Information Engineering, Luoyang Institute of Science and Technology, Luoyang, China
  • 2State Key Laboratory of Bioelectronics; School of Biological Science & Medical Engineering, Southeast University, Nanjing, China

Introduction: Pocillopora damicornis, a key species of stony corals, has been the subject of considerable scientific study. However, the cellular composition of P. damicornis and the roles of these cells in endosymbiosis and biomineralization remain elusive. The development of single-cell technology has provided new opportunities for researching the cellular and molecular mechanisms underlying symbiosis and mineralization. Nevertheless, the stringent environmental requirements, the complexity of the cellular components, and the paucity of high-quality reference genomes of P. damicornis have posed significant challenges for single-cell transcriptome research.

Methods: In this study, we quantified the transcriptomic expression of P. damicornis by aligning its single-cell transcriptome (scRNA-seq) data to multiple species, including Stylophora pistillata, P. damicornis, and Pocillopora verrucosa. We determined the cell types of P. damicornis by comparing its cluster-specific genes with the published cell type-specific genes of S.pistillata and conducted gene function and enrichment analyses.

Results: Unsupervised clustering analysis yielded the identification of ten distinct cell populations, including epidermis cells, gastrodermis cells, algae-hosting cells, calicoblast, cnidocytes, and immune cells. In addition, we identified 53 genes that were highly similar to known sequences in the symbiotic zooxanthellae. These genes were mainly expressed in four different cell populations, corresponding to active symbiotic populations.

Conclusion: This study identified cell types closely associated with symbiosis and calcification in P. damicornis, along with their marker genes, which are consistent with the findings in S. pistillata. These results offer insights into the cellular functions and symbiotic mechanisms of P. damicornis.

Introduction

Coral reefs have been likened to tropical rainforests in the context of the oceans. These ecosystems provide a habitat for fish and play a role in the biodiversity of the ocean (Hughes et al., 2017). Coral reefs are primarily formed by the accumulation of calcium carbonate secreted by reef-building corals over extended periods of time. The survival and reef-building activities of stony corals are inseparable from the important contributions of the algae that live in symbiosis with their cells, known as Symbiodinium (Wilson and Hilferty, 1931). In the event of a disruption to the external environment, such as an increase in temperature, the symbiont undergoes disintegration, a process referred to as coral bleaching (Odum and Odum, 1955; Muscatine, 1967; Moberg and Folke, 1999; Costanza et al., 2014; Huang et al., 2019; Reimer et al., 2019). The ongoing threat to the coral ecosystem, precipitated by global warming and marine environmental pollution, has led to coral bleaching as a primary factor in the degradation of coral reefs (Swain et al., 2016). The study of the physiological characteristics and symbiotic mechanisms of corals is of great significance to the preservation of coral ecosystems. The primary methodologies employed in coral research encompass the use of isotope labeling (Clayton and Lasker, 1984), and metabolite chromatography analysis (Biel et al., 2007; Gates et al., 1995). As coral genome data continues to improve (Shinzato et al., 2011; Polato et al., 2011; Levy et al., 2021), transcriptomic analysis has become a staple of coral research. In comparison with conventional transcriptional analysis, single-cell RNA sequencing (scRNA-seq) has been shown to identify the transcriptional profile of individual cells and to identify multiple types of cells in Stylophora pistillata, including immune cells, gastrodermis cells, neurons cells, alga-hosting cells, and calicoblasts (Mansour et al., 2016; Hayward et al., 2011; Reyes-Bermudez et al., 2009; Bay et al., 2009). Nevertheless, corals have evolved over 500 million years, and there are substantial disparities between their genomes. Consequently, the analysis of single-celled organisms belonging to other stony coral families remains a significant area of research (Reimer et al., 2019).

To this end, the present study employs P. damicornis to elucidate its cell specificity. Both P. damicornis and S. pistillata are classified under the monophyletic group of the Pocilloporidae family. The family Pocilloporidae, belonging to the order Scheractinia(stony corals),is predominantly distributed in the Indo-Pacific Ocean. The family includes several genera, notably Pocillopra and Stylophora (Schmidt-Roach et al., 2014). P. damicornis exhibits a broad distributed across tropical and subtropical reefs, with notable populations in the Indo-Pacific and Eastern regions (Schmidt-Roach et al., 2014) (Barbosa et al., 2016). This study sampled cells from different parts of adult P. damicornis while retaining a small number of symbiotic algae genes to facilitate the identification of symbiotic cells.

Materials and methods

Collection and cultivation of samples

The biological samples utilized in this study were obtained from the Xisha Islands, situated in the South China Sea (15°40’–17°10’ N, 111°–113°E). The coral samples were cultivated in a coral tank specifically designed to replicate their natural habitat conditions. All P. damicornis specimens were maintained in a RedSea® tank (redsea575, Red Sea Aquatics Ltd., London, UK) at 26°C and 1.025 salinity (Red Sea Aquatics Ltd). The temperature was maintained at 26°C and the salinity at 1.025 (as per Red Sea Aquatics Ltd. specifications). The coral cultivation system was equipped with three coral lamps (AI®, Red Sea Aquatics Ltd), a protein skimmer (regal250s, Reef Octopus), a water chiller (tk1000, TECO Ltd, Port Louis, Mauritius), two wave-making devices (VorTechTM MP40, EcoTech Marine Ltd), and a calcium reactor (Calreact 200, Reef Octopus).

Single-cell suspension preparation, library construction, and sequencing

To prepare the single ‐cell suspension, we first diluted a 0.02 M PBS (Solarbio, P1010‐2L, Beijing, China)to obtain a 0.035 M PBS solution with a pH of 8.0. The sample was thoroughly rinsed with this solution, and the rinsing process was repeated three times.

Next, we prepared a digestion solution by mixing a 0.1% trypsin (Sigma-Aldrich, T1426-500MG, St. Louis, MO, USA) with the 0.035 M PBS solution. The sample was then placed in this digestion solution at room temperature for 2 - 3 hours. Subsequently, the solution was filtered through a filter with a pore size of 40 μm Fisherbrand™ Sterile Cell Strainers (Thermo Fisher Scientific, 22‐363‐547, Waltham, MA, USA) and centrifuged at a centrifugal force of 350 G for 5‐6 minutes. The enzymatic reaction was terminated by adding a stop solution containing 0.4 mg of BSA.

We employed the trypan blue staining method. An automatic cell counter (TC20™, Bio‐Rad Laboratories, Inc, Hercules, CA, USA) was used to add a 0.4% trypan blue (Invitrogen, T10282, Carlsbad, CA, USA) to the cell suspension at a 1:1 ratio. To ensure the optimal single - cell capture effect, we defined the qualified criteria as a cell survival rate greater than 60% and a cell concentration ranging from 1×105 to 1×107 cells/mL. After washing and resuspension, the cell suspension was diluted to a concentration of 800 - 1200 cells/μL. The prepared cell suspension was then sent to Tianjin Novogene Company for library construction(10x Genomics, Inc, Pleasanton, CA, USA) and quality control. Finally, sequencing was carried out on the HiSeq PE150 platform(Illumina, Inc, San Diego, CA, USA).

Bioinformatics analysis

Reference genome construction and quantification

Due to the current paucity of high-quality coral genome data and the presence of numerous gap regions, the species with relatively high-quality genome data in the family Pocilloporidae, suborder Astrocoeniina, order Scleractinia on NCBI currently include S. pistillata, P. damicornis, and Pocillopora verrucosa. This study utilizes these three species as references, in conjunction with our laboratory’s P. damicornis full-length transcriptome data (Han et al., 2022), employing the CellRanger (V8.0.1) mkref function to construct a reference genome and the count function for gene expression quantification.

Identification of symbiotic zooxanthellae genes

To identify symbiotic zooxanthellae genes, we utilized the genome of S. pilosum (NCBI accession: GCA_964212085.1) as a reference. Through CellRanger count analysis to detect symbiodinium homologous genes in coral host cell.

Gene function annotation

The gene functions were annotated based on the previously published P. damicornis full-length transcriptome (Han et al., 2022) using Pfam (protein family; Finn et al., 2016); KOG/COG (EuKaryotic Orthologous Groups; Clusters of Orthologous Groups of proteins; Tatusov et al., 2003); KEGG (Kyoto Encyclopedia of Genes and Genomes; Kanehisa et al., 2004); and GO (Gene Ontology; Ashburner et al., 2000).

Cluster and marker gene identification

An unsupervised clustering algorithm with a resolution parameter set at 0.1 was employed for cell clustering. For the analysis of differential expression, we established thresholds of a minimum expression percentage (min.pct) of ≥ 0.01 and a log fold change (logfc.threshold) of > 0.25. Dimensionality reduction was conducted utilizing the Uniform Manifold Approximation and Projection (UMAP) algorithm, which facilitates the low-dimensional visualization of high-dimensional data.

Cell type identification

Presently, the extant literature on marker genes for coral cell types is scant, with the study on S. pistillata constituting a rare resource. The present study employs a three-step approach to address this paucity of knowledge. First, the significantly upregulated genes that are uniquely expressed in each cluster identified in the P. damicornis dataset will be determined. Second, these cluster-specific marker genes are then cross-referenced with the cell type-specific markers reported for S. pistillata. Finally, combining the enrichment analysis result of the overlapping marker genes will allow for the inference of the corresponding cell types present in P. damicornis.

Function enrichment analysis

Functional enrichment analysis was conducted using the Molecular Signatures Database (MSigDB) framework (Liberzon et al., 2015). The analysis pipeline incorporated multiple resources, including GO terms for biological processes, molecular functions, and cellular components; KEGG pathways for metabolic and signaling networks; KOG for evolutionary relationships and Pfam for protein domain information. The statistical significance of the enrichment patterns was assessed using Fisher’s exact test, with a significance threshold of P < 0.05. Enriched terms were prioritized based on their statistical significance (−log10[P-value]). The resulting functional categories were then visualized using the ggplot2 package, with particular emphasis on biological processes and molecular pathways relevant to coral biology.

Result

Single-cell genome mapping

In this study, single-cell RNA sequencing (scRNA-seq) analysis was performed on five reference genomesfrom the family Pocilloporidae: S. pistillata, P. damicornis, P. verrucosa, and P. damicornis full-length transcriptome data, and the full-length transcriptome data of S. pistillata and P. damicornis. Among the species examined, P. verrucosa exhibited the highest number of expressed genes (32,915), followed by P. damicornis (23,077), S. pistillata (22,408), and P. damicornis (14,402) (see Table 1 for details).

Table 1
www.frontiersin.org

Table 1. Comparative analysis of reference genomes and integrated datasets used for single-cell transcriptome alignment.

We observed samples-specific variations in the average number of genes expressed per cell. P. verrucosa demonstrated the highest average expression levels, with 87 genes per cell. In contrast, P. damicornis, S. pistillata, and P. damicornis exhibited average expressions of 75, 74, and 74 genes per cell, respectively. Integration of transcriptomic data from all five species yielded a comprehensive dataset comprising 53,661 genes, with an average of 260 genes expressed per cell across the integrated dataset (Table 1).

Clustering and differential expression analysis

UMAP visualization revealed ten distinct clusters in two-dimensional space (Figure 1A). Differential expression analysis was performed to characterize the molecular signatures of these clusters, with significance criteria set as: Absolute value of log2 fold change > 0.25, expression in > 1% of cells, and P < 0.05. The differential expression results of each cluster are shown in Table 2. The specific differentially expressed genes for each cluster are in Supplementary Table S1.

Figure 1
www.frontiersin.org

Figure 1. A study of the coral symbiotic system via single-cell transcriptome analysis. (A) Utilizing UMAP dimensionality reduction of single-cell transcriptome data, ten transcriptionally distinct cell subpopulations were identified. The utilization of distinct colors and numerical identifiers enables the clear representation of these diverse cell clusters. (B) Annotation of major cell types on the UMAP plot, including epidermis cells, gastrodermis cells, alga hosting cells, calicoblast, cnidocytes, and immune cells. (C) The expression distribution of algal symbiont genes on the UMAP plot is represented by different colors, with different colors representing different expression levels. (D) A schematic diagram of the main cell types identified is shown, with their spatial relationships within the coral symbiotic system illustrated.

Table 2
www.frontiersin.org

Table 2. Cell cluster statistics and marker gene identification.

Symbiotic cell identification

Among the 53 identified Symbiodinium genes, 25 were functionally annotated (Supplementary Tables S3, S4). These genes were primarily highly expressed in clusters 4, 5, 8, and 9 (Figure 1C). These genes are mainly involved in energy metabolism, biosynthesis, stress response (including HSP family), cytoskeleton regulation, and protein synthesis and modification.

Cell type identification

Cross-species comparison of cell-type markers between S. pistillata and P. damicornis revealed distinct cell population correspondences. Based on shared marker genes, we identified six major cell types: cluster 0 as epidermis cells, cluster 1 as gastrodermis cells, cluster 2 as algal host cells, cluster 3 as calicoblasts, cluster 6 as cnidocytes, and cluster 7 as immune cells(Figure 1B; Supplementary Table S2). This cell-type assignment was supported by overlapping marker genes between the two species (Supplementary Table S1, S8), providing molecular evidence for conserved cell populations across coral species.

The ectoderm predominantly comprises of epidermis cells (Cluster 0), which form the outermost protective barrier and contain cnidocytes (Cluster 6) for defense and predation. The endoderm, on the other hand, contains gastrodermis cells (Cluster 1) and various symbiotic cell subpopulations (Clusters 2, 4, 5, 8, and 9), which are responsible for nutrient absorption and maintaining the symbiotic relationship with zooxanthellae. The calicoblast of the basal ectoderm (Cluster 3) is responsible for bone formation, while the immune cells distributed in the mesoderm (Cluster 7) are responsible for immune defense (Figures 1B, D; Supplementary Table S3).

Differentially expressed gene enrichment analysis

Gene set enrichment analysis was performed using MSigDB against KEGG and GO databases. Additionally, transcripts were categorized based on GO, KEGG, KOG/COG, and Pfam annotations, followed by Fisher’s exact test (P < 0.05 considered significant; Supplementary Table S4).

Epidermis cells

GO analysis of upregulated genes in coral epidermis cells revealed enrichment in transcription elongation, DNA binding, and ribosome biogenesis, suggesting active gene expression and protein synthesis regulation (Mayfield et al., 2014). KEGG analysis identified enriched pathways in metabolism and signal transduction (Rädecker et al., 2021). KOG and Pfam analyses highlighted specific regulatory elements, including BTG1/TOB anti-proliferative factors, Mafk bZIP transcription factor, and basic leucine zipper domains, which may modulate epidermal cell functions (Li et al., 2018; Supplementary Figure S1).

Gastrodermis cells

Functional enrichment analysis of upregulated genes in gastrodermis cells revealed multiple digestive-related pathways (Figure 2). These included catalytic activity regulation and cysteine peptidase activity, indicating active enzyme secretion (e.g., pepsinogen); lipid transport and lipoprotein pathways; and N-acetyltransferase activity for substrate modification. Signal transduction analysis showed enrichment in the PI3K-Akt pathway, while structural features were highlighted by enrichment in cell adhesion and actin cytoskeleton components, reflecting their epithelial barrier function (Raz-Bahat et al., 2017). The enrichment of lysosomal and protein digestion pathways further supports their role in protein processing (Coffey and De Duve, 1968; Moya et al., 2012).

Figure 2
www.frontiersin.org

Figure 2. Enrichment analysis of genes and proteins in coral gastrodermis cells. The figure shows the functional enrichment analysis results of genes and proteins categorized by biological process (BP), molecular function (MF), cellular component (CC), and KEGG pathways, as well as the enrichment fold changes of specific genes and proteins. The x-axis represents the negative logarithm transformation of the p-value (-logPV), with higher values indicating higher enrichment levels. “-logPV” represents the negative logarithm (base 10) transformation of the p-value.

Calicoblast

GO functional enrichment analysis of coral calicoblast revealed significant pathways in energy metabolism, ion transport, and cytoskeletal regulation. Biological processes showed enrichment in ATP synthesis coupled proton transport, mitochondrial electron transport, and anion transport. Molecular functions highlighted cytochrome c oxidase activity, NADH dehydrogenase activity, and ion transport mechanisms. Cellular component analysis identified enrichment in mitochondrial structures and ion transport complexes. KEGG pathway analysis revealed enrichment in oxidative phosphorylation and related metabolic pathways (Supplementary Figure S2). KOG and Pfam analyses identified key components including carbonic anhydrases, cytoskeletal proteins, ion transporters (Na+/K+-ATPase), and regulatory domains (FerB, calponin homology, tRNA synthetase classes).

Cnidocyte

GO functional enrichment analysis of coral cnidocytes revealed significant enrichment in protein deubiquitination, calcium ion homeostasis, and protein processing pathways. Molecular function analysis identified enrichment in various enzymatic activities (metallocarboxypeptidase, ubiquitin-specific protease), calcium channels, and neurotransmitter transport. Cellular component analysis highlighted the DNA-dependent RNA polymerase III complex, suggesting active transcriptional regulation.

KEGG pathway analysis showed enrichment in protein export, endoplasmic reticulum processing, and mucin-type O-glycosylation biosynthesis, potentially involved in toxin modification (Supplementary Figure S3). KOG classification revealed enrichment in signal transduction, cytoskeleton organization, and protein modification. Pfam analysis identified domains crucial for cnidocyte function, including EF-hand calcium-binding, dynein heavy chain, myosin head, and toxin-related domains.

These findings indicate that cnidocytes possess specialized features for toxin synthesis and secretion, including protein modification systems, calcium signaling pathways, and dynamic cytoskeletal regulation (Jouiaei et al., 2015; Özbek, 2011; Beckmann and Özbek, 2012; Rachamim et al., 2015).

Immune cells

GO analysis revealed enrichment in protein deubiquitination, catalytic activity regulation, DNA binding, and exosome complex functions. KEGG pathways highlighted immune-related processes, including lysosomal and platelet activation pathways. KOG and Pfam analyses identified regulatory proteins such as E3 ubiquitin ligases and interferon factors, suggesting roles in coral immune regulation (Supplementary Figure S4).

Discussion

As the cornerstone of the marine ecosystem, the survival of coral depends on the functional collaboration between cell groups (Levy et al., 2021). We revealed the molecular characteristics of certain cell populations in P. damicornis and their synergistic mechanisms in symbiosis, calcification and immunity through single-cell transcriptome analysis.

Functional collaboration network between cell populations

Epithelial cells act as the first barrier to the external environment, and they quickly respond to stress by highly expressing HSP90AB1 and Jun genes (Cziesielski et al., 2018).The cytokines they secrete may regulate the metabolic activity of adjacent gastric epithelial cells through paracrine signals. The significant upregulation of ribosomal proteins (RPL13, RPS20) and collagen (COL1A2, COL5A1) in gastric epithelial cells indicates that they are not only responsible for nutrient absorption and transport, but may also provide structural support for the symbiotic algae through extracellular matrix (ECM) remodeling (Mass et al., 2014). It is worth noting that the genes for lipid metabolism (ALDH3A2, ACSL4) in symbiotic cells are highly expressed together with the mitochondrial oxidative phosphorylation pathway (such as COX6B1, NDUFA4) in calicoblast (Drake et al., 2020). This may suggest that photosynthetic products are directly supplied via lipid carriers to ATP required for calcification.

The role of multiple cell types in symbiosis

The maintenance of the coral-alga symbiosis depends on the coordinated action of multiple cell types (Levy et al., 2021). A lipid metabolism network comprising ALDH3A2, APOD and ACSL4 in the symbiotic cells supports nutrient exchange, while the expression of ATP6V1G1 maintains an acidic pH suitable for algal photosynthesis (Barott et al., 2015). Gastric epithelial cells play a key role in regulating nutrient transport and cell proliferation by upregulating ribosomal proteins and signaling molecules such as MST1R, YAP1 and JAG2 (Rosic et al., 2014).

Multicellular population synergies for environmental adaptation

We found high co-expression of mesoheat shock protein (HSP90AB1) in epidermal cells with oxidative stress genes (e.g. SOD2) in symbiotic cells (Cziesielski et al., 2018). This may suggest the activation of our transcellular antioxidant network. Meanwhile the high expression of thioredoxin (Trx) in calcified cells not only supports mitochondrial energy metabolism, but may also protect calcified frontiers from oxidative damage by scavenging reactive oxygen species (ROS) (Drake et al., 2020; Ramos-Silva et al., 2013; Guo et al., 2021).This multilayered defense mechanism may be the molecular basis of coral response to climate change (DeSalvo et al., 2010; Barshis et al., 2013).

Co-regulation of the calcification process

High expression of thioredoxin in calcified cells provides redox balance and energy support for biomineralization (Drake et al., 2020). Gastrodermal cells produce ATP through ETS-mediated aerobic respiration to provide energy support for the calcification process (Tambutté et al., 2011). Meanwhile, the expression of collagen proteins such as COL1A2, COL3A1, and COL5A1 in gastrodermal cells provides an extracellular matrix scaffold for the calcification process (Mass et al., 2014), ARHGEF17 and ACTB-mediated cytoskeleton promotes transport of calcified material (Zoccola et al., 2015).

Multi-layered immune defense system

The immune cells and the cnidocyte together form the defense line of the coral (Palmer, 2018; Sun et al., 2024). Immune cells express a variety of stress-related proteins that are involved in the basal defense (Rachamim et al., 2015). Cnidocyte regulate the synthesis of nematocystin through the expression of PDIA3 and RCN1 In terms of tissue defense, both immune cells and cnidocyte express specific functional genes such as FGF1 and IRF transcription factors in immune cells and ADAMTS9 and CMH1 in cnidocytes (Libro et al., 2013). The toxin and release mechanisms of cnidocyte exhibit dynamic regulatory mechanisms (Columbus-Shenkar et al., 2018).

Conclusions

In this study, we revealed the division of labor among multiple cell populations in corals in the maintenance of physiological symbiosis through single-cell transcriptome analysis. Epidermal, gastrodermal, and calicoblast collectively maintain coral symbiosis, calcification processes, and environmental adaptation through the synergy of metabolic networks, signaling pathways, and structural support. These findings provide new perspectives for understanding the cellular molecular basis of coral physiological functions, as well as a theoretical basis for predicting and improving coral adaptations to the environment.

Limitations and prospects

Although the present study reveals the functional collaboration network among coral cell populations, there are still some shortcomings: firstly, the single-cell transcriptome analysis reflects the gene expression profiles at a specific point in time, which makes it difficult to capture the spatial and temporal characteristics of the dynamic cell-cell interactions; secondly, many of the hypothesized intercellular signaling pathways still lack functional validation. In the future, we can: 1. analyze and elucidate the dynamics of cellular interactions by combining with time series; 2. validate the key signaling pathways by combining with in vitro culture and other experiments; 3. explore the plasticity of cellular collaborative networks in response to environmental stresses; 4. verify the spatial location of coral cells by in situ hybridization experiments.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Author contributions

WD: Data curation, Methodology, Software, Supervision, Validation, Writing – original draft. SZ: Data curation, Methodology, Software, Writing – original draft, Writing – review & editing. XZ: Methodology, Software, Visualization, Writing – review & editing. PW: Data curation, Visualization, Writing – review & editing. BD: Software, Visualization, Writing – review & editing. ZH: Data curation, Funding acquisition, Writing – review & editing. YL: Methodology, Software, Validation, Writing – original draft. HC: Conceptualization, Formal Analysis, Investigation, Project administration, Resources, Supervision, Writing – review & editing. ZL: Conceptualization, Formal Analysis, Funding acquisition, Investigation, Project administration, Resources, Supervision, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Science and Technology Project of Henan Province (242102210192, 242102210127, 242102210008).

Acknowledgments

This research is supported by the Big Data Computing Center of Southeast University and HongHe Artificial Intelligence Technology (Luoyang) Co.,Ltd., which provide computing power support.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2025.1551763/full#supplementary-material

References

Ashburner M., Ball C. A., Blake J. A., Botstein D., Butler H., Cherry J. M., et al. (2000). Gene ontology: Tool for the unification of biology. Nat. Genet. 25 (1), 25–29. doi: 10.1038/75556

PubMed Abstract | Crossref Full Text | Google Scholar

Barbosa M., Valentão P., Andrade P. B. (2016). Biologically active oxylipins from enzymatic and nonenzymatic routes in macroalgae. Mar. Drugs 14 (1), 23. doi: 10.3390/md14010023

PubMed Abstract | Crossref Full Text | Google Scholar

Barott K. L., Venn A. A., Perez S. O., Tambutteeé S., Tresguerres M., Somero G. N. (2015). Coral host cells acidify symbiotic algal microenvironment to promote photosynthesis. Proc. Natl. Acad. Sci. U.S.A. 112 (2), 607–612. doi: 10.1073/pnas.1413483112

PubMed Abstract | Crossref Full Text | Google Scholar

Barshis D. J., Ladner J. T., Oliver T. A., Seneca F. O., Traylor-Knowles N., Palumbi S. R. (2013). Genomic basis for coral resilience to climate change. Proc. Natl. Acad. Sci. U.S.A. 110 (4), 1387–1392. doi: 10.1073/pnas.1210224110

PubMed Abstract | Crossref Full Text | Google Scholar

Bay L. K., Nielsen H. B., Jarmer H., Seneca F., van Oppen M. J. H. (2009). Transcriptomic variation in a coral reveals pathways of clonal organization. Mar. Genomics 2 (2), 119–125. doi: 10.1016/j.margen.2009.07.004

PubMed Abstract | Crossref Full Text | Google Scholar

Beckmann A., Özbek S. (2012). The nematocyst: A molecular map of the cnidarian stinging organelle. Int. J. Dev. Biol. 56 (6-8), 577–582. doi: 10.1387/ijdb.113472ab

PubMed Abstract | Crossref Full Text | Google Scholar

Biel K. Y., Gates R. D., Muscatine L. (2007). Effects of free amino acids on the photosynthetic carbon metabolism of symbiotic dinoflagellates. Russian J. Plant Physiol. 54, 171–183. doi: 10.1134/S1021443707020033

Crossref Full Text | Google Scholar

Clayton W. S., Lasker H. R. (1984). Host feeding regime and zooxanthellal photosynthesis in the anemone, Aiptasia pallida (Verrill). Biol. Bull. 167 (3), 590–600. doi: 10.2307/1541412

PubMed Abstract | Crossref Full Text | Google Scholar

Coffey J. W., De Duve C. (1968). Digestive activity of lysosomes. I. The digestion of proteins by extracts of rat liver lysosomes. J. Biol. Chem. 243 (12), 3255–3263. doi: 10.1016/S0021-9258(18)93301-6

Crossref Full Text | Google Scholar

Columbus-Shenkar Y. Y., Sachkova M. Y., Macrander J., Fridrich A., Modepalli V., Reitzel A. M., et al. (2018). Dynamics of venom composition across a complex life cycle. Elife 7, e35014. doi: 10.7554/eLife.35014

PubMed Abstract | Crossref Full Text | Google Scholar

Costanza R., de Groot R., Sutton P., van der Ploeg S., Anderson S. J., Kubiszewski I., et al. (2014). Changes in the global value of ecosystem services. Global Environ. Change 26, 152–158. doi: 10.1016/j.gloenvcha.2014.04.002

Crossref Full Text | Google Scholar

Cziesielski M. J., Liew Y. J., Cui G., Schmidt-Roach S., Campana S., Marondedze C., et al. (2018). Multi-omics analysis of thermal stress response in a zooxanthellate cnidarian reveals the importance of associating with thermotolerant symbionts. Proc. R. Soc. B.: Biol. Sci. 285 (1877), 20172654. doi: 10.1098/rspb.2017.2654

PubMed Abstract | Crossref Full Text | Google Scholar

DeSalvo M. K., Sunagawa S., Voolstra C. R., Medina M. (2010). Transcriptomic responses to heat stress and bleaching in the elkhorn coral Acropora palmata. Mar. Ecol. Prog. Ser. 402, 97–113. doi: 10.3354/meps08372

Crossref Full Text | Google Scholar

Drake J. L., Mass T., Stolarski J., Von Euw S., van de Schootbrugge B., Falkowski P. G. (2020). How corals made rocks through the ages. Glob. Chang. Biol. 26 (1), 31–53. doi: 10.1111/gcb.14912

PubMed Abstract | Crossref Full Text | Google Scholar

Finn R. D., Coggill P., Eberhardt R. Y., Eddy S. R., Mistry J., Mitchell A. L., et al. (2016). The Pfam protein families database: Towards a more sustainable future. Nucleic Acids Res. 44 (D1), D279–D285. doi: 10.1093/nar/gkv1344

PubMed Abstract | Crossref Full Text | Google Scholar

Gates R. D., Hoegh-Guldberg O., Mcfall-Ngai M. J., Bil’ K. Y., Muscatine L. (1995). Free amino acids exhibit anthozoan “host factor” activity: They induce the release of photosynthate from symbiotic dinoflagellates in vitro. Proc. Natl. Acad. Sci. U.S.A. 92 (16), 7430–7434. doi: 10.1073/pnas.92.16.7430

PubMed Abstract | Crossref Full Text | Google Scholar

Guo Z., Liao X., Han T., Chen J., He C., Lu Z. (2021). Full-length transcriptomics reveal the gene expression profiles of reef-building coral pocillopora damicornis and symbiont zooxanthellae. Diversity. 13 (11), 543. doi: 10.3390/d13110543

Crossref Full Text | Google Scholar

Han T., Liao X., Zhu Y., Liu Y., Lu N., Li Y., et al. (2022). Full-length transcriptome maps of reef-building coral illuminate the molecular basis of calcification, symbiosis, and circadian genes. Int. J. Mol. Sci. 23 (19), 11135. doi: 10.3390/ijms231911135

PubMed Abstract | Crossref Full Text | Google Scholar

Hayward D. C., Hetherington S., Behm C. A., Grasso L. C., Forêt S., Miller D. J., et al. (2011). Differential gene expression at coral settlement and metamorphosis - A subtractive hybridization study. PloS One 6 (10), e26411. doi: 10.1371/journal.pone.0026411

PubMed Abstract | Crossref Full Text | Google Scholar

Huang C., Leng D., Sun S., Zhang X. D. (2019). Re-analysis of the coral Acropora digitifera transcriptome reveals a complex lncRNAs-mRNAs interaction network implicated in Symbiodinium infection. BMC Genomics 20, 1–15. doi: 10.1186/s12864-019-5429-3

PubMed Abstract | Crossref Full Text | Google Scholar

Hughes T. P., Barnes M. L., Bellwood D. R., Cinner J. E., Cumming G. S., Jackson J. B. C., et al. (2017). Coral reefs in the anthropocene. Nature 546 (7656), 82–90. doi: 10.1038/nature22901

PubMed Abstract | Crossref Full Text | Google Scholar

Jouiaei M., Sunagar K., Federman Gross A., Scheib H., Alewood P. F., Moran Y., et al. (2015). Evolution of an ancient venom: Recognition of a novel family of cnidarian toxins and the common evolutionary origin of sodium and potassium neurotoxins in sea anemone. Mol. Biol. Evol. 32 (6), 1598–1610. doi: 10.1093/molbev/msv050

PubMed Abstract | Crossref Full Text | Google Scholar

Kanehisa M., Goto S., Kawashima S., Okuno Y., Hattori M. (2004). The KEGG resource for deciphering the genome. Nucleic Acids Res. 32 (suppl_1), D277–D280. doi: 10.1093/nar/gkh063

PubMed Abstract | Crossref Full Text | Google Scholar

Levy S., Elek A., Grau-Bové X., Menéndez-Bravo S., Iglesias M., Tanay A., et al. (2021). A stony coral cell atlas illuminates the molecular and cellular basis of coral symbiosis, calcification, and immunity. Cell 184 (11), 2973–2938. doi: 10.1016/j.cell.2021.04.005

PubMed Abstract | Crossref Full Text | Google Scholar

Li D., Xiao L., Ge Y., Fu Y., Zhang W., Cao H., et al. (2018). High expression of Tob1 indicates poor survival outcome and promotes tumor progression via a Wnt positive feedback loop in colon cancer. Mol. Cancer 17, 1–6. doi: 10.1186/s12943-018-0907-9

PubMed Abstract | Crossref Full Text | Google Scholar

Liberzon A., Birger C., Thorvaldsdóttir H., Ghandi M., Mesirov J. P., Tamayo P. (2015). The molecular signatures database hallmark gene set collection. Cell Syst. 1 (6), 417–425. doi: 10.1016/j.cels.2015.12.004

PubMed Abstract | Crossref Full Text | Google Scholar

Libro S., Kaluziak S. T., Vollmer S. V. (2013). RNA-seq profiles of immune related genes in the staghorn coral Acropora cervicornis Infected with white band disease. PloS One 8 (11), e81821. doi: 10.1371/journal.pone.0081821

PubMed Abstract | Crossref Full Text | Google Scholar

Mansour T. A., Rosenthal J. J. C., Brown C. T., Roberson L. M. (2016). Transcriptome of the Caribbean stony coral Porites asteroides from three developmental stages. Gigascience 5 (1), s13742-016. doi: 10.1186/s13742-016-0138-1

PubMed Abstract | Crossref Full Text | Google Scholar

Mass T., Drake J. L., Peters E. C., Jiang W., Falkowski P. G. (2014). Immunolocalization of skeletal matrix proteins in tissue and mineral of the coral Stylophora pistillata. Proc. Natl. Acad. Sci. U.S.A. 111 (35), 12728–12733. doi: 10.1073/pnas.1408621111

PubMed Abstract | Crossref Full Text | Google Scholar

Mayfield A. B., Wang Y., Chen C. S., Lin C. Y., Chen S. H. (2014). Compartment-specific transcriptomics in a reef-building coral exposed to elevated temperatures. Mol. Ecol. 23 (23), 5816–5830. doi: 10.1111/mec.12982

PubMed Abstract | Crossref Full Text | Google Scholar

Moberg F., Folke C. (1999). Ecological goods and services of coral reef ecosystems. Ecol. Econ. 29 (2), 215–233. doi: 10.1016/S0921-8009(99)00009-9

Crossref Full Text | Google Scholar

Moya A., Huisman L., Ball E. E., Hayward D. C., Grasso L. C., Chua C. M., et al. (2012). Whole Transcriptome Analysis of the Coral Acropora millepora Reveals Complex Responses to CO 2-driven Acidification during the Initiation of Calcification. Mol. Ecol. 21, 2440–2454. doi: 10.1111/j.1365-294X.2012.05554.x

PubMed Abstract | Crossref Full Text | Google Scholar

Muscatine L. (1967). Glycerol excretion by symbiotic algae from corals and Tridacna and its control by the host. Sci. (1979). 156 (3774), 516–519. doi: 10.1126/science.156.3774.516

PubMed Abstract | Crossref Full Text | Google Scholar

Odum H. T., Odum E. P. (1955). Trophic structure and productivity of a windward coral reef community on eniwetok atoll. Ecol. Monogr. 25 (3), 291–320. doi: 10.2307/1943285

Crossref Full Text | Google Scholar

Özbek S. (2011). The cnidarian nematocyst: A miniature extracellular matrix within a secretory vesicle. Protoplasma 248 (4), 635–640. doi: 10.1007/s00709-010-0219-4

PubMed Abstract | Crossref Full Text | Google Scholar

Palmer C. V. (2018). Immunity and the coral crisis. Commun. Biol. 1 (1), 91. doi: 10.1038/s42003-018-0097-4

PubMed Abstract | Crossref Full Text | Google Scholar

Polato N. R., Vera J. C., Baums I. B. (2011). Gene discovery in the threatened elkhorn coral: 454 sequencing of the acropora palmata transcriptome. PloS One 6 (12), e28634. doi: 10.1371/journal.pone.0028634

PubMed Abstract | Crossref Full Text | Google Scholar

Rachamim T., Morgenstern D., Aharonovich D., Brekhman V., Lotan T., Sher D. (2015). The dynamically evolving nematocyst content of an Anthozoan, a scyphozoan, and a hydrozoan. Mol. Biol. Evol. 32 (3), 740–753. doi: 10.1093/molbev/msu335

PubMed Abstract | Crossref Full Text | Google Scholar

Rädecker N., Pogoreutz C., Gegner H. M., Cárdenas A., Roth F., Bougoure J., et al. (2021). Heat stress destabilizes symbiotic nutrient cycling in corals. Proc. Natl. Acad. Sci. U.S.A. 118 (5), e2022653118. doi: 10.1073/pnas.2022653118

PubMed Abstract | Crossref Full Text | Google Scholar

Ramos-Silva P., Kaandorp J., Huisman L., Marie B., Zanella-Cléon I., Guichard N., et al. (2013). The skeletal proteome of the coral acropora millepora: The evolution of calcification by co-option and domain shuffling. Mol. Biol. Evol. 30 (9), 2099–2112. doi: 10.1093/molbev/mst109

PubMed Abstract | Crossref Full Text | Google Scholar

Raz-Bahat M., Douek J., Moiseeva E., Peters E. C., Rinkevich B. (2017). The digestive system of the stony coral Stylophora pistillata. Cell Tissue Res. 368, 311–323. doi: 10.1007/s00441-016-2555-y

PubMed Abstract | Crossref Full Text | Google Scholar

Reimer J. D., Kise H., Wee H. B., Lee C. L., Soong K. (2019). Crown-of-thorns starfish outbreak at oceanic Dongsha Atoll in the northern South China Sea. Mar. Biodivers. 49, 2495–2497. doi: 10.1007/s12526-019-01021-2

Crossref Full Text | Google Scholar

Reyes-Bermudez A., DeSalvo M. K., Voolstra C. R., Sunagawa S., Szmant A. M., Iglesias-Prieto R., et al. (2009). Gene expression microarray analysis encompassing metamorphosis and the onset of calcification in the scleractinian coral Montastraea faveolata. Mar. Genomics 2 (3-4), 149–159. doi: 10.1016/j.margen.2009.07.002

PubMed Abstract | Crossref Full Text | Google Scholar

Rosic N., Kaniewska P., Chan C. K. K., Ling E. Y. S., Edwards D., Dove S., et al. (2014). Early transcriptional changes in the reef-building coral Acropora aspera in response to thermal and nutrient stress. BMC Genomics 15, 1–17. doi: 10.1186/1471-2164-15-1052

PubMed Abstract | Crossref Full Text | Google Scholar

Schmidt-Roach S., Miller K. J., Lundgren P., Andreakis N. (2014). With eyes wide open: A revision of species within and closely related to the Pocillopora damicornis species complex (Scleractinia; Pocilloporidae) using morphology and genetics. Zool. J. Linn. Soc. 170 (1), 1–33. doi: 10.1111/zoj.12092

Crossref Full Text | Google Scholar

Shinzato C., Shoguchi E., Kawashima T., Hamada M., Hisata K., Tanaka M., et al. (2011). Using the Acropora digitifera genome to understand coral responses to environmental change. Nature 476 (7360), 320–323. doi: 10.1038/nature10249

PubMed Abstract | Crossref Full Text | Google Scholar

Sun Y., Lan Y., Rädecker N., Sheng H., Diaz-Pulido G., Qian P. Y., et al. (2024). Gene expression of Pocillopora damicornis coral larvae in response to acidification and ocean warming. BMC Genomic Data 25 (1), 28. doi: 10.1186/s12863-024-01211-3

PubMed Abstract | Crossref Full Text | Google Scholar

Swain T. D., Vega-Perkins J. B., Oestreich W. K., Triebold C., DuBois E., Henss J., et al. (2016). Coral bleaching response index: a new tool to standardize and compare susceptibility to thermal bleaching. Glob. Chang. Biol. 22 (7), 2475–2488. doi: 10.1111/gcb.13276

PubMed Abstract | Crossref Full Text | Google Scholar

Tambutté S., Holcomb M., Ferrier-Pagès C., Reynaud S., Tambutté É., Zoccola D., et al. (2011). Coral biomineralization: From the gene to the environment. J. Exp. Mar. Biol. Ecol. 408 (1-2), 58–78. doi: 10.1016/j.jembe.2011.07.026

Crossref Full Text | Google Scholar

Tatusov R. L., Fedorova N. D., Jackson J. D., Jacobs A. R., Kiryutin B., Koonin E. V., et al. (2003). The COG database: an updated version includes eukaryotes. BMC Bioinf. 4, 1–14. doi: 10.1186/1471-2105-4-41

PubMed Abstract | Crossref Full Text | Google Scholar

Wilson E. B., Hilferty M. M. (1931). The distribution of chi-square. Proc. Natl. Acad. Sci. 17 (12), 684–688. doi: 10.1073/pnas.17.12.684

PubMed Abstract | Crossref Full Text | Google Scholar

Zoccola D., Ganot P., Bertucci A., Caminiti-Segonds N., Techer N., Voolstra C. R., et al. (2015). Bicarbonate transporters in corals point towards a key step in the evolution of cnidarian calcification. Sci. Rep. 5 (1), 9983. doi: 10.1038/srep09983

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: Pocillopora damicornis, ScRNA-seq, endosymbiosis, gastrodermis, epidermis, cnidocyte, calicoblast

Citation: Diao W, Zhang S, Zhu X, Wu P, Du B, Han Z, Liu Y, He C and Lu Z (2025) Single-cell transcriptome sequencing reveals the cell populations and functions associated with stony coral Pocillopora damicornis. Front. Mar. Sci. 12:1551763. doi: 10.3389/fmars.2025.1551763

Received: 26 December 2024; Accepted: 10 February 2025;
Published: 28 February 2025.

Edited by:

Ana Varela Coelho, Universidade Nova de Lisboa, Portugal

Reviewed by:

Pedro Martinez, University of Barcelona, Spain
Xiaopeng Yu, Guangxi University, China

Copyright © 2025 Diao, Zhang, Zhu, Wu, Du, Han, Liu, He and Lu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yunqing Liu, bGl1eXFAbGl0LmVkdS5jbg==; Chunpeng He, Y3BoZUBzZXUuZWR1LmNu; Zuhong Lu, emhsdUBzZXUuZWR1LmNu

These authors have contributed equally to this work and share first authorship

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Research integrity at Frontiers

Man ultramarathon runner in the mountains he trains at sunset

94% of researchers rate our articles as excellent or good

Learn more about the work of our research integrity team to safeguard the quality of each article we publish.


Find out more