SciELO - Scientific Electronic Library Online

 
vol.43 número1Contribución del hidroperiodo y la microtopografía a la regeneración natural de Avicennia germinans en un sitio de restauración ecológica índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay artículos similaresSimilares en SciELO

Compartir


Ciencias marinas

versión impresa ISSN 0185-3880

Cienc. mar vol.43 no.1 Ensenada mar. 2017  Epub 27-Mayo-2021

https://doi.org/10.7773/cm.v43i1.2692 

Articles

Population structure of the Pacific angel shark (Squatina californica) along the northwestern coast of Mexico based on the mitochondrial DNA control region

Estructura poblacional del tiburón ángel del Pacífico (Squatina californica) a lo largo de la costa noroccidental de México con base en la regíon control del ADN mitocondrial

Sergio Ramírez-Amaro1 

Dení Ramírez-Macías2 

Ricardo Vázquez-Juárez3  * 

Sergio Flores-Ramírez4 

Felipe Galván-Magaña5 

Jesús N Gutiérrez-Rivera3 

1 Universitat de les Illes Balears, Laboratori de Genètica, Carretera de Valldemossa km 7.5, 07122, Palma de Mallorca, Spain.

2Tiburón Ballena México de ConCIENCIA México, Manatí 4, no. 802, colonia Esperanza III, CP 23090, La Paz, Baja California Sur, México.

3 Centro de Investigaciones Biológicas del Noroeste, Mar Bermejo 195, colonia Playa Palo de Santa Rita Sur, CP 23096, La Paz, Baja California Sur, México.

4 Departamento de Biología Marina, Universidad Autónoma de Baja California Sur, Apartado Postal 12-B, CP 23081, La Paz, Baja California Sur, México.

5 Instituto Politécnico Nacional (IPN), Centro Interdisciplinario de Ciencias Marinas, Av. IPN s/n, colonia Playa Palo de Santa Rita, CP 23096, La Paz, Baja California Sur, México.


Abstract

Population structure has direct consequences on species evolution, local adaptation, and the capacity of a metapopulation to adapt to climate change. In this study we assess the genetic structure and demographic history of the Pacific angel shark, Squatina californica, along the northwestern coast of Mexico using the mitochondrial control region. Results showed high levels of genetic diversity and reveal genetic differentiation between the samples from the Pacific coast of the Baja California Peninsula and those from the Gulf of California. These genetic differences are consistent with ontogenetic differences previously reported for S. californica off the northwestern coast of Mexico. Our findings can be explained by the low dispersal ability and habitat preference of the angel shark and by the complex geomorphology off the northwestern coast of Mexico. Contrasting historical demographic patterns were observed: population expansion in the gulf and population stability, with a declining trend, in the Pacific; however, more studies are required to corroborate demographic patterns in both populations. Our results point to 2 distinct management units of the Pacific angel shark in the study area, and this may have direct implications for the management and conservation of this species in Mexico.

Key words: demersal shark; ecological clusters; genetic differentiation; Gulf of California; Pacific coast

Resumen

La estructura poblacional tiene consecuencias directas sobre la evolución de las especies, la adaptación local y la capacidad de una metapoblación para adaptarse al cambio climático. En el presente estudio, se analizó la estructura genética y la historia demográfica del tiburón ángel del Pacífico, Squatina californica, a lo largo de la costa noroccidental de México con base en la región control del ADN mitocondrial. Los resultados mostraron elevados niveles de diversidad genética y diferenciación genética entre las muestras del Pacífico frente a la península de Baja California y aquellas del golfo de California. Estas diferencias genéticas son consistentes con las diferencias ontogenéticas reportadas previamente para S. californica en la costa noroccidental de México. Los resultados pueden ser explicados por la baja capacidad de dispersión y la preferencia de hábitat de S. californica y por la compleja geomorfología frente a la costa noroccidental de México. Además, se observaron diferentes patrones demográficos históricos: la población del golfo exhibió una tendencia a la expansión y la población del Pacífico exhibió estabilidad con tendencia a disminuir; sin embargo, es necesario realizar más estudios para corroborar los patrones demográficos detectados en ambas poblaciones. Los resultados sugieren la presencia de 2 diferentes unidades para la gestión pesquera del tiburón ángel del Pacífico en el área de estudio, lo cual tiene implicaciones directas para la conservación y gestión de esta especie en México.

Palabras clave: tiburón demersal; grupos ecológicos; diferenciación genética; golfo de California; costa del Pacífico

Introduction

