RAD-Seq-derived SNPs reveal no local population structure in the commercially important deep-sea queen snapper (Etelis oculatus) in Puerto Rico

Abstract

The queen snapper (Etelis oculatus Valenciennes in Cuvier & Valenciennes, 1828) is a deep-sea snapper whose commercial importance continues to increase in the US Caribbean. However, little is known about the biology and ecology of this species. In this study, the presence of a fine-scale population structure and genetic diversity of queen snapper from Puerto Rico was assessed through 16,188 SNPs derived from the Restriction site Associated DNA Sequencing (RAD-Seq) technique. Summary statistics estimated low genetic diversity (HO = 0.333–0.264) and did not reveal population differentiation within our samples (F_ST = − 0.001–0.025). Principal component analysis and a model-based clustering method did not detect a fine-scale subpopulation structure among sampling sites, however, there was genetic variability within regions and sites. Our results have revealed comparable genetic and dispersal patterns to those observed in other shallow-water snapper species in Puerto Rico waters. It is crucial to further enhance our understanding of the ecological and biological aspect of the queen snapper to effectively manage and conserve this species as fishing pressure has been extended to deep water species in the US Caribbean.

Special Topic: Ecology & Environmental Biology.

The online version contains supplementary material available at https://doi.org/10.1007/s42995-025-00289-7.

Keywords

Etelis oculatus / Genetic diversity / Caribbean Fisheries / Genomics

References

[1] Allen GR (1985). FAO Species Catalogue. Vol. 6. Snappers of the world. An annotated and illustrated catalogue of lutjanid species known to date. FAO Fish Synop 125:208 Rome: FAO

[2] AndrewsKR, WilliamsAJ, Fernandez-SilvaI, NewmanSJ, CopusJM, WakefieldCB, RandallJE, BowenBW. Phylogeny of deepwater snappers (Genus Etelis) reveals a cryptic species pair in the Indo-Pacific and Pleistocene invasion of the Atlantic. Mol Phylogenet Evol, 2016, 100: 361-371.

[3] AndrewsKR, CopusJM, WilcoxC, WilliamsAJ, NewmanSJ, WakefieldCB, BowenBW. Range-Wide population structure of 3 deepwater eteline snappers across the Indo-Pacific basin. J Hered, 2020, 111: 471-485.

[4] BacoAR, EtterRJ, RibeiroPA, Von der HeydenS, BeerliP, KinlanBP. A synthesis of genetic connectivity in deep-sea fauna and implications for marine reserve design. Mol Ecol, 2016, 25: 3276-3298.

[5] BaumsIB, ParisCB, CherubinLM. A bio-oceanographic filter to larval dispersal in a reef-building coral. Limnol Oceanogr, 2006, 51: 1969-1981.

[6] BeltránDM, SchizasNV, AppeldoornRS, PradaC. Effective dispersal of Caribbean reef fish is smaller than current spacing among marine protected areas. Sci Rep, 2017, 7: 1-10.

[7] BiggsCR, NemethRS. Spatial and temporal movement patterns of two snapper species at a multi-species spawning aggregation. Mar Ecol Prog Ser, 2016, 558: 129-142.

[8] BorsEK, HerreraS, MorrisJAJr, ShanTM. Population genomics of rapidly invading lionfish in the Caribbean reveal signals of range expansion in the absence of spatial population structure. Ecol Evol, 2019, 9: 3306-3320.

[9] Bryan-BrownDN, BrownCJ, HughesJM, ConnollyRM. Patterns and trends in marine population connectivity research. Mar Ecol Prog Ser, 2017, 585: 243-256.

[10] CarsonEW, SaillantEA, RenshawMA, CummingsNJ, GoldJR. Population structure, long-term connectivity, and effective size of mutton snapper (Lutjanus analis) in the Caribbean Sea and Florida Keys. Fish Bull, 2011, 109: 416-428

[11] CatchenJ, HohenlohePA, BasshamS, AmoresA, CreskoWA. Stacks: an analysis tool set for population genomics. Mol Ecol, 2013, 22: 3124-3140.

