Population structure of long-snout seahorse Hippocampus reidi in Southwestern Atlantic and implications for management

Abstract

Hippocampus reidi represents the most abundant species of the genus Hippocampus along the Brazilian coast. Despite being charismatic, the species is globally threatened due to habitat degradation and commercial exploration, especially in Brazil, which is the leader in exportation and consumption of the species. Through mitochondrial (cytochrome b and control region) and nuclear (1st intron S7) data, the current study investigates the variation and genetic structure of H. reidi along the Brazilian coast, from Pará to Santa Catarina states. The mitochondrial data indicate the presence of two lineages: (1) North/Northeast and (2) South/Southeast, which was partially recovered by nuclear data. This scenario could be related to temperature differences and circulation patterns of the Brazil and North-Brazil currents, which define these groups into biogeographic sub-provinces. The lineages occur in sympatry in Bahia state, which can be explained by the occurrence of secondary contact during the last glacial maximum. Despite presenting two lineages, for management and conservation, three units are indicated: (1) North/Northeast, (2) Bahia, and (3) South/Southeast. The North/Northeast unit proved to be more vulnerable, presenting the lowest genetic diversity indices, representing a priority for future conservation actions.

Keywords: Genetic diversity; Gene flow; Management units; Marine conservation; Secondary contact

Resumo

Hippocampus reidi representa a espécie mais abundante do gênero na costa brasileira. Apesar de carismáticos, encontram-se globalmente ameaçados devido à degradação de habitat e intensa exploração comercial, especialmente no Brasil, líder na exportação e consumo da espécie. Através de dados mitocondriais (citocromo b e região controle) e nucleares (1st íntron S7), este estudo investigou a variação e estrutura genética de H. reidi em toda a costa brasileira, do estado do Pará até Santa Catarina. Os dados mitocondriais indicam a existência de duas linhagens de H. reidi na costa brasileira: (1) Norte/Nordeste e (2) Sudeste/Sul, padrão parcialmente recuperado pelos dados nucleares. Este cenário pode ser explicado por diferenças na temperatura e padrões de circulação das correntes do Brasil e Norte do Brasil, que definem estes grupos como subprovíncias biogeográficas. As linhagens ocorrem em simpatria no estado da Bahia, o que pode ser explicado pela ocorrência de contato secundário durante o último glacial máximo. Apesar de apresentar duas linhagens, para fins de manejo e conservação, são indicadas três unidades: (1) Norte/Nordeste, (2) Bahia, e (3) Sudeste/Sul. A unidade do Norte/Nordeste mostrou-se a mais vulnerável devido aos baixos índices de diversidade genética apresentados e representa uma prioridade para futuras ações de conservação.

Palavras chave: Contato secundário; Conservação marinha; Diversidade genética; Fluxo gênico; Unidades de Manejo

INTRODUCTION

Anthropic pressures, such as pollution, habitat degradation and loss, climatic changes, and, especially, harvesting, are the major agents of marine defaunation, reducing the effective population size (Pan et al., 2013; Pinsky, Palumbi, 2014; Allendorf et al., 2014; McCauley et al., 2015; Martínez-Candelas et al., 2020). Since small populations are more susceptible to genetic drift, compromising genetic diversity and evolutive potential, which can lead several species to extinction, understanding the genetic diversity distribution along the adaptive landscape is crucial for successful management plans (Allendorf et al., 2014; Cadrin et al., 2014; Cadrin, 2020).

Due to the supposed absence of clear gene flow barriers in marine environments, the management of marine species is generally based on panmixia (Cowen et al., 2006). In addition, a large number of marine species have high dispersion capacity, through both active, such as adult migrations, and passive pathways, during the larval phase, which is strongly associated with genetic homogeneity (Palumbi, 1992; Selkoe et al., 2014, 2016). However, ocean currents, convergence zones, and oceanic gyres, as well as differences in temperature, salinity, philopatry, and historical phenomena, can promote isolation, reducing the gene flow, and leading to diversification events (Grant, Bowen, 1998; Cowen et al., 2000; Machado-Schiaffino et al., 2010; Laurrabaquio-A et al., 2019; Lehnert et al., 2019; Chen et al., 2020; Faria et al., 2020; McKeown et al., 2020; Torrado et al., 2020; Zhao et al., 2020; Sadeghi et al., 2021). To better characterize different management units and identify barriers between species and populations, molecular tools are effective (e.g., Haney et al., 2010; da Silva et al., 2016; Healey et al., 2018; Azpelicueta et al., 2019; Jacobina et al., 2020; Klanten et al., 2020; Neves et al., 2020; Andrade et al., 2021; Labrador et al., 2021; Zarei et al., 2021), especially when using both mitochondrial and nuclear data, to more closely establish the species history, obtaining more robust results.

In this context, the genus Hippocampus (Syngnathidae), commonly known as seahorses, stands out. The relevance of the current study is related to observed decreases in several seahorse populations, especially linked to overfishing and habitat loss (Foster, Vincent, 2004). Currently, most seahorse species are listed in the International Union for Conservation of Nature (IUCN) red list, and the entire Hippocampus genus is listed in Appendix II of the Convention on the International Trade in Endangered Species of Wild Fauna and Flora (CITES). Until 2014, the occurrence of two species was confirmed along the Brazilian coast: H. reidi Ginsburg, 1933, and H. erectus Perry, 1810. However, using both molecular and morphological methods, Silveira et al. (2014) provided evidence of the presence of H. patagonicus Piacentino & Luzzatto, 2004. These three species are classified as Vulnerable (VU) in the Brazilian Red List of the Brazilian Ministry of Environment (ICMBio, 2018), and the planning for the species conservation was recently established at a meeting of experts promoted by IUCN in the Northeastern Brazil (RAT, 2024, participation).

Among the Hippocampus species, the long-snout seahorse H. reidi is the most abundant species, inhabiting coastal and shallow waters, in the Atlantic Ocean, from North Carolina (United States), Gulf of Mexico and Caribbean Sea, to Southeast Brazil (Lourie et al., 1999; Musick et al., 2000; Hercos, Giarrizzo, 2007; Silveira, 2011; Silveira et al., 2014). As it is the most common seahorse in Brazilian estuaries, H. reidi is a tourist attraction, which includes direct interactions with the humans (Ternes et al., 2016). Some life history strategies make the long-snout seahorse extremely vulnerable to these practices, as well as to harvesting and habitat loss, such as the formation of stable reproductive pairs, low mobility, patchy distribution, small home ranges, low reproductive rates, and a long period of parental care. For these reasons, this species is globally classified as Near Threatened (NT) by IUCN (Foster, Vincent, 2004; Oliveira, Pollom, 2017).

In the Brazilian Northeast, nuclear genetics of the Inter Simple Sequence Repeat (ISSR) showed a structure pattern between Maracaípe (Pernambuco) and Jericoacora (Ceará) (Montes et al., 2018). In addition, Carmo et al. (2022) showed a stable population structure between two bays of Rio de Janeiro state, and no differences between seasons. However, despite the species being charismatic, crucial information for conservation plans (e.g., population/genetic structure, genetic diversity, connectivity, and number of management units) remains poorly explored, especially in Brazil, which is considered the largest seahorse exporter for the international aquarium trade in Latin America as well as the largest consumer market (Baum, Vincent, 2005; Rosa, 2005),

Thus, the present study aimed to characterize the H. reidi genetic structure in the Southwest Atlanticthrough a multiloci approach, using both mitochondrial (Cytochrome b and control region) and nuclear (first intron of the ribosomal protein S7) data. More specifically, we asked how many H. reidi management units exist along the Brazilian coast and if there are any priority areas for conservation, using genetic diversity, population genetics, and demographic parameters. These informations can and will be directly used for the management plan of long-snout seahorses in Brazil.

MATERIAL AND METHODS

Sample collection and molecular procedures. In total, 362 tissue samples of H. reidi were obtained from trade (dried specimens; Fig. S1) or collected from individuals found in situ, along the Brazilian Atlantic coast, in estuarine and reef environments, from Pará to Santa Catarina states (Fig. 1), between 2012 and 2015. For the in situ sampling, the non-destructive methods of collecting the dorsal fin (fin-clipping) or dermal filaments were used (Lourie, 2003; Planas et al., 2007). As it is a Near threatened (NT) species, vouchers were not collected. All individuals were photographed (Fig. 2) and returned to the same location. For trade seahorses, tissue was removed from the base of the tail for DNA analysis. All tissue samples were stored in 96% ethanol and kept at -20 ºC.All details about the samples, including geographical coordinates, can be found in supplementary material Tab. S2.

The total genomic DNA was extracted using the DNeasy Blood and Tissue (Qiagen®) kit, following the protocol suggested by the manufacturer. The extracted DNA was visualized in 1% electrophoresis gel and stained with GelredTM, and posteriorly quantified using a nano spectrophotometer Nanodrop 2000 (Thermo Scientific).

