The complete mitochondrial genome of the sycamore lace bug Corythucha ciliata (Hemiptera: Tingidae)

The complete mitochondrial genome of the sycamore lace bug, Corythucha ciliata, was sequenced in this study. It represents the first sequenced mitogenome of family Tingidae in Heteroptera. The mitogenome of C. ciliata is 15,257 bp and contains 37 genes including 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes and a large non-coding region. Gene arrangement, nucleotide content, codon usage, and amino acid composition and asymmetry indicate a high degree of conservation with six other species of Cimicomorpha. The 13 PCGs initiated with ATN as the start codon and terminated with TAA, TA or T as stop codon. The evolutionary rate of each PCG was different, among which ATP8 showed the highest rate while ATP6 indicated the lowest rate. The 22 tRNAs genes apparently fold into a typical cloverleaf structure; however, the anticodon (TTC) of trnSer (AGN) differs from other Heteropteran insects. Secondary structure modeling of rRNA genes revealed similarity to other insects, except for two incomplete helices (H1648 and H2735) in lrRNA. The predicted secondary structure of lrRNA indicates 45 helices in six domains, whereas srRNA has 27 helices in three domains. Three potential stem–loops and two tandem repeats (–TCTAAT–) were identified in the A+T-rich region. Phylogenetic analysis indicated that C. ciliata is a sister group to other Heteroptera species based on analysis of the 13 PCGs.     

Complete mitochondrial genome sequence of bighead croaker Collichthys niveatus (Perciformes, Sciaenidae): a mitogenomic perspective on the phylogenetic relationships of Pseudosciaeniae

The monophyly and phylogenetic relationships of Pseudosciaeniae have long been controversial. Here we describe the mitochondrial genome (mitogenome) sequence of Collichthys niveatus. It is a circular double-stranded DNA molecule of 16,450 base pairs (bp) in length with a standard set of 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs), 13 protein-coding genes as well as a non-coding control region. The mitogenome of C. niveatus shared common features with those of other bony fishes in terms of gene arrangement, base composition, and tRNA structures. The C. niveatus mitogenome exhibited pronounced strand-specific asymmetry in nucleotide composition, which was also reflected in the codon usage of genes oriented in opposite directions. Contrary to the typical structure of the control region, the central conserved blocks (CSB-D, -E, and -F) could not be detected in C. niveatus mitogenome. Phylogenetic analysis based on whole mitogenome sequences provided strong support for the monophyly of Pseudosciaeniae, and sister-group relationships of C. niveatus + Collichthys lucidus and Larimichthys crocea + Larimichthys polyactis, which was consistent with the traditional taxonomy. Unexpected divergence was found in two C. niveatus mitogenomes and several hypotheses were proposed to explain this observation including misidentification and introgressive hybridization between C. niveatus and L. polyactis, and polyphyletic origin of C. niveatus. We considered species misidentification to be the main hypothesis. However, additional data is essential to test these proposed hypotheses.     

Skeletal anatomy of the Late Cretaceous shark,Squalicorax (Neoselachii: Anacoracidae)

The paleobiology of the Cretaceous neoselachian shark,Squalicorax, has largely been based on isolated teeth. We examined partial and nearly complete skeletons of three species ofSqualicorax, S. falcatus (Aoassiz),S. kaupi (Agassiz), andS. pristodontus (Agassiz), that were collected from the U.S.A. These specimens suggest that the total body length (TL) ofS. falcatus typically measured 1.8–2.0 m, and probably did not exceed 3 m. Moderatesized individuals ofS. kaupi andS. pristodontus perhaps measured about 3 m TL. AlthoughS. pristodontus was the largest form among the three species examined, this taxon possessed a set of large jaws (with large but fewer teeth) relative to its body size compared toS. falcatus orS. kaupi. This suggests that tooth size is not an accurate indicator of the TL if one compares oneSqualicorax species to another. Neurocranial features suggest that the vision ofSqualicorax was not as acute as that of a contemporaneous macrophagous lamniform shark,Cretoxyrhina mantelli (Agassiz) , but olfaction ofSqualicorax may have been better thanC. mantelli. The morphology of placoid scales suggests thatSqualicorax was capable of fast swimming. New skeletal data support the view that the feeding dynamics ofSqualicorax was similar to the modern tiger shark (Galeocerdo Müller & Henle). The present data do not allow for exact ordinal placement, but, contrary to some previous interpretations,Squalicorax can be excluded from the Hexanchiformes and Orectolobiformes. The taxon should more appropriately be placed within the Lamniformes or Carcharhiniformes.      

