Comparative study of genetic variation at 15 STR loci in three isolated populations of the Bosnian mountain area

Fifteen autosomal STR loci (D3S1358, TH01, D21S11, D18S51, Penta E, D5S818, D13S317, D7S820, D16S539, CSF1PO, Penta D, VWA, D8S1179, TPOX, and FGA) were studied in three geographically close but isolated populations from the Bosnian mountain area. The three villages are Bobovica, Dejcici, and Lukomir. DNA was obtained from 83 individuals, and the allele frequencies and genetic diversity among the three sample groups were compared. In addition, seven of the STR loci (CSF1PO, D13S317, D3S1358, D5S818, D7S820, FGA, TH01) were used in a comparative population analysis of the Bjelasnica-Treskavica region and the Adriatic islands of Brac, Hvar, and Korcula. Although the sample sizes are relatively small, the observed variation within any of the small isolated populations is high and comparable to less isolated groups. In addition, even though the populations are geographically isolated, the STR data are similar among the populations. The most significant frequency differences were observed at the TH01 locus. Although the specific allele distributions in any untyped population cannot be determined a priori, we find support for a high degree of diversity for the STR loci in most populations. In addition, the multiple locus profile is highly informative not only for various population studies but also for forensic studies, even when specific population data are not available.     

Determination of combined sibship indices “gray zone” using 15 STR loci for central Bosnian human population

In our previous population studies of Bosnia and Herzegovina human population, we have used autosomal STR, Y-STR, and X-STR loci, as well as Y-chromosome NRY biallelic markers. All obtained results were included in Bosnian referent database. In order of future development of applied population molecular genetics researches of Bosnia and Herzegovina human population, we have examined the effectiveness of 15 STR loci system in determination of sibship by using 15 STR loci and calculating different cut-off points of combined sibship indices (CSI) and distribution of sharing alleles. From the perspective of its application, it is very difficult and complicated to establish strict CSI cut-off values for determination of the doubtless sibship. High statistically significant difference between the means of CSI values and in distribution of alleles sharing in siblings and non-siblings was noticed (P < 0.0001). After constructing the “gray zone”, only one false positive result was found in three CSI cut-off levels with the highest percent of determined sibship/non-sibship at the CSI = 0.067, confirming its practical benefit. Concerning the distribution of sharing alleles, it is recommended as an informative estimator for its usage within Bosnia and Herzegovina human population.     

Genetic characterization of grey wolves (<Emphasis Type="Italic">Canis lupus</Emphasis> L. <Emphasis Type="Italic">1758</Emphasis>) from Bosnia and Herzegovina: implications for conservation

The grey wolves of Bosnia and Herzegovina form a subpopulation of the Dinaric-Balkan wolf population and represent one of Europe’s least studied wolf populations. Since the Dinaric-Balkan population is a valuable source of genetic diversity for neighboring populations, comprehensive assessments are warranted. We aimed to determine the genetic variability and structure of the grey wolf population from Bosnia and Herzegovina, as well as estimate levels of gene flow and inbreeding and evaluate genetic signals of a bottleneck. To do this, we analyzed the variability of eighteen microsatellite loci. We found moderately high genetic heterozygosity for wolves from Bosnia and Herzegovina, as described for other Dinaric-Balkan wolf populations. We reveal weak genetic structuring with two genetic clusters identified. Wolves from the eastern part of the region formed a relatively distinct cluster, whereas individuals in the second cluster overlapped quite considerably with admixed individuals. Despite the signal of genetic structure being weak, clustering of individuals from the eastern part of the country extended through all analyses. Thus, this cluster could be considered a separate management unit, perhaps requiring specific conservation attention.     

Effects of forced migrations on spatial and temporal changes in genetic structure of two local populations in Bosnia