FIGURE 1 |
Sampling locations of Hippocampus reidi along the Brazilian coast. PA = Pará; PI = Piauí; CE = Ceará; RN = Rio Grande do Norte; PB = Paraíba; PE = Pernambuco; AL = Alagoas; BA = Bahia; ES = Espírito Santo; RJ = Rio de Janeiro; SP = São Paulo; SC = Santa Catarina.

FIGURE 2 |
Specimens of the long-snout seahorse Hippocampus reidi collected along the Brazilian coast. A. Pará, B. Piauí, C. Ceará, D. Rio Grande do Norte, E. Paraíba, F. Pernambuco, G. Alagoas, H. Bahia, I. Espírito Santo, J. Rio de Janeiro, K. São Paulo, L. Santa Catarina.

Three regions were amplified via PCR: the Cytochrome b gene (Cytb) and the control region (CR), representing the mitochondrial genome, and the first intron of the ribosomal protein S7 (S7), representing the nuclear genome. The Cytb fragment was amplified using shf2 and shr2 primers (Lourie et al., 2005), following the cycle presented by the authors. The CR region was amplified using the HCAL2 and HCAH2 primers (Teske et al., 2003) following the cycle presented by the authors. For each region, the PCR reactions were carried out in a total volume of 25 µL, using: 12.5 µl of 2X Taq Master Mix (Vivantis®) (1.25U of Polymerase Taq, 1X of buffer, 0.2 mM of dNTPs and 1.5 mM of MgCl2), 0.5 µl of MgCl2 (50 µM), 1.0 µl of each primer (10 µM), 2.5 µl of DNA (2-10 ng/µl), and 8 µl of ultrapure water. For the S7, the universal primers S7RPEX1F and S7RPEX2R were used (Chow, Hazama, 1998), following the amplification protocol described by Teske et al. (2004): 94 °C for 5 min, followed by 35 cycles of 30 s at 94 °C, 1 min at 60 °C, 1 min at 72 °C, and a final extension for 10 min at 72 ºC. The PCR reaction was carried out in a total volume of 25 µL using: 12.5 µl of 2X Taq Master Mix (Vivantis®) (11.25U of Polymerase Taq, 1X of buffer, 0.2 mM of dNTPs and 1.5 mM of MgCl2), 0.5 µl of MgCl2, 0.5 µl of each primer (10 µM/µL), 3 µl of DNA (2–10 ng/ µl), and 8 µl of ultrapure water.

The PCR products were purified using the ExoSAP-IT kit (QIAquick® PCR Purification Kit), following the protocol suggested by the manufacturer. The sequencing was carried out in a forward direction using the Bigdye Terminator v. 3.1 cycle Sequencing Ready Reaction kit (Applied Biosystems), in an automatic sequencer ABI 3500-Applied Biosystems. All sequences were deposited in the GenBank database under accession codes PQ134117 - PQ134478 (Cytb), PQ127407 - PQ127768 (CR), and PQ130575 - PQ131134 (S7).

Data analysis. All sequences were edited and aligned using the ClustalW algorithm (Thompson et al., 1994) implemented in BioEdit Sequence Alignment Editor v. 7.0. (Hall, 1999). The mitochondrial fragments (Cytb and CR) were concatenated as a single non-recombinant marker and, from now on, will be called mtDNA. Due to the presence of polymorphic sites in diploid regions, such as intron S7, the alleles were reconstructed with the PHASE v. 2.1 tool (Stephens et al., 2004) implemented in DnaSP v. 6.0 (Librado, Rozas, 2009), using default parameters and considering only allelic states with probabilities higher than 70% (Stephens et al., 2004). The analyses described below were performed for each marker (mtDNA and intron S7) separately. The genetic diversity indices [number of haplotypes (H) and polymorphic sites (S), private haplotypes (%Hp), haplotype (h) and nucleotide (π) diversity] was assessed by sample sites and cluster/genetic groups identified by other analysis, using the software Arlequin v. 3.5 (Excoffier, Lischer, 2010).

The population structure of H. reidi was tested through three approaches. First, the genealogical relationship between the haplotypes and their sample sites was investigated through a haplotype network using the TCS method in PopART (Clement et al., 2002; Leigh, Bryant, 2015). Second, the population structure of H. reidi was tested by the Bayesian Analysis of Population Structure – BAPS v. 6.0 (BAPS; Corander, Marttinen, 2006; Corander et al., 2008), firstly using a genetic mixture analysis with the sequences, followed by a population admixture analysis with a total of 10,000 interactions.

Lastly, to directly associate the geographic information with the genetic structure, the Geneland package (Guillot et al., 2005) on the R platform (http://www.R-project.org) was used. Due to limitations in the analysis regarding the sample size, only 300 individuals were included (Tab. S2). The number of groups (k) analyzed was 1–12 groups, with 9 independent runs, 1 million Markov chain interactions (MCMC), and a Thinning value = 1000.

The phylogenetic relationships of H. reidi were reconstructed through Bayesian Inference topologies in MrBayes v. 3.2.6 (Ronquist, Huelsenbeck, 2003). The nucleotide substitution model was estimated in jModelTest v. 2.1.7 (Darriba et al., 2012) under the Akaike Information Criterion [HKY+I (mtDNA) and HKY+G (intron S7)]. GenBank sequences of H. erectus [NC_022722.1 and KF557652.1 (mtDNA); KX646492.1 and KX646493.1 nuDNA)], H. trimaculatus Leach, 1814 [MF579378.1 and MF579379.1 (nuDNA)], Syngnathus typhle Linnaeus, 1758 [NC_030279.1 and KU925872.1 (mtDNA)], S. schlegeli Kaup, 1853 [AP012318.1 and KP861226.1 (mtDNA)] and S. temminckii Kaup, 1856 [AY277308.1 (nuDNA)] were used as external groups. Each database was analyzed with a burn-in of 10% and 10 million MCMC.

The genetic differentiation was assessed through the pairwise FST, among all sample sites and genetic clusters in Arlequin v. 3.5 (Excoffier, Lischer, 2010) using 1,000 permutations (p < 0.05). The Analysis of Molecular Variance (AMOVA) tested different hypotheses based on the population structure results for mitochondrial data: (a) Null hypotheses of one single group (Brazil); (b) Two groups, without the Bahia samples: (1) North/Northeast (NNE – PA, PI, CE, RN, PB, PE, AL) and (2) South/Southeast (SES – ES, RJ, SP, SC); (c) Two groups, with Bahia samples in the NNE group: (1) NNE + BA and (2) SES; (d) Two groups, with Bahia samples in the SES group: (1) NNE and (2) SES+BA; (e) Three groups: NNE (PA, PI, CE, RN, PB, PE, AL), (2) BA and (3) SES (ES, RJ, SP, SC).

Oscillations in population size were investigated through the Bayesian Skyline Plot analysis (BSP; Drummond et al., 2005) in Beast v. 2.6.4 (Bouckaert et al., 2014). Based on the population structure results, three groups were defined: (1) NNE, (2) BA, and (3) SES. For each mitochondrial marker, the nucleotide substitution model was estimated in jModelTest v. 2.1.7 (Darriba et al., 2012) under the Akaike Information Criterion, selecting the HKY model, except for control region data from NNE (JC model) and Cytb from BA (HKY+G model). In Beauti, both markers were linked. The mutational rates used for calibration were 1 x 10-8 per site per year (Mobley et al., 2010) and 5 x 10-8 per site per year (Bowen et al., 2006), for Cytb and CR, respectively. Three independent runs of 10 million (NNE and SSE groups) and 15 million (BA group) MCMC were performed with a burn-in of 25%. The log and tree files were combined using the LogCombiner tool implemented on Beast v. 2.4.6, and the parameters convergence (ESS>200; Effective Sample Size) was checked on Tracer v. 1.7.1 (Rambaut et al., 2014). In addition, the traditional neutrality tests [Fu’s FS (Fu, 1997) and Tajima’s D (Tajima, 1989)] in Arlequin v. 3.5. (Excoffier, Lischer, 2010), and the Mismatch Distribution analysis in DNAsp v. 5.1 (Rogers, Harpending, 1992; Librado, Rozas, 2009) were performed as complementary approaches.

Voucher specimens. As the long-snout seahorse Hippocampus reidi is globally classified as Near Threatened (NT), and Vulnerable (VU) in Brazil, vouchers were not collected. Thus, all individuals were photographed and, after collecting tissue samples, they were returned to the same location.

RESULTS

For mitochondrial data, concatenated fragments of 1,162 bp from 362 H. reidi individuals defined 69 haplotypes, of which 54 are private to some sample sites. The genetic diversity ranged from 0.308 (PB) to 0.952 (BA and RJ), for haplotype, and from 0.0007 (RN and PB) to 0.0059 (BA), for nucleotide diversity (Tab. 1).