Understanding the spatiotemporal patterns of gene flow among isolated populations has long been a major focus in ecology (Spaet et al. 2015). Sharks are likely to display different patterns of population structure. Low levels of intraspecific genetic differentiation have been detected mainly in highly migrant pelagic species (e.g., Rhincodon typus) across large ocean basins (Vignaud et al. 2014), while high levels of population structure have been detected mainly in demersal species (e.g., Stegostoma fasciatum) on relatively small geographic scales (Dudgeon et al. 2009). Absence of population structure was detected for the deep-water Centroscymnus coelolepis (Veríssimo et al. 2011), whereas almost global panmixia was found for Cetorhinus maximus (Hoelzel et al. 2006). These population patterns are commonly promoted by geographic barriers or dispersal, or can be influenced by environmental gradients throughout continuous areas (Spaet et al. 2015).

Knowing the processes that affect genetic diversity and population structure is crucial for the conservation of exploited species (Arenas et al. 2012). Overfishing can decrease genetic variation through bottlenecks, cause loss of genetically unique stocks, select for early sexual maturation, and ultimately reduce the adaptability of wild populations (Hutchings and Reynolds 2004). Understanding population structure and connectivity of marine resources is an essential prerequisite for devising effective fishery management strategies in at least 2 critical ways: (i) the delineation of an appropriate spatial scale for management and (ii) the specification of subareas, nested within broader management units, that must be protected from exploitation (Fogarty and Botsford 2007, Hilário et al. 2015).

The Pacific angel shark, Squatina californica Ayres 1859, is a demersal shark found off the coast of North America, from Alaska to the Gulf of California, where it is associated with soft-bottom habitats (Compagno et al. 2005). It is a lecithotrophic viviparous shark that has low fecundity (10 embryos per year), with gestation lasting 10 months (Compagno et al. 2005, Romero-Caicedo et al. 2016). This species is frequently targeted by the artisanal elasmobranch fishery along the norwestern coast of Mexico (Villavicencio-Garayzar 1996, Ramírez-Amaro et al. 2013) . The population trend remains unknown, but reported landing data (2003-2015) from Baja California (Mexico) suggest population decline (Ramírez-Amaro et al. 2013). The Pacific angel shark is categorized as Near Threatened in the International Union for Conservation of Nature Red List of Endangered Species, with a suspected population decline approaching 30% over 3 generations (Cailliet et al. 2016). Currently, no information on stock or structure has been documented for this species.

Assessing the differences in the population structure of a fishery resource is necessary to understand the extent to which populations are genetically isolated from each other and need different management plans to protect the species (Nielsen and Powers 1995). It is important to select a locus or loci with a relative high mutation rate in order to detect sufficient polymorphism in population-level studies (Heist et al. 1996). The mitochondrial DNA genome includes a small non-coding region known as the mitochondrial control region (mtCR), which serves as the origin of replication for the mitochondrial genome and is usually more variable than other coding genes (Avise 1994). Several shark population studies have used the mtCR to infer processes responsible for observed patterns of spatial genetic differentiation and to identify fishery management and conservation units (see Dudgeon et al. 2012 and references therein). In the present study, we assess the levels of genetic diversity, population structure, and demographic history of the Pacific angel shark along the northwestern coast of Mexico based on mtCR.

Materials and methods

DNA extraction, amplification, and sequencing

Sampling was conducted at 2 localities where artisanal fishermen often land elasmobranchs: the port of Ensenada (27 samples) on the western coast of the Baja California Peninsula and El Portugués fishing camp (11 samples) on the eastern coast (Fig. 1). A small piece of muscular tissue was excised from each individual and stored at 4 ºC in absolute ethanol. Genomic DNA was extracted by a standard phenol-chloroform method (Aljanabi and Martinez 1997). The polymerase chain reaction (PCR) was carried out to amplify mtCR using universal primers: cytb-1040 (5 -GGN CAA CCM GTY GAA CAA CC-3 ) and 12 rev 326 (5 -ACT CGT ATA ACG CGG TGG CT-3 ) (J Sandoval-Castillo pers. comm.). Cycling conditions were as follows: an initial denaturation step at 96 ºC for 5 min; 35 cycles for 60 s each at 94 ºC, 60 s at 60 ºC, and 60 s at 72 ºC; and a final extension step for 10 min at 72 ºC. PCR products were purified using the QIAquick PCR Purification Kit (QIAGEN, Hilden, Germany). Both heavy and light strands were sequenced (Arizona Research Laboratories, Tucson, AZ) using direct sequencing reactions with a 3730XL DNA Analyzer (Applied Biosystems, Carlsbad, CA) and a BigDye Terminator Cycle Sequencing Kit (Applied Biosystems, Foster City, CA). The resulting DNA sequences were imported into ChromasPro v2.0.1 (Technelysium, Brisbane, Australia) and checked for quality and accuracy of nucleotide base assignment, and multiple sequence alignments were obtained using Clustal X implemented in MEGA v3.1 (Kumar et al. 2008).

Figure 1 Location of the 2 fishery localities (black circles) on the Baja California Peninsula (Mexico) where samples of Squatina californica were collected: the port of Ensenada and El Portugués fishing camp. 

