New perspectives on the genetic structure of dotted gizzard shad (Konosirus punctatus) based on RAD-seq

Abstract

To maintain, develop and rationally utilize marine organisms, understanding their genetic structure and habitat adaptation pattern is necessary. Konosirus punctatus, which is a commercial fish species inhabiting the Indo-west Pacific Ocean, has shown an obvious annual global capture and aquaculture production decline due to climate changes and human activities. In the present study, restriction-site associated DNA sequencing (RAD-seq) was used to describe its genome-wide single nucleotide polymorphisms panel (SNPs). Among 146 individuals collected at nine locations scattered in China, Korea and Japan, a set of 632,090 SNPs were identified. Population genetic analysis showed that K. punctatus individuals were divided into two significant genetic clusters. Meanwhile, potential genetic differentiation between northern and southern population of K. punctatus was found. Treemix results indicated that gene flow existed among sampling locations of K. punctatus, especially from southern Japan to others. Moreover, candidate genes associated with habitat adaptations of K. punctatus were identified, which are involved in diverse physiological processes of K. punctatus including growth and development (e.g., KIDINS220, PAN3), substance metabolism (e.g., PGM5) and immune response (e.g., VAV3, CCT7, HSPA12B). Our findings may aid in understanding the possible mechanisms for the population genetic structure and local adaptation of K. punctatus, which is beneficial to establish the management and conservation units of K. punctatus, guiding the rational use of resources, with reference significance for a profound understanding of the adaptative mechanisms of other marine organisms to the environment.

Keywords

Population structure / Konosirus punctatus / SNP / RAD-seq / Local adaptation

References

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

[] Assis J, Tyberghein L, Bosch S, Verbruggen H, Serrão EA, De Clerck O. Bio-ORACLE v2. 0: extending marine data layers for bioclimatic modelling. Glob Ecol Biogeogr, 2018, 27: 277-284.

[] Bai Y, Qu A, Liu Y, Chen X, Wang J, Zhao J, Ke Q, Chen L, Chi H, Gong H. Integrative analysis of GWAS and transcriptome reveals p53 signaling pathway mediates resistance to visceral white-nodules disease in large yellow croaker. Fish Shellfish Immunol, 2022, 130: 350-358.

[] Billington N, Hebert PD. Mitochondrial DNA diversity in fishes and its implications for introductions. Can J Fish Aquat Sci, 1991, 48: 80-94.

[] Campbell D, Duchesne P, Bernatchez L. AFLP utility for population assignment studies: analytical investigation and empirical comparison with microsatellites. Mol Ecol, 2003, 12: 1979-1991.

[] Cano JM, Shikano T, Kuparinen A, Merilä J. Genetic differentiation, effective population size and gene flow in marine fishes: implications for stock management. J Integr Field Biol, 2008, 5:1-10.

[] Cantalapiedra CP, 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, 2021, 38: 5825-5829.

[] Catchen J, Hohenlohe PA, Bassham S, Amores A, Cresko WA. Stacks: an analysis tool set for population genomics. Mol Ecol, 2013, 22: 3124-3140.

[] Conover D, Clarke L, Munch S, Wagner G. Spatial and temporal scales of adaptive differentiation in marine fishes and the implications for conservation. J Fish Biol, 2006, 69: 21-47.

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

[] Devine B, Munday P, Jones G. Rising CO₂ concentrations affect settlement behaviour of larval damselfishes. Coral Reefs, 2012, 31: 229-238.

[] Dhaliwal A. DNA extraction and purification. Mater Methods, 2013, 3: 191.

[] Dong J, Wen J. Structure and function of small molecule GTPases. Chem Life, 2005, 25:422-424.

[] Eaton DA, Ree RH. Inferring phylogeny and introgression using RADseq data: an example from flowering plants (Pedicularis: Orobanchaceae). Syst Biol, 2013, 62: 689-706.

[] Etter PD, Bassham S, Hohenlohe PA, Johnson EA, Cresko WA. SNP discovery and genotyping for evolutionary genetics using RAD sequencing. Methods Mol Biol, 2011, 772: 157-178.

[] Fernández MV, Heras S, Maltagliati F, Turco A, Roldán MI. Genetic structure in the blue and red shrimp Aristeus antennatus and the role played by hydrographical and oceanographical barriers. Mar Ecol Prog Ser, 2011, 421: 163-171.

[] Fitak RR. OptM: estimating the optimal number of migration edges on population trees using Treemix. Biol Methods Protoc, 2021, 6: bpab017.