Ten qualitative traits were observed in two isolated human populations in Bosnia for prewar and postwar periods. Due to recent war in Bosnia and Herzegovina two million citizens was forced to migrate. Dramatic migration effects occurred in the investigated populations. Population characteristics found were geographic isolation, high levels of marital distance “zero”, patrilocality and propagation isolation. Relative recessive allele frequencies of ten qualitative traits were calculated. Average heterozigosity, genetic differentiation, total heterozigosity, within populations heterozigosity, genetic distance and relative measures based on results of genetic distance analysis were observed. The total number of observed individuals was 1875. We have observed 1209 individuals in prewar period (Bijela population 731 and Memici population 478) and 666 individuals in the postwar period (Bijela population 248 and Memici population 478). Results ofFst andGst analysis showed changed degree of genetic differentiation across observed loci. Significant changes in differentiation were recorded for crooked little finger and PTC tasting. Modest changes of gene differentiation were recorded for ear lobe type and thumb proximal extensibility. Genetic distance analisis between Bijela and Memici has lower value after then before the war, but genetic distance between Bijela-prewar and Bijela-postwar has higher value, as well as Memici-prewar and Memici-postwar genetic distance. Genetic distance between Bijela-prewar and Memiciprewar, as well as Bijela-postwar and Memici-postwar showed similarity. Constructed dendrogram based on results of genetic distance analyses indicated two clusters groups (1. Bijela-prewar, Memici-prewar; 2. Bijela-postwar, Memici-postwar). Changed genetic differentiation and results of genetic distance analyses indicated possibility of significant impact of forced migration in genetic structure of observed populations.     

Ancient DNA analysis of human neolithic remains found in northeastern Siberia

We successfully extracted DNA from a bone sample of a Neolithic skeleton (dated 3,600 +/- 60 years BP) excavated in northeastern Yakutia (east Siberia). Ancient DNA was analyzed by autosomal STRs (short tandem repeats) and by sequencing of the hypervariable region I (HV1) of the mitochondrial DNA (mtDNA) control region. The STR profile, the mitochondrial haplotype, and the haplogroup determined were compared with those of modern Eurasian and Native American populations. The results showed the affinity of this ancient skeleton with both east Siberian/Asian and Native American populations.     

Genetic variation and island biogreography: Microsatellite and mitochondrial DNA variation in island populations of the Australian bush rat, Rattus fuscipes greyii

To understand the impact of various factors on the maintenance of genetic variation in natural populations, we need to focus on situations where at least some of these factors are removed or controlled. In this study, we used highly variable, presumably neutral, microsatellite and mtDNA markers to assess the nature of genetic variation in 14 island and two mainland populations of the Australian bush rat, where there is no migration between islands. Thus we are controlling for selection and gene flow. Both marker sets revealed low levels of diversity within the small island populations and extreme differentiation between populations. For six microsatellite loci, all of the small island populations had less genetic variation than the mainland populations; reduction in allelic diversity was more pronounced than loss of heterozygosity. Kangaroo Island, the large island population, had similar levels of diversity to the mainland populations. A 442 base pair (bp) section of the mtDNA control region was screened for variation by outgroup heteroduplex analysis/temperature gradient gel electrophoresis (OHA/TGGE). Only three of the 13 small island populations showed haplotypic diversity: Gambier (2), Waldegrave (2), and Eyere (3). The level of haplotypic diversity in the small island populations was similar to that on the mainland, most likely reflecting a recent population bottleneck on the mainland. In contrast, Kangaroo Island had 9 mtDNA haplotypes. The dominant factor influencing genetic diversity on the islands was island size. No correlation was detected between genetic diversity and the time since isolation or distance form the mainland. The combination of genetic drift within and complete isolation among the small island populations has resulted in rapid and extreme population divergence. Population pair-wise comparisons of allele frequency distributions showed significant differences for all populations for all loci (F _st = 0.11–0.84, R _st = 0.07–0.99). For the mtDNA control region, 92.6% of variation was apportioned between populations; only the Pearson islands shared a haplotype. Mantel tests of pair-wise genetic distance with pair-wise geographic distance showed no significant geographical clustering of haplotypes. However, population substructuring was detected within populations where sampling was conducted over a broader geographical range, as indicated by departures from Hardy-Weinberg equilibrium. Thus substructuring in the ancestral population cannot be ruled out. The dominant evolutionary forces on the islands, after the initial founder event, are stochastic population processes such as genetic drift and mutation. This revised version was published online in July 2006 with corrections to the Cover Date.     

