Genomic evidence for demographic fluctuations, genetic burdens and adaptive divergence in fourfinger threadfin Eleutheronema rhadinum

Abstract

Declining populations and bottlenecks lead to the accumulation of deleterious mutations in fish populations. These processes also trigger genetic purging, which is a key genetic factor in reducing the deleterious burdens and increasing population viability. However, there is a lack of empirical evidence on the interaction between demographic history and the genome-wide pattern of deleterious variations. Here, we generated genome resequencing data of Eleutheronema rhadinum from China and Thailand, representing the major distribution of the species’ southern regions. E. rhadinum had exceptionally low genome-wide variability and experienced dramatic population expansions followed by continuous declines. The geographical divergence, which occurred ~ 23,000 years ago, shaped different demographic trajectories and generated different regional patterns of deleterious mutations in China and Thailand populations. Several lines of evidence revealed that this geographical pattern of deleterious mutation was driven by the purging of highly deleterious mutations. We showed that purifying selection had inbreeding-associated fitness costs and was more efficient against missense mutations in the Thailand population, which had the lowest genetic burden of homozygous deleterious mutations. Multiple evolutionarily conserved protein domains were disrupted by the loss-of-function mutations, posing a high probability of gene functionality elimination. Moreover, thermal and salinity genes (Trpm3, Nek4, Gtf2f2, Cldn14) were identified in genomic divergence regions of E. rhadinum among China and Thailand populations. Our findings highlight the importance of demographic history factors shaping the geographical patterns of deleterious mutations. The results serve to deepen our understanding of the adaptive evolution and divergence of E. rhadinum with implications for other marine fish.

References

[1] Akbarzadeh A, Günther OP, Houde AL, Li S, Ming TJ, Jeffries KM, Hinch SG, Miller KM. Developing specific molecular biomarkers for thermal stress in salmonids. BMC Genom, 2018, 19: 1-28.

[2] Alexa A, Rahnenfuhrer J. topGO: enrichment analysis for gene ontology. R Package Vers, 2010, 2: 2010

[3] Alexander DH, Novembre J, Lange K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res, 2009, 19: 1655-1664.

[4] Andrews S (2010) FastQC: a quality control tool for high throughput sequence data. Available online at: http://www.bioinformatics.babraham.ac.uk/projects/fastqc.

[5] Barrio AM, Lamichhaney S, Fan G, Rafati N, Pettersson M, Zhang H, Dainat J, Ekman D, Höppner M, Jern P (2016). The genetic basis for ecological adaptation of the Atlantic herring revealed by genome sequencing. eLife 5:e12081

[6] Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics, 2014, 30: 2114-2120.

[7] Bradbury IR, Hubert S, Higgins B, Borza T, Bowman S, Paterson IG, Snelgrove PV, Morris CJ, Gregory RS, Hardie DC, Hutchings JA, Ruzzante DE, Taggart CT, Bentzen P. Parallel adaptive evolution of Atlantic cod on both sides of the Atlantic Ocean in response to temperature. Proc Biol Sci, 2010, 277: 3725-3734

[8] Bui P, Kelly SP. Claudin-6, -10d and -10e contribute to seawater acclimation in the euryhaline puffer fish Tetraodon nigroviridis. J Exp Biol, 2014, 217: 1758-1767

[9] Bui P, Bagherie-Lachidan M, Kelly S. Cortisol differentially alters claudin isoform mRNA abundance in a cultured gill epithelium from puffer fish (Tetraodon nigroviridis). Mol Cell Endocrinol, 2010, 317: 120-126.

[10] Ceballos FC, Joshi PK, Clark DW, Ramsay M, Wilson JF. Runs of homozygosity: windows into population history and trait architecture. Nat Rev Genet, 2018, 19: 220-234.

[11] Charlesworth B. Effective population size and patterns of molecular evolution and variation. Nat Rev Genet, 2009, 10: 195-205.

[12] Cingolani P, Platts A, Wang LL, Coon M, Nguyen T, Wang L, Land SJ, Lu X, Ruden DM. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly, 2012, 6: 80-92.