Analysis of mtCR

Molecular diversity indices for the mtCR fragment, including the number of polymorphic sites, nucleotide composition, haplotype diversity (h), nucleotide diversity, and nucleotide divergence, were estimated for each sample using Arlequin v3.1 (Excoffier et al. 2005) and DnaSP v5.0 (Librado and Rozas 2009).

The level of differentiation among populations was assessed by estimating the pairwise F ST values for mtCR, based on the mean number of pairwise differences among sequences, using Arlequin. The hierarchical genetic structure of populations was assessed for the mtCR by an analysis of molecular variance (AMOVA) implemented in Arlequin. Statistical significance was estimated using 10,000 permutations. The Bayesian phylogeographic and ecological clustering (BPEC) method (Manolopoulou et al. 2011, Manolopoulou and Emerson 2012) was used to detect genetically and geographically distinct population clusters. This test was implemented in R 3.3.1 (R Development Core Team). Different tree and migration event scenarios were explored through the Markov chain Monte Carlo procedure to estimate the geographically clustered structure. Migration events were assumed to occur when a haplotype (with or without a mutation from its parent haplotype) migrated to a new cluster. The inputs for the BPEC analysis were 107 iterations, 2 dimensions, 1 parsimony relaxation parameter, and a maximum of 3 migrations.

Network and phylogenetic analyses

The relationships among the mtCR haplotypes were represented with a median-joining network (Bandelt et al. 1999) that was constructed using DNA Alignment v1.2 and PopART v3.0 (AWC 2007).

Phylogenetic analysis was performed by Bayesian inference (Huelsenbeck and Ronquist 2001, Huelsenbeck et al. 2001). The Bayesian phylogeny was estimated using the GTR + I + G model as the substitution model selected with the Bayesian information criterion implemented in MrModeltest v2.3 (Nylander 2004). Bayesian tree inferences were performed using MrBayes v3.1.2 (Ronquist and Huelsenbeck 2003). We used Markov chain Monte Carlo sampling with 4 chains, which ran simultaneously for 107 iterations with trees sampled every 10,000 generations. The first 25% of the trees were discarded as burn-in and the remaining trees were used to construct a 70% majority-rule concensus phylogeny. For phylogenetic reconstruction Squatina formosa was used as an outgroup (GenBank accession number: KM084865.1).

Demographic history

Two methods were applied to the mtCR data set to explore the demographic history of the Pacific angel shark samples that showed genetic differentation. First, evidence for population expansion was evaluated with Tajima’s D (Tajima 1989) and Fu’s F S (Fu 1997) statistics using Arlequin. Second, to further explore the possibility of demographic changes, a mismatch distribution analysis (distribution of pairwise differences) (Rogers and Harpending 1992) was performed for both geographical groups in Arlequin. Mismatch analysis models assume that any population growth or decline will reveal a genetic signature (Rogers and Harpending 1992, Schneider and Excoffier 1999).

Results

Genetic analysis

A fragment of 980 bp from the mtCR was sequenced for all samples (GenBank accession numbers: KP792162to KP792190). The nucleotide composition of this fragment consisted of 21.59% cytosine, 35.11% thymine, 29.66% adenine, and 13.64% guanine. Altogether, 29 haplotypes were found (Table 1) . Overall haplotype and nucleotide sequence diversities were 0.977 and 0.013, respectively. A total of 19 and 10 haplotypes occurred at the Gulf of California (El Portugués) and Pacific (Ensenada) sites, respectively. The most common haplotypes, GCF5 and GCF7, occurred at the gulf site. There were no shared haplotypes between sampling areas.

Table 1 Genetic diversity indices based on the mitochondrial control region sequences of Squatina californica from the Gulf of California and the Pacific coast of the Baja California Peninsula: n, number of samples; PS, number of polymorphic sites; SN, number of substitutions; Tr, number of transitions; Tv, number of transversions; NH, number of haplotypes found; h, haplotype diversity; π, nucleotide diversity. The asterisk indicates statistical significance at P < 0.05. 

Sampling area n PS SN Tr Tv NH h π D 1 F S 2
Gulf of California 27 52 54 18 36 19 0.982 0.008 -1.978* -2.985*
Pacific Ocean 11 90 90 17 73 10 0.957 0.055 -1.875* -2.377*
Total 38 127 134 33 101 29 0.977 0.013

1 Tajima (1989); 2Fu (1997)

The results of the molecular variance and pairwise analyses were concordant. The null hypothesis of global panmixia across the study area was rejected because there was significant genetic heterogeneity between sampling areas (Table 2). The pairwise F ST value was significant (F ST = 0.123, P < 0.01) for the pairwise comparison between the Pacific and Gulf of California samples. This genetic differentiation is also consistent with the number of clusters identified by BPEC. Based on the mtCR data set, BPEC revealed the highest probability (P = 0.97) at 2 clusters, indicating this as the best estimate of the true number of genetic clusters, which clearly correspond to the gulf and Pacific areas (Fig. 2).