Mitochondrial and nuclear genetic relationships among Pacific Island and Asian populations.

Mitochondrial and autosomal short tandem-repeat (STR) genetic distances among 28 Pacific Island and Asian populations are significantly correlated (r=.25, P<.01) but describe distinct patterns of relationships. Maternally inherited-mtDNA data suggest that Remote Oceanic Islanders originated in island Southeast Asia. In contrast, biparental STR data reveal substantial genetic affinities between Remote Oceanic Islanders and Near Oceanic populations from highland Papua New Guinea and Australia. The low correlation between maternal and biparental genetic markers from the same individuals may reflect differences in genome-effective population sizes or in sex-biased gene flow. To explore these possibilities, we have examined genetic diversity, gene flow, and correlations among genetic, linguistic, and geographic distances within four sets of populations representing potential geographic and cultural spheres of interaction. GST estimates (a measure of genetic differentiation inversely proportional to gene flow) from mtDNA sequences vary between 0.13 and 0.39 and are typically five times greater than GST estimates from STR loci (0.05-0.08). Significant correlations (r>.5, P<.05) between maternal genetic and linguistic distances are coincident with high mtDNA GST estimates (>0.38). Thus, genetic and linguistic distances may coevolve, and their correspondence may be preserved under conditions of genetic isolation. A significant correlation (r=.65, P<.01) between biparental genetic and geographic distances is coincident with a low STR GST estimate (0.05), indicating that isolation by distance is observed under conditions of high nuclear-gene flow. These results are consistent with an initial settlement of Remote Oceania from island Southeast Asia and with extensive postcolonization male-biased gene flow with Near Oceania.     

Phylogenetic Relationships Among Fragmented Asian Black Bear (Ursus Thibetanus) Populations in Western Japan

The number of Asian black bears (Ursus thibetanus) in Japan has been reduced and their habitats fragmented and isolated because of human activities. Our previous study examining microsatellite DNA loci revealed significant genetic differentiation among four local populations in the western part of Honshu. Here, an approximate 700-bp nucleotide sequence of mitochondrial DNA (mtDNA) control region was analysed in 119 bears to infer the evolutionary history of these populations. Thirteen variable sites and variation in the number of Ts at a T-repeat site were observed among the analysed sequences, which defined 20 mtDNA haplotypes with the average sequence divergence of 0.0051 (SD = 0.00001). The observed haplotype frequencies differed significantly among the four populations. Phylogeographic analysis of the haplotypes suggested that black bears in this region have gone through two different colonisation histories, since the observed haplotypes belonged to two major monophyletic lineages and the lineages were distributed with an apparent border. The spatial genetic structure revealed by using mtDNA was different from that observed using microsatellite DNA markers, probably due to female philopatry and male-biased dispersal. Since nuclear genetic diversity will be lost in the three western populations because of the small population size and genetic isolation, their habitats need to be preserved, and these four populations should be linked to each other by corridors to promote gene flow from the easternmost population with higher nuclear genetic diversity.     

Mitochondrial sequence variation in the Guahibo Amerindian population from Venezuela

New data were obtained on mitochondrial DNA (mtDNA) from Guahibo from Venezuela, a group so far not studied using molecular data. A population sample (n = 59) was analyzed for mtDNA variation in two control-region hypervariable segments (HV1 and HV2) by sequencing. The presence or absence of a 9-bp polymorphism in the COII/tRNA(Lys) region was studied by direct amplification and electrophoretic identification. Thirty-eight variable sites were detected in regions HV1 and HV2, defining 26 mtDNA lineages; 23.7% of these were present in a single individual. The 9-bp deletion was found in 3.39% of individuals. Nucleotide and haplotype diversities were relatively high compared with other New World populations. The identified sequence haplotypes were classified into four major haplogroups (A-D) according to previous studies, with high frequencies for A (47.46%) and C (49.15%), low frequency for B (3.39%), and an absence of D.     

