in

The de novo transcriptome of the freshwater copepod Cyclops abyssorum tatricus reveals high-elevation adaptation


Abstract

Copepods, small aquatic crustaceans, are one of the most abundant zooplankton in the world. These animals play a critical ecological role in aquatic ecosystems such as oceans, streams, or, as in this study, alpine lakes. In these ecosystems, copepods have adapted to stressful and dynamic environments, a process which can be investigated via comparative transcriptomics. An assembled transcriptome is a pre-requisite to applying transcriptomic tools in physiological research, but there are few transcriptome assemblies available for copepods. To address this gap, we assembled a de novo transcriptome of Cyclops abyssorum tatricus by merging Pacific Bioscience long reads from copepods collected in two lakes at two different time points. The final assembly consisted of 52,521 contigs with a BUSCO score of 80.7%. We annotated a total of 26,255 (49.99%) protein sequences using the eggNOG database. Gene ontology analyses revealed that most gene annotations were involved in cellular processes and signaling (34.61%). Comparisons with two other copepod species showed that the transcriptome assembly of C. abyssorum tatricus is enriched for cold acclimation genes, consistent with its long-term adaptation to cold water environments. This de novo transcriptome will enable comparative transcriptomic studies in this species, allowing us to investigate physiological adaptations to alpine environments.

Data availability

All elements of this study are available to the general public via the NCBI BioProject PRJNA1377012 database (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA1377012/). This unannotated transcriptome shotgun assembly project has been deposited at DDBJ/EMBL/GenBank under the accession number GLLG00000000. The version described in this paper is the first version, GLLG01000000, while the two raw PacBio datasets are accessible via the accession numbers SAMN53755002 (run number: SRR36346890 in .bam or SRR37062719 in .fastq) and SAMN53755001 (run number: SRR36346891 in .bam or SRR37062720 in .fastq). No custom code has been used in the creation of the dataset for this article.

