Discovery of microsatellites and gene functional annotation based on the transcriptome of Oulastrea crispata

Lin Liuqing, Yang Zhao, Zhang Han, Zhang Die, Shi Tuo, Mao Shuoqian, Ding Peng, Zheng Xinqing, Zhou Xijie, He Congying

Journal of Marine Sciences ›› 2026, Vol. 44 ›› Issue (3) : 66-75.

PDF(2354 KB)
PDF(2354 KB)
Journal of Marine Sciences ›› 2026, Vol. 44 ›› Issue (3) : 66-75. DOI: 10.27049/j.cnki.1001-909X.2025089

Discovery of microsatellites and gene functional annotation based on the transcriptome of Oulastrea crispata

Author information +
History +

Abstract

Oulastrea crispata is a widely distributed reef-building coral with strong adaptability to various environmental stressors. However, existing studies lack systematic genetic background analysis of this species as a host organism to elucidate its underlying adaptive mechanisms. In this study, coral bleaching was induced by high-temperature treatment to obtain coral tissues with low symbiotic algal background. Transcriptome sequencing and bioinformatics analysis were then performed based on the Illumina HiSeqTM high-throughput sequencing platform. Following data quality assessment, EST-SSR loci were systematically identified, and EST-SSR molecular markers were developed. The results showed that the percentage of bases with Q30 scores (96.32%) of all samples met the high-quality sequencing standard. A total of 594 447 unigenes were assembled, with a transcript N50 length of 1 699 bp, indicating good transcriptome data quality. A total of 101 034 EST-SSR loci were identified from 87 979 unigenes, among which mono- and tri-nucleotide repeats were the dominant motif types, and the EST-SSR repeat number showed a decreasing trend with increasing motif length. Multi-database joint functional annotation revealed significant enrichment of EST-SSR loci in functional genes involved in environmental information processing, suggesting that EST-SSR marker-associated genes were closely linked to the core physiological mechanisms underlying coral environmental stress response. The EST-SSR marker system constructed in this study can serve as a high-quality candidate genetic marker library, provide not only important molecular tools for deciphering the environmental adaptation potential and evolutionary mechanisms of Oulastrea crispata populations, but also scientific references for coral population genetic analysis, cross-latitude adaptability research, and coral reef resource conservation and restoration.

Key words

Oulastrea crispata / EST-SSR / adaptation / genetic connectivity / coral reef restoration / functional annotation / transcriptome

Cite this article

Download Citations
Lin Liuqing , Yang Zhao , Zhang Han , et al . Discovery of microsatellites and gene functional annotation based on the transcriptome of Oulastrea crispata[J]. Journal of Marine Sciences. 2026, 44(3): 66-75 https://doi.org/10.27049/j.cnki.1001-909X.2025089

References