Genetic variation of European grayling (Thymallus thymallus) populations in the Western Balkans

In order to elucidate genetic composition of European grayling (Thymallus thymallus) populations in the Western Balkans, the partial mitochondrial DNA (mtDNA) control region was sequenced and 12 microsatellite loci genotyped in 14 populations originating from tributaries of the Adriatic and Danube drainages. Eleven mtDNA haplotypes were found, one confined to the Adriatic clade, one to the Alpine group and the rest to the ‘Balkan’ grayling phylogenetic clade. Haplotypes from the Balkan clade were confined to the Danube drainage and constituted two groups: northern group with haplotypes found in the Slovenian part of the Danube drainage, and southern group, consisting from Bosnia–Herzegovina and Montenegro. Substantial genetic distance between northern and southern groups of haplotypes (0.75–1.8%) and well supported divisions within the northern group indicate very structured grayling population within the studied Danube basin that most probably did not evolve due to vicariance but rather as a consequence of multiple colonization waves that might have occurred during the Pleistocene. Furthermore, genetic distance of ~4% between Adriatic and Danube populations’ haplotypes, suggest that their separation occurred in mid-Pliocene. These findings imply a complex colonization pattern of the Western Balkans drainages. Microsatellite data also confirm high genetic diversity in Western Balkans populations of grayling (on average 7.5 alleles per microsatellite locus and H _exp 0.58). Limited stocking activities were detected based on microsatellites and mtDNA data. Regarding current knowledge of grayling phylogeography appropriate management strategies were proposed to preserve unique, autochthonous grayling populations in Western Balkan.     

STR polymorphism of mtDNA D-loop in rhesus macaques of Bangladesh

Molecular variation of mitochondrial DNA (mtDNA) was investigated for rhesus macaques (Macaca mulatta) of Bangladesh. A partial sequence (583–599 bp) of mtDNA containing the second variable region of the D-loop was compared for 39 individuals from five localities in the country. A total of seven haplotypes were detected with substitutional or insertion/deletion mutations. They contained a unique polymorphism of pentanucleotide STRs (short tandem repeats). There were at least four different length types, from two to five repeats of the unit nucleotide. One site of substitution and one site of single nucleotide insertion/deletion were also involved in the polymorphism. The mutation hot spots of the STR polymorphism were located between the first and second conserved sequence blocks (CSB1 and CSB2), as observed previously in some other mammals. The geographical distribution of the STR polymorphism revealed local differences; the northeastern population was polymorphic with three STR haplotypes, but other local populations were simply monomorphic with a single STR haplotype. Molecular phylogenetic analysis with reported sequences from outside Bangladesh indicated a low substitution diversity of mtDNA in Bangladesh. Clustering results suggested a close relationship to India and divergence from Laos and China.     

Overview of human population-genetic studies in Bosnia and Herzegovina during the last three centuries: history and prospective

Modern Bosnia and Herzegovina is a multi-ethnic and multi-religion country, with a very stormy history. Certain archaeological findings indicate continuous population of its territory since the Paleolithic. In time, vast number of different factors jointly influenced fascinating diversity of local human populations. A great number of small, more or less isolated, indigenous populations, make this area quite attractive for population-genetic surveys of different levels and approaches. Austro-Hungarian military physicians conducted the very first known bio-anthropological analyses of Bosnia-Herzegovina population at the end of the 19th century. Thus, the first step towards resolving the genetic structures of local B&H human populations was made. The studies that followed (conducted throughout most of the 20th century) were primarily based on the observation of various phenotypic traits. This stage was followed by the examination of various cytogenetic and fundamental DNA based molecular markers. The efforts undertaken over the last three centuries revealed "human genetic treasure" in Bosnia and Herzegovina. However, even now, after all the studies that were conducted, many interesting features remain to be discovered and described within the existing local human populations.     

Population data for 17 short tandem repeat loci on Y chromosome in northern Croatia