Identification and characterization of polymorphic microsatellite loci in the blue shark Prionace glauca, and cross-amplification in other shark species

Two to 14 alleles were found to be segregating per locus (mean 5·2), with observed and expected heterozygosities ranging from 0·08 to 0·78 and 0·08 to 0·94, respectively. Cross-amplification of six of these microsatellite loci indicated that they are also polymorphic in three species of Carcharhiniformes and two species of Lamniformes. The newly developed primers reported here constitute a useful tool for genetic population analyses on Prionace glauca and, potentially, other related species.     

Length-weight relationships of eight elasmobranch species captured along the Coromandel coast of Tamil Nadu,Eastern Indian Ocean

This study reports the length-weight relationships (LWR) for eight elasmobranch species; four shark species (Carcharhiniformes and Orectolobiformes), two species of rays (Myliobatiformes) and two species of guitarfishes (Rhinopristiformes) from the Coromandel Coast of Tamil Nadu, Eastern Indian Ocean. The specimens of these elasmobranch species were collected from the catches of various gears like trawl net, bottom set and drift gill net, hook & line at Chennai—Royapuram Fishing Harbour, Cuddalore Fishing Harbourand Nagapattinam Fishing Harbour situated along this coastal region fortnightly during June 2019–March 2020. The values of the parameter ‘b’ remained within the expected range of 2.5–3.5. Length- weight / Disc-width-weight relationships showed good fit with r² values varying from 0.8923 to 0.9869. This study also reports a new maximum TL length (L_max)for a shark species (Chiloscyllium griseum).     

Shark tales: a molecular species-level phylogeny of sharks (Selachimorpha, Chondrichthyes)

Sharks are a diverse and ecologically important group, including some of the ocean's largest predatory animals. Sharks are also commercially important, with many species suffering overexploitation and facing extinction. However, despite a long evolutionary history, commercial, and conservation importance, phylogenetic relationships within the sharks are poorly understood. To date, most studies have either focused on smaller clades within sharks, or sampled taxa sparsely across the group. A more detailed species-level phylogeny will offer further insights into shark taxonomy, provide a tool for comparative analyses, as well as facilitating phylogenetic estimates of conservation priorities. We used four mitochondrial and one nuclear gene to investigate the phylogenetic relationships of 229 species (all eight Orders and 31 families) of sharks, more than quadrupling the number of taxon sampled in any prior study. The resulting Bayesian phylogenetic hypothesis agrees with prior studies on the major relationships of the sharks phylogeny; however, on those relationships that have proven more controversial, it differs in several aspects from the most recent molecular studies. The phylogeny supports the division of sharks into two major groups, the Galeomorphii and Squalimorphii, rejecting the hypnosqualean hypothesis that places batoids within sharks. Within the squalimorphs the orders Hexanchiformes, Squatiniformes, Squaliformes, and Pristiophoriformes are broadly monophyletic, with minor exceptions apparently due to missing data. Similarly, within Galeomorphs, the orders Heterodontiformes, Lamniformes, Carcharhiniformes, and Orectolobiformes are broadly monophyletic, with a couple of species 'misplaced'. In contrast, many of the currently recognized shark families are not monophyletic according to our results. Our phylogeny offers some of the first clarification of the relationships among families of the order Squaliformes, a group that has thus far received relatively little phylogenetic attention. Our results suggest that the genus Echinorhinus is not a squaliform, but rather related to the saw sharks, a hypothesis that might be supported by both groups sharing 'spiny' snouts. In sum, our results offer the most detailed species-level phylogeny of sharks to date and a tool for comparative analyses.     

Feeding habits of the whale shark (<Emphasis Type="Italic">Rhincodon typus</Emphasis>) inferred by fatty acid profiles in the northern Mexican Caribbean