[1]
Birkeland C. Coral reefs in the Anthropocene[M]. Dordrecht: Springer, 2015: 1-15.
[2]
Hoegh-Guldberg O, Mumby P J, Hooten A J, et al. Coral reefs under rapid climate change and ocean acidification[J]. Science, 2007, 318(5857): 1737-1742.
Atmospheric carbon dioxide concentration is expected to exceed 500 parts per million and global temperatures to rise by at least 2 degrees C by 2050 to 2100, values that significantly exceed those of at least the past 420,000 years during which most extant marine organisms evolved. Under conditions expected in the 21st century, global warming and ocean acidification will compromise carbonate accretion, with corals becoming increasingly rare on reef systems. The result will be less diverse reef communities and carbonate reef structures that fail to be maintained. Climate change also exacerbates local stresses from declining water quality and overexploitation of key species, driving reefs increasingly toward the tipping point for functional collapse. This review presents future scenarios for coral reefs that predict increasingly serious consequences for reef-associated fisheries, tourism, coastal protection, and people. As the International Year of the Reef 2008 begins, scaled-up management intervention and decisive action on global emissions are required if the loss of coral-dominated ecosystems is to be avoided.
[3]
Caballero A, García-Dorado A. Allelic diversity and its implications for the rate of adaptation[J]. Genetics, 2013, 195(4): 1373-1384.
[4]
Chen C A, Lam K K, Nakano Y, et al. A stable association of the stress-tolerant zooxanthellae, Symbiodinium clade D, with the low-temperature-tolerant coral, Oulastrea crispata (Scleractinia: Faviidae) in subtropical non-reefal coral communities[J]. Zoological Studies, 2003, 42(4): 540-550.
[5]
Ng T Y, Ang P. Low symbiont diversity as a potential adaptive strategy in a marginal non-reefal environment: A case study of corals in Hong Kong[J]. Coral Reefs, 2016, 35(3): 941-957.
[6]
Alvarado-Cerón V, Muñiz-Castillo A I, León-Pech M G, et al. A decade of population genetics studies of scleractinian corals: A systematic review[J]. Marine Environmental Research, 2023, 183: 105781.
[7]
Huang W, Li M, Yu K F, et al. Genetic diversity and large-scale connectivity of the scleractinian coral Porites lutea in the South China Sea[J]. Coral Reefs, 2018, 37(4): 1259-1271.
[8]
Tóth G, Gáspári Z, Jurka J. Microsatellites in different eukaryotic genomes: Survey and analysis[J]. Genome Research, 2000, 10(7): 967-981.
We examined the abundance of microsatellites with repeated unit lengths of 1-6 base pairs in several eukaryotic taxonomic groups: primates, rodents, other mammals, nonmammalian vertebrates, arthropods, Caenorhabditis elegans, plants, yeast, and other fungi. Distribution of simple sequence repeats was compared between exons, introns, and intergenic regions. Tri- and hexanucleotide repeats prevail in protein-coding exons of all taxa, whereas the dependence of repeat abundance on the length of the repeated unit shows a very different pattern as well as taxon-specific variation in intergenic regions and introns. Although it is known that coding and noncoding regions differ significantly in their microsatellite distribution, in addition we could demonstrate characteristic differences between intergenic regions and introns. We observed striking relative abundance of (CCG)(n)*(CGG)(n) trinucleotide repeats in intergenic regions of all vertebrates, in contrast to the almost complete lack of this motif from introns. Taxon-specific variation could also be detected in the frequency distributions of simple sequence motifs. Our results suggest that strand-slippage theories alone are insufficient to explain microsatellite distribution in the genome as a whole. Other possible factors contributing to the observed divergence are discussed.
[9]
李福宇, 陈丹丹, 李元超, 等. 澄黄滨珊瑚微卫星标记的开发[J]. 基因组学与应用生物学, 2021, 40(S2): 2513-2521.
Li F Y, Chen D D, Li Y C, et al. Development of microsatellite markers for the hermatypic coral Porites lutea[J]. Genomics and Applied Biology, 2021, 40(S2): 2513-2521.
[10]
Wu Q, Huang W, Chen B, et al. Genetic structure of Turbinaria peltata in the northern South China Sea suggest insufficient genetic adaptability of relatively high-latitude scleractinian corals to environment stress[J]. Science of the Total Environment, 2021, 775: 145775.
[11]
Song Y, Jiang X, Zhang M, et al. Differences of EST-SSR and genomic-SSR markers in assessing genetic diversity in poplar[J]. Forestry Studies in China, 2012, 14(1): 1-7.
[12]
罗燕秋, 黄雯, 余克服, 等. 基于高通量测序的鹿角杯形珊瑚微卫星标记的开发[J]. 广西大学学报:自然科学版, 2020, 45(4): 958-965.
Luo Y Q, Huang W, Yu K F, et al. Development of microsatellite in Pocillopora damicornis using high-throughput sequencing[J]. Journal of Guangxi University: Natural Science Edition, 2020, 45(4): 958-965.
[13]
Dwivedi A, Suthar K P, Hamid R, et al. Exploitation of novel drought responsive EST-SSR markers in tetraploid cotton (Gossypium hirsutum L.)[J]. Gene Reports, 2025, 38: 102097.
[14]
盛岩, 郑蔚虹, 裴克全, 等. 微卫星标记在种群生物学研究中的应用[J]. 植物生态学报, 2002, 26(S1): 119-126.
Sheng Y, Zheng W H, Pei K Q, et al. Applications of microsatellites in population biology[J]. Chinese Journal of Plant Ecology, 2002, 26(S1): 119-126.
[15]
Jia S W, Li Y C, Chen S Q, et al. Microsatellite markers for Monitipora digitata designed using restriction-site associated DNA sequencing[J]. Frontiers in Marine Science, 2023, 10: 1019419.
[16]
王沛政, 盖如新, 夏利栋, 等. 两种软珊瑚EST-SSR通用性引物的研发[J]. 海南热带海洋学院学报, 2023, 30(2): 17-24.
Wang P Z, Gai R X, Xia L D, et al. Development and characterization of EST-SSR markers for Sinularia ceramensis and Sarcophyton ehrenbergi[J]. Journal of Hainan Tropical Ocean University, 2023, 30(2): 17-24.
[17]
Liu S Y V, Dai C F, Fan T Y, et al. Cloning and characterization of microsatellite loci in a Gorgonian coral, Junceella juncea (Anthozoa; Octocorallia; Ellisellidae) and its application in clonal genotyping[J]. Marine Biotechnology, 2005, 7(1): 26-32.
[18]
Wang S, Zhang L, Matz M. Microsatellite characterization and marker development from public EST and WGS databases in the reef-building coral Acropora millepora (Cnidaria, Anthozoa, Scleractinia)[J]. Journal of Heredity, 2009, 100(3): 329-337.
Mining for microsatellites (also called simple sequence repeats [SSRs]) in public sequence databases of a common Indo-Pacific coral Acropora millepora identified 191 SSRs from 10 258 expressed sequence tag (EST) and 618 SSRs from 14 625 whole-genome shotgun (WGS) sequences. In contrast to other animals, trinucleotide repeats, rather than dinucleotide repeats, are dominant in the WGS-SSRs, and AAT is the most frequent trinucleotide motif in EST-SSRs. We successfully developed 40 polymorphic markers from EST-SSRs and WGS-SSRs. Both EST- and WGS-SSRs show high levels of polymorphism within corals from the same reef patch. Interestingly, markers WGS079 and WGS227 revealed SSR duplications in a few individuals, suggesting recent duplication events. Genotypic linkage disequilibrium was identified in 5 pairs of SSR markers, which will be invaluable for high-resolution studies of genetic admixture in natural populations of A. millepora. Transferability analysis showed that 25 of these markers can be successfully amplified in one of the most ubiquitous Indo-Pacific corals Acropora hyacinthus. The marker collection reported here is the largest ever developed for any reef-building coral. It holds great potential for addressing coral reef connectivity across the Indo-Pacific with an unprecedented precision, especially taking into account the cross-species transferability of a substantial number of markers.
[19]
Studer R A, Dessailly B H, Orengo C A. Residue mutations and their impact on protein structure and function: Detecting beneficial and pathogenic changes[J]. The Biochemical Journal, 2013, 449(3): 581-594.
[20]
Dreisigacker S, Zhang P, Warburton M L, et al. SSR and pedigree analyses of genetic diversity among CIMMYT wheat lines targeted to different megaenvironments[J]. Crop Science, 2004, 44(2): 381-388.
Improved bread wheat (Triticum aestivum L.) cultivars for diverse agroecological environments are important for success in the effort to increase food production. In the 1980s, CIMMYT introduced the megaenvironment (ME) concept to breed wheats specifically adapted to different areas. Our objective was to analyze the genetic diversity among 68 advanced CIMMYT wheat lines targeted to different MEs by using 99 simple sequence repeats (SSRs) and the coefficient of parentage (COP). The average number of alleles detected was higher for the 47 genomic SSRs (5.4) than for the 52 SSRs derived from expressed sequence tags (EST) (3.3), but gene diversity between MEs was similar for both types of markers. No significant differences among the five MEs were observed for the means of SSR‐based genetic similarities (GS), calculated as 1 − Rogers' distance, and COP values. Both measures showed a low correlation (r = 0.43). High levels of genetic diversity were found within the germplasm targeted to each ME. However, principle coordinate analysis based on modified Rogers' distances did not separate the genotypes according to their targeted MEs. We conclude that presence of a single core germplasm can reflect large phenotypic differences. A sufficient number of diverse breeding lines for each ME is required because MEs generally combine various production areas. SSRs represent a powerful tool to quantify genetic diversity in wheat, but genotypic differentiation for adaptation to specific MEs in the CIMMYT program could not be proven.
PDF(2354 KB)

Accesses

Citation

Detail

Sections
Recommended

/

〈 〉