[] Foll M, Beaumont MA, Gaggiotti O. An approximate Bayesian computation approach to overcome biases that arise when using amplified fragment length polymorphism markers to study population structure. Genetics, 2008, 179: 927-939.

[] Forester BR, Lasky JR, Wagner HH, Urban DL. Comparing methods for detecting multilocus adaptation with multivariate genotype–environment associations. Mol Ecol, 2018, 27: 2215-2233.

[] Forsgren E, Dupont S, Jutfelt F, Amundsen T. Elevated CO₂ affects embryonic development and larval phototaxis in a temperate marine fish. Ecol Evol, 2013, 3: 3637-3646.

[] Gibbard P, Kolfschoten T. Gradstein FM, Ogg JG, Smith AG. The pleistocene and holocene epochs. A geologic time scale 2004, 2005 Cambridge Cambridge University Press 441-452.

[] Govindaraju DR. Variation in gene flow levels among predominantly self-pollinated plants. J Evol Biol, 1989, 2: 173-181.

[] Gustafsson R, Eckhard U, Ye W, Enbody ED, Pettersson M, Jemth P, Andersson L, Selmer M. Structure and characterization of phosphoglucomutase 5 from Atlantic and Baltic herring—an inactive enzyme with intact substrate binding. Biomolecules, 2020, 10: 1631.

[] Gwak WS, Lee YD, Nakayama K. Population structure and sequence differentiation in the mitochondrial DNA control region of gizzard shad Konosirus punctatus in Korea and Japan. Ichthyol Res, 2015, 62: 379-385.

[] Hamilton SL, Regetz J, Warner RR. Postsettlement survival linked to larval life in a marine fish. Proc Natl Acad Sci, 2008, 105: 1561-1566.

[] Han G, Huang W. Pacific decadal oscillation and sea level variability in the Bohai, Yellow, and East China seas. J Phys Oceanogr, 2008, 38: 2772-2783.

[] Han ZQ, Gao TX, Yanagimoto T, Sakurai Y. Deep phylogeographic break among white croaker Pennahia argentata (Sciaenidae, Perciformes) populations in North-western Pacific. Fish Sci, 2008, 74: 770-780.

[] Hellberg ME, Burton RS, Neigel JE, Palumbi SR. Genetic assessment of connectivity among marine populations. Bull Mar Sci, 2002, 70:273-290.

[] Howe K, Clark MD, Torroja CF, Torrance J, Berthelot C, Muffato M, Collins JE, Humphray S, McLaren K, Matthews L. The zebrafish reference genome sequence and its relationship to the human genome. Nature, 2013, 496: 498-503.

[] Jeffery NW, DiBacco C, Van Wyngaarden M, Hamilton LC, Stanley RRE, Bernier R, FitzGerald J, Matheson K, McKenzie CH, Ravindran PN, Beiko R, Bradbury IR. RAD sequencing reveals genomewide differentiation between independent invasions of the European green crab (Carcinus maenas) in the Northwest Atlantic. Ecol Evol, 2017, 7: 2513-2524.

[] Johansen J, Jones G. Increasing ocean temperature reduces the metabolic performance and swimming ability of coral reef damselfishes. Glob Change Biol, 2011, 17: 2971-2979.

[] Kincade PW, Cooper MD. Development and distribution of immunoglobulin-containing cells in the chicken: an immunofluorescent analysis using purified antibodies to µ, γ and light chains. J Immunol, 1971, 106: 371-382.

[] Kong L, Kawasaki M, Kuroda K, Kohno H, Fujita K. Spawning characteristics of the konoshiro gizzard shad in Tokyo and Sagami Bays, central Japan. Fish Sci, 2004, 70: 116-122.

[] Korneliussen TS, Albrechtsen A, Nielsen R. ANGSD: analysis of next generation sequencing data. BMC Bioinform, 2014, 15: 1-13.

[] Krück NC, Innes DI, Ovenden JR. New SNP s for population genetic analysis reveal possible cryptic speciation of eastern Australian sea mullet (Mugil cephalus). Mol Ecol Resour, 2013, 13: 715-725.

[] Lasky JR, Des Marais DL, McKAY JK, Richards JH, Juenger TE, Keitt TH. Characterizing genomic variation of Arabidopsis thaliana: the roles of geography and climate. Mol Ecol, 2012, 21: 5512-5529.

[] Lecomte F, Grant WS, Dodson J, Rodríguez-Sánchez R, Bowen B. Living with uncertainty: genetic imprints of climate shifts in East Pacific anchovy (Engraulis mordax) and sardine (Sardinops sagax). Mol Ecol, 2004, 13: 2169-2182.