References

  1. Sommer, U. et al. Beyond the plankton ecology group (PEG) model: mechanisms driving plankton succession. Annu. Rev. Ecol. Evol. Syst. 43, 429–448. https://doi.org/10.1146/annurev-ecolsys-110411-160251 (2012).

    Google Scholar 

  2. Whitman, D. & Ananthakrishnan, T. Phenotypic plasticity of insects. Phenotypic Plast. insects: Mech. consequences. 1–63. https://doi.org/10.1201/b10201 (2009).

  3. Lauritano, C., Procaccini, G. & Ianora, A. Gene expression patterns and stress response in marine copepods. Mar. Environ. Res. 76, 22–31. https://doi.org/10.1016/j.marenvres.2011.09.015 (2012).

    Google Scholar 

  4. DeBiasse, M. B. & Kelly, M. W. Plastic and evolved responses to global change: what can we learn from comparative transcriptomics? J. Hered. 107, 71–81. https://doi.org/10.1093/jhered/esv073 (2016).

    Google Scholar 

  5. Adrian, R. et al. Lakes as sentinels of climate change. Limnol. Oceanogr. 54, 2283–2297. https://doi.org/10.4319/lo.2009.54.6_part_2.2283 (2009).

    Google Scholar 

  6. Tartarotti, B. et al. Distribution and UV protection strategies of zooplankton in clear and glacier-fed alpine lakes. Sci. Rep. 7, 4487. https://doi.org/10.1038/s41598-017-04836-w (2017).

    Google Scholar 

  7. Conesa, A. et al. A survey of best practices for RNA-seq data analysis. Genome Biol. 17, 13. https://doi.org/10.1186/s13059-016-0881-8 (2016).

    Google Scholar 

  8. Kelly, M., DeBiasse, M. B., Villela, V., Roberts, H. & Cecola, C. F. Adaptation to climate change: trade-offs among responses to multiple stressors in an intertidal crustacean. Evol. Appl. 9, 1147–1155. https://doi.org/10.1111/eva.12394 (2016).

    Google Scholar 

  9. Schoville, S. D., Barreto, F. S., Moy, G. W., Wolff, A. & Burton, R. S. Investigating the molecular basis of local adaptation to thermal stress: population differences in gene expression across the transcriptome of the copepod Tigriopus californicus. BMC Evol. Biol. 12, 170 (2012).

    Google Scholar 

  10. Roncalli, V., Cieslak, M. C. & Lenz, P. H. Transcriptomic responses of the calanoid copepod Calanus finmarchicus to the saxitoxin producing dinoflagellate Alexandrium fundyense. Sci. Rep. 6, 25708. https://doi.org/10.1038/srep25708 (2016).

    Google Scholar 

  11. Bron, J. E. et al. Observing copepods through a genomic lens. Front. Zool. 8, 22. https://doi.org/10.1186/1742-9994-8-22 (2011).

    Google Scholar 

  12. Jørgensen, T. S. et al. The whole genome sequence and mRNA transcriptome of the tropical cyclopoid copepod Apocyclops royi. G3 9, 1295–1302; (2019). https://doi.org/10.1534/g3.119.400085

  13. Kelly, M. W., Pankey, S. M., DeBiasse, M. B. & Plachetzki, D. C. Adaptation to heat stress reduces phenotypic and transcriptional plasticity in a marine copepod. Funct. Ecol. 31, 398–406. https://doi.org/10.1111/1365-2435.12725 (2017).

    Google Scholar 

  14. Barreto, F. S. et al. Genomic signatures of mitonuclear coevolution across populations of Tigriopus californicus. Nat. Ecol. Evol. 2, 1250–1257. https://doi.org/10.1038/s41559-018-0588-1 (2018).

    Google Scholar 

  15. Du, Z., Gelembiuk, G., Moss, W., Tritt, A. & Lee, C. E. The genome architecture of the copepod Eurytemora carolleeae – the highly invasive atlantic clade of the Eurytemora affinis species complex. Genom Proteom. Bioinform. 22. https://doi.org/10.1093/gpbjnl/qzae066 (2024).

  16. Rhoads, A. & Au, K. F. PacBio sequencing and its applications. Genom Proteom. Bioinform. 13, 278–289. https://doi.org/10.1016/j.gpb.2015.08.002 (2015).

    Google Scholar 

  17. BioBam & Bioinformatics OmicsBox. Bioinformatics Made Easy ( (BioBam Bioinformatics, 2019).

  18. Götz, S. et al. High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 36, 3420–3435. https://doi.org/10.1093/nar/gkn176 (2008).

    Google Scholar 

  19. Tatusov, R., Koonin, E. & Lipman, D. A genomic perspective on protein families. Science 278, 631–637. https://doi.org/10.1126/science.278.5338.631 (1997).

    Google Scholar 

  20. Sun, J. et al. OrthoVenn3: an integrated platform for exploring and visualizing orthologous data across genomes. Nucleic Acids Res. 51, W397–W403. (2023). https://doi.org/10.1093/nar/gkad313 (2023). 

    Google Scholar 

  21. Chen, F., Mackey, A. J., Stoeckert, C. J. & Roos, D. S. OrthoMCL-DB: querying a comprehensive multi-species collection of ortholog groups. Nucleic Acids Res. 34, D363–D368. https://doi.org/10.1093/nar/gkj123 (2006).

    Google Scholar 

  22. Xu, L. et al. OrthoVenn2: a web server for whole-genome comparison and annotation of orthologous clusters across multiple species. Nucleic Acids Res. 47, W52–W58. https://doi.org/10.1093/nar/gkz333 (2019).

    Google Scholar 

  23. Hartline, D. K. et al. De novo transcriptomes of six calanoid copepods (Crustacea): a resource for the discovery of novel genes. Sci. Data. 10, 242. https://doi.org/10.1038/s41597-023-02130-1 (2023).

    Google Scholar 

  24. Lee, B. Y. et al. RNA-seq based whole transcriptome analysis of the cyclopoid copepod Paracyclopina nana focusing on xenobiotics metabolism. Comp. Biochem. Physiol. D. 15, 12–19. https://doi.org/10.1016/j.cbd.2015.04.002 (2015).

    Google Scholar 

  25. Dahms, H. U. et al. Potential of the small cyclopoid copepod Paracyclopina nana as an invertebrate model for ecotoxicity testing. Aquat. Toxicol. 180, 282–294. https://doi.org/10.1016/j.aquatox.2016.10.013 (2016).

    Google Scholar 

  26. Drak Alsibai, K. et al. High positive correlations between ANRIL and p16-CDKN2A/p15-CDKN2B/p14-ARF Gene cluster overexpression in multi-tumor types suggest deregulated activation of an ANRIL-ARF bidirectional promoter. Non-coding RNA. 5. https://doi.org/10.3390/ncrna5030044 (2019).

  27. Schroeder, A. et al. The RIN: an RNA integrity number for assigning integrity values to RNA measurements. BMC Mol. Biol. 7, 3. https://doi.org/10.1186/1471-2199-7-3 (2006).

    Google Scholar 

  28. Kang, S. et al. The genome of the Antarctic-endemic copepod, Tigriopus kingsejongensis. GigaScience 6, 1–9; (2017). https://doi.org/10.1093/gigascience/giw010

  29. Chow, S., Ueno, Y., Toyokawa, M., Oohara, I. & Takeyama, H. Preliminary analysis of length and GC content variation in the ribosomal first internal transcribed spacer (ITS1) of marine animals. Mar. Biotech. 11, 301–306. https://doi.org/10.1007/s10126-008-9153-2 (2009).

    Google Scholar 

  30. Jørgensen, T. S. et al. The genome and mRNA transcriptome of the cosmopolitan calanoid copepod Acartia tonsa Dana improve the understanding of copepod genome size evolution. Genome Biol. Evol. 11, 1440–1450. https://doi.org/10.1093/gbe/evz067 (2019).

    Google Scholar 

  31. Shao, C. et al. The enormous repetitive Antarctic krill genome reveals environmental adaptations and population insights. Cell 186, 1279–1294. https://doi.org/10.1016/j.cell.2023.02.005 (2023). e19.

    Google Scholar 

  32. Bernot, J. P., Boxshall, G. A. & Crandall, K. A. A synthesis tree of the Copepoda: integrating phylogenetic and taxonomic data reveals multiple origins of parasitism. PeerJ 9, e12034. https://doi.org/10.7717/peerj.12034 (2021).

    Google Scholar 

  33. Wyngaard, G. A. & Rasch, E. M. Patterns of genome size in the copepoda. Hydrobiologia 417, 43–56. https://doi.org/10.1023/A:1003855322358 (2000).

    Google Scholar 

  34. Ma, X. & Johnson, K. B. Comparative phototaxis of calanoid and harpacticoid copepods. Mar. Biol. 164 https://doi.org/10.1007/s00227-016-3054-0 (2017).

  35. van Duren, L. A. & Videler, J. J. Escape from viscosity: the kinematics and hydrodynamics of copepod foraging and escape swimming. J. Exp. Biol. 206, 269–279. https://doi.org/10.1242/jeb.00079 (2003).

    Google Scholar 

  36. Tartarotti, B. et al. Seasonal plasticity in photoprotection modulates UV-induced hsp gene expression in copepods from a clear lake. Limnol. Oceanogr. 63, 1579–1592. https://doi.org/10.1002/lno.10793 (2018).

    Google Scholar 

  37. Isoseq v3 ((Pacific Biosciences of California, 2022)).

  38. Fu, L., Niu, B., Zhu, Z., Wu, S. & Li, W. CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinformatics 28, 3150–3152. https://doi.org/10.1093/bioinformatics/bts565 (2012).

    Google Scholar 

  39. Li, W. & Godzik, A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 22, 1658–1659. https://doi.org/10.1093/bioinformatics/btl158 (2006).

    Google Scholar 

  40. Raghavan, V., Kraft, L., Mesny, F. & Rigerte, L. A simple guide to de novo transcriptome assembly and annotation. Brief. Bioinform. 23. https://doi.org/10.1093/bib/bbab563 (2022).

  41. Shen, W., Le, S., Li, Y., Hu, F. & SeqKit A cross-platform and ultrafast toolkit for FASTA/Q file manipulation. PloS One. 11, e0163962. https://doi.org/10.1371/journal.pone.0163962 (2016).

    Google Scholar 

  42. Seppey, M., Manni, M., Zdobnov, E. M. & BUSCO Assessing genome assembly and annotation completeness. Meth. Mol. Biol. 227–245; (1962). https://doi.org/10.1007/978-1-4939-9173-0_14 (2019).

  43. Manni, M., Berkeley, M. R., Seppey, M., Simão, F. A. & Zdobnov, E. M. BUSCO Update: novel and streamlined workflows along with broader and deeper phylogenetic coverage for scoring of eukaryotic, prokaryotic, and viral genomes. Mol. Biol. Evol. 38, 4647–4654. https://doi.org/10.1093/molbev/msab199 (2021).

    Google Scholar 

  44. Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. Basic local alignment search tool. J. Mol. Biol. 215, 403–410. https://doi.org/10.1016/S0022-2836(05)80360-2 (1990).

    Google Scholar 

  45. Haas, B. TransDecoderThe Broad Institute, Github, (2012).

  46. Cantalapiedra, C. P., Hernández-Plaza, A., Letunic, I., Bork, P. & Huerta-Cepas, J. eggNOG-mapper v2: Functional annotation, orthology assignments, and domain prediction at the metagenomic scale. Mol. Biol. Evol. 38, 5825–5829. https://doi.org/10.1093/molbev/msab293 (2021).

    Google Scholar 

  47. Huerta-Cepas, J. et al. eggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses. Nucleic Acids Res. 47, D309–D314. https://doi.org/10.1093/nar/gky1085 (2019).

    Google Scholar 

  48. Gomes, S. I. L., Roca, C. P., Scott-Fordsmand, J. J. & Amorim, M. J. B. Identifying conserved UV exposure genes and mechanisms. Sci. Rep. 8, 8605. https://doi.org/10.1038/s41598-018-26865-9 (2018).

    Google Scholar 

  49. Bahrman, N. et al. Identification of genes differentially expressed in response to cold in Pisum sativum Using RNA sequencing analyses. Plants 8 https://doi.org/10.3390/plants8080288 (2019).

Download references

Acknowledgements

We would like to thank Paulina Kalita, Ethan Salino and Eric Weniger for their help in the field. We are grateful to Rujuta Vaidya and Abdul Ada for their support on the statistical analysis.

Funding

This research was funded by the Austrian Science Fund (FWF) [grant DOI https://doi.org/10.55776/P35886 to B.T].

Author information

Authors and Affiliations

Authors

Contributions

A.P. collected and identified specimens, extracted RNA, performed bioinformatic and statistical analyses and interpretation, curated the data for publication in the repository, and wrote the manuscript. M.K. provided statistical support and reviewed and edited the manuscript. B.T. designed and supervised the study, acquired the funding and made critical manuscript revisions. All authors approved of the final manuscript.

Corresponding authors

Correspondence to
Placide Ambre or Tartarotti Barbara.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Below is the link to the electronic supplementary material.

Supplementary Material 1 (download PDF )

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Cite this article

Ambre, P., Morgan, K. & Barbara, T. The de novo transcriptome of the freshwater copepod Cyclops abyssorum tatricus reveals high-elevation adaptation.
Sci Rep (2026). https://doi.org/10.1038/s41598-026-46084-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1038/s41598-026-46084-x

Keywords

  • Zooplankton
  • Long-read PacBio sequencing
  • Gene ontology
  • Alpine lakes


Source: Ecology - nature.com

Bacillus megaterium strain KGA3 increases saline–alkaline tolerance of maize by recruiting keystone taxa in rhizosphere soil

Antarctic Sanctuary: fishing effort responses to an international MPA in the Southern Ocean

Back to Top