[12] ChhatreVE, EmersonKJ. StrAuto: automation and parallelization of STRUCTURE analysis BMC. Bioinformatics, 2017, 18: 1-5

[13] CowenRK, SponaugleS. Larval dispersal and marine population connectivity. Annu Rev Mar Sci, 2009, 1: 443-466.

[14] CowenRK, ParisCB, SrinivasanA. Scaling of connectivity in marine populations. Science, 2006, 311: 522-527.

[15] Crabtree RE (2019) Comprehensive Fishery Management Plan for the Puerto Rico Exclusive Economic Zone and Environmental Assessment. https://repository.library.noaa.gov/view/noaa/45269

[16] CrawfordDL, OleksiakMF. Ecological population genomics in the marine environment. Brief Funct Genomics, 2016, 15: 342-351.

[17] Cummings NJ, Matos-Caraballo D (2003) Summarized reported commercial landings in Puerto Rico from 1969–2001 with specific notes on the silk snapper landings category. Sustainable Fisheries Division Contribution No. SFD 2003-0023 and Caribbean Deepwater SEDAR Data Workshop Report 4 SEDAR https://sedarweb.org/documents/s4dw05-summarized-reported-commercial-landings-in-puerto-rico-from-1969-2001-with-specific-notes-on-the-silk-snapper-landing-category/

[18] D’AlessandroEK, SponaugleS, SerafyJE. Larval ecology of a suite of snappers (family: Lutjanidae) in the Straits of Florida, western Atlantic Ocean. Mar Ecol Prog Ser, 2010, 410: 159-175.

[19] de Oliveira LeisM, Barragán-PaladinesMJ, SaldañaA, BishopD, JinJH, KerežiV, ChuenpagdeeRSalasS, Barragán-PaladinesMJ, ChuenpagdeeR. Overview of small-scale fisheries in Latin America and the Caribbean: challenges and prospects. Viability and sustainability of small-scale fisheries in Latin America and the Caribbean, 2019MARE Publication. 1547.

[20] EvannoG, RegnautS, GoudetJ. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol, 2005, 14: 2611-2620.

[21] FollM, GaggiottiO. A genome-scan method to identify selected loci appropriate for both dominant and codominant markers: a Bayesian perspective. Genetics, 2008, 180: 977-993.

[22] Galindo-CortesG, Jiménez-BadilloL, MeinersCSalasS, Barragán-PaladinesMJ, ChuenpagdeeR. Moving from stock assessment to fisheries management in Mexico: the finfish fisheries from the southern Gulf of Mexico and Caribbean Sea. Viability and sustainability of small-scale fisheries in Latin America and the Caribbean, 2019MARE Publication. 243263.

[23] García-Sais JR, Castro R, Sabater-Clavell J, Carlo M, Esteves-Amador RF, Bruckner AW, Fairey L (2007) Characterization of benthic habitats and associated reef communities at Bajo de Sico seamount, Mona Passage, Puerto Rico final report. https://repository.library.noaa.gov/view/noaa/562

[24] GarySF, FoxAD, BiastochA, RobertsJM, CunninghamSA. Larval behaviour, dispersal and population connectivity in the deep sea. Sci Rep, 2020, 1010675.

[25] GobertB, GuillouA, MurrayP, BerthouP, TurciosMDO, LopezE, Jérôme-HuetPL, DiazN, GervainP. Biology of queen snapper (Etelis oculatus: Lutjanidae) in the Caribbean. Fish Bull, 2005, 103: 417-425

[26] Goudet J, Jombart T (2015) hierfstat: estimation and tests of hierarchical F-statistics. R package version 0.04-22, 10

[27] HellbergME. Gene flow and isolation among populations of marine animals. Ann Rev Ecol Evol Syst, 2009, 40291.

[28] HernándezCM, ParisCB, VazAC, JonesBT, KellnerJB, RichardsonDE, SponaugleS, CowenRK, LlopizJK. Diverse patterns of larval coral reef fish vertical distribution and consequences for dispersal and connectivity. Coral Reefs, 2023, 42: 453-465.