[] Lee T. Age composition and reproductive period of the shad, Konoserus punctztus, in Cheonsu Bay. J Korean Oceanogr Soc, 1983, 18:161 168 (in Korean)

[] Legendre P, Legendre L. Ordination in reduced space. Dev Environ Model, 2012, 24:425-520.

[] Letunic I, Bork P. Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res, 2021, 49: W293-W296.

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

[] Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R. The sequence alignment/map format and SAMtools. Bioinformatics, 2009, 25: 2078-2079.

[] Li M, Shi S, Wang M, Lun Z, Yang T. Complete mitochondrial genome of the dotted gizzard shad Konosirus punctatus (Teleostei, Clupeidae). Mitochondrial DNA, 2012, 23: 295-297.

[] Liu B (2017) Population genetic structure and local adaptation of the small yellow croaker (Larimichthys polyactis) and Japanese eel (Anguilla japonica). Dissertation, Institute of Oceanology, Chinese Academy of Sciences

[] Liu X, Fu YX. Stairway Plot 2: demographic history inference with folded SNP frequency spectra. Genome Biol, 2020, 21:1-9.

[] Liu JX, Gao TX, Yokogawa K, Zhang YP. Differential population structuring and demographic history of two closely related fish species, Japanese sea bass (Lateolabrax japonicus) and spotted sea bass (Lateolabrax maculatus) in Northwestern Pacific. Mol Phylogenet Evol, 2006, 39: 799-811.

[] Liu JX, Gao TX, Zhuang ZM, Jin XS, Yokogawa K, Zhang YP. Late Pleistocene differentiation and subsequent population expansion of two closely related fish species, Japanese anchovy (Engraulis japonicus) and Australian anchovy (Engraulis australis). Mol Phylogenet Evol, 2006, 40: 712-723.

[] Liu JX, Gao TX, Wu SF, Zhang YP. Pleistocene isolation in the Northwestern Pacific marginal seas and limited dispersal in a marine fish, Chelon haematocheilus (Temminck & Schlegel, 1845). Mol Ecol, 2007, 16: 275-288.

[] Liu M, Wang Y, Gao T, Yanagimoto T, Zhuang Z, Sakurai Y. Morphological comparison of Pacific herring Clupea pallasii in China and Japan. J Ocean Univ China, 2007, 37:131-136.

[] Liu B, Zhang K, Zhu K, Shafi M, Gong L, Jiang L, Liu L, Muhammad F, Lü Z. Population genetics of Konosirus punctatus in Chinese coastal waters inferred from two mtDNA genes (COI and Cytb). Front Mar Sci, 2020, 7: 534.

[] Liu BJ, Zhang K, Zhang SF, Liu YF, Li JS, Peng Y, Jin X, Wang YP, Zheng SX, Gong L, Liu LQ, Lyu ZM. Chromosome-level genome assembly of the dotted gizzard shad (Konosirus punctatus) provides insights into its adaptive evolution. Zool Res, 2022, 43: 217.

[] Liu D, Gui L, Zhu Y, Xu C, Zhou W, Li M. Chromosome-level assembly of male Opsariichthys bidens genome provides insights into the regulation of the gnrh signaling pathway and genome evolution. Biology, 2022, 11: 1500.

[] Lopes AF, Morais P, Pimentel M, Rosa R, Munday PL, Gonçalves EJ, Faria AM. Behavioural lateralization and shoaling cohesion of fish larvae altered under ocean acidification. Mar Biol, 2016, 163: 1-10.

[] Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. Embnet J, 2011, 17: 10-12.

[] Maulvault AL, Santos LH, Paula JR, Camacho C, Pissarra V, Fogaça F, Barbosa V, Alves R, Ferreira PP, Barceló D. Differential behavioural responses to venlafaxine exposure route, warming and acidification in juvenile fish (Argyrosomus regius). Sci Total Environ, 2018, 634: 1136-1147.

[] McCormack JE, Hird SM, Zellmer AJ, Carstens BC, Brumfield RT. Applications of next-generation sequencing to phylogeography and phylogenetics. Mol Phylogenet Evol, 2013, 66: 526-538.

[] 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.

[] Meirmans PG, Hedrick PW. Assessing population structure: FST and related measures. Mol Ecol Resour, 2011, 11: 5-18.

[] Metz CW, Bridges CB. Incompatibility of mutant races in Drosophila. Proc Natl Acad Sci, 1917, 3: 673-678.

