ORIGINAL RESEARCH article

Front. Mar. Sci., 27 September 2022

Sec. Marine Molecular Biology and Ecology

Volume 9 - 2022 | https://doi.org/10.3389/fmars.2022.991391

Genome of elegance coral Catalaphyllia jardinei (Euphylliidae)

  • 1. School of Life Sciences, Simon F.S. Li Marine Science Laboratory, State Key Laboratory of Agrobiotechnology, The Chinese University of Hong Kong, Hong Kong, Hong Kong SAR, China

  • 2. Dovetail Genomics, California, CA, United States

  • 3. School of Biological Science, The University of Hong Kong, Hong Kong, Hong Kong SAR, China

  • 4. Department of Biology, Queen’s University, Kingston, ON, Canada

  • 5. School of Life Sciences, The Chinese University of Hong Kong, Hong Kong, Hong Kong SAR, China

  • 6. Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou), Guangzhou, China

  • 7. Department of Ocean Science and Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou), The Hong Kong University of Science and Technology, Hong Kong, Hong Kong SAR, China

  • 8. Department of Biology, Hong Kong Baptist University, Hong Kong, Hong Kong SAR, China

  • 9. CAS and Shandong Province Key Laboratory of Experimental Marine Biology, Center for Ocean Mega-Science, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, China

Abstract

Coral reefs are under stress throughout the world. To better understand the molecular mechanisms underlying coral biology and their genomic evolution, here we sequenced the genome and transcriptomes of elegance coral Catalaphyllia jardinei (Euphylliidae). This monotypic genus stony coral is widespread but rare, being found across the Indo-West Pacific, from the northern Indian Ocean, Australia, Philippines, to the South China Sea. Due to its popularity among aquarium hobbyists, it is an overexploited species collected in large quantities from the wild for aquarium trade. The assembled genome is ~ 651.3 Mb in total length and of high physical contiguity with a scaffold N50 size of 28.9 Mb. The gene copy numbers of abiotic stress regulator (heat shock protein family genes) and neuropeptides (GLWamide, GRFamide, PRGamide and HIRamide) are similar to other sequenced anthozoans, and we have also identified the first set of sesquiterpenoid biosynthetic pathway genes in coral. Sequencing of small RNAs allows us to identify 35 microRNAs in C. jardinei and update the number of conserved microRNAs in cnidarians. This study established a foundation for further investigation into the roles of sesquiterpenoids and microRNAs in development of coral and understand their responses to climate change. Due to the easiness to culture C. jardinei in reef tanks and the established resources in this study, we propose this species be adopted as a new laboratory model in environmental and ecological experiments aiming to understand coral biology and responses to environmental stressors.

Introduction

The elegance coral Catalaphyllia jardinei was first described as Pectinia jardinei found in the shallow water on the reefs at the Warrior Reef (latitude 9° 45’ S), Thursday Island (latitude 10° 35’ S), and Albany Pass (latitude 10° 45’ S) (). Based on the morphological and color identification, the species was later determined to be a new genus under Euphylliidae (; ). Given that its skeleton is homomorphic to that of Euphyllia, and in recognition of Catala’s contribution to tropical marine fauna, Pectinia was changed to Catalaphyllia (). Today, the genus Catalaphyllia is considered monospecific with C. jardinei as its only member ().

Coral reef ecosystems are important in providing a favorable habitat for its associated organisms (). Throughout the world coral reefs are under stresses that include habitat destruction (), overexploitation (; ), pollution (; ), diseases (), predation (), and climate change (; ; ; ; ; ). The reef-building coral C. jardinei has a broad geographical distribution and can be found both in the Indian Ocean and Pacific Ocean (), including Australia, Indonesia, Japan, Maldives, Malaysia, Papua New Guinea, Philippines, Seychelles and Vietnam (). Nevertheless, given its small population size and overexploitation, it is considered by the IUCN Red List as having a “Vulnerable” status ().