[29] HerreraS, WatanabeH, ShankTM. Evolutionary and biogeographical patterns of barnacles from deep-sea hydrothermal vents. Mol Ecol, 2015, 24: 673-689.

[30] HollenbeckCM, PortnoyDS, SaillantE, GoldJR. Population structure of red snapper (Lutjanus campechanus) in US waters of the western Atlantic Ocean and the northeastern Gulf of Mexico. Fish Res, 2015, 172: 17-25.

[31] JombartT. adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics, 2008, 24: 1403-1405.

[32] KadisonE, BrandtM, NemethR, MartensJ, BlondeauJ, SmithT. Abundance of commercially important reef fish indicates different levels of over-exploitation across shelves of the US Virgin Islands. PLoS ONE, 2017, 12. e0180063

[33] KamvarZN, TabimaJF, GrünwaldNJ. Poppr: an R package for genetic analysis of populations with clonal, partially clonal, and/or sexual reproduction. PeerJ, 2014, 2. e281

[34] KnausBJ, GrünwaldNJ. vcfr: a package to manipulate and visualize variant call format data in R. Mol Ecol Resour, 2017, 17: 44-53.

[35] KoolJT, MoilanenA, TremlEA. Population connectivity: recent advances and new perspectives. Landscape Ecol, 2013, 28: 165-185.

[36] KopelmanNM, MayzelJ, JakobssonM, RosenbergNA, MayroseI. Clumpak: a program for identifying clustering modes and packaging population structure inferences across K. Mol Ecol Resour, 2015, 15: 1179-1191.

[37] LeisJM, LeeK. Larval development in the lutjanid subfamily Etelinae (Pisces): the genera Aphareus, Aprion, Etelis and Pristipomoides. Bull Mar Sci, 1994, 551: 46-125

[38] LiYL, LiuJX. StructureSelector: a web based software to select and visualize the optimal number of clusters using multiple methods. Mol Ecol Resour, 2018, 18: 176-177.

[39] Matos-CaraballoD. Overview of Puerto Rico’s small-scale fisheries statistics, 1994–1997. Proc Gulf Caribb Fish Inst, 2000, 55: 103-141

[40] Matos-CaraballoD. Lessons learned from the Puerto Rico’s commercial fishery, 1988–2008. Proc Gulf Caribb Fish Inst, 2009, 61: 123-129

[41] Matos-Caraballo D (2012) Puerto Rico/NMFS Cooperative Fisheries Statistics Program April 2007–September 2012. Final report to NOAA:NA07NMF4340039. https://estadisticas.pr/files/Inventario/publicaciones/DRNA_CommercialFisheriesStatistics_2007-2011_0.pdf

[42] MeirmansPG. Genodive version 30: easy-to-use software for the analysis of genetic data of diploids and polyploids. Mol Ecol Resour, 2020, 20: 1126-1131.

[43] MussmannSM, DouglasMR, ChafinTK, DouglasME. BA3-SNPs: contemporary migration reconfigured in BayesAss for next-generation sequence data. Methods Ecol Evol, 2019, 10: 1808-1813.

[44] NazarenoAG, BemmelsJB, DickCW, LohmannLG. Minimum sample sizes for population genomics: an empirical study from an Amazonian plant species. Mol Ecol Resour, 2017, 6: 1136-1147.

[45] Neuwirth E, Brewer RC (2014) ColorBrewer palettes. R package version 1:4

[46] NielsenEE, Hemmer-HansenJ, LarsenPF, BekkevoldD. Population genomics of marine fishes: identifying adaptive variation in space and time. Mol Ecol, 2009, 18: 3128-3150.

[47] NOAA (2011) National Marine Fisheries Service 2010 report to Congress on the status of U.S. Fisheries., National Oceanographic and Atmospheric Administration. https://repository.library.noaa.gov/view/noaa/15614