Table 2 Variance of pairwise genetic distances and estimated nucleotide divergence (K) based on the mitochondrial control region sequences of Squatina californica from waters off northwestern Mexico. 

Source of variation Variane % total P F ST K (%)
Among populations 3.77 44.51 0.00 0.4451 2.13
Within populations 4.70 55.49

Figure 2 Results from Bayesian phylogeographic and ecological clustering analysis based on the mitochondrial control region sequences of Squatina californica from the Pacific coast of the Baja California Peninsula and the Gulf of California. The simulated contour areas are centered for each population cluster, and the shaded areas show the radius of 50% concentration contours around it. 

The Bayesian phylogenetic tree for the S. californica samples collected off northwestern Mexico showed 2 clusters (Fig. 3), the G cluster comprising the gulf samples and the P cluster comprising the Pacific samples. Both clusters were reflected in the haplotype network because no haplotypes were shared between the Pacific and the gulf samples (Fig. 4).

Figure 3 Bayesian phylogenetic reconstruction based on the mitochondrial control region sequences of Squatina californica. Posteriori probabilities are indicated in the nodes. Squatina formosa (SF) was used as the outgroup. P, sample from the Pacific coast of the Baja California Peninsula; G, sample from the Gulf of California. Colors show the clusters according to the Bayesian phylogeographic and ecological clustering results. 

Figure 4 Haplotype network based on the mitochondrial control region sequences of Squatina californica. Black circles indicate haplotypes from the Pacific coast of the Baja California Peninsula (PCM) and blue circles indicate haplotypes from the Gulf of California (GCF). Relative area of circles reflects sample size and the red nodes represent possible haplotype ancestors. Concerning bars represent mutation steps. 

Demographic history

The unimodal mismatch distribution analysis (Fig. 5a), together with the negative values of Tajima’s D and Fu’s F S statistics (Table 1), indicated a recent population expansion scenario for the gulf population. In contrast, the mismatch distribution analysis showed a multimodal distribution for the Pacific population, suggesting demographic stability (Fig. 5b).

Figure 5 Pairwise mismatch distribution of relative effective population size of Squatina californica based on the mitochondrial control region sequences of samples from the (a) Gulf of California and (b) Pacific coast of the Baja California Peninsula. The solid line represents expected frequency and the dotted line represents observed frequency. 

Discussion

The present study provides a first approach to learning the genetic diversity, population structure, and demographic history of S. californica off northwestern Mexico using mtCR. The mtCR sequence of this shark species is shorter than that observed in most cartilaginous fishes (~1,200 bp, Stoner et al. 2003). Overall haplotypic diversity was high, mainly due to the presence of many unique haplotypes. These high diversity values are similar to those reported for Pacific angel sharks in the Gulf of California (h = 0.97; Grijalva-Chon et al. 2002). High haplotype and nucleotide diversity could likely be due to high mutation rates commonly recorded in the mtCR (Randi 2000).

Given the significant levels of genetic heterogeneity and differentiation and the lack of shared haplotypes between the Pacific and gulf samples, we suggest that there has not been sufficient gene flow between both areas. In fact, this population division was observed in the global phylogenetic analysis of angel shark (genus Squatina) samples performed by Stelbrink et al. (2009). The regional genetic variation observed in our study could be linked to ontogenetic differences in S. californica between study areas. Length at sexual maturity (LSM) and maximum total length (MTL) are lower in individuals from the gulf area (LSM = 70.0-77.7 cm; MTL = 88-99 cm) than in individuals from the Pacific area (LSM = 99 cm; MTL = 117 cm) (Villavicencio-Garayzar 1996, Romero-Caicedo et al. 2016). Morphological differences and genetic differentiation are both in agreement with the population structure of this species in the study area.

The genetic differentiation observed in the Pacific angel shark can be attributed to the low dispersal ability of this species because of its life history traits, particularly its preference for coastal and benthic habitats and the lack of pelagic development stages such as the larval stage (Compagno 2005). In general, angel shark species have limited ability for sustained swimming due to their morphological and anatomical features, which limit swim power and suggest low dispersal potential and short migration events (Pittenger 1984). Taking into account the characteristics of S. californica, we consider that the submarine canyons constitute a natural barrier for its dispersal. Geomorphology off northwestern Mexico is complex, and there are numerous submarine canyons and great depths (at the mouth of the Gulf of California, the ocean basin drops to over 3,000 m) (Gutiérrez et al. 2004). Furthermore, the continental shelf along the Pacific coast of Baja California is typically narrow (usually <20 km), with great depths close to it (Allen and Horn 2006). The effect of deep marine basins as barriers to dispersal has also been suggested for the Pacific angel shark in waters between the northern and southern Channel Islands (California) (Gaida 1997) and even for other demersal elasmobranchs such as the thornback ray (Raja clavata, Chevolot et al. 2006) and the nurse shark (Ginglymostoma cirratum, Karl et al. 2012).