[13] Danecek P, Auton A, Abecasis G, Albers CA, Banks E, DePristo MA, Handsaker RE, Lunter G, Marth GT, Sherry ST, et al.. The variant call format and VCFtools. Bioinformatics, 2011, 27: 2156-2158.

[14] Facciolo RM, Crudo M, Giusi G, Canonaco M. GABAergic influences on ORX receptor-dependent abnormal motor behaviors and neurodegenerative events in fish. Toxicol Appl Pharmacol, 2010, 243: 77-86.

[15] Fay JC, Wyckoff GJ, Wu CI. Positive and negative selection on the human genome. Genetics, 2001, 158: 1227-1234.

[16] Feng C, Li X, Sha H, Luo X, Zou G, Liang H. Comparative transcriptome analysis provides novel insights into the molecular mechanism of the silver carp (Hypophthalmichthys molitrix) brain in response to hypoxia stress. Comp Biochem Physiol D, 2022, 41 100951

[17] Frankham R, Bradshaw CJ, Brook BW. Genetics in conservation management: revised recommendations for the 50/500 rules, Red List criteria and population viability analyses. Biol Conserv, 2014, 170: 56-63.

[18] Gao G, Nome T, Pearse DE, Moen T, Naish KA, Thorgaard GH, Lien S, Palti Y. New single nucleotide polymorphism database for rainbow trout generated through whole genome resequencing. Front Genet, 2018, 9: 147.

[19] Glémin S. How are deleterious mutations purged? Drift versus nonrandom mating. Evolution, 2003, 57: 2678-2687

[20] Grossen C, Guillaume F, Keller LF, Croll D. Purging of highly deleterious mutations through severe bottlenecks in Alpine ibex. Nat Commun, 2020, 11: 1001.

[21] Grummer JA, Beheregaray LB, Bernatchez L, Hand BK, Luikart G, Narum SR, Taylor EB. Aquatic landscape genomics and environmental effects on genetic variation. Trends Ecol Evol, 2019, 34: 641-654.

[22] Han ZQ, Guo XY, Liu Q, Liu SS, Zhang ZX, Xiao SJ, Gao TX. Whole-genome resequencing of Japanese whiting (Sillago japonica) provide insights into local adaptations. Zool Res, 2021, 42: 548.

[23] Horne JB, Momigliano P, Welch DJ, Newman SJ, Van Herwerden L. Searching for common threads in threadfins: phylogeography of Australian polynemids in space and time. Mar Ecol Prog Ser, 2012, 449: 263-276.

[24] Jaureguiberry P, Titeux N, Wiemers M, Bowler DE, Coscieme L, Golden AS, Guerra CA, Jacob U, Takahashi Y, Settele J, Díaz S, Molnár Z, Purvis A (2022) The direct drivers of recent global anthropogenic biodiversity loss. Sci Adv 8: eabm9982

[25] Jones P, Binns D, Chang H-Y, Fraser M, Li W, McAnulla C, McWilliam H, Maslen J, Mitchell A, Nuka G. InterProScan 5: genome-scale protein function classification. Bioinformatics, 2014, 30: 1236-1240.

[26] Keller LF, Waller DM. Inbreeding effects in wild populations. Trends Ecol Evol, 2002, 17: 230-241.

[27] Khan A, Patel K, Shukla H, Viswanathan A, van der Valk T, Borthakur U, Nigam P, Zachariah A, Jhala YV, Kardos M. Genomic evidence for inbreeding depression and purging of deleterious genetic variation in Indian tigers. Proc Natl Acad Sci USA, 2021, 118. e2023018118

[28] Knutsen H, Catarino D, Rogers L, Sodeland M, Mattingsdal M, Jahnke M, Hutchings JA, Mellerud I, Espeland SH, Johanneson K. Combining population genomics with demographic analyses highlights habitat patchiness and larval dispersal as determinants of connectivity in coastal fish species. Mol Ecol, 2022, 31: 2562-2577.

[29] Li H, Durbin R. Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics, 2010, 26: 589-595.