TABLE 1 |
Molecular parameters of Hippocampus reidi populations and lineages/groups along the Brazilian coast. **In nuclear marker, the N number corresponds to alleles; PA = Pará; PI = Piauí; CE = Ceará; RN = Rio Grande do Norte; PB = Paraíba; PE = Pernambuco; AL = Alagoas; BA = Bahia; ES = Espírito Santo; RJ = Rio de Janeiro; SP = São Paulo; SC = Santa Catarina; NNE = Lineage formed by Northeast samples, except Bahia samples; BA = Group formed by Bahia samples; SSA = Lineage formed by South and Southeast samples, except Bahia samples; N = Sample size; H = Haplotype number; Hp = Private haplotypes (percentage in parenthesis); S = Polymorphic sites; h = Haplotype diversity; π = Nucleotide diversity; *significant P values (p < 0.05).

Marker	Location	N	H	Hp	S	h	π	Fs de Fu	D de Tajima
mtDNA	PA PI CE RN PB PE AL BA ES RJ SP SC NNE BA SSE Total	9 30 41 30 30 32 40 47 40 25 7 31 212 47 103 362	3 7 8 7 5 8 6 28 6 14 6 11 17 28 30 69	0 1 (1.45%) 1 (1.45%) 1 (1.45%) 2 (2.9%) 1 (1.45%) 3 (4.35%) 23 (33.3%) 4 (5.8%) 9 (13.04%) 1 (1.45%) 8 (11.6%) 15 (21.7%) 23 (33.3%) 25 (36.2%) 54 (78.2%)	3 15 8 7 8 7 10 45 17 28 18 21 29 45 54 94	0.639 0.462 0.685 0.462 0.308 0.688 0.359 0.952 0.391 0.92 0.952 0.843 0.516 0.952 0.871 0.814	0.0012 0.0011 0.0013 0.0007 0.0007 0.0012 0.0007 0.0059 0.0015 0.0034 0.0049 0.0025 0.00098 0.0059 0.0029 0.0048	1.176 -1.826 -1.896 -3.527* -1.309 -2.607 -2.063 -11.742* 0.2927 -4.903* -0.9 -2.287 -9.9* -11.742* -15.976* -24.735*	0.794 -2.21* -0.65 -1.598* -1.876* -0.658 -1.996* -1.148 -1.868* -1.784* -1.313 -1.538 -2.158* -1.148 -2.135* -1.835*
S7** (nuDNA)	PA PI CE RN PB PE AL BA ES RJ SP SC Total	14 60 76 60 58 80 78 80 82 42 12 40 682	6 5 5 4 4 5 5 5 3 7 6 3 11	0 0 0 0 0 0 1 0 0 0 0 0 1	6 6 6 6 6 6 6 5 4 4 4 4 6	0.747 0.648 0.607 0.633 0.67 0.634 0.649 0.634 0.576 0.684 0.879 0.668 0.684	0.0043 0.0044 0.0042 0.0044 0.0046 0.0043 0.0044 0.0034 0.0034 0.0037 0.0033 0.0032 0.0042	-0.681 2.93 3.053 4.23 4.412 3.289 3.441 2.187 4.998 -0.0304 -1.748 3.608 1.234	0.655 1.898 1.791 1.905 2.089 1.948 2.101 1.611 2.448 2.552 1.029 1.815 2.742

The intron S7 data recovered the alleles from 341 H. reidi individuals, of which 120 was homozygote and 221 was heterozygote, with fragments of 520 bp, defining 11 haplotypes. The genetic diversity ranged from 0.576 (ES) to 0.879 (SP), for haplotype, and form 0.0032 (SC) to 0.0046 (PB), for nucleotide diversity (Tab. 1).

The mitochondrial haplotype network recovered two lineages separated by 6 mutational steps. In general, these groups correspond to the North/Northeast (NNE) and South/Southeast (SES) regions, and the Bahia state presented both lineages. Furthermore, one haplotype from Piauí state grouped with the SES samples (Fig. 3A). The individual results of each mitochondrial marker recovered a similar pattern and can be found in Fig. S3. The nuclear data from intron S7 did not recover any population structure pattern, with all haplotypes being shared among all sample sites, except for the Hap6, exclusive from Alagoas state (Fig. 3B).

FIGURE 3 |
Haplotype networks based on the TCS method generated in PopART software of Hippocampus reidi for mitochondrial data (A), and nuclear data (B). The circles represent the haplotypes and different colors represent the sampling locations. Lines between haplotypes represent the mutation steps and black circles are missing or unidentified haplotypes. PA = Pará; PI = Piauí; CE = Ceará; RN = Rio Grande do Norte; PB = Paraíba; PE = Pernambuco; AL = Alagoas; BA = Bahia; ES = Espírito Santo; RJ = Rio de Janeiro; SP = São Paulo; SC = Santa Catarina.

FIGURE 4 |
Bayesian Analysis of Population Structure BAPS of Hippocampus reidi. A. Mitochondrial, B. Nuclear data. PA = Pará; PI = Piauí; CE = Ceará; RN = Rio Grande do Norte; PB = Paraíba; PE = Pernambuco; AL = Alagoas; BA = Bahia; ES = Espírito Santo; RJ = Rio de Janeiro; SP = São Paulo; SC = Santa Catarina.

The BAPS analysis identified 3 genetic profiles (k = 3; p = 1) for the mtDNA data. The green and blue/red profiles are almost exclusively of the sample sites from the NNE and SES groups, respectively. The BA group presented the 3 genetic profiles (Fig. 4A). The nuclear data recovered 8 genetic profiles (k = 8; p = 1), shared, in general, among all sample sites, in different frequencies, except by the blue/green, which is exclusive of the samples from the NNE group, and pink, which is exclusive of SES samples (Fig. 4B). These three genetic profiles are present in the BA group.

For the mtDNA, the Geneland data showed a higher probability of K = 4: Cluster 1 – Bahia, Cluster 2 –Espírito Santo, Cluster 3 –Rio de Janeiro, São Paulo, and Santa Catarina; and Cluster 4 –Pará, Piauí, Ceará, Rio Grande do Norte, Paraíba, Pernambuco, and Alagoas (Figs. 5A–E). The nuclear data recovered a similar result, with a higher probability of K = 3: Cluster 1 –Bahia and Espírito Santo, Cluster 2 –Rio de Janeiro, São Paulo, and Santa Catarina, Cluster 3 –Pará, Piauí, Ceará, Rio Grande do Norte, Paraíba, Pernambuco, and Alagoas (Figs. 5F–I).

The mitochondrial Bayesian topology sustains H. reidi from the Brazilian coast as a monophyletic group, also recovering the two lineages identified by the haplotype network. The monophyletic reciprocity between these groups presented a high branch support (posterior probability > 0.9; Fig. S4). The nuclear data failed to recover any clades (Fig. S5).

FIGURE 5 |
Posterior probability maps generated by Geneland analysis for mitochondrial (A–E) and nuclear (F–I) data of Hippocampus reidi. White tones indicate a greater probability of the samples (black circles) forming a particular population. A. Probability graph of densities obtained for the possible ‘K’ genetic populations for mithocondrial data, B. Cluster 1: Bahia, C. Cluster 2: Espírito Santo, D. Cluster 3: Rio de Janeiro, São Paulo, Santa Catarina, E. Cluster 4:Pará, Piauí, Ceará, Rio Grande do Norte, Paraíba, Pernambuco, Alagoas, F. Probability graph of densities obtained for the possible ‘K’ genetic populations for nuclear data, G. Cluster 1:Bahia, Espírito Santo, H. Cluster 2: Rio de Janeiro, São Paulo, Santa Catarina, I. Cluster 3:Pará, Piauí, Ceará, Rio Grande do Norte, Paraíba, Pernambuco, Alagoas.

For mitochondrial data, the significant pairwise FST ranged from moderate [PE vs. PI (FST = 0.045)] to very high [ES vs. AL (FST = 0.862)]. The ES samples presented high and significant values in all comparisons, ranging from 0.237, when compared to BA samples, to 0.862 when compared to AL samples (Fig. 6A). In addition, the BA samples presented high and significant values in all comparisons, except when compared to SP samples (FST = 0.053; p > 0.05) (Fig. 6A). The individual results of each mitochondrial marker recovered a similar pattern and can be found in Tab. S6. The nuclear data revealed a similar scenario, with significant pairwise FST ranging from low [CE vs. AL (FST = 0.001; p > 0.05)] to high [BA vs. SC (FST = 0.198; p > 0.05)]. The ES and BA samples presented moderate-high values in all comparisons, except when compared to SP and CE, respectively (Fig. 6B).

FIGURE 6 |
Heat map of pairwise FST values among each sample site of Hippocampus reidi represented in both x and y axis. A. Mitochondrial, B. Nuclear data. The asterisk represents significant values (p < 0.05). PA = Pará; PI = Piauí; CE = Ceará; RN = Rio Grande do Norte; PB = Paraíba; PE = Pernambuco; AL = Alagoas; BA = Bahia; ES = Espírito Santo; RJ = Rio de Janeiro; SP = São Paulo; SC = Santa Catarina.