Whale shark (Rhincodon typus, Smith, 1828) is an endangered species with anthropogenic pressures due to increasing demand of encounter tourism activities. Research efforts to identify management and conservation strategies for this species are needed. The Northern Mexican Caribbean is one of the most important feeding aggregation sites of whale sharks worldwide. In this study, Mexican Caribbean whale shark feeding habits are assessed by means of fatty acid (FA) signature analysis, a biochemical non-destructive technique widely applied in trophic ecology studies. Sub-dermal tissue biopsies of 68 whale sharks and samples of their potential prey (zooplankton) were collected during 2010 and 2011 in two areas with high R. typus abundance. Zooplankton samples (n?=?17) were divided in two categories: mixed zooplankton (several groups of zooplankton) and fish eggs (> 95% of sample components were fish eggs). FA profiles of whale shark tissue sampled between years showed significant variability; while there was no intraspecific differences in FA signature related to sex, size and location. FA profiles of whale sharks and their potential prey were dominated by saturated fatty acids (SFA). R. typus FA signature was significantly different from that of mixed zooplankton; on the other hand, whale shark and fish egg FA profiles formed groups with overlapping values and registered high levels of oleic acid. PUFA average ω3/ ω6 ratio on whale shark FA profiles for both years was below 1. Arachidonic acid (ARA) percentage was higher in whale shark biopsies (13.2% in 2010, 6.8% in 2011) compared to values observed in fish eggs (2.0%) and mixed zooplankton (1.4%). Similarity between FA profiles of whale sharks and fish eggs, low levels of bacterial FA found in R. typus biopsies, as well as whale shark feeding behavior observations in the study area, suggest that R. typus is feeding mainly on surface zooplankton in Mexican Caribbean; however, elevated ARA percentages in whale shark samples may indicate that this species has complementary feeding sources, such as demersal zooplankton, which has been reported in other aggregation sites. Results obtained contribute to the knowledge of the whale shark trophic ecology in the area, but are inconclusive. Further studies are recommended to evaluate whale shark FA profiles from different tissues (muscle or blood); also, broader information is needed about zooplankton FA signature in the study area.     

Complete mitochondrial genome of the stonefly Cryptoperla stilifera Sivec (Plecoptera: Peltoperlidae) and the phylogeny of Polyneopteran insects

We present the complete mitogenome of a stonefly, Cryptoperla stilifera Sivec (Plecoptera; Peltoperlidae). The mitogenome was a circular molecule consisting of 15,633 nucleotides, 37 genes and a A + T-rich region. C. stilifera mitogenome was similar to Pteronarcys princeps mitogenome (Plecoptera; Pteronarcyidae). All transfer RNA genes (tRNAs) had typical cloverleaf secondary structures except for trnSer (AGN), where the stem-loop structure of the dihydrouridine (DHU) arm was missing. The A + T-rich region of C. stilifera had two stem-loops and each had two interlink. Three conserved sequence blocks (CSBs) were present in the A + T-rich regions of C. stilifera, Peltoperla tarteri and Peltoperla arcuata. Moreover, many polynucleotide stretches (Poly N, N = A, T and C) in the A + T-rich region of C. stilifera Phylogenetic relationships of Polyneopteran species were constructed based on the nucleotide sequences of 13 protein coding genes (PCGs). Both maximum likelihood (ML) and Bayesian inference (BI) analyses supported Grylloblattodea as the sister group to Plecoptera + Dermaptera and Embiidina and Phasmatodea as sister groups.     

Molecular phylogeny of elasmobranchs inferred from mitochondrial and nuclear markers

The elasmobranchs (sharks, rays and skates) being the extant survivors of one of the earliest offshoots of the vertebrate evolutionary tree are good model organisms to study the primitive vertebrate conditions. They play a significant role in maintaining the ecological balance and have high economic value. Due to over-exploitation and illegal fishing worldwide, the elasmobranch stocks are being decimated at an alarming rate. Appropriate management measures are necessary for restoring depleted elasmobranch stocks. One approach for restoring stocks is implementation of conservation measures and these measures can be formulated effectively by knowing the evolutionary relationship among the elasmobranchs. In this study, a total of 30 species were chosen for molecular phylogeny studies using mitochondrial cytochrome c oxidase subunit I, 12S ribosomal RNA gene and nuclear Internal Transcribed Spacer 2. Among different genes, the combined dataset of COI and 12S rRNA resulted in a well resolved tree topology with significant bootstrap/posterior probabilities values. The results supported the reciprocal monophyly of sharks and batoids. Within Galeomorphii, Heterodontiformes (bullhead sharks) formed as a sister group to Lamniformes (mackerel sharks): Orectolobiformes (carpet sharks) and to Carcharhiniformes (ground sharks). Within batoids, the Myliobatiformes formed a monophyly group while Pristiformes (sawfishes) and Rhinobatiformes (guitar fishes) formed a sister group to all other batoids.     