[] Miller MR, Dunham JP, Amores A, Cresko WA, Johnson EA. Rapid and cost-effective polymorphism identification and genotyping using restriction site associated DNA (RAD) markers. Genome Res, 2007, 17: 240-248.

[] Mohr OL. Character changes caused by mutation of an entire region of a chromosome in Drosophila. Genetics, 1919, 4: 275.

[] Molt S, Bührdel JB, Yakovlev S, Schein P, Orfanos Z, Kirfel G, Winter L, Wiche G, van der Ven PF, Rottbauer W. Aciculin interacts with filamin C and Xin and is essential for myofibril assembly, remodeling and maintenance. J Cell Sci, 2014, 127: 3578-3592.

[] Munday PL, Dixson DL, Donelson JM, Jones GP, Pratchett MS, Devitsina GV, Døving KB. Ocean acidification impairs olfactory discrimination and homing ability of a marine fish. Proc Natl Acad Sci, 2009, 106: 1848-1852.

[] Myoung SH, Kim JK. Genetic diversity and population structure of the gizzard shad, Konosirus punctatus (Clupeidae, Pisces), in Korean waters based on mitochondrial DNA control region sequences. Genes Genomics, 2014, 36: 591-598.

[] Nielsen EE, Hemmer-Hansen J, Poulsen NA, Loeschcke V, Moen T, Johansen T, Mittelholzer C, Taranger G-L, Ogden R, Carvalho GR. Genomic signatures of local directional selection in a high gene flow marine organism; the Atlantic cod (Gadus morhua). BMC Evol Biol, 2009, 9: 1-11.

[] Okamura O, Yamada U (1986) Fishes of the East China Sea and the Yellow Sea. Contrib Seikai Reg Fish Res Lab 422

[] Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, McGlinn D, Minchin PR, O’hara R, Simpson G, Solymos P. The vegan package. Community Ecol Package, 2007, 10:719

[] Okubo K, Nagahama Y. Structural and functional evolution of gonadotropin-releasing hormone in vertebrates. Acta Physiol, 2008, 193: 3-15.

[] Ouborg NJ, Pertoldi C, Loeschcke V, Bijlsma RK, Hedrick PW. Conservation genetics in transition to conservation genomics. Trends Genet, 2010, 26: 177-187.

[] Pankhurst NW, Munday PL. Effects of climate change on fish reproduction and early life history stages. Mar Freshw Res, 2011, 62: 1015-1026.

[] Park SY, Kim JY, Yoon JM. Geographic variations and DNA polymorphisms in gizzard-shad (Konosirus punctatus). Korean J Ichthyol, 2006, 18:300-310.

[] Peng Y, Li J, Zhang K, Liu Y, Li X, Zhang H, Gong L, Liu L, Lü Z, Liu B. Identification of a large dataset of SNPs in hair-fin anchovy (Setipinna tenuifilis) based on RAD-seq. Anim Genet, 2021, 52: 371-374.

[] Pickrell J, Pritchard J, Tang H. Inference of population splits and mixtures from genome-wide allele frequency data. PLoS Genet, 2012, 8.

[] Pörtner HO, Knust R. Climate change affects marine fishes through the oxygen limitation of thermal tolerance. Science, 2007, 315: 95-97.

[] 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.

[] Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics, 2010, 26: 841-842.

[] Retief JD. Phylogenetic analysis using PHYLIP. Bioinformatics methods and protocols, 2000 Berlin Springer 243-258.

[] Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol, 1987, 4: 406-425.

[] Sauer DJ, Heidinger BJ, Kittilson JD, Lackmann AR, Clark ME. No evidence of physiological declines with age in an extremely long-lived fish. Sci Rep, 2021, 11: 1-10.

[] Shanks AL, Grantham BA, Carr MH. Propagule dispersal distance and the size and spacing of marine reserves. Ecol Appl, 2003, 13: 159-169.

[] Skoda AM, Simovic D, Karin V, Kardum V, Vranic S, Serman L. The role of the hedgehog signaling pathway in cancer: a comprehensive review. Bosn J Basic Med Sci, 2018, 18: 8.

[] Song N, Gao T, Ying Y, Yanagimoto T, Han Z. Is the Kuroshio Current a strong barrier for the dispersal of the gizzard shad (Konosirus punctatus) in the East China Sea?. Mar Freshw Res, 2016, 68: 810-820.

[] Steinig EJ, Neuditschko M, Khatkar MS, Raadsma HW, Zenger KR. netview p: a network visualization tool to unravel complex population structure using genome-wide SNPs. Mol Ecol Resour, 2016, 16: 216-227.