Considering the null hypothesis, the AMOVA presented a high and significant FST value for mitochondrial (FST = 0.59; p < 0.05), and low for nuclear data (FST = 0.04; p < 0.05) (Tab. 2). In the different scenarios tested for mitochondrial data, based on the results of the population structure analyses, the highest genetic differentiation between groups (FCT) was found when Bahia (BA) samples were excluded, considering the North/Northeast (NNE) and South/Southeast (SSE) as distinct groups (FCT = 0.79; p < 0.05). Similar FCT values were found considering the following groups: NNE, BA, ES, SES (FCT = 0.69; p < 0.05), NNE, SES+BA (FCT = 0.67; p < 0.05), and NNE, BA+ES, SES (FCT = 0.67; p < 0.05) (Tab. 2).

The demographic analysis of BSP indicates expansion events at 5 and 20 thousand years ago for the NNE and BA groups, respectively (Fig. 7). The SES group seems to have suffered a recent contraction in population size (ca. 5 thousand years ago; Fig. 7).

The neutrality tests were negative and simultaneously significant only for mitochondrial data for total BR samples, RN and RJ (Tab. 1). The mismatch distribution analysis for both NNE and SES presented a unimodal curve, and for BA a bimodal curve (Fig. S7).

TABLE 2 |
AMOVA of Hippocampus reidi testing different hypothesis by mitochondrial data. *Significant p values (p < 0.05). BR: represents all sample sites in one single group; NNE represents PA, PI, CE, RN, PB, PE, AL populations; SES represents ES, RJ, SP, SC populations; BA represents Bahia population; ES represents Espírito Santo population.

Hypothesis		BR		1)NNE 2)SES	1)NNE+BA 2)SES	1)NNE 2)SES+BA	1)NNE 2)BA 3)SES	1)NNE 2)BA 3)ES 4)SES	1)NNE 2)BA+ES 3)SES
Variation Source (%)	Between groups Between populations Within populations	nuDNA - 4.02 95.98	mtDNA - 59.98 40.02	78.66 2.95 18.39	60.54 12.76 26.7	66.96 5.43 27.6	67.16 3.08 29.76	69.07 0.28 30.65	66.54 3.32 30.15
Fixation indices	FSC FST FCT	- 0.04* -	- 0.59* -	0.14* 0.82* 0.79*	0.32* 0.73* 0.61*	0.16* 0.72* 0.67*	0.094* 0.7* 0.67*	0.009* 0.7* 0.69*	0.099* 0.7* 0.67*

FIGURE 7 |
Bayesian Skyline plot reconstructions for mitochondrial data of Hippocampus reidi. The x axis represents the time in years, and the y axis the effective population size (Ne). A. NNE (North/Northeast), B. BA (Bahia), C. SES (South/Southeast).

DISCUSSION

Population structure, genetic diversity and demographic parameters ofH. reidi. The concatenated mitochondrial data reveal the presence of two H. reidi lineages in Brazilian coast. The first one represents the North and Northeast (NNE) samples and the other represents the South and Southeast (SES) samples. In addition, these lineages occur in sympatry in Bahia state, which seems to be a contact zone between them. Despite detecting some signals of structure (e.g., population differentiation in BA and ES population by pairwise FST and BAPS; three genetic clusters by Geneland), the nuclear data did not recover the two lineages identified by mitochondrial data.

These incongruences could be related to mutational rate differences between the molecular markers, since the mtDNA has higher mutational rates than nuDNA (Zink, Barrowclough, 2008; Calcagnotto, 2012; Toews, Brelsford, 2012). In addition, the subtle genetic differentiation in nuDNA could be influenced by male movements, acting as gene flow mediators, allowing a higher miscegenation between NNE and SES (Hellberg et al., 2002; Murray et al., 2017; Green et al., 2018; Day et al., 2019; Roycroft et al., 2019). However, more accurate data about the movement patterns of H. reidi are necessary in order to explain why males juveniles possibly disperse more, as suggested by nuDNA data observed herein.

The genetic structure in two lineages is reinforced by several mtDNA analyses (haplotype network, Bayesian topology, BAPS, and pairwise FST). These lineages appear to be distributed in three groups: (1) North/Northeast (except the BA samples), representing Lineage A, (b) South/Southeast, representing Lineage B, and (c) Bahia (BA), representing a mixture of both lineages, suggesting being a contact zone between them. This hypothesis is supported by AMOVA (FCT = 0.67)and Geneland results. However, it is important to highlight that the BA samples are genetically closer to the South/Southeast samples, presenting lowest values of pairwise FST. The absence of shared haplotypes between NNE and SSE may indicate a reduced gene flow.

This structure pattern can be partially explained by the Isolation by Distance (IBD) hypothesis, since a positive correlation between the genetic differentiation and geographic distance was found (unpublished data; Defavari, 2016). The IBD can be related to long-snout seahorse life history strategies, such as the closer relationship with the estuarine/mangrove environments (Lourie et al., 1999). However, two geographically closer populations presented a moderate and significant pairwise FST (PA vs. PI; FST = 0.134), while populations separated by 2,000 km presented a negative FST (PE vs. PA; FST = -0.023). In this way, the genetic similarity between populations from the same group can be explained by the high dispersion of H. reidi during the larval or juvenile phases as previously observed (Foster, Vincent, 2004; Lourie et al., 2004). The absence of clear physical gene flow barriers in the marine environment could also facilitates this phenomenon (Cowen et al., 2000, 2006), and the broad-scalehomogeneity occurs in several marine species with planktotrophic larvae (Hellberg et al., 2002). Thus, this pattern can be used to explain the absence of H. reidi genetic structure in several areas along the Brazilian coast.

The NNE and SSE represent different biogeographic sub-provinces of the Brazilian province (Pinheiro et al., 2018), and their sample sites represent different marine Ecoregions (Spalding et al., 2007). This geographic division seems to be deeply related with the split of the South Equatorial current in the North-Brazil and Brazil currents, species composition, and the São Francisco River mouth. It is important to highlight that, despite being the most common seahorse species along the Brazilian coast, presenting tolerance to soft changes in the salinity levels (Tseng et al., 2020), abrupt changes in salinity levels due to discharge of freshwater can affect the survival of the H. reidi individuals (da Hora et al., 2016) and can act as gene flow barrier.In addition, while the NNE presents warmer waters, the SES presents colder waters, and differences in temperature can reduce the gene flow, allowing local adaptation and isolation (Santos et al., 2003; Cunha et al., 2014). Thus, these temperatures gradients could be able to explain the genetic pattern found in H. reidi.

For both mitochondrial and nuclear data, the major pairwise differentiation was related to two sample sites: Espírito Santo (ES) and Bahia (BA). The ES presented only two shared haplotypes by mtDNA data [with BA (Hap32) and with BA and SP (Hap27)], revealing possible gene flow loss. Despite presenting significant FST values when compared with the SES samples, the ES samples is genetically closer to this group. These samples were collected below the Doce River mouth, and freshwater discharge as this one can limit the seahorses’ movements in NNE direction, as argued above. Additionally, in ES, the continental platform is narrow and contains a submersed mountain chain, the Vitória-Trindade (VTC), which could have acted as a glacial refuge during the Pleistocene (Pinheiro et al., 2015), and is being associated with genetic differentiation of other marine species (e.g., Santos et al., 2006; Pinheiro et al., 2015; dos Santos Freitas et al., 2017; Nauer et al., 2019; Souza et al., 2019). Thus, these features could favor the reduction of the gene flow between SES and NNE longsnout seahorses.

The BA samples presented high differentiation levels by pairwise FST, including when compared to NNE group. Of 28 haplotypes, only five are shared. However, it is genetically closer to SES group. Pinheiro et al. (2018) grouped the Bahia state into the same sub-province that SES populations, and Spalding et al. (2007) considered the BA state as a different marine Ecoregion, grouped with the ES into Eastern Brazil, which can explain the genetic similarity found between them. Despite that, BA presented both H. reidi lineages, reinforcing the idea of a contact zone mentioned above.

This sympatry can be explained by the demographic expansion that occurred after the Last Glacial Maximum (LGM; ca. 15 thousand years ago). The use of different refugia followed by gene flow during the glacial (sea level retraction) and interglacial (sea level expansion) cycles, respectively, has been associated with marine species diversification (e.g., Chen et al., 2020; Neves et al., 2020). These phenomena can be potentialized in areas with a narrow continental shelf (Dolby et al., 2016, 2018), such as Bahia, which presents the narrowest continental shelf on the Brazilian coast (ca. 14 km wide; Dominguez et al., 2012). Thus, after the LGM, the two H. reidi lineages may have had secondary contact in the Bahia coast, which is reinforced by the bimodal pattern of the mismatch distribution analysis, suggesting two episodes of population expansion.

The high degrees of genetic diversity were similar to those found for other Hippocampus species (Lourie et al., 2004, 2005; Goswami et al., 2009; Panithanarak et al., 2010; Saarman et al., 2010). Despite the highest sample size, the NNE group presented the lowest diversity level, and signals of population size contractions were identified by mismatch distribution analysis. However, the BSP did not recover any contraction events; on the contrary, indicated a recent populational expansion (ca. 5 thousand years ago). The SES group presented an opposite scenario. Despite the high genetic diversity, the BSP revealed a recent populational contraction (ca. 5 thousand years ago), after a long expansion period. These events could be related to anthropic actions, such as trade, tourism, and habitat degradation and loss.