[48] Ojeda-SerranoE, AppeldoornRS, Ruiz-ValentinI. Reef fish spawning aggregations of the Puerto Rican shelf. Proc Gulf Caribb Fish Inst, 2007, 59: 467-474

[49] OkuyamaJ, ShishidouH, HayashibaraT. Post-release horizontal and vertical behavior and philopatry of deepwater longtail red snapper Etelis coruscans around a bank. Fish Sci, 2019, 85: 361-368.

[50] OleksiakMF, RajoraOPOleksiakMF, RajoraOP. Marine population genomics: challenges and opportunities. Population genomics: marine organisms, 2020Springer Nature. 335.

[51] OverlyKE, ShervetteVR. Caribbean deepwater snappers: application of the bomb radiocarbon age estimation validation in understanding aspects of ecology and life history. PLoS One, 2023, 18. e0295650

[52] PalumbiSR. Population genetics, demographic connectivity, and the design of marine reserves. Ecol Appl, 2003, 13: 146-158.

[53] ParadisE, ClaudeJ, StrimmerK. APE: analyses of phylogenetics and evolution in R language. Bioinformatics, 2004, 20: 289-290.

[54] PetersonBK, WeberJN, KayEH, FisherHS, HoekstraHE. Double digest RADseq: an inexpensive method for de novo SNP discovery and genotyping in model and non-model species. PLoS One, 2012, 7. e37135

[55] PritchardJK, StephensM, DonnellyP. Inference of population structure using multilocus genotype data. Genetics, 2000, 155: 945-959.

[56] PurcellS, NealeB, Todd-BrownK, ThomasL, FerreiraMA, BenderD, MallerJ, SklarP, De BakkerPI, DalyMJ, ShamPC. PLINK: a tool set for whole-genome association and population-based linkage analyses. Ame J Hum Genet, 2007, 81: 559-575.

[57] QuWM, LiangN, WuZK, ZhaoYG, ChuD. Minimum sample sizes for invasion genomics: empirical investigation in an invasive whitefly. Ecol Evol, 2020, 10: 38-49.

[58] Rivera-ColónAG, CatchenJVerdeC, GiordanoD. Population genomics analysis with RAD, reprised: stacks 2. Marine genomics. Methods in molecular biology, 2022, New York. Humana. 99149. 2498

[59] Rosado-NicOdJ, HoganJD, Lara-ArenasJH, Rosas-LuisR, CarrilloL, Villegas-SánchezCA. Gene flow between subpopulations of gray snapper (Lutjanus griseus) from the Caribbean and Gulf of Mexico. PeerJ, 2020, 8. e8485

[60] Rosario A, Rojas J, Piñeiro E, Figuerola M, Peña N, Torres W (2006) Completion Report to National Marine Fisheries Service Reproductive Cycle of Queen Snapper (Etelis oculatus) and the Wenchman (Pristipomoides macrophthalmus). Mayaguez, Puerto Rico, Fisheries Research Laboratory, 31 pp

[61] SaillantEA, RenshawMA, CummingsNJ, GoldJR. Conservation genetics and management of yellowtail snapper, Ocyurus chrysurus, in the US Caribbean and South Florida. Fish Manage Ecol, 2012, 4: 301-312.

[62] Salas S, Chuenpagdee R, Charles A, Seijo, JC (2011) Coastal fisheries of Latin America and the Caribbean region: issues and trends. In: Salas S, Chuenpagdee R, Charles A, Seijo JC (eds) Coastal fisheries of Latin America and the Caribbean. FAO Fisheries and Aquaculture Technical Paper No 544 Rome, FAO, pp 1–12

[63] Schärer-UmpierreMT, Mateos-MolinaD, AppeldoornR, BejaranoI, Hernández-DelgadoEA, NemethRS, Valdés-PizziniM, SmithTB. Marine managed areas and associated fisheries in the US Caribbean. Adv Mar Biol, 2014, 69: 129-152.

[64] SEDAR26 (2011) SEDAR 26 Stock assessment report U.S. Caribbean Queen Snapper. Charleston, SC, USA, South Atlantic Fishery Management Council. https://sedarweb.org/documents/sedar-26-caribbean-queen-snapper-final-stock-assessment-report/