The Baja California Peninsula has been reported as an important biogeographic barrier for several species of elasmobranchs such as the shovelnose guitarfish (Rhinobatos productus, Sandoval-Castillo et al. 2004), California butterfly ray (Gymnura marmorata, Smith et al. 2009), golden cownose ray (Rhinoptera steindachneri, Sandoval-Castillo and Rocha-Olivares 2011), and banded guitarfish (Zapteryx exasperata , Castillo-Páez et al. 2014). The dynamic history of the area as a transition zone (Álvarez-Borrego 1983) could have produced multiple events of population fragmentation, likely influencing the genetic differentiation in elasmobranch populations (Sandoval-Castillo and Beheregaray 2015).

Reproductive behavior may also promote the observed genetic differentiation of Pacific angel shark populations. Philopatric behavior reduces the possibility of gene flow among breeding areas (Dudgeon et al. 2012). Female reproductive philopatry to nursery areas has been recorded for several elasmobranchs, such as the zebra shark (Stegostoma fasciatum, Dudgeon et al. 2009), thornback ray (R. clavata, Chevolot et al. 2006), and small-spotted catshark (Scyliorhinus canicula, Kousteni et al. 2015), and it even has been associated with the genetic variation in the banded guitarfish (Z. exasperata, Castillo-Páez et al. 2014) between the Pacific and gulf areas. However, no tagging or biological marker studies are available to determine whether S. californica returns to the same breeding areas every reproductive season. In this sense, tagging studies are required to assess reproductive philopatry.

In the present study, we found private or unique haplotypes for both sampling areas probably due to the small sample size (Kalinowski 2004). An adequate sampling scheme is needed to properly analyze biological questions regarding geographic genetic structure (Mona et al. 2014). Small sample sizes can overestimate population division and fail to yield statistical significance in the pairwise AMOVA (Sinclair and Hobbs 2008). Because of this and given the genetic variation, we consider that fine-scale sampling and larger sample sizes are required to detect intermediate populations.

The demographic analyses suggest a sudden growth for the gulf population, which could be associated with the formation of the Gulf of California. The formation of the Gulf of California (Pleistocene, Jacobs et al. 2004) provided new habitat and propitious conditions. The Pacific population showed patterns of population stability, but the mismatch distribution showed a ragged distribution (Fig. 5a), indicating a population in disequilibrium that is less likely to expand. A similar distribution has been documented for the Guadalupe fur seal (Arctocephalus townsendi, Weber et al. 2004), suggesting a substantial loss of genetic variability associated with a recent population bottleneck. Divergence time studies using more molecular markers are nevertheless required to adequately assess the demographic history of both Pacific angel shark populations.

Our results have implications for the design of management and conservation policies because they suggest 2 populations or stocks for the Pacific angel shark off northwestern Mexico. The adoption of genetic and evolutionary criteria in the management of natural resources has led to the recognition of management units, which represent functionally independent populations, or groups of populations, characterized by low levels of gene flow (Moritz 1994). In this sense, we suggest that both populations should be considered as separate management units. Future stock assessment should be done separately for each population. In Mexico, however, tools for legislation (e.g., NOM- 029-PESC) do not account for population or stock structure or population dynamics of shark resources. An inadequate management strategy can cause a stock to collapse or lead to local extirpation with ecological and economic consequences (Dulvy and Forrest 2010).

In particular, demographic analyses for the Pacific population of S. californica suggests population stability but with tendency to decline, likely due to the effect of fisheries. Reported landing data (2003-2015) from Baja California suggest that Pacific angel shark landings declined by >99% (Cailliet et al. 2016). This shark species has been heavily exploited off the western coast of the Baja California Peninsula (e.g., San Ignacio Lagoon and Magdalena Bay) over the past 3 decades (mainly in the 1980s and 1990s; Villavicencio-Garayzar and Abitia-Cardenas 1994, Leet et al. 2001). A recent study of elasmobranch fisheries along the western coast of the Baja California Peninsula (Ramírez-Amaro et al. 2013) found that Pacific angel sharks were not reported in previously known landing areas. It is therefore necessary to pay attention to the current and future effects of fishing on the population of S. californica off the western coast of the Baja California Peninsula.

Acknowledgments