Conservation implications. Solve taxonomic uncertainties, identify management units, and investigate the genetic diversity are crucial steps to management success. Here, we discuss some H. reidi conservation issues considering the presence of three management units, despite the presence of two lineages. This evidence justifies their separate management, setting different protection and sustainable actions. However, for sustainable use of the longsnout seahorse, we suggest that periods of non-harvesting be established during the months of October to February, which are the reproductive peaks of the species in Brazil.

Management Unit I consist of the NNE group, represented by North and Northeast populations, ranging from Pará to Alagoas states. This unit contains only the H. reidi Lineage A and is characterized by low genetic diversity (mtDNA), absence of genetic structure, shared haplotypes, and a low-moderate pairwise FST, suggesting a gene flow between the populations. This unit is considered the most vulnerable due to both low genetic diversity and possible contraction in population size by mismatch distribution analysis. Thus, we recommend that new conservation units should be created, or existing units should be amplified, to allow gene flow between the populations, avoiding the erosion of genetic variation. Furthermore, since H. reidi is traditionally used in the Brazilian Northeast as ornamental fishes, supervision efforts should be concentrated in trade and transport.

Management Unit II consists of the BA group, made up of specimens from the Bahia coast, representing a contact zone between Lineages A and B. This unit is characterized by the highest genetic diversity levels (mtDNA). In addition, the highest differentiation levels in all comparisons [except when compared to São Paulo (mtDNA and nuDNA) and Espírito Santo (nuDNA)] suggest a gene flow loss. The BA unit represents a priority area for conservation due to the high genetic levels and sympatry of the two lineages, representing a significant gene pool portion of H. reidi from the Brazilian coast. Thus, we suggest the integral protection of the longsnout seahorse along the Bahia coast, without allowing sustainable uses, and the supervision efforts should be concentrated in trade and transport.

Management Unit III consists of the SSE group, represented by South and Southeast, ranging from Espírito Santo to Santa Catarina states. This unit contains only the H. reidi Lineage B. Although the AMOVA hypothesis that considered the ES as a distinct group showed a high and significant differentiation between groups (FCT), the shared Lineage and exclusive genetic profiles with the others SES populations indicates a collaborative management.Contradicting the recent population contraction observed by the BSP analysis, the SES unit is characterized by high genetic diversity and the absence of an interpopulation structure, despite presenting some high pairwise FST values. Thus, this unit is in a reasonable conservation state, and the existing preservation efforts seem to be effective. Nevertheless, we reinforce the need to maintain the existing conservation units along the SES coast, especially due to the presence of another seahorse species, the H. patagonicus (unpublished data; Defavari, 2016).

ACKNOWLEDGEMENTS

This study was funded by the Fundação Grupo Boticário (Grant Number: 0964_20122) to ILR, and by the Fundação de Amparo à Ciência e Tecnologia do Estado de Pernambuco (Grant Number: 12/2010) to RAT. MCGQB and GRD are thankful to the CNPq and CAPES for the scholarships.

REFERENCES

• Allendorf FW, Berry O, Ryman N So long to genetic diversity, and thanks for all the fish. Mol Ecol. 2014; 23(1):23–25. https://doi.org/10.1111/mec.12574
  » https://doi.org/10.1111/mec.12574
• Andrade FRS, Afonso AS, Hazin FHV, Mendonça FF, Torres RA Population genetics reveals global and regional history of the apex predator Galeocerdo cuvier (Carcharhiniformes) with comments on mitigating shark attacks in north-eastern Brazil. Mar Ecol. 2021; 42(2):e12640. https://doi.org/10.1111/maec.12640
  » https://doi.org/10.1111/maec.12640
• Azpelicueta MLM, Delpiani SM, Cione AL, Oliveira C, Marceniuk AP, Díaz De Astarloa JM Morphology and molecular evidence support the validity of Pogonias courbina (Lacepède, 1803) (Teleostei: Sciaenidae), with a redescription and neotype designation. PLoS ONE. 2019; 14(6):e0216280. https://doi.org/10.1371/journal.pone.0216280
  » https://doi.org/10.1371/journal.pone.0216280
• Baum JK, Vincent ACJ Magnitude and inferred impacts of the seahorse trade in Latin America. Environ Conserv. 2005; 32(4):305–19. https://doi.org/10.1017/S0376892905002481
  » https://doi.org/10.1017/S0376892905002481
• Bouckaert R, Heled J, Kühnert D, Vaughan T, Wu CH, Xie D, Suchard MA, Rambaut A, Drummond AJ BEAST 2: a software platform for Bayesian evolutionary analysis. PLOS Comput Biol. 2014; 10(4):e1003537. https://doi.org/10.1371/journal.pcbi.1003537
  » https://doi.org/10.1371/journal.pcbi.1003537
• Bowen BW, Muss A, Rocha LA, Grant WS Shallow mtDNA coalescence in Atlantic pygmy angelfishes (genus Centropyge) indicates a recent invasion from the Indian Ocean. J Hered. 2006; 97(1):1–12. https://doi.org/10.1093/jhered/esj006
  » https://doi.org/10.1093/jhered/esj006
• Cadrin SX, Karr LA, Mariani S Stock identification methods: an overview. In: Cadrin SX, Karr LA, Mariani S, editors. Stock identification methods. Academic Press; 2014. p.1–05.
• Cadrin SX Defining spatial structure for fishery stock assessment. Fish Res. 2020; 221:105397. https://doi.org/10.1016/j.fishres.2019.105397
  » https://doi.org/10.1016/j.fishres.2019.105397
• Calcagnotto D Taxas de evolução e relógios moleculares. In: Matioli SR, Fernandes FMC, editors. Biologia molecular e evolução. Holos Editora, Ribeirão Preto, São Paulo; 2012. p.61–73.
• Carmo TF, Santos LN, Bertoncini AA, Freret-Meurer NV Population structure of the seahorse Hippocampus reidi in two Brazilian estuaries. Ocean Coast Res. 2022; 70:e22009. http://doi.org/10.1590/2675-2824070.21016tfdc
  » http://doi.org/10.1590/2675-2824070.21016tfdc
• Chen W, Li C, Chen F, Li Y, Yang J, Li J, Li X Phylogeographic analyses of a migratory freshwater fish (Megalobrama terminalis) reveal a shallow genetic structure and pronounced effects of sea-level changes. Gene. 2020; 737:144478. https://doi.org/10.1016/j.gene.2020.144478
  » https://doi.org/10.1016/j.gene.2020.144478
• Chow S, Hazama K Universal primers for S7 ribosomal protein gene introns in fish. Mol Ecol. 1998; 7(9):1255–56.
• Clement MJ, Snell Q, Walker P, Posada D, Crandall KA TCS: estimating gene genealogies. IPDPS. 2002; 311(4):1110–16. Available from: http://www.hicomb.org/HiCOMB2003/papers/HICOMB2002-03.pdf
  » http://www.hicomb.org/HiCOMB2003/papers/HICOMB2002-03.pdf
• Corander J, Marttinen P Bayesian identification of admixture events using multilocus molecular markers. Mol Ecol. 2006; 15(10):2833–43. https://doi.org/10.1111/j.1365-294X.2006.02994.x
  » https://doi.org/10.1111/j.1365-294X.2006.02994.x
• Corander J, Marttinen P, Sirén J, Tang J Enhanced Bayesian modelling in BAPS software for learning genetic structures of populations. BMC Bioinformatics. 2008; 9:1–14. https://doi.org/10.1186/1471-2105-9-539
  » https://doi.org/10.1186/1471-2105-9-539
• Cowen RK, Lwiza KM, Sponaugle S, Paris CB, Olson DB Connectivity of marine populations: open or closed? Science. 2000; 287(5454):857–59. https://doi.org/10.1126/science.287.5454.857
  » https://doi.org/10.1126/science.287.5454.857
• Cowen RK, Paris CB, Srinivasan A Scaling of connectivity in marine populations. Science. 2006; 311(5760):522–27. https://doi.org/10.1126/science.112203
  » https://doi.org/10.1126/science.112203
• Cunha IMC, Souza AS, Dias Jr. EA, Amorim KDJ, Soares RX, Costa GWWF, García-Machado E, Galetti Jr. PM, Molina WF Genetic multipartitions based on d-loop sequences and chromosomal patterns in brown Chromis, Chromis multilineata (Pomacentridae), in the western Atlantic. BioMed Res Int. 2014; 2014. https://doi.org/10.1155/2014/254698
  » https://doi.org/10.1155/2014/254698
• Darriba D, Taboada GL, Doallo R, Posada D jModelTest 2: more models, new heuristics and parallel computing. Nat Methods. 2012; 9:772.
• Day J, Clark JA, Williamson JE, Brown C, Gillings M Population genetic analyses reveal female reproductive philopatry in the oviparous Port Jackson shark. Mar Freshwater Res. 2019; 70(7):986–94. https://doi.org/10.1071/MF18255
  » https://doi.org/10.1071/MF18255