[65] SEDAR4 (2003) Data Report of SEDAR 4, Atlantic and Caribbean Deepwater Snapper-Grouper, Caribbean Species. Charleston, SC, USA, South Atlantic Fishery Management Council. https://sedarweb.org/documents/sedar-4-caribbean-data-workshop-report/

[66] SelkoeKA, AloiaCC, CrandallED, IaccheiM, LigginsL, PuritzJB, von der HeydenS, ToonenRJ. A decade of seascape genetics: contributions to basic and applied marine connectivity. Mar Ecol Prog Ser, 2016, 554: 1-9.

[67] ShermanKD, ParisJR, KingRA, MooreKA, DahlgrenCP, KnowlesLC, StumpK, TylerCR, StevensJR. RAD-seq analysis and in situ monitoring of nassau grouper reveal fine-scale population structure and origins of aggregating fish. Front Mar Sci, 2020, 7157.

[68] SwearerSE, CaselleJE, LeaDW, WarnerRR. Larval retention and recruitment in an island population of a coral-reef fish. Nature, 1999, 402: 799-802.

[69] TaylorMS, HellbergME. Comparative phylogeography in a genus of coral reef fishes: biogeographic and genetic concordance in the Caribbean. Mol Ecol, 2006, 15: 695-707.

[70] Tonioli FC, Agar JJ (2011) Synopsis of Puerto Rican commercial fisheries. NOAA Technical Memorandum NMFS-SEFSC-622:69. https://repository.library.noaa.gov/view/noaa/24765

[71] Tovar VerbaJ, StowA, BeinB, PenninoMG, LopesPF, FerreiraBP, MortierM, Queiroz LimaSM, PereiraRJ. Low population genetic structure is consistent with high habitat connectivity in a commercially important fish species (Lutjanus jocu). Mar Biol, 2023, 1705.

[72] VineeshN, ShihabI, AkhileshKV, SajeelaKA, MukthaM, GopalakrishnanA. Is Lutjanus argentimaculatus genetically connected along the Arabian Sea and Bay of Bengal? A study using microsatellite markers. Aquaculture, 2023, 562738863e.

[73] WeirBS, CockerhamCC. Estimating F-statistics for the analysis of population structure. Evolution, 1984, 1: 1358-1370

[74] WeistP, JentoftS, TørresenOK, SchadeFM, PampoulieC, KrummeU, HanelR. The role of genomic signatures of directional selection and demographic history in the population structure of a marine teleost with high gene flow. Ecol Evol, 2022, 12. e9602

[75] Williams SM, Prada C, Beltrán DM (2022) Characterization of Prey Diversity of the Commercially Important Queen Snapper (Cartucho), Etelis oculatus.https://www.caribbeanfmc.com/Briefing_Books/178th_Meeting/CFMC_finalreport_v10_april14_2.pdf

[76] WilliamsSM. Prey diversity in the deep ocean: metabarcoding feeding ecology of the commercially important queen snapper in the US Caribbean. Front Mar Sci, 2024, 111409336.

[77] WillisSC, HollenbeckCM, PuritzJB, PortnoyDS. Genetic recruitment patterns are patchy and spatiotemporally unpredictable in a deep-water snapper (Lutjanus vivanus) sampled in fished and protected areas of western Puerto Rico. Conserv Genet, 2022, 3: 435-447.

[78] WilsonGA, RannalaB. Bayesian inference of recent migration rates using multilocus genotypes. Genetics, 2003, 163: 1177-1191.

[79] Xu H (2022) Oceanographic pathways and hydrodynamic connectivity between marine protected areas in the US Virgin islands and eastern Puerto Rico. Ph.D. Dissertation, University of Puerto Rico at Mayagüez https://hdl.handle.net/20.500.11801/2989

Funding

Museo Nacional de Ciencias Naturales (MNCN)

RIGHTS & PERMISSIONS

The Author(s)