Development and microstructure of tooth histotypes in the blue shark,Prionace glauca (Carcharhiniformes: Carcharhinidae) and the great white shark,Carcharodon carcharias (Lamniformes: Lamnidae)

Elasmobranchs exhibit two distinct arrangements of mineralized tissues in the teeth that are known as orthodont and osteodont histotypes. Traditionally, it has been said that orthodont teeth maintain a pulp cavity throughout tooth development whereas osteodont teeth are filled with osteodentine and lack a pulp cavity when fully developed. We used light microscopy, scanning electron microscopy, and high‐resolution micro‐computed tomography to compare the structure and development of elasmobranch teeth representing the two histotypes. As an example of the orthodont histotype, we studied teeth of the blue shark, Prionace glauca (Carcharhiniformes: Carcharhinidae). For the osteodont histotype, we studied teeth of the great white shark, Carcharodon carcharias (Lamniformes: Lamnidae). We document similarities and differences in tooth development and the microstructure of tissues in these two species and review the history of definitions and interpretations of elasmobranch tooth histotypes. We discuss a possible correlation between tooth histotype and tooth replacement and review the history of histotype differentiation in sharks. We find that contrary to a long held misconception, there is no orthodentine in the osteodont teeth of C. carcharias. J. Morphol. 276:797–817, 2015. © 2015 Wiley Periodicals, Inc.     

Ensemble modeling of the potential distribution of the whale shark in the Atlantic Ocean

The whale shark (Rhincodon typus) is an endangered marine fish species which can be adversely affected by the fishing activities of the industrial purse seine fleet targeting tropical tuna. Tuna tend to aggregate around all types of floating objects, including whale sharks. We analyzed and modeled the spatial distribution and environmental preferences of whale sharks based on the presence and absence data from fishing observations in the Atlantic Ocean. We used a thorough multialgorithm analysis, based on a new presence–absence dataset, and endeavored to follow the most recent recommendations on best practices in species distribution modeling. First, we selected a subset of relevant variables using a generalized linear model that addressed multicollinearity, statistical errors, and information criteria. We then used the selected variables to build a model ensemble including 19 different algorithms. After eliminating models with insufficient performance, we assessed the potential distribution of whale sharks using the mean of the predictions of the selected models. We also assessed the variance among the predictions of different algorithms, in order to identify areas with the highest model consensus. The results show that several coastal regions and warm shallow currents, such as the Gulf Stream and the Canary and Benguela currents, are the most suitable areas for whale sharks under current environmental conditions. Future environmental projections for the Atlantic Ocean suggest that some of the suitable regions will shift northward, but current concentration areas will continue to be suitable for whale shark, although with less productivity, which could have negative consequences for conservation of the species. We discuss the implications of these predictions for the conservation and management of this charismatic marine species.     

Population genetic structure of Earth's largest fish, the whale shark (Rhincodon typus)

Large pelagic vertebrates pose special conservation challenges because their movements generally exceed the boundaries of any single jurisdiction. To assess the population structure of whale sharks (Rhincodon typus), we sequenced complete mitochondrial DNA control regions from individuals collected across a global distribution. We observed 51 single site polymorphisms and 8 regions with indels comprising 44 haplotypes in 70 individuals, with high haplotype (h = 0.974 +/- 0.008) and nucleotide diversity (pi = 0.011 +/- 0.006). The control region has the largest length variation yet reported for an elasmobranch (1143-1332 bp). Phylogenetic analyses reveal no geographical clustering of lineages and the most common haplotype was distributed globally. The absence of population structure across the Indian and Pacific basins indicates that oceanic expanses and land barriers in Southeast Asia are not impediments to whale shark dispersal. We did, however, find significant haplotype frequency differences (AMOVA, Phi(ST) = 0.107, P < 0.001) principally between the Atlantic and Indo-Pacific populations. In contrast to other recent surveys of globally distributed sharks, we find much less population subdivision and no evidence for cryptic evolutionary partitions. Discovery of the mating and pupping areas of whale sharks is key to further population genetic studies. The global pattern of shared haplotypes in whale sharks provides a compelling argument for development of broad international approaches for management and conservation of Earth's largest fish.     