• Dolby GA, Hechinger R, Ellingson RA, Findley LT, Lorda J, Jacobs DK Sea-level driven glacial-age refugia and post-glacial mixing on subtropical coasts, a palaeohabitat and genetic study. Proc Royal So B. 2016; 283(1843):20161571. https://doi.org/10.1098/rspb.2016.1571
  » https://doi.org/10.1098/rspb.2016.1571
• Dolby GA, Ellingson RA, Findley LT, Jacobs DK How sea level change mediates genetic divergence in coastal species across regions with varying tectonic and sediment processes. Mol Ecol. 2018; 27(4):994–1011. https://doi.org/10.1111/mec.14487
  » https://doi.org/10.1111/mec.14487
• Dominguez JML, Ramos JMF, Rebouças RC, Nunes AS, Melo LD A plataforma continental do município de Salvador: geologia, usos múltiplos e recursos minerais. In: Barbosa JSF, Mascarenhas JF, Corrêa Gomes LC, editors. Geologia da Bahia vol. 2. Companhia Baiana de Pesquisa Mineral, Salvador, Bahia; 2012. p.427–96.
• Drummond AJ, Rambaut A, Shapiro BE, Pybus OG Bayesian coalescent inference of past population dynamics from molecular sequences. Mol Biol Evol. 2005; 22(5):1185–92. https://doi.org/10.1093/molbev/msi103
  » https://doi.org/10.1093/molbev/msi103
• Excoffier L, Lischer HEL Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010; 10(3):564–67. https://doi.org/10.1111/j.1755-0998.2010.02847.x
  » https://doi.org/10.1111/j.1755-0998.2010.02847.x
• Faria DM, Silva Jr. JM, Pires Costa L, Rezende Paiva S, Marino CL, Rollo Jr. MM, Cazerta Farro AP Low mtDNA diversity in a highly differentiated population of spinner dolphins (Stenella longirostris) from the Fernando de Noronha Archipelago, Brazil. PLoS ONE. 2020; 15(4):e0230660. https://doi.org/10.1371/journal.pone.0230660
  » https://doi.org/10.1371/journal.pone.0230660
• Foster SJ, Vincent ACJ Life history and ecology of seahorses: implications for conservation and management. J Fish Biol. 2004; 65(1):1–61. https://doi.org/10.1111/j.0022-1112.2004.00429.x
  » https://doi.org/10.1111/j.0022-1112.2004.00429.x
• Fu Yun-Xin Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics. 1997; 147(2):915–25. https://doi.org/10.1093/genetics/147.2.915
  » https://doi.org/10.1093/genetics/147.2.915
• Goswami M, Thangaraj K, Chaudhary BK, Bhaskar LVSK, Gopalakrishnan A, Joshi MB, Singh L, Lakra WS Genetic heterogeneity in the Indian stocks of seahorse (italic)Hippocampus kuda(/italic) and (italic)Hippocampus trimaculatus(/italic) inferred from mtDNA cytochrome b gene. Hydrobiologia. 2009; 621:213–21. https://doi.org/10.1007/s10750-008-9642-3
  » https://doi.org/10.1007/s10750-008-9642-3
• Grant WAS, Bowen BW Shallow population histories in deep evolutionary lineages of marine fishes: insights from sardines and anchovies and lessons for conservation. J Hered. 1998; 89(5):415–26. https://doi.org/10.1093/jhered/89.5.415
  » https://doi.org/10.1093/jhered/89.5.415
• Green ME, Anastasi BR, Hobbs JPA, Feldheim K, McAuley R, Peverell S, Stapley J, Johnson G, Appleyard SA, White WT, Simpfendorfer CA, van Herwerden L Mixed-marker approach suggests maternal philopatry and sex-biased behaviours of narrow sawfish Anoxypristis cuspidata Endanger Species Res. 2018; 37:45–54. https://doi.org/10.3354/esr00912
  » https://doi.org/10.3354/esr00912
• Guillot G, Mortier F, Estoup A GENELAND: a computer package for landscape genetics. Mol Ecol Notes. 2005; 5(3):712–15. https://doi.org/10.1111/j.1471-8286.2005.01031.x
  » https://doi.org/10.1111/j.1471-8286.2005.01031.x
• Hall TA BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acid S. 1999; 41(41):95–98.
• Haney RA, Silliman BR, Rand DM Effects of selection and mutation on mitochondrial variation and inferences of historical population expansion in a Caribbean reef fish. Mol Phylogenet Evol. 2010; 57(2):821–28. https://doi.org/10.1016/j.ympev.2010.07.014
  » https://doi.org/10.1016/j.ympev.2010.07.014
• Healey AJ, Mckeown NJ, Taylor AL, Provan J, Sauer W, Gouws G, Shaw PW Cryptic species and parallel genetic structuring in Lethrinid fish: Implications for conservation and management in the southwest Indian Ocean. Ecol Evol. 2018; 8(4):2182–95. https://doi.org/10.1002/ece3.3775
  » https://doi.org/10.1002/ece3.3775
• Hellberg ME, Burton RS, Neigel JE, Palumbi SR Genetic assessment of connectivity among marine populations. B Mar Sci. 2002; 70(1):273–90.
• Hercos AP, Giarrizzo T Pisces, Syngnathidae, Hippocampus reidi: filling distribution gaps. Check List. 2007; 3(4):287–90. https://doi.org/10.15560/3.4.287
  » https://doi.org/10.15560/3.4.287
• da Hora MSC, Joyeux J-C, Rodrigues RV, Sousa-Santos LP, Gomes LC, Tsuzuki MY Tolerance and growth of the longsnout seahorse Hippocampus reidi at different salinities. Aquaculture. 2016; 463:1–06. https://doi.org/10.1016/j.aquaculture.2016.05.003
  » https://doi.org/10.1016/j.aquaculture.2016.05.003
• Instituto Chico Mendes de Conservação da Biodiversidade (ICMBio). Livro Vermelho da Fauna Brasileira Ameaçada de Extinção: Volume VI - Peixes. In: Instituto Chico Mendes de Conservação da Biodiversidade (Org.). Livro Vermelho da Fauna Brasileira Ameaçada de Extinção. Brasília: ICMBio; 2018.
• Jacobina UP, Torres RA, Affonso PRAM, Santos EV, Calado LL, Bitencourt JA DNA barcoding reveals cryptic diversity and peculiar phylogeographic patterns in mojarras (Perciformes: Gerreidae) from the Caribbean and South-western Atlantic. J Mar Biol Assoc UK. 2020; 100(2):277–83. https://doi.org/10.1017/S0025315419001206
  » https://doi.org/10.1017/S0025315419001206
• Klanten OS, Gaither MR, Greaves S, Mills K, O’keeffe K, Turnbull J, Booth DJ Genomic and morphological evidence of distinct populations in the endemic common (weedy) seadragon Phyllopteryx taeniolatus (Syngnathidae) along the east coast of Australia. PLoS ONE. 2020; 15(12):e0243446. https://doi.org/10.1371/journal.pone.0243446
  » https://doi.org/10.1371/journal.pone.0243446
• Labrador K, Agmata A, Palermo JD, Ravago-Gotanco R, Pante MJ Mitochondrial DNA reveals genetically structured haplogroups of Bali sardinella Sardinella lemuru in Philippine waters. Reg Stud Mar Sci. 2021; 41:101588. https://doi.org/10.1016/j.rsma.2020.101588
  » https://doi.org/10.1016/j.rsma.2020.101588
• Laurrabaquio-A NS, Islas-Villanueva V, Adams DH, Uribe-Alcocer M, Alvarado-Bremer J.R, Díaz-Jaimes P Genetic evidence for regional philopatry of the bull shark Carcharhinus leucas, to nursery areas in estuaries of the Gulf of Mexico and western North Atlantic Ocean. Fish Res. 2019; 209:67–74. https://doi.org/10.1016/j.fishres.2018.09.013
  » https://doi.org/10.1016/j.fishres.2018.09.013
• Lehnert SJ, Dibacco C, Van Wyngaarden M, Jeffery NW, Lowen JB, Sylvester EV, Bradbury IR Fine-scale temperature-associated genetic structure between inshore and offshore populations of sea scallop Placopecten magellanicus Heredity. 2019; 122(1):69–80. https://doi.org/10.1038/s41437-018-0087-9
  » https://doi.org/10.1038/s41437-018-0087-9
• Leigh JW, Bryant D POPART: full-feature software for haplotype network construction. Methods Ecol Evol. 2015; 6(9):1110–16. https://doi.org/10.1111/2041-210X.12410
  » https://doi.org/10.1111/2041-210X.12410
• Librado P, Rozas J DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009; 25(11):1451–52. https://doi.org/10.1093/bioinformatics/btp187
  » https://doi.org/10.1093/bioinformatics/btp187