The authors acknowledge funding support from the National Council for Science and Technology (CONACYT, Mexico). FGM thanks IPN for the granted fellowships (COFFA, EDI). We thank Leonardo Castillo Géniz (INAPESCA), Daniel Cartamil (Scripps Institution of Oceanography), IEMANYA OCEANICA A.C., and the CICIMAR elasmobranch team for collecting the samples. Hever Latisnére provided assistance with laboratory analyses. Adrián Munguía helped with sequencing procedures. Gerardo González Barba provided helpful comments. Two anonymous reviewers and the editor are also greatly acknowledged for their constructive comments that improved the quality of the manuscript.

References

Aljanabi SM, Martinez I. 1997. Universal and rapid salt-extraction of high quality genomic DNA for PCR-based techniques. Nucleic Acids Res. 25(22): 4692-4693. [ Links ]

Allen LG, Horn MH. 2006. The Ecology of Marine Fishes. University of California Press, Berkeley, CA, 670 pp. [ Links ]

Álvarez-Borrego J. 1983. Gulf of California. In: Ketchum BH (ed.), Estuaries and Enclosed Seas. Ecosystems of the World. Elsevier Publishing, Oxford, pp. 427-449. [ Links ]

Arenas M, Ray N, Currat M, Excoffier L. 2012. Consequences of range contractions and range shifts on molecular diversity. Mol. Biol. Evol. 29: 207-218. [ Links ]

Avise JC. 1994. Molecular Markers, Natural History and Evolution. Chapman and Hall, Springer, New York, 511 pp. [ Links ]

[AWC] Allan Wilson Centre. 2007. PopART (Population Analysis with Reticulate Trees). Free Software Foundation, Inc. [Accessed 21 February 2016] http://popart.otago.ac.nz/index.shtmlLinks ]

Bandelt HJ, Forster P, Rol A. 1999. Median-Joining networks for inferring intraespecific phylogenies. Mol. Biol. Evol. 16: 37-48. [ Links ]

Cailliet GM, Chabot CL, Nehmens MC, Carlisle AB. 2016. Squatina californica. The IUCN Red List of Threatened Species 2016: e.T39328A80671059. [Accessed 24 July 2016] http://dx.doi.org/10.2305/IUCN.UK.2016-2.RLTS.T39328A80671059.en [ Links ]

Castillo-Páez A, Sosa-Nishizaki O, Sandoval-Castillo J, Galván-Magaña F, Rocha-Olivares A. 2014. Strong population structure and shallow mitochondrial phylogeny in the banded guitarfish, Zapteryx exasperata (Jordan and Gilbert 1880), from the northern Mexican Pacific. J. Hered. 105(1): 91-100. http://dx.doi.org/10.1093/jhered/est067 [ Links ]

Chevolot M, Hoarau G, Rijnsdorp AD, Stam WT, Olsen JL. 2006. Phylogeography and population structure of thornback rays (Raja clavata L. Rajidae). Mol. Ecol. 15: 3693-3705. [ Links ]

Compagno L, Dando M, Fowler S. 2005. Sharks of the World. Princeton University Press, 496 pp. [ Links ]

Dudgeon CL, Blower DC, Broderick D, Giles JL, Holmes BJ, Kashiwagi T, Krück NC, Morgan JAT, Tillett BJ, Ovenden JR. 2012. A review of the application of molecular genetics for fisheries management and conservation of sharks and rays. J. Fish Biol. 80(5): 1789-843. http://dx.doi.org/10.1111/j.1095-8649.2012.03265.x [ Links ]

Dudgeon CL, Broderick D, Ovenden JR. 2009. IUCN classification zones concord with, but underestimate, the population genetic structure of the zebra shark Stegostoma fasciatum in the Indo-West Pacific. Mol. Ecol. 18(2): 248-261. http://dx.doi.org/10.1111/j.1365-294X.2008.04025.x [ Links ]

Dulvy NK, Forrest RE. 2010. Life histories, population dynamics, and extinction risks in Chondrichthyans. In: Carrier JC, Musick JA, Heithaus MR (eds.), Sharks and their Relatives II: Biodiversity, Adaptive Physiology, and Conservation. CRC Press, Boca Raton, FL, pp. 639−679. [ Links ]

Excoffier L, Laval G, Schneider S. 2005. ARLEQUIN v. 3.1: An integrated software package for population genetics data analysis. Evol. Bioinform. Online 1: 47-50. [ Links ]

Fu YX. 1997. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics 147(2): 915-925. [ Links ]

Gaida IH. 1997. Population structure of the Pacific angel shark, Squatina californica (Squatiniformes: Squatinidae), around the California Channel Islands. Copeia 4: 738-744. [ Links ]

Grijalva-Chon JM, Kaichi A, Numachi K. 2002. Homogeneidad genética en tiburón angelito (Squatina californica) del Golfo de California, evidencia por análisis PCR-RFLP de la región control del ADN mitocondrial. Ciencia y Mar 17: 37-42. [ Links ]

