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
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).
Whitman, D. & Ananthakrishnan, T. Phenotypic plasticity of insects. Phenotypic Plast. insects: Mech. consequences. 1–63. https://doi.org/10.1201/b10201 (2009).
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).
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).
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).
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).
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).
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).
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).
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).
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).
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
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).
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).
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).
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).
BioBam & Bioinformatics OmicsBox. Bioinformatics Made Easy ( (BioBam Bioinformatics, 2019).
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).
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).
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).
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).
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).
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).
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).
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).
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).
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).
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
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).
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).
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.
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).
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).
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).
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).
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).
Isoseq v3 ((Pacific Biosciences of California, 2022)).
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).
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).
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).
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).
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).
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).
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).
Haas, B. TransDecoderThe Broad Institute, Github, (2012).
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).
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).
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).
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).
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
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
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
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