• Lourie SA, Vincent AC, Hall HJ Seahorses: an identification guide to the world’s species and their conservation. Project Seahorse: London; 1999.
• Lourie S Fin-clipping procedure for seahorses. Project Seahorse Technical Bulletin No. 3, Version1.1. Project Seahorse, Fisheries Centre, University of British Columbia; 2003.
• Lourie SA, Foster SJ, Cooper EW, Vincent AC A guide to the identification of seahorses. Project Seahorse and TRAFFIC North America: Washington D.C.: University of British Columbia and World Wildlife Fund; 2004.
• Lourie SA, Green DM, Vincent ACJ Dispersal, habitat differences, and comparative phylogeography of Southeast Asian seahorses (Syngnathidae: Hippocampus). Mol Ecol. 2005; 14(4):1073–94. https://doi.org/10.1111/j.1365-294X.2005.02464.x
  » https://doi.org/10.1111/j.1365-294X.2005.02464.x
• Machado-Schiaffino G, Juanes F, Garcia-Vazquez E Introgressive hybridization in North American hakes after secondary contact. Mol Phylogenet Evol. 2010; 55(2):552–58. https://doi.org/10.1016/j.ympev.2010.01.034
  » https://doi.org/10.1016/j.ympev.2010.01.034
• Martínez-Candelas IA, Pérez-Jiménez JC, Espinoza-Tenorio A, McClenachan L, Méndez-Loeza I Use of historical data to assess changes in the vulnerability of sharks. Fish Res. 2020; 226:105526. https://doi.org/10.1016/j.fishres.2020.105526
  » https://doi.org/10.1016/j.fishres.2020.105526
• Mccauley DJ, Pinsky ML, Palumbi SR, Estes JA, Joyce FH, Warner RR Marine defaunation: animal loss in the global ocean. Science. 2015; 347(6219):1255641. https://doi.org/10.1126/science.1255641
  » https://doi.org/10.1126/science.1255641
• Mckeown NJ, Gwilliam MP, Healey AJ, Skujina I, Potts WM, Sauer WH, Shaw PW Deep phylogeographic structure may indicate cryptic species within the Sparid genus Spondyliosoma J Fish Biol. 2020; 96(6):1434–43. https://doi.org/10.1111/jfb.14316
  » https://doi.org/10.1111/jfb.14316
• Mobley KB, Small CM, Jue NK, Jones AG Population structure of the dusky pipefish (Syngnathus floridae) from the Atlantic and Gulf of Mexico, as revealed by mitochondrial DNA and microsatellite analyses. J Biogeogr. 2010; 37(7):1363–77. https://doi.org/10.1111/j.1365-2699.2010.02288.x
  » https://doi.org/10.1111/j.1365-2699.2010.02288.x
• Montes MA, Cardoso MLV, Neves CHCB, Garcia ACL, Silva JC, Silveira RB Genetic diversity and populational structure of the seahorse Hippocampus reidi (Syngnathidae) in north-eastern Brazil: a conservationist approach. Aquat Conserv Mar Freshw Ecosyst. 2018; 28(5):1114–22. https://doi.org/10.1002/aqc.2919
  » https://doi.org/10.1002/aqc.2919
• Murray T, Cowley P, Childs AR, Bennett R Philopatry and dispersal of juvenile leervis Lichia amia (Teleostei: Carangidae) tagged in a warm-temperate South African estuary. Afr J Mar Sci. 2017; 39(1):59–68. https://doi.org/10.2989/1814232X.2017.1303401
  » https://doi.org/10.2989/1814232X.2017.1303401
• Musick JA, Harbin MM, Berkeley SA, Burgess GH, Eklund AM, Findley L, Gilmore RG, Golden JT, Ha DS, Huntsman GR, McGovern JC, Sedberry GR, Parker SJ, Poss SG, Sala E, Schmidt TW, Weeks H, Wright SG Marine, estuarine, and diadromous fish stocks at risk of extinction in North America (Exclusive of Pacific Salmonids). Fisheries. 2000; 25(11):6–30. https://doi.org/10.1577/1548-8446(2000)025<0006:MEADFS>2.0.CO;2
  » https://doi.org/10.1577/1548-8446(2000)025<0006:MEADFS>2.0.CO;2
• Nauer F, Deluqui Gurgel CF, Ayres-Ostrock LM, Plastino EM, Oliveira MC Phylogeography of the Hypnea musciformis species complex (Gigartinales, Rhodophyta) with the recognition of cryptic species in the western Atlantic Ocean. J Phycol. 2019; 55(3):676–87. https://doi.org/10.1111/jpy.12848
  » https://doi.org/10.1111/jpy.12848
• Neves A, Vieira AR, Sequeira V, Paiva RB, Gordo LS, Paulo OS Highly regional population structure of Spondyliosoma cantharus depicted by nuclear and mitochondrial DNA data. Sci Rep. 2020; 10:4063. https://doi.org/10.1038/s41598-020-61050-x
  » https://doi.org/10.1038/s41598-020-61050-x
• Oliveira T, Pollom R Hippocampus reidi. The IUCN Red List of Threatened Species. 2017. Available from: https://dx.doi.org/10.2305/IUCN.UK.20173.RLTS.T10082A17025021.en
  » https://dx.doi.org/10.2305/IUCN.UK.20173.RLTS.T10082A17025021.en
• Palumbi SR Marine speciation on a small planet. Trends Ecol Evol. 1992; 7(4):114–18. https://doi.org/10.1016/0169-5347(92)90144-Z
  » https://doi.org/10.1016/0169-5347(92)90144-Z
• Pan J, Marcova MA, Bazzini SM, Vallina MV, Marco SG Coastal marine biodiversity challenges and threats. In: Arias AH, Menendez MC, editors. Marine ecology in a changing world. Boca Raton, FL: CRC Press; 2013. p.43–67.
• Panithanarak T, Karuwancharoen R, Na-Nakorn U, Nguyen TTT Population genetics of the spotted seahorse Hippocampus kuda in Thai Waters: implications for conservation. Zool Stud. 2010; 49(4):564–76. https://zoolstud.sinica.edu.tw/Journals/49.4/564.pdf
  » https://zoolstud.sinica.edu.tw/Journals/49.4/564.pdf
• Planas M, Vilar A, López A, Bouza C Técnicas de identificação individual de caballitos de mar (Hippocampus guttulatus) para uma correta gestão de reproductores em cautividad: uso de collares e amostras não invasivas para análise genética. Anais do XI Congreso Nacional de Acuicultura. 2007; p.203–06.
• Pinheiro HT, Mazzei E, Moura RL, Amado-Filho GM, Carvalho-Filho A, Braga AC, Costa PAS, Ferreira BP, Ferreira CEL, Floeter SR, Francini-Filho RB, Gasparini JL, Macieira RM, Martins AS, Olavo G, Pimentel CR, Rocha LA, Sazima I, Simon T, Teixeira JB, Xavier LB, Joyeux JC Fish biodiversity of the Vitória-Trindade Seamount Chain, southwestern Atlantic: an updated database. PLoS ONE. 2015; 10(3):e0118180. https://doi.org/10.1371/journal.pone.0118180
  » https://doi.org/10.1371/journal.pone.0118180
• Pinheiro HT, Rocha LA, Macieira RM, Carvalho-Filho A, Anderson AB, Bender MG, DiDario F, Ferreira CEL, Figueiredo-Filho J, Francini-Filho R, Gasparini JL, Joyeux JC, Luiz OJ, Mincarone MM, Moura RL, Nunes JACC, Quimbayo JP, Rosa RS, Sampaio CLS, Sazima I, Simon T, Villa-Nova DA, Floeter SR South-western Atlantic reef fishes: zoogeographical patterns and ecological drivers reveal a secondary biodiversity centre in the Atlantic Ocean. Divers Distrib. 2018; 24(7):951–65. https://doi.org/10.1111/ddi.12729
  » https://doi.org/10.1111/ddi.12729
• Pinsky ML, Palumbi SR Meta-analysis reveals lower genetic diversity in overfished populations. Mol Ecol. 2014; 23(1):29–39. https://doi.org/10.1111/mec.12509
  » https://doi.org/10.1111/mec.12509
• Rogers AR, Harpending H Population growth makes waves in the distribution of pairwise genetic differences. Mol Biol Evol. 1992; 9(3):552–69. https://doi.org/10.1093/oxfordjournals.molbev.a040727
  » https://doi.org/10.1093/oxfordjournals.molbev.a040727
• Ronquist F, Huelsenbeck JP MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003; 19(12):1572–74. https://doi.org/10.1093/bioinformatics/btg180
  » https://doi.org/10.1093/bioinformatics/btg180
• Rosa IL National Report Brazil. In: Bruckner AW, Fields JD, Daves N, editors. Proceedings of the international workshop on CITES implementation for seahorse conservation and trade, Mazatlán, Sinaloa: Mexico; 2005. Available from: https://repository.library.noaa.gov/view/noaa/456
  » https://repository.library.noaa.gov/view/noaa/456
• Roycroft EJ, Le Port A, Lavery SD Population structure and male-biased dispersal in the short-tail stingray Bathytoshia brevicaudata (Myliobatoidei: Dasyatidae). Conserv Genet. 2019; 20(4):717–28. https://doi.org/10.1007/s10592-019-01167-3
  » https://doi.org/10.1007/s10592-019-01167-3