Nuclear and mitochondrial DNA reveals isolation of imperilled grey nurse shark populations (Carcharias taurus)

Loss of sharks and other upper-trophic marine predators has sparked worldwide concern for the stability of ocean ecosystems. The grey nurse (ragged-tooth or sand tiger) shark ( Carcharias taurus ) is Vulnerable on a global scale, Critically Endangered in Australia and presumed extinct in parts of its historical range. We used 193 muscle and fin samples collected from six extant populations to assess global mtDNA and microsatellite diversity and the degree of global population genetic structure. Control region mtDNA diversity was low in every population, and two populations (eastern Australia and Japan) contained only a single mtDNA haplotype. Genetic signatures of recent losses of genetic variation were not yet apparent at microsatellite loci, indicating that this low mtDNA variation is not a result of anthropogenic population declines. Population differentiation was substantial between each population pair except Brazil and South Africa, F _ST values ranged from 0.050 to 0.699 and 0.100 to 1.00 for microsatellite and mitochondrial data respectively. Bayesian analysis clearly partitioned individuals into five of the populations from which they were sampled. Our data imply a low frequency of immigrant exchange among each of these regions and we suggest that each be recognized as a distinct evolutionary significant unit. In contrast to pelagic species such as whale shark and white shark that may cross ocean basins and where cooperative international efforts are necessary for conservation, grey nurse shark, like many coastal species, need to be managed regionally.     

Complete mitochondrial genome of the atlas moth, Attacus atlas (Lepidoptera: Saturniidae) and the phylogenetic relationship of Saturniidae species

Mitochondrial genome (mitogenome) can provide information for genomic structure as well as for phylogenetic analysis and evolutionary biology. In this study, we present the complete mitogenome of the atlas moth, Attacus atlas (Lepidoptera: Saturniidae), a well-known silk-producing and ornamental insect with the largest wing surface area of all moths. The mitogenome of A. atlas is a circular molecule of 15,282 bp long, and its nucleotide composition shows heavily biased towards As and Ts, accounting for 79.30%. This genome comprises 13 protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and an A + T-rich region. It is of note that this genome exhibits a slightly positive AT skew, which is different from the other known Saturniidae species. All PCGs are initiated by ATN codons, except for COI with CGA instead. Only six PCGs use a common stop codon of TAA or TAG, whereas the remaining seven use an incomplete termination codon T or TA. All tRNAs have the typical clover-leaf structure, with an exception for tRNA^Ser(AGN). The A. atlas A + T-rich region contains non-repetitive sequences, but harbors several features common to the Bombycoidea insects. The phylogenetic relationships based on Maximum Likelihood method provide a well-supported outline of Saturniidae, which is in accordance with the traditional morphological classification and recent molecular works.     

Multilocus sequence evaluation for differentiating species of the trematode Family Gastrothylacidae,with a note on the utility of mitochondrial COI motifs in species identification

Amphistomiasis, a neglected trematode infectious disease of ruminants, is caused by numerous species of amphistomes belonging to six families under the Superfamily Paramphistomoidea. In the present study, four frequently used DNA markers, viz. nuclear ribosomal 28S (D1–D3 regions), 18S and ITS2 and mitochondrial COI genes, as well as sequence motifs from these genes were evaluated for their utility in species characterization of members of the amphistomes' Family Gastrothylacidae commonly prevailing in Northeast India. In sequence and phylogenetic analyses the COI gene turned out to be the most useful marker in identifying the gastrothylacid species, with the exception of Gastrothylax crumenifer, which showed a high degree of intraspecific variations among its isolates. The sequence analysis data also showed the ITS2 region to be effective for interspecies characterization, though the 28S and 18S genes were found unsuitable for the purpose. On the other hand, sequence motif analysis data revealed the motifs from the COI gene to be highly conserved and specific for their target species which allowed accurate in silico identification of the gastrothylacid species irrespective of their intraspecific differences. We propose the use of COI motifs generated in the study as a potential tool for identification of these species.     

Mitogenomic insights into a recently described and rarely observed killer whale morphotype