Gutiérrez OQ, Marinone SG, Parés-Sierra A. 2004. Lagrangian surface circulation in the Gulf of California from a 3D numerical model. Deep-Sea Res. II 51: 659-672. [ Links ]

Fogarty MJ, Botsford LW. 2007. Population connectivity and spatial management of marine fisheries. Oceanography 20: 112-123. [ Links ]

Heist J, Musick JA, Graves JE. 1996. Genetic population structure of the shortfin mako (Isurus oxyrinchus) inferred from restriction fragment length polymorphism analysis of mitochondrial DNA. Can. J. Fish. Aquat. Sci. 53(3): 583-588. [ Links ]

Hilário A, Metaxas A, Gaudron SM, Howell KL, Mercier A, Mestre NC, Ross RE, Thurnherr AM, Young C. 2015. Estimating dispersal distance in the deep sea: Challenges and applications to marine reserves. Front. Mar. Sci. 2: 1-14. [ Links ]

Hoelzel AR, Shivji MS, Magnussen J, Francis MP. 2006. Low worldwide genetic diversity in the basking shark (Cetorhinus maximus). Biol. Lett. 2: 639-642. [ Links ]

Huelsenbeck JP, Ronquist F. 2001. MrBayes: Bayesian estimation of phylogenetic tree. Bioinformatics 17(8): 754-755. https://doi.org/10.1093/bioinformatics/17.8.754 [ Links ]

Huelsenbeck JP, Ronquist F, Nielsen R, Bollback JP. 2001. Bayesian inference of phylogeny and its impact on evolutionary biology. Science 294: 2310-2314. [ Links ]

Hutchings JA, Reynolds JD. 2004. Marine fish population collapses: Consequences for recovery and extinction risk. Bioscience 54: 297-309. [ Links ]

Jacobs DK, Haney TA, Louie KD. 2004. Genes, diversity, and geologic process on the Pacific coast. Annu. Rev. Earth Planet. Sci. 32: 601-652. [ Links ]

Kalinowski ST. 2004. Counting alleles with rarefaction: Private alleles and hierarchical sampling designs. Conserv. Genet. 5: 539-543. [ Links ]

Karl SA, Castro ALF, Garla RC. 2012. Population genetics of the nurse shark (Ginglymostoma cirratum) in the western Atlantic. Mar. Biol. 159: 489-498. [ Links ]

Kousteni V, Kasapidis P, Kotoulas G, Megalofonou P. 2015. Strong population genetic structure and contrasting demographic histories for the small-spotted catshark (Scyliorhinus canicula) in the Mediterranean Sea. Heredity 114(3): 333-343. http://dx.doi.org/10.1038/hdy.2014.107 [ Links ]

Kumar S, Dudley J, Nei M, Tamura K. 2008. MEGA: biologist-centric software for evolutionary analysis of DNA and protein sequences. Brief. Bioinform. 9: 299-306. [ Links ]

Leet WS, Dewees CM, Klingbeil R, Larson EJ. 2001. California’s living marine resources: A status report. The Resources Agency, California Department of Fish and Game, Sacramento, CA, 585 pp. [ Links ]

Librado P, Rozas J. 2009. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25: 1451-1452. [ Links ]

Manolopoulou I, Emerson BC. 2012. Phylogeographic ancestral inference using the coalescent model on haplotype trees. J. Comp. Biol. 19: 745-755. [ Links ]

Manolopoulou I, Legarreta L, Emerson BC, Brooks S, Tavaré S. 2011. A Bayesian approach to phylogeographic clustering. Interface Focus 1: 909-921. [ Links ]

Mona S, Ray N, Arenas M, Excoffier L. 2014. Genetic consequences of habitat fragmentation during a range expansion. Heredity 11: 291-299. [ Links ]

Moritz C. 1994. Defining ‘evolutionarily significant units’ for conservation. Trends Ecol. Evol. 9(10): 373-375. http://dx.doi.org/10.1016/0169-5347(94)90057-4 [ Links ]

Nielsen JL, Powers DA. 1995. Evolution and the aquatic ecosystem: Defining unique units in population conservation. Am. Fish. Soc. Symp. 435 pp. [ Links ]

Nylander JAA. 2004. MrModeltest v2. Program distributed by the author. Evolutionary Biology Center, Uppsala University. [ Links ]

Pittenger G. 1984. Movements, distributions, feeding, and growth of the Pacific angel sharks, Squatina californica, off Santa Barbara, California. Copeia 4: 987-994. [ Links ]

Ramírez-Amaro SR, Cartamil D, Galván-Magaña F, González-Barba G, Graham JB, Carrera-Fernández M, Escobar-Sánchez O, Sosa-Nishizaki O, Rochin-Alamillo A. 2013. The artisanal elasmobranch fishery of the Pacific coast of Baja California Sur, México, management implications. Sci. Mar. 77(3): 473-487. [ Links ]