[30] Liao PC, Havanond S, Huang S. Phylogeography of Ceriops tagal (Rhizophoraceae) in Southeast Asia: the land barrier of the Malay Peninsula has caused population differentiation between the Indian Ocean and South China Sea. Conserv Genet, 2007, 8: 89-98.

[31] Liu Z, Ji J, Zheng D, Su L, Peng T. Calpain-2 protects against heat stress-induced cardiomyocyte apoptosis and heart dysfunction by blocking p38 mitogen-activated protein kinase activation. J Cell Physiol, 2019, 234: 10761-10770.

[32] Liu J, Liu Y, Li Q, Lu Y. Heat shock protein 70 and Cathepsin B genes are involved in the thermal tolerance of Aphis gossypii. Pest Manag Sci, 2023, 79: 2075-2086.

[33] Lynch M, Conery J, Burger R. Mutation accumulation and the extinction of small populations. Am Nat, 1995, 146: 489-518.

[34] Massullo P, Sumoza-Toledo A, Bhagat H, Partida-Sánchez S. TRPM channels, calcium and redox sensors during innate immune responses. Semin Cell Dev Biol, 2006, 17: 654-666.

[35] McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res, 2010, 20: 1297-1303.

[36] Motomura H, Iwatsuki Y, Kimura S, Yoshino T. Revision of the Indo-West Pacific polynemid fish genus Eleutheronema (Teleostei: Perciformes). Ichthyol Res, 2002, 49: 47-61.

[37] Newman S, Pember M, Rome B, Mitsopoulos G, Skepper C, Allsop Q, Saunders T, Ballagh A, Van Herwerden L, Garrett R. Stock structure of blue threadfin Eleutheronema tetradactylum across northern Australia as inferred from stable isotopes in sagittal otolith carbonate. Fish Manag Ecol, 2011, 18: 246-257.

[38] Papakostas S, Vasemaegi A, Vähä JP, Himberg M, Peil L, Primmer CR. A proteomics approach reveals divergent molecular responses to salinity in populations of European whitefish (Coregonus lavaretus). Mol Ecol, 2012, 21: 3516-3530.

[39] Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, De Bakker PI, Daly MJ. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet, 2007, 81: 559-575.

[40] Robinson JA, Brown C, Kim BY, Lohmueller KE, Wayne RK. Purging of strongly deleterious mutations explains long-term persistence and absence of inbreeding depression in island foxes. Curr Biol, 2018, 28: 3487-3494.

[41] Rosen D. The relationships and taxonomic position of the halfbreaks, killifishes, siversides, and their relatives. Bull Mus Nat Hist, 1964, 27: 219-267

[42] Samphan P. Abundance and species composition of Demersal fish, with descriptions of dominance species in Moo-Kho Bulon, Satun Province, Thailand. Int J Agric Technol, 2016, 12: 63-73

[43] Sun X, Xu D, Lou B, Zhang T, Xin J, Guo Y, Ma S. Genetic diversity and population structure of Eleutheronema rhadinum in the East and South China Seas revealed in mitochondrial COI sequences. Chin J Oceanol Limn, 2013, 31: 1276-1283.

[44] Terhorst J, Kamm JA, Song YS. Robust and scalable inference of population history from hundreds of unphased whole genomes. Nat Genet, 2017, 49: 303-309.

[45] Tiffin P, Ross-Ibarra J. Advances and limits of using population genetics to understand local adaptation. Trends Ecol Evol, 2014, 29: 673-680.

[46] Van Der Valk T, Díez-del-Molino D, Marques-Bonet T, Guschanski K, Dalén L. Historical genomes reveal the genomic consequences of recent population decline in eastern gorillas. Curr Biol, 2019, 29: 165-170.

[47] Vaser R, Adusumalli S, Leng SN, Sikic M, Ng PC. SIFT missense predictions for genomes. Nat Protoc, 2016, 11: 1-9.

[48] Velmurugan BK, Chan CR, Weng CF. Innate-immune responses of tilapia (Oreochromis mossambicus) exposure to acute cold stress. J Cell Physiol, 2019, 234: 16125-16135.