Identifying evolutionary divergent taxonomic units, e.g. species and subspecies, is important for conservation and evolutionary biology. The ‘type D’ killer whale, Orcinus orca, is a rarely observed morphotype with a pelagic, circumpolar subantarctic distribution, making dedicated research and therefore taxonomic study extremely difficult to date. In this study, we used DNA target enrichment hybridisation capture coupled to high throughput sequencing, to obtain the first DNA sequence from the only known museum specimen of this recently described morphotype. The high coverage, complete mitogenome sequence was compared to a previously published global dataset of 139 individuals, indicating that this type is highly divergent to all previously genetically sequenced killer whale forms. The estimated divergence time (390,000 years ago) from its most recent common ancestor with other extant killer whale lineages was the second oldest split within the killer whale phylogeny. This study provides the first genetic support of type D potentially being a distinct subspecies or species of killer whale, although further samples are needed to identify whether there is monophyly of mitogenome sequences and whether nuclear DNA also indicates reproductive isolation. These findings also highlight the value of natural history museum collections and new technologies to investigate the taxonomy of rare, cryptic or difficult to access species.     

Characteristics of the mitochondrial genome of Rana omeimontis and related species in Ranidae: Gene rearrangements and phylogenetic relationships

The Omei wood frog (Rana omeimontis), endemic to central China, belongs to the family Ranidae. In this study, we achieved detail knowledge about the mitogenome of the species. The length of the genome is 20,120 bp, including 13 protein‐coding genes (PCGs), 22 tRNA genes, two rRNA genes, and a noncoding control region. Similar to other amphibians, we found that only nine genes (ND6 and eight tRNA genes) are encoded on the light strand (L) and other genes on the heavy strand (H). Totally, The base composition of the mitochondrial genome included 27.29% A, 28.85% T, 28.87% C, and 15.00% G, respectively. The control regions among the Rana species were found to exhibit rich genetic variability and A + T content. R. omeimontis was clustered together with R. chaochiaoensis in phylogenetic tree. Compared to R. amurensis and R. kunyuensi, it was more closely related to R. chaochiaoensis, and a new way of gene rearrangement (ND6‐trnE‐Cytb‐D‐loop‐trnL2 (CUN)‐ND5‐D‐loop) was also found in the mitogenome of R. amurensis and R. kunyuensi. Our results about the mitochondrial genome of R. omeimontis will contribute to the future studies on phylogenetic relationship and the taxonomic status of Rana and related Ranidae species.     

Seasonal trends in whale shark Rhincodon typus sightings in an established tourism site in the Gulf of California,Mexico

The number of individual whale shark Rhincodon typus sightings registered October 2015–March 2018 in Bahia de La Paz, Mexico, ranged from 73 to 129 animals per season (total 1662), with the majority of them identified as males. Density plot maps showed high correlations in sightings between sample seasons and revealed potential hotspot areas. Our study provides an essential baseline of information for the management of whale shark conservation and the associated touristic activities in the region.     

Phylogeography and conservation of the bull shark (Carcharhinus leucas) inferred from mitochondrial and microsatellite DNA

The bull shark (Carcharhinus leucas) is a widely distributed, large coastal shark species known to travel long distances. These characteristics, coupled with the species?? long life span and late age of maturity, would lead one to predict significant global genetic exchange among bull shark populations. By contrast, data show localized depletion in some areas of large coastal shark fisheries, indicating some geographic isolation may exist. We examined genetic variation in the control region of mitochondrial DNA and at five nuclear microsatellite loci in bull sharks sampled from the western Atlantic to investigate the degree of population subdivision. The average per sample haplotype and nucleotide diversity in the mtDNA (0.51 ± 0.26 and 0.12% ± 0.12, respectively) and expected heterozygosity (0.84) in the microsatellite loci contrast sharply in having lower and higher values (respectively) relative to many other shark species. Significant structure exists between the Brazilian and all northern populations at the mtDNA control region (pairwise ??_ST > 0.8, P < 0.001), but not at the nuclear microsatellite loci. Adjacent northern populations show weak to no genetic differentiation for both markers. These results are congruent with restricted maternal gene flow between populations caused by female site fidelity to nursery areas. We estimate the current effective population size to be around 160,000 and 221,000 individuals for the southern and northern Atlantic populations, respectively. The philopatric habits and the relatively low levels of mtDNA genetic diversity observed in bull sharks must be considered in the conservation of this species. Our results indicate that effective bull shark management strategies will require local, regional, and international attention and cooperation.