Human Y-short tandem repeats (STRs) are tandem repeat arrays of two to seven base pair units on non-recombining region (NRY) of the human Y chromosome. Studies on Y-STR are interesting in both population genetics and forensics. The aim of this study was to investigate the population genetic properties of 17 STR loci on Y chromosome in the northern Croatia region. We carried out a statistical analysis of the data from previously performed genetic analysis collected during routine forensic work by the Forensic Science Centre “Ivan Vu?eti?”. A total of 220 unrelated healthy men from northern Croatia were selected for the purpose of this study. Genomic DNA was extracted using Chelex procedure from FTA^® cards. Y-chromosomal STRs were determined using the AmpFISTR Yfiler PCR amplification kit. The haplotype frequencies were determined by direct counting and analyzed using Arlequin 3.1 and analysis of molecular variance calculated with the Y chromosome haplotype reference database online analysis tool. A total of 210 haplotypes were identified, 200 of which were unique. Total haplotype diversity was 0.995. Locus diversity varied from 0.331 for DYS392 to 0.783 for DYS385 locus. Allele frequencies diversity was 0.662. Discrimination capacity was 95.7%. The use of European minimal haplotype set indicated the most resemblance of this population to the Croatian capital of Zagreb, with modest resemblance to Bosnia and Herzegovina, Serbia and Hungary. This article provides the first overview of the Y chromosome STR variability in northern Croatia, thus providing the referent point for any future forensic and genetic epidemiology efforts in this region.     

Microsatellite analysis of the natterjack toad (Bufo calamita) in Denmark: populations are islands in a fragmented landscape

The European natterjack toad (Bufo calamita) has declined rapidly in recent years, primarily due to loss of habitat, and in Denmark it is estimated that 50% of the isolated populations are lost each decade. To efficiently manage and conserve this species and its genetic diversity, knowledge of the genetic structure is crucial. Based on nine polymorphic microsatellite loci, the genetic diversity, genetic structure and gene flow were investigated at 12 sites representing 5–10% of the natterjack toad localities presently known in Denmark. The expected heterozygosity (H _E) within each locality was generally low (range: 0.18–0.43). Further analyses failed to significantly correlate genetic diversity with population size, degree of isolation and increasing northern latitude, indicating a more complex combination of factors in determining the present genetic profile. Genetic differentiation was high (overall θ = 0.29) and analyses based on a Bayesian clustering method revealed that the dataset constituted 11 genetic clusters, defining nearly all sampling sites as distinct populations. Contemporary gene flow among populations was undetectable in nearly all cases, and the failure to detect a pattern of isolation by distance within major regions supported this apparent lack of a gene flow continuum. Indications of a genetic bottleneck were found in three populations. The analyses suggest that the remaining Bufo calamita populations in Denmark are genetically isolated, and represent independent units in a highly fragmented gene pool. Future conservation management of this species is discussed in light of these results. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Conserving the endangered Mexican fishing bat (<Emphasis Type="Italic">Myotis vivesi</Emphasis>): genetic variation indicates extensive gene flow among islands in the Gulf of California

The endangered Mexican fishing bat, Myotis vivesi, appears to have suffered widespread extinction and population decline on islands throughout the Gulf of California, largely due to predation by introduced cats and rats. To restore populations of fishing bats and other native species, conservation efforts have focused on eradicating introduced vertebrates from several Gulf islands. These efforts assume that individuals from existing populations will recolonize islands and that continued dispersal will help sustain vulnerable populations thereafter. However, the extent of inter-island dispersal in fishing bats is unknown. In this study we analyzed patterns of genetic variation to gauge the extent of gene flow and, thus, potential dispersal among islands. DNA was sampled from 257 fishing bats on 11 Gulf islands (separated by ca. 6–685 km of open water), and individuals were genotyped at six microsatellite loci and haplotyped at a 282 bp fragment of the mtDNA control region. With microsatellites, we found weak population genetic structure and a pattern of isolation by distance, while with mtDNA we found strong structure but no isolation by distance. Our results indicate that island subpopulations separated by large expanses of open water are nonetheless capable of maintaining high genetic diversity and high rates of gene flow. Unfortunately, little is known about the spatial patterns of dispersal or mating system of fishing bats, and these behavioral factors, in particular female philopatry, might reduce the probability of the species recolonizing Gulf islands.     