• Sadeghi R, Esmaeili HR, Zarei F, Reichenbacher B Matrilineal evidence for genetic structure and Late Pleistocene demographic expansion of the Ornate goby Istigobius ornatus (Teleostei: Gobiidae) in the Persian Gulf and Oman Sea. Mar Ecol. 2021; 42(1):e12629. https://doi.org/10.1111/maec.12629
  » https://doi.org/10.1111/maec.12629
• Saarman NP, Louie KD, Hamilton H Genetic differentiation across eastern Pacific oceanographic barriers in the threatened seahorse Hippocampus ingens Conserv Genet. 2010; 11(5):1989–2000. https://doi.org/10.1007/s10592-010-0092-x
  » https://doi.org/10.1007/s10592-010-0092-x
• Santos S, Hrbek T, Farias IP, Schneider H, Sampaio I Population genetic structuring of the king weakfish, Macrodon ancylodon (Sciaenidae), in Atlantic coastal waters of South America: deep genetic divergence without morphological change. Mol Ecol. 2006; 15(14):4361–73. https://doi.org/10.1111/j.1365-294X.2006.03108.x
  » https://doi.org/10.1111/j.1365-294X.2006.03108.x
• Santos S, Schneider H, Sampaio I Genetic differentiation of Macrodon ancylodon (Sciaenidae, Perciformes) populations in Atlantic coastal waters of South America as revealed by mtDNA analysis. Genet Mol Biol. 2003; 26(2):151–61. https://doi.org/10.1590/S1415-47572003000200008
  » https://doi.org/10.1590/S1415-47572003000200008
• dos Santos Freitas A, Silva R, Sampaio I, Schneider H The mitochondrial control region reveals genetic structure in southern kingcroaker populations on the coast of the Southwestern Atlantic. Fish Res. 2017; 191:87–94. https://doi.org/10.1016/j.fishres.2017.03.008
  » https://doi.org/10.1016/j.fishres.2017.03.008
• Selkoe KA, Gaggiotti OE, Bowen BW, Toonen RJ Emergent patterns of population genetic structure for a coral reef community. Mol Ecol. 2014; 23(12):3064–79. https://doi.org/10.1111/mec.12804
  » https://doi.org/10.1111/mec.12804
• Selkoe KA, Aloia CC, Crandall ED, Iacchei M, Liggins L, Puritz JB, Toonen RJ A decade of seascape genetics: contributions to basic and applied marine connectivity. Mar Ecol Prog Ser. 2016; 554:1–19. https://doi.org/10.3354/meps11792
  » https://doi.org/10.3354/meps11792
• da Silva R, Sampaio I, Schneider H, Gomes G Lack of spatial subdivision for the snapper Lutjanus purpureus (Lutjanidae-Perciformes) from Southwest Atlantic based on multi-locus analyses. PLoS ONE. 2016; 11(8):e0161617. https://doi.org/10.1371/journal.pone.0161617
  » https://doi.org/10.1371/journal.pone.0161617
• Silveira RB Registros de cavalos-marinhos (Syngnathidae: Hippocampus) ao longo da costa Brasileira. Oecol Aust. 2011; 15(2):316–325.
• Silveira RB, Siccha-Ramirez R, Silva JRS, Oliveira C Morphological and molecular evidence for the occurrence of three Hippocampus species (Teleostei: Syngnathidae) in Brazil. Zootaxa. 2014; 3861(4):317–32. http://dx.doi.org/10.11646/zootaxa.3861.4.2
  » http://dx.doi.org/10.11646/zootaxa.3861.4.2
• Souza ASD, Dias Júnior EA, Perez MF, Cioffi MB, Bertollo LAC, Garcia-Machado E, Vallinoto MNS, Galetti Jr. PM, Molina WF Phylogeography and historical demography of two sympatric Atlantic snappers: Lutjanus analis and L. jocu. Front Mar Sci. 2019; 6:545. https://doi.org/10.3389/fmars.2019.00545
  » https://doi.org/10.3389/fmars.2019.00545
• Spalding MD, Fox HE, Allen GR, Davidson N, Ferdaña ZA, Finlayson M, Halpern BS, Jorge MA, Lombana A, Lourie SA, Martin KD, McManus E, Molnar J, Recchia CA, Robertson J Marine ecoregions of the world: a bioregionalization of coastal and shelf areas. BioScience. 2007; 57(7):573–83. https://doi.org/10.1641/B570707
  » https://doi.org/10.1641/B570707
• Stephens AM, Smith NJ, Donnelly P, Stephens CM, Li CFN Documentation for PHASE, version 2.1. 2004. Available from: https://www.animalgenome.org/bioinfo/resources/manuals/PHASE.pdf
  » https://www.animalgenome.org/bioinfo/resources/manuals/PHASE.pdf
• Spalding MD, Fox HE, Allen GR, Davidson N, Ferdaña ZA, Finlayson M, Halpern BS, Jorge MA, Lombana A, Lourie SA, Martin KD, McManus E, Molnar J, Recchia CA, Robertson J Marine ecoregions of the world: a bioregionalization of coastal and shelf areas. BioScience. 2007; 57(7):573–83. https://doi.org/10.1641/B570707
  » https://doi.org/10.1641/B570707
• Ternes MLF, Gerhardinger LC, Schiavetti A Seahorses in focus: local ecological knowledge of seahorse-watching operators in a tropical estuary. J Ethnobiol Ethnomed. 2016; 12(52):1–12. https://doi.org/10.1186/s13002-016-0125-8
  » https://doi.org/10.1186/s13002-016-0125-8
• Teske PR, Cherry MI, Matthee CA Population genetics of the endangered Knysna seahorse, Hippocampus capensis. Mol Ecol. 2003; 12(7):1703–15. https://doi.org/10.1046/j.1365-294X.2003.01852.x
  » https://doi.org/10.1046/j.1365-294X.2003.01852.x
• Teske PR, Cherry MI, Matthee CA The evolutionary history of seahorses (Syngnathidae: Hippocampus): molecular data suggest a West Pacific origin and two invasions of the Atlantic Ocean. Mol Phylogenet Evol. 2004; 30(2):273–86. https://doi.org/10.1016/S1055-7903(03)00214-8
  » https://doi.org/10.1016/S1055-7903(03)00214-8
• Thompson JD, Higgins DG, Gibson TJ CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994; 22(22):4673–80. https://doi.org/10.1093/nar/22.22.4673
  » https://doi.org/10.1093/nar/22.22.4673
• Torrado H, Carreras C, Raventos N, Macpherson E, Pascual M Individual-based population genomics reveal different drivers of adaptation in sympatric fish. Sci Rep-UK. 2020; 10:12683. https://doi.org/10.1038/s41598-020-69160-2
  » https://doi.org/10.1038/s41598-020-69160-2
• Toews DPL, Brelsford A The biogeography of mitochondrial and nuclear discordance in animals. Mol Ecol. 2012; 21(16):3907–30. https://doi.org/10.1111/j.1365-294X.2012.05664.x
  » https://doi.org/10.1111/j.1365-294X.2012.05664.x
• Tseng CC, Chien JH, Chu TW, Cheng AC, Shiu YL, Han TW, Liu CH Effects of food type, temperature and salinity on the growth performance and antioxidant status of the longsnout seahorse, Hippocampus reidi. Aquac. Res. 2020; 51(11):4793–804. https://doi.org/10.1111/are.14826
  » https://doi.org/10.1111/are.14826
• Zarei F, Esmaeili HR, Abbasi K, Sayyadzadeh G, Eagderi S, Coad BW Genealogical concordance, comparative species delimitation, and the specific status of the Caspian pipefish Syngnathus caspius (Teleostei: Syngnathidae). Mar Ecol. 2021; 42(1):e12624. https://doi.org/10.1111/maec.12624
  » https://doi.org/10.1111/maec.12624
• Zhao L, Shan B, Song N, Gao T Genetic diversity and population structure of Acanthopagrus schlegelii inferred from mtDNA sequences. Reg Stud Mar Sci. 2020; 41:101532. https://doi.org/10.1016/j.rsma.2020.101532
  » https://doi.org/10.1016/j.rsma.2020.101532
• Zink RM, Barrowclough GF Mitochondrial DNA under siege in avian phylogeography. Mol Ecol. 2008; 17(9):2107–21. https://doi.org/10.1111/j.1365-294X.2008.03737.x
  » https://doi.org/10.1111/j.1365-294X.2008.03737.x

ADDITIONAL NOTES

• HOW TO CITE THIS ARTICLE
  Queiroz-Brito MCG, Defavari GR, Rosa IL, Torres RA. Population structure of long-snout seahorse Hippocampus reidi in Southwestern Atlantic and implications for management. Neotrop Ichthyol. 2024; 22(3):e240027. https://doi.org/10.1590/1982-0224-2024-0027

• Edited-by
  Osmar Luiz

Publication Dates

• Publication in this collection
  23 Sept 2024
• Date of issue
  2024

History

• Received
  21 Mar 2024
• Accepted
  18 July 2024