[] Sutachan JJ, Chao MV, Ninan I. Regulation of inhibitory neurotransmission by the scaffolding protein ankyrin repeat-rich membrane spanning/kinase D-interacting substrate of 220 kDa. J Neurosci Res, 2010, 88: 3447-3456.

[] Team RC R: a language and environment for statistical computing. Computing, 2013, 1:12-21.

[] Ulc A, Gottschling C, Schäfer I, Wegrzyn D, Sv L, Luft V, Reinhard J, Faissner A. Involvement of the guanine nucleotide exchange factor Vav3 in central nervous system development and plasticity. Biol Chem, 2017, 398: 663-675.

[] Vousden KH, Lu X. Live or let die: the cell's response to p53. Nat Rev Cancer, 2002, 2: 594-604.

[] Waits L, Sork V. Landscape genetics. Mol Ecol, 2010, 19: 3489-3835.

[] Wampler PJ, Rediske RR, Molla AR. Using ArcMap, Google Earth, and global positioning systems to select and locate random households in rural Haiti. Int J Health Geogr, 2013, 12: 1-8.

[] Wang P. Response of Western Pacific marginal seas to glacial cycles: paleoceanographic and sedimentological features. Mar Geol, 1999, 156: 5-39.

[] Weng AP, Ferrando AA, Lee W, Morris JP IV Silverman LB, Sanchez-Irizarry C, Blacklow SC, Look AT, Aster JC. Activating mutations of NOTCH1 in human T cell acute lymphoblastic leukemia. Science, 2004, 306: 269-271.

[] Wu L, Qin Z, Liu H, Lin L, Ye J, Li J. Recent advances on phagocytic B cells in teleost fish. Front Immunol, 2020, 11: 824.

[] Yang M, Wu J, Zhang L, Huang Y, Zi Y, Zhang Y. Effects of miR-143 on invasion and migration of B cell lymphoma cells through regulating expression of PAN3 protein. Chin J Immunol, 2017, 33:1774

[] Yao S, Bo Z, Xuewu G, Jun W, Qisheng T. Feeding, growth and ecological conversion efficiency of Konosirus punctatus—field gastric content method in vivo. Mar Fish Res, 1999, 20:12 16 (in Chinese)

[] Ying Y (2011) Population genetic structure and molecular phylogeography of Konosirus punctatus and Sardinella zunasi. Dissertation, Ocean University of China

[] Yokogawa K. Morphological and genetic differences between Japanese and Chinese Red Ark Shell (Scapharca broughtonii). Fish Sci, 1997, 63: 332-337.

[] Yoo DG, Lee GS, Kim GY, Kang NK, Yi BY, Kim YJ, Chun JH, Kong GS. Seismic stratigraphy and depositional history of late Quaternary deposits in a tide-dominated setting: an example from the eastern Yellow Sea. Mar Pet Geol, 2016, 73: 212-227.

[] Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS, 2012, 16: 284-287.

[] Zhang S (2001) Fauna sinica: Osteichthyes, Acipenseriformes, Elopiformes, Clupeiformes, Gonorhynchiformes. Science Press

[] Zhang BD, Xue DX, Wang J, Li YL, Liu BJ, Liu JX. Development and preliminary evaluation of a genomewide single nucleotide polymorphisms resource generated by RAD-seq for the small yellow croaker (Larimichthys polyactis). Mol Ecol Resour, 2016, 16: 755-768.

[] Zhang BD, Xue DX, Li YL, Liu JX. RAD genotyping reveals fine-scale population structure and provides evidence for adaptive differentiation in a commercially important fish from the northwestern Pacific Ocean. PeerJ, 2019, 7.

[] Zhao Q (2018) A study on optimizing markduplicate in genome sequencing pipeline. In: Proceedings of the 2018 5th international conference on bioinformatics research and applications, pp 8–15

[] Zhou Y (2011) Function research of Kidins220/ARMS on the mouse embrvonic stem cells. Dissertation, Zhejiang University

[] Zhou H, Qian J, Li C, Li J, Zhang X, Ding Z, Gao X, Han Z, Cheng Y, Liu L. Attenuation of cardiac dysfunction by HSPA12B in endotoxin-induced sepsis in mice through a PI3K-dependent mechanism. Cardiovasc Res, 2010, 89: 109-118.

[] Zulliger D, Schnyder E, Gugerli F. Are adaptive loci transferable across genomes of related species? Outlier and environmental association analyses in Alpine Brassicaceae species. Mol Ecol, 2013, 22: 1626-1639.