[49] Visudtiphole V, Watthanasurorot A, Klinbunga S, Menasveta P, Kirtikara K. Molecular characterization of Calreticulin: a biomarker for temperature stress responses of the giant tiger shrimp Penaeus monodon. Aquaculture, 2010, 308: S100-S108.

[50] Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nuleic Acids Res, 2010, 38: e164.

[51] Willi Y, Kristensen TN, Sgrò CM, Weeks AR, Ørsted M, Hoffmann AA. Conservation genetics as a management tool: the five best-supported paradigms to assist the management of threatened species. Proc Natl Acad Sci USA, 2022, 119. e2105076119

[52] Wu Q, Zheng P, Hu Y, Wei F. Genome-scale analysis of demographic history and adaptive selection. Protein Cell, 2014, 5: 99-112.

[53] Wurster CM, Bird MI, Bull ID, Creed F, Bryant C, Dungait JA, Paz V. Forest contraction in north equatorial Southeast Asia during the Last Glacial Period. Proc Natl Acad Sci USA, 2010, 107: 15508-15511.

[54] Xiao J, Wang WX. Linking HIF oxygen-sensing system diversity to hypoxia fitness in Eleutheronema: molecular characterization and transcriptional response to hypoxia exposure. Sci Total Environ, 2024, 909. 168646

[55] Xiao J, Lyu S, Iqbal TH, Hajisamae S, Tsim KW, Wang WX. Molecular phylogenetic and morphometric analysis of population structure and demography of endangered threadfin fish Eleutheronema from Indo-Pacific waters. Sci Rep, 2022, 12: 3455.

[56] Xiao J, Tsim KW, Hajisamae S, Wang WX. Chromosome-level genome and population genomics provide novel insights into adaptive divergence in allopatric Eleutheronema tetradactylum. Int J Biol Macromol, 2023, 244. 125299

[57] Xie C, Mao X, Huang J, Ding Y, Wu J, Dong S, Kong L, Gao G, Li CY, Wei L. KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nuleic Acids Res, 2011, 39: W316-W322.

[58] Xu S, Zhao L, Xiao S, Gao T. Whole genome resequencing data for three rockfish species of Sebastes. Sci Data, 2019, 6: 97.

[59] Xuan Z, Wang X, Hajisamae S, Tsim KW, Yang J, Wang WX. Otolith microchemistry reveals different environmental histories for two endangered fourfinger threadfin species. Mar Ecol Prog Ser, 2022, 700: 161-178.

[60] Xue Y, Prado-Martinez J, Sudmant PH, Narasimhan V, Ayub Q, Szpak M, Frandsen P, Chen Y, Yngvadottir B, Cooper DN, et al.. Mountain gorilla genomes reveal the impact of long-term population decline and inbreeding. Science, 2015, 348: 242-245.

[61] Yang J, Lee SH, Goddard ME, Visscher PM. GCTA: a tool for genome-wide complex trait analysis. Am J Hum Genet, 2011, 88: 76-82.

[62] Zhang X, Wen H, Wang H, Ren Y, Zhao J, Li Y. RNA-Seq analysis of salinity stress–responsive transcriptome in the liver of spotted sea bass (Lateolabrax maculatus). PLoS ONE, 2017, 12. e0173238

[63] Zhang C, Dong SS, Xu JY, He WM, Yang TL. PopLDdecay: a fast and effective tool for linkage disequilibrium decay analysis based on variant call format files. Bioinformatics, 2019, 35: 1786-1788.

[64] Zhang Y, Nie H, Yan X. Transient receptor potential (TRP) channels in the Manila clam (Ruditapes philippinarum): characterization and expression patterns of the TRP gene family under heat stress in Manila clams based on genome-wide identification. Gene, 2023, 854. 147112

[65] Zhong L, Wang M, Li D, Tang S, Chen X. Mitochondrial genome of Eleutheronema rhadinum with an additional non-coding region and novel insights into the phylogenetics. Front Mar Sci, 2021, 8. 746598

RIGHTS & PERMISSIONS

The Author(s)