Randi E. 2000. Mitochondrial DNA. In: Baker AJ (ed.), Molecular Methods in Ecology. Blackwell Science, Malden, MA, pp. 136-167. [ Links ]

Rogers AR, Harpending H. 1992. Population growth makes waves in the distribution of pairwise differences. Mol. Biol. Evol. 9: 552-569. [ Links ]

Romero-Caicedo AF, Hernández-Herrera A, Carrera-Fernández M, Galván-Magaña F. 2016. Reproductive parameters of the Pacific angel shark, Squatina californica (Selachii: Squatinidae), in the northern Mexican Pacific. J. Fish Biol. 88: 1430-1440. [ Links ]

Ronquist F, Huelsenbeck JP. 2003. Mr Bayes 3 v3.1.2: Bayesian phylogenetic inference under mixed models. Bioinformatics 19: 1572-1574. [ Links ]

Sandoval-Castillo J, Beheregaray LB. 2015. Metapopulation structure informs conservation management in a heavily exploited coastal shark (Mustelus henlei). Mar. Ecol. Prog. Ser. 533: 191-293. [ Links ]

Sandoval-Castillo J, Rocha-Olivares A. 2011. Deep mitochondrial divergence in Baja California populations of an aquilopelagic elasmobranch: The golden cownose ray. J. Hered. 102(3): 269-274. [ Links ]

Sandoval-Castillo J, Rocha-Olivares A, Villavicencio-Garayzar C, Balart E. 2004. Cryptic isolation of Gulf of California shovelnose guitarfish evidenced by mitochondrial DNA. Mar. Biol. 145: 983-988. [ Links ]

Schneider S, Excoffier L. 1999. Estimation of demographic parameters from the distribution of pairwise differences when the mutation rates vary among sites: Application to human mitochondrial DNA. Genetics 152: 1079-1089. [ Links ]

Sinclair EA, Hobbs R. 2008. Sample size effects on estimates of population genetic structure: Implications for ecological restoration. Restor. Ecol. 17: 837-844. [ Links ]

Smith WD, Bizzarro JJ, Richards VP, Nielsen J, Márquez-Farías F, Shivji MS. 2009. Morphometric convergence and molecular divergence: The taxonomic status and evolutionary history of Gymnura crebripunctata and Gymnura marmorata in the eastern Pacific Ocean. J. Fish Biol. 75: 761-783. [ Links ]

Spaet JLY, Jabado RW, Henderson AC, Moore ABM, Berumen ML. 2015. Population genetics of four heavily exploited shark species around the Arabian Peninsula. Ecol. Evol. 5: 2317-2332. [ Links ]

Stelbrink B, von Rintelen T, Cliff G, Kriwet J. 2009. Molecular systematics and global phylogeography of angel shark (genus Squatina). Mol. Phylogenet. Evol. 2: 395-404. [ Links ]

Stoner DS, Grandy JM, Priede KA, Quattro JM. 2003. Amplification primer for the mitochondrial control region and sixth intron of the nuclear-encoded lactate dehydrogenase A gene in elasmobranch fishes. Conserv. Genet. 4: 805-808. [ Links ]

Tajima F. 1989. Statistical methods to test for nucleotide mutation hypothesis by DNA polymorphism. Genetics 123: 585-595. [ Links ]

Veríssimo A, McDowell JR, Graves JE. 2011. Population structure of a deep-water squaloid shark, the Portuguese dogfish (Centroscymnus coelolepis). ICES J. Mar. Sci. 68: 555-563. [ Links ]

Vignaud TM, Maynard JA, Leblois R, Meekan MG, Vázquez-Juárez R, Ramírez-Macías D, Pierce S, Rowat D, Berumen ML, Beeravolu C, Baksay S, Planes S. 2014. Genetic structure of populations of whale sharks among ocean basins and evidence for their historic rise and recent decline. Mol. Ecol. 23: 2590-2601. [ Links ]

Villavicencio-Garayzar CJ. 1996. Aspectos poblacionales del angelito, Squatina californica, Ayres, en Baja California, México. Rev. Invest. Cient. Ser. Cienc. Mar. UABCS 1: 15-21. [ Links ]

Villavicencio-Garayzar CJ, Abitia-Cardenas LA. 1994. Elasmobranquios de Bahía Magdalena y Laguna San Ignacio, Baja California Sur, México. Rev. Invest. Cient. Ser. Cienc. Mar. UABCS 5: 63-67. [ Links ]

Weber DS, Stewart BS, Lehman N. 2004. Genetic consequences of a severe population bottleneck in the Guadalupe fur seal (Arctocephalus townsendi). J. Hered. 95: 144-153. [ Links ]

Received: November 01, 2016; Accepted: February 01, 2017

* Corresponding author. E-mail: rvazquez04@cibnor.mx

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License