Three reciprocally monophyletic mtDNA lineages elucidate the taxonomic status of Grant’s gazelles

The intraspecific phylogeography of Grant’s gazelles Nanger granti was assessed with mitochondrial DNA control region sequences. Samples of 177 individuals from 17 Kenyan and Tanzanian populations were analysed. Three highly divergent, reciprocally monophyletic lineages were found, with among group net nucleotide distances of 8–12%. The three lineages—notata, granti and petersii—grouped populations according to their geographic origin, encompassing populations in the north, southwest, and east, respectively. The mtDNA lineages reflected distinct evolutionary trajectories, and the data are discussed in reference to the four currently recognised subspecies. We suggest Grant’s gazelles be raised to the superspecies Nanger (granti) comprising three taxonomic units corresponding to the three mtDNA lineages. There was no evidence of gene flow between the notata and granti lineages, despite their geographic proximity, suggesting reproductive isolation. These constitute evolutionary significant units within the adaptive evolutionary framework. Due to its restricted geographic distribution and genetic and morphological distinctiveness, we suggest the petersii lineage be raised to the species Nanger (granti) petersii within the Grant’s gazelles superspecies.     

Identification of conservation units in the European Mergus merganser based on nuclear and mitochondrial DNA markers

The conservation status of small breeding areas of the Goosander (Mergus merganser merganser) in Central Europe is unclear. Geographic isolation of these areas suggests restricted gene flow to and from large North-European populations. On the other hand, migrating Goosanders from northern Europe join the Central European breeding population for wintering. To evaluate the conservation status of the small breeding areas we assessed the genetic structure of M. merganser populations in Europe by examining two nuclear marker systems (microsatellites and Single Nucleotide Polymorphisms, SNP) and mitochondrial (mtDNA) control region sequence variation for Goosanders in 11 sampling areas representing three of five distinct breeding areas and two subspecies (M. m. merganser and M. m. americanus). Overall population differentiation estimates including both subspecies were high, both based on mtDNA () and nuclear markers (θ _ST = 0.219; 95% CI 0.088–0.398, SNP and microsatellites combined). Within Europe, mtDNA revealed a strong overall () and significant pairwise population differentiation between almost all comparisons. In contrast, both nuclear marker systems combined revealed only a small overall genetic differentiation (θ _ST = 0.022; 95% CI 0.003–0.041). The strong genetic differentiation based on female-inherited mtDNA but not on biparentally inherited nuclear markers can be explained by sex-biased dispersal and strong female philopatry. Therefore, small breeding areas in Europe are endangered despite large male-mediated gene-flow, because when these populations decline, only males—but due to strong philopatry not females—can be efficiently supplemented by migration from the large North European populations. We therefore propose to manage the small breeding areas independently and to strengthen conservation efforts for this species in Central Europe.     

Genetic diversity and structuring of the grey wolf population from the Central Balkans based on mitochondrial DNA variation

The Dinaric-Balkan grey wolf population used to be at a border between the large remaining Eastern European populations and the largely eradicated Western European populations. During the last few decades we have witnessed the Western European wolf population recovery. Substantial genetic variation has previously been reported in the Balkan wolf population, but rigorous genetic characterization has not been done for its central parts. The aims of this research were to determine genetic diversity based on mtDNA sequence variability, to infer possible population structuring, to find genetic signals of population expansions or bottlenecks and to evaluate phylogenetic position of the grey wolf population from the Central Balkans. Six haplotypes were detected, of which three have only been found in the Balkan region. These haplotypes belong to both haplogroups previously determined in Europe. Based on our mtDNA sequence analyses, the Dinaric-Balkan wolf population is vertically differentiated into “western” (Croatia/Bosnia and Herzegovina) and “eastern” (Serbia/Macedonia) subpopulations. None of the results support assumption of population expansion. Instead, significantly positive values for Tajima's D and Fu's Fs may suggest recent population bottleneck. Obtained data may be helpful in observation to which extent gene pool from the Balkans contribute to newly founded populations in Western Europe.     

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.