Since the late 1980s, the popularity of keeping live corals at home has increased rapidly (). Large-scale culture of C. jardinei is a current business in Indonesia and corals are exported globally (). In addition to its appearance as described in its common name, C. jardinei is a popular coral species due to its easiness to culture ex situ. Naturally found in deep flats (greater than 35 meters) with sandy/silty bottoms, lagoons, and nearshore mudflats (), C. jardinei can tolerate and survive in fluctuating or less favorable environmental conditions (; http://www.edgeofexistence.org/species/elegance-coral/), making it a potential new model in the laboratory to study coral biology and their responses to environmental stressors.

To facilitate the use of C. jardinei to study the molecular mechanisms underlying coral biology, this study sequenced and provided the first genome and transcriptomic resources of this coral species.

Materials and methods

Sample collection and genome sequencing

Adult Catalaphyllia jardinei (Euphylliidae) (Figure 1A) were purchased from mainland China (original source from Indonesia) and Hong Kong suppliers (original source from Australia), and were cultured in an aquarium tank with circulating artificial seawater (salinity: 30-35 ppt) at around 23 ± 1°C and other equipment to simulate pristine waters at The Chinese University of Hong Kong. Corals were allowed to acclimatize under the laboratory conditions for 15 days prior to be used for DNA and RNA extraction.

Figure 1

Sample for genome sequencing originated from tentacles from a single individual (Figure 1A). Genomic DNA (gDNA) was extracted using the QIAamp DNA minikit (Qiagen, Hilden, Germany) following the manufacturer’s instructions. Extracted gDNA was subjected to quality control using a Nanodrop spectrophotometer (Thermo Scientific) and gel electrophoresis. Qualified samples were sent to Novogene and Dovetail Genomics for library preparation and sequencing. The resulting library was sequenced on an Illumina HiSeq X platform to produce 2 × 150 paired-end sequences. The length-weighted mean molecule length was 17.9 Kb, and details of the sequencing data were highlighted in Supplementary Information S1.

Dovetail Omni-C library preparation and sequencing

For each Dovetail Omni-C library, chromatin was fixed with formaldehyde and extracted. Fixed chromatin was digested with DNase I, and chromatin ends were repaired and ligated to a biotinylated bridge adapter followed by proximity ligation of adapter containing ends. After proximity ligation, crosslinks were reversed and the DNA was purified. Purified DNA was treated to remove biotin that was not internal to ligated fragments. Sequencing libraries were generated using NEBNext Ultra enzymes and Illumina-compatible adapters. Biotin-containing fragments were isolated using streptavidin beads before PCR enrichment of each library. The library was sequenced on an Illumina HiSeqX platform to produce 196.9 millions of 150 bp reads, and the raw data were listed in Supplementary Information S1.

Transcriptome sequencing

Total RNA, including messenger RNA and small RNA, from tentacles from single individual was isolated using mirVana™ miRNA Isolation Kit (Ambion) according to the manufacturer’s instructions, and subjected to quality control using a Nanodrop spectrophotometer (Thermo Scientific), gel electrophoresis, and Agilent 2100 Bioanalyzer (Agilent RNA 6000 Nano Kit). Qualified samples underwent library construction and 150 bp paired-end (PE) sequencing at Novogene; polyA-selected RNA-Sequencing libraries were prepared using TruSeq RNA Sample Prep Kit v2. Insert sizes and library concentrations of final libraries were determined using an Agilent 2100 bioanalyzer instrument (Agilent DNA 1000 Reagents) and real-time quantitative PCR (TaqMan Probe) respectively. Small RNA was submitted to BGI Hong Kong for HiSeq Small RNA library construction and 50 bp single-end (SE) sequencing. Details of the sequencing data can be found in Supplementary Information S1.

Genome assembly

K-mers of the HiseqXten 150PE (Chromium WGS) were counted using jellyfish version 2.2.5 () with k-mers = 21, and estimation of genome size, repeat content, and heterozygosity were analyzed based on a k-mer-based statistical approach in the GenomeScope2 webtool (, Supplementary Information S3). Kraken2 () (k2_standard_20210517 (kraken2 standard database) from https://benlangmead.github.io/aws-indexes/k2, which contains archaea, bacteria, viral, plasmid, human1, UniVec_Core) was used to mark and remove the contamination from the gDNA and transcriptome reads, 3.75% and 0.29% of contaminated reads were removed respectively. Processed Chromium WGS reads were used to construct a de novo assembly using Supernova (v 2.1.1) with default parameters (raw coverage = 78.86x), which generates phased, whole-genome de novo assemblies from a Chromium-prepared library () (https://support.10xgenomics.com/de-novo-assembly/software/pipelines/latest/using/running). The supernova output pseudohap assembly were used to generate the dedupe assembly by dedupe.sh (version 19Feb, 2020) of BBMap (Bushnell B. - sourceforge.net/projects/bbmap/), and then scaffolded using Dovetail Omni-C library reads by HiRise, a software pipeline designed specifically for using proximity ligation data to scaffold genome assemblies (). Dovetail Omni-C library sequences were aligned to the draft input assembly using bwa (https://github.com/lh3/bwa). The separations of Dovetail Omni-C read pairs mapped within draft scaffolds were analyzed by HiRise to produce a likelihood model for genomic distance between read pairs, and the model was used to identify and break putative misjoins, to score prospective joins, and make joins above a threshold. To further check for potential symbiotic contamination, 11 symbiodinium genomes were downloaded from NCBI (Supplementary Information S2) and BLASTN () were carried out on the assembled genomes (e-value 1e-3). Scaffolds with alignment lengths >= 500 bp were further searched on nt database with BLASTN, and no symbiotic contamination was found in the final assembly genome.

Gene model prediction

Gene models were predicted as described in the jellyfish genomes (). Briefly, the gene models were trained and predicted using funannotate (v1.8.9, https://github.com/nextgenusfs/funannotate) () with the following parameters: “–protein_evidence uniprot_sprot.fasta –genemark_mode ET –busco_seed_species fly –optimize_augustus –busco_db arthropoda –organism other –max_intronlen 350000”. In the “funannotate predict” step, GeneMark-ES v4.32 () was used for ab initio gene prediction and Augustus () was trained using BUSCO. The gene models from the prediction sources “Augustus”, high-quality Augustus predictions (HiQ) with exon evidence, “pasa” ()_, “snap” (), “GlimmerHMM” (), and “GeneMark” were passed to Evidence Modeler (EVM Weights: {‘GeneMark’: 1, ‘HiQ’: 2, ‘pasa’: 6, ‘proteins’: 1, ‘Augustus’: 1, ‘GlimmerHMM’: 1, ‘snap’: 1, ‘transcripts’: 1}) to generate the final annotation files. Finally, followed by PASA to update the EVM consensus predictions, and add UTR annotations and models for alternatively spliced isoforms. Protein-coding genes were searched with BLASTP against the nr and SwissProt databases by diamond (v0.9.24) () with parameters “–more-sensitive –evalue 1e-3”, and mapped by HISAT2 (version 2.1.0) () with transcriptome reads. Gene models with no similarity to any known proteins and no mRNA support were removed from the final version. The gene models were screened by CGAT () and only the longest genes were retained.

Synteny analysis

Orthologues and orthogroups in C. jardinei and other 33 cnidarians proteomes were identified and inferred using OrthoFinder v2.5.2 () with “-M msa” activated and default values to other parameters. To symbolise the gene families, the longest protein of each gene was taken as the representative in orthologues analysis. The species tree was constructed with fasttree using STAG algorithm and Hydra vulgaris was identified as the best outgroup species to construct the rooted tree using STRIDE method by OrthoFinder (; ). Tree supporting value were the proportion of species tree derived from gene tree supporting each bipartition. The species tree was then converted into divergence tree using r8s (). The tree structure and annotations were visualised using R package ‘ggtree’ (; ) in R v4.1.2 (). Orthologues anchored by mutual best Diamond blastp hits (e-value cut-off 0.001) between C. jardinei and other 33 cnidarians as mentioned above were used in macrosynteny analyses. Oxford synteny plot were generated following previously described methods (; ). Detailed analytical scripts and materials were distributed at https://github.com/xieyichun50/coral_genome_Catalaphyllia_jardinei.

microRNA annotation

Adaptor sequences were trimmed from small RNA sequencing reads and Phred quality score less than 20 were removed. Processed reads of length within 18 bp and 27 bp were then mapped to the genomes using mapper.pl module of the mirDeep2 package (). To identify known miRNA, the predicted coral microRNA hairpins were compared against metazoan microRNA precursor sequences from miRBase () using BLASTN (e-value <1e-2). For miRNAs with no significant sequence similarity to any of the miRNAs in miRBase were then checked manually. Novel microRNAs were defined when they fulfilled the criteria of microRNAs (). The expression levels of different arms of a miRNA were calculated based on the number of sequencing reads mapped to the respective arm region in the predicted miRNA hairpin by bowtie. In addition to the above, the precursor sequences of microRNAs from other corals in previous studies (; ) were also used to carry out BLASTN searches for identification of any missed microRNA annotations.

Annotation of sesquiterpenoid pathway and neuropeptide genes

17 gene sequences involved in the sesquiterpenoid biosynthetic pathway retrieved from Nematostella vectensis and 4 gene sequences of neuropeptide (GLWamide, GRFamide, PRGamide and HIRamide) sequences from Hydra vulgaris were used for searching against the gene models and genome using BLASTP and TBLASTN, respectively. Putatively identified orthologues with a threshold of E-value equals to 10E-3 were tested by reciprocal searches in the NCBI nr database using BLASTP. For neuropeptide preprohormones, amino acid sequences were translated manually and aligned with MEGA 7.0 (), while the signaling sequences were analysed with SignalP 3.0 (). The cleavage sites were predicted by software ProP 1.0 () and checked manually accordingly to the criteria suggested by . The potential mature neuropeptides were aligned using MEGA 7.0 ().

Results and discussion

High quality genome

Here, we presented a high-quality chromosome-level genome assembly of C. jardinei (2n = 28) (; ; ), the assembled genome size was 651 Mb, which is compatible with the kmer 21 estimate of a genome size of 685.5 Mb (Figure 1, Supplementary Information S3). 69.64% of the genomic sequences were contained on 14 pseudomolecules (Supplementary Information S1.3). The BUSCO score for complete genes is 88.05% (BUSCO version 5.1.3, metazoa_odb10), with complete and single-copy BUSCOs (S) being 84.28%, complete and duplicated BUSCOs (D) being 3.77%, fragmented BUSCOs (F) being 6.60% and Missing BUSCOs (M) being 5.35% (Figure 1C). Its scaffold N50 size (28.9Mb) represents the highest physical contiguity among the published coral genomes (Table 1). Using the transcriptomes, we predicted a total of 44,407 gene models including 37,003 protein coding genes and 7,404 tRNA genes, in the C. jardinei genome, which is within the range of 21,369 to 39,160 for other published coral genomes (Table 1).

Table 1

SpeciesAssembled genome size (bp)N50Number of ProteinsSum of Amino AcidsSum of Exons (bp)Sum of Introns (bp)Sum of gene region (bp)References
Acropora acuminata397,781,7771,004,88822,30611,844,58566,338,114234,162,220234,529,491
Acropora awi430,482,9021,088,40722,65311,755,03065,880,915249,482,354241,506,016
Acropora cytherea430,309,1691,083,94623,36312,256,65767,552,835247,246,641245,571,636
Acropora digitifera415,842,4891,856,31222,32610,509,35962,321,781215,573,311233,203,300
Acropora echinata409,669,9001,812,70121,94311,333,55165,943,777246,278,370236,046,259
Acropora florida442,860,187751,37623,85712,087,36864,685,340240,937,980244,285,300
Acropora gemmifera405,217,1651,134,58122,24711,698,57665,230,873237,013,391233,312,317
Acropora hyacinthus452,727,5921,562,59223,14712,513,95169,389,302254,370,978250,891,995
Acropora intermedia416,883,045577,31223,34311,998,37764,476,112232,793,965237,988,598
Acropora microphthalma387,940,0751,050,19622,61811,614,44365,390,591234,712,549235,057,382
Acropora millepora475,381,25319,840,54328,18812,836,11749,237,008158,255,732207,492,740
Acropora muricata420,750,014574,62723,64612,079,124237,931,003242,349,835
Acropora nasuta419,444,1681,045,28923,31912,115,88467,461,458248,696,816245,388,094
Acropora selago393,159,869657,17223,11511,696,41665,591,025233,912,084234,017,974
Acropora tenuis406,871,8331,160,22023,11811,917,12968,612,749192,068,808237,497,938
Acropora yongei443,890,7703,019,48223,50012,103,78367,790,062248,352,452249,216,954
Astreopora myriophthalma381,883,8221,590,75128,71113,287,97642,124,462116,058,587155,181,959
Catalaphyllia jardinei651,363,12328,947,99937,00313,815,76659,489,494221,415,763282,642,405This study
Dendronephthya gigantea286,150,6281,445,52328,74117,717,91874,944,700345,593,618147,908,393
Fungia sp.606,319,501323,14938,20915,470,32156,333,027194,127,360250,460,387
Galaxea fascicularis334,165,88087,93322,4189,704,88342,730,289132,589,881175,320,170
Goniastrea aspera764,857,003518,94935,90115,491,30651,803,637238,014,146289,817,783
Montipora cactus652,728,006898,51121,98213,478,39950,412,387311,983,610337,284,991
Montipora capitata614,509,607185,53736,69113,448,12740,454,043222,885,069148,531,395
Montipora efflorescens643,312,9411,132,31621,36913,130,86649,077,197308,033,328331,889,023
Orbicella faveolata485,548,9391,162,44632,58716,554,99595,196,829556,036,559266,375,485
Pachyseris speciosa984,430,632766,64939,16017,175,30772,350,345264,758,179337,108,524
Pocillopora damicornis234,350,878326,13325,18313,735,37963,957,019274,910,693151,175,419
Pocillopora verrucosa380,505,698333,69627,43917,115,89066,501,999140,069,492208,011,644
Porites lutea552,020,673660,70831,12614,351,12356,249,572199,362,641255,612,213
Stylophora pistillata400,120,318457,45333,25218,275,23886,199,075513,138,303223,994,894

Statistics of published coral genomes.

Underline, the statistics of C. jardinei genome

Synteny analysis

To clarify the syntenic relationship and gene linkage blocks between C. jardinei and other cnidarians, we generated the Oxford synteny plots of species pairs using the orthologues anchored by mutual best BLAST hit (e-value 1e-3) (Figure 2, Supplementary Information S4). Global gene orthogroup inference identified 26,260 orthogroups. The phylogenetic tree was constructed using 1,029 single-copy orthogroups presented in 31 out of 34 species, where single-copy indicated only one gene was assigned to specific orthogroup for each species. The chromosomal organisations are conserved in most soft corals, and each pseudo-chromosome in C. jardinei was associated with a different chromosome in another soft coral Acropora millipora. Further comparisons of C. jardinei against the more distant species, such as the sea anemone (Nematostella vectensis) and the true jellyfish (Rhopilema esculentum), displayed the one-to-one relationship on most synteny blocks, and a few one-to-many blocks (Supplementary Figures 13). Evidence on breakage, fusion or large-scale duplication was rarely found in these investigated cnidarians, suggesting that the genomes of these cnidarians have not undergone massive evolutionary changes.

Figure 2

Homeobox genes

Homeobox genes are important developmental genes in animals and have been used as markers to understand genomic changes in animal evolution (e.g. , , ; ; ; ; ). A total of 93 ANTP-class homeobox genes have been identified in the C. jardinei genome, including clusters of Hox-like and NK-like genes, and a ParaHox gene cluster (Figure 3; Supplementary Information S5). The data here further strengthens the idea that if distinct Hox and NK genes clusters are associated with ectodermal and mesodermal patterning, it could be a bilaterian-lineage specific event (; ), and the three gene ParaHox gene cluster could only be retained in the medusozoans than anthozoans among the extant cnidarians.

Figure 3

Heat shock protein family genes

Corals experiences abiotic stresses such as temperature, salinity, and oxygen concentration, and heat shock proteins (HSP) can assist and reconstruct the structure of damaged proteins under stresses (; ; ; Supplementary Figure 4). In general, HSP can be divided into six subfamilies, including HSPE (GroEL, cpn10), small heat shock proteins, HSP40, HSP70, HSP90, and HSP110, and HSP70 has been shown as chaperone proteins participated in response to heat stress in reef-building scleractinian corals (), (). In C. jardinei, we identified members of six HSP subfamilies, and their numbers were similar to those identified in the stony coral Acropora millepora chromosomal-level genome (Supplementary Figure 5, Supplementary Information S6). The identified heat shock protein gene families will be useful for future investigation of their roles in C. jardinei under different environmental conditions.

Neuropeptide genes

Neuropeptides play essential roles as endocrine factors and neurotransmitters, and are involved in the cellular differentiation and development of cnidarians (). Here, Scleractinia conserved neuropeptides GLWamide, GRFamide, PRGamide, HIRamide, RWamide, RPamide and pQITRFamide were identified in C. jardinei (Figures 4A, B, Supplementary Information S7). A total of 6 mature peptides were identified from the preprohormone sequences of GLWamide in C. jardinei, and sequence alignment showed a conserved C-terminal GXW motif contained in these sequences (Figures 4C). Two of the GLWamides had “-GIW” replacement and such mutations were also identified in hard coral Euphyllia ancora () (Figure 4C, Supplementary Information S8). For GRFamide preprohormones, 16 mature peptides were predicted. Among the identified PRGamide genes in C. jardinei, a single preprohormone gene was found to yield 46 different mature peptides (Figure 4B, Supplementary Information S8). For HIRamides, it was found that there was no glutamine residue identified from the predicted HIRamides in C. jardinei, where this phenomenon is also observable in other cnidarian HIRamides (; Supplementary Information S8). For RWamide and RPamide, only 1 and 2 mature peptide(s) are yielding from the precursor respectively. Meanwhile, the precursor of pQITRFamide could be identified from the genome, but there was no expression shown in the sequenced transcriptomes, thus the number and form of mature peptides could not be identified (Figures 4B, C). The identified neuropeptides will be useful for future investigation of their developmental roles in C. jardinei.

Figure 4

Sesquiterpenoid biosynthetic pathway genes

Sesquiterpenoid hormones, such as juvenile hormone, regulate the development, metamorphosis, and reproduction in insects, and have long been thought to be confined to the arthropods (e.g. ; ). A recent study has unexpectedly identified the sesquiterpenoid hormonal system in jellyfishes, suggesting that the sesquiterpenoid hormone system is established at the bilaterian-cnidarian ancestor (). To test whether coral species contain a sesquiterpenoid hormone system, we searched the sesquiterpenoid biosynthetic pathway genes in C. jardinei genome (Figure 4D) and identified the mevalonate pathway, aldehyde dehydrogenase (ALDH), farnesyl transferase (FNT), Ste 24 endopeptidase (ZMPSTE24), prenyl protein peptidase (RCE1), isoprenylcysteine carboxymethyl transferase (ICMT) and prenylcysteine oxidase (PCYOXIL) in a similar manner as the jellyfishes (Figure 4D, Supplementary Information S9). This result suggests that coral also has the sesquiterpenoid hormone system, and the functional roles of these hormones remain to be seriously tested in different cnidarians.

microRNA gain and loss in corals

Unlike the situation in bilaterians which contain lots of conserved microRNAs across different lineages, the cnidarians and bilaterians are well known to share only one conserved microRNA miR-100 (). Information of small RNA sequences among different cnidarians reveals limited number of conserved microRNAs among different cnidarian lineages (; ). Here, we sequenced the small RNA from C. jardinei, and annotated a total of 35 microRNAs in this species (Figure 5, Supplementary Information S10). In comparison with other published cnidarian microRNAs, in congruent to previous studies (e.g. ), there were more lineage-specific microRNAs than conserved microRNAs among different cnidarian lineages. We also identified conserved microRNAs at different key nodes in the cnidarian phylogeny, including miR-2022 and miR-2030 that are conserved across all cnidarians, six other microRNAs (miR-2023, miR-2025, miR-2036, miR-2037, miR-2050, miR-9425) that are conserved across all anthozoans, and 4 and 1 other microRNAs that are conserved in the Actiniaria and Scleractinia, respectively (Figure 5). Our findings established the foundation for further investigation of the different roles of these conserved microRNAs in cnidarians.

Figure 5

Conclusion

In the present study, we provided a high-quality genome assembly and transcriptome resources of the elegance coral Catalaphyllia jardinei (Euphylliidae). Abiotic stress regulator, hormonal gene families, and microRNAs were annotated, which will serve as a foundation for further investigation into their roles in development and responses to climate change. Due to the easiness to culture C. jardinei in reef tanks and the established resources in this study, we propose this species to be adopted as a new laboratory model.

Funding

This work was supported by the TUYF Charitable Trust, Hong Kong Research Grant Council Collaborative Research Grant (C4015-20EF, C7013-19G), General Research Fund (14100919, 14100420), the Open Collaborative Research Fund and Operation Fund from the Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou) (HKB L20200008), Center for Ocean Research in Hong Kong and Macau (CORE) (CORE is a joint research center for ocean research between QNLM and HKUST), and the Direct Grant of The Chinese University of Hong Kong (4053489, 4053547). YY was supported by the PhD studentship of The Chinese University of Hong Kong.

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.

Statements

Data availability statement

The final chromosome assembly was submitted to NCBI Assembly under accessionnumber JAEMPF000000000 in NCBI. The raw reads generated in this study havebeen deposited to the NCBI database under the BioProject accessions:PRJNA687327, the genome annotation files were deposited in the Figshare (https://Emms, D. M., & Kelly, S. (2018). STAG:species tree inference from all genes. BioRxiv, 267914.doi.org/10.6084/m9.figshare.19619700).

Author contributions

JHLH conceived and supervised the study. YY performed the extraction of DNA and RNA, and analyses of homeobox gene, heat shock protein, and microRNA. WN, JH, and TS performed genome assembly. WN performed gene model prediction, and transcriptomic analyses of mRNA and sRNA. WLS performed the neuropeptide analyses. YY and WLS performed the sesquiterpenoid analyses. YX performed synteny analyses. HYY and YY maintained the animal culture. DMB, WGB, TFC, APYC, KFL, PYQ, JWQ, BT, FX, JHLH obtained the funding. YY, WN, WLS, YX, JHLH wrote the first draft of the manuscript. All authors approve the final version of the manuscript.

Conflict of interest

JHa and TS are employees of Dovetail Genomics, but have no roles in experimental design.

The remaining 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.

Supplementary material

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

Supplementary Information 1

Genome and transcriptome sequencing data statistics.

Supplementary Information 2

Symbiodinium genomes.

Supplementary Information 3

Reports of GenomeScope2.

Supplementary Information 4

Species used in orthologues searches.

Supplementary Information 5

ANTP-class homeobox genes.

Supplementary Information 6

Heat shock proteins.

Supplementary Information 7

Anthozoan conserved neuropeptides

Supplementary Information 8

Mature neuropeptides.

Supplementary Information 9

Sesquiterpenoid biosynthetic pathway genes.

Supplementary Information 10

MicroRNAs.

References

Summary

Keywords

coral, genome, Catalaphyllia, microRNA, sesquiterpenoid

Citation

Yu Y, Nong W, So WL, Xie Y, Yip HY, Haimovitz J, Swale T, Baker DM, Bendena WG, Chan TF, Chui APY, Lau KF, Qian P-Y, Qiu J-W, Thibodeau B, Xu F and Hui JHL (2022) Genome of elegance coral Catalaphyllia jardinei (Euphylliidae). Front. Mar. Sci. 9:991391. doi: 10.3389/fmars.2022.991391

Received

11 July 2022

Accepted

29 August 2022

Published

27 September 2022

Volume

9 - 2022

Edited by

Yong Wang, Tsinghua University, China

Reviewed by

Thomas Lund Koch, University of Copenhagen, Denmark; Miriam Ruocco, Anton Dohrn Zoological Station, Italy

Updates

Copyright

*Correspondence: Jerome H. L. Hui,

†These authors have contributed equally to this work

This article was submitted to Marine Molecular Biology and Ecology, a section of the journal Frontiers in Marine Science

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics