Global population genetic dynamics of a highly migratory,apex predator shark

Knowledge of genetic connectivity dynamics in the world's large‐bodied, highly migratory, apex predator sharks across their global ranges is limited. One such species, the tiger shark (Galeocerdo cuvier), occurs worldwide in warm temperate and tropical waters, uses remarkably diverse habitats (nearshore to pelagic) and possesses a generalist diet that can structure marine ecosystems through top‐down processes. We investigated the phylogeography and the global population structure of this exploited, phylogenetically enigmatic shark by using 10 nuclear microsatellites (n = 380) and sequences from the mitochondrial control region (CR, n = 340) and cytochrome oxidase I gene (n = 100). All three marker classes showed the genetic differentiation between tiger sharks from the western Atlantic and Indo‐Pacific ocean basins (microsatellite F_ST > 0.129; CR Φ_ST > 0.497), the presence of North vs. southwestern Atlantic differentiation and the isolation of tiger sharks sampled from Hawaii from other surveyed locations. Furthermore, mitochondrial DNA revealed high levels of intraocean basin matrilineal population structure, suggesting female philopatry and sex‐biased gene flow. Coalescent‐ and genetic distance‐based estimates of divergence from CR sequences were largely congruent (d_corr = 0.0015–0.0050), indicating a separation of Indo‐Pacific and western Atlantic tiger sharks <1 million years ago. Mitochondrial haplotype relationships suggested that the western South Atlantic Ocean was likely a historical connection for interocean basin linkages via the dispersal around South Africa. Together, the results reveal unexpectedly high levels of population structure in a highly migratory, behaviourally generalist, cosmopolitan ocean predator, calling for management and conservation on smaller‐than‐anticipated spatial scales.     

Population structure,connectivity, and demographic history of an apex marine predator,the bull shark Carcharhinus leucas

Knowledge of population structure, connectivity, and effective population size remains limited for many marine apex predators, including the bull shark Carcharhinus leucas. This large‐bodied coastal shark is distributed worldwide in warm temperate and tropical waters, and uses estuaries and rivers as nurseries. As an apex predator, the bull shark likely plays a vital ecological role within marine food webs, but is at risk due to inshore habitat degradation and various fishing pressures. We investigated the bull shark's global population structure and demographic history by analyzing the genetic diversity of 370 individuals from 11 different locations using 25 microsatellite loci and three mitochondrial genes (CR, nd4, and cytb). Both types of markers revealed clustering between sharks from the Western Atlantic and those from the Western Pacific and the Western Indian Ocean, with no contemporary gene flow. Microsatellite data suggested low differentiation between the Western Indian Ocean and the Western Pacific, but substantial differentiation was found using mitochondrial DNA. Integrating information from both types of markers and using Bayesian computation with a random forest procedure (ABC‐RF), this discordance was found to be due to a complete lack of contemporary gene flow. High genetic connectivity was found both within the Western Indian Ocean and within the Western Pacific. In conclusion, these results suggest important structuring of bull shark populations globally with important gene flow occurring along coastlines, highlighting the need for management and conservation plans on regional scales rather than oceanic basin scale.     

COMPUTER PROGRAMS: onesamp: a program to estimate effective population size using approximate Bayesian computation

The estimation of effective population size from one sample of genotypes has been problematic because most estimators have been proven imprecise or biased. We developed a web-based program, onesamp that uses approximate Bayesian computation to estimate effective population size from a sample of microsatellite genotypes. onesamp requires an input file of sampled individuals' microsatellite genotypes along with information about several sampling and biological parameters. onesamp provides an estimate of effective population size, along with 95% credible limits. We illustrate the use of onesamp with an example data set from a re-introduced population of ibex Capra ibex.     

Ontogeny of head and caudal fin shape of an apex marine predator: The tiger shark (Galeocerdo cuvier)

How morphology changes with size can have profound effects on the life history and ecology of an animal. For apex predators that can impact higher level ecosystem processes, such changes may have consequences for other species. Tiger sharks (Galeocerdo cuvier) are an apex predator in tropical seas, and, as adults, are highly migratory. However, little is known about ontogenetic changes in their body form, especially in relation to two aspects of shape that influence locomotion (caudal fin) and feeding (head shape). We captured digital images of the heads and caudal fins of live tiger sharks from Southern Florida and the Bahamas ranging in body size (hence age), and quantified shape of each using elliptical Fourier analysis. This revealed changes in the shape of the head and caudal fin of tiger sharks across ontogeny. Smaller juvenile tiger sharks show an asymmetrical tail with the dorsal (upper) lobe being substantially larger than the ventral (lower) lobe, and transition to more symmetrical tail in larger adults, although the upper lobe remains relatively larger in adults. The heads of juvenile tiger sharks are more conical, which transition to relatively broader heads over ontogeny. We interpret these changes as a result of two ecological transitions. First, adult tiger sharks can undertake extensive migrations and a more symmetrical tail could be more efficient for swimming longer distances, although we did not test this possibility. Second, adult tiger sharks expand their diet to consume larger and more diverse prey with age (turtles, mammals, and elasmobranchs), which requires substantially greater bite area and force to process. In contrast, juvenile tiger sharks consume smaller prey, such as fishes, crustaceans, and invertebrates. Our data reveal significant morphological shifts in an apex predator, which could have effects for other species that tiger sharks consume and interact with. J. Morphol. 277:556–564, 2016. © 2016 Wiley Periodicals, Inc.     

Genetic effective size of a wild primate population: influence of current and historical demography

A comprehensive assessment of the determinants of effective population size (N(e)) requires estimates of variance in lifetime reproductive success and past changes in census numbers. For natural populations, such information can be best obtained by combining longitudinal data on individual life histories and genetic marker-based inferences of demographic history. Independent estimates of the variance effective size (N(ev), obtained from life-history data) and the inbreeding effective size (N((eI), obtained from genetic data) provide a means of disentangling the effects of current and historical demography. The purpose of this study was to assess the demographic determinants of N(e) in one of the most intensively studied natural populations of a vertebrate species: the population of savannah baboons (Papio cynocephalus) in the Amboseli Basin, southern Kenya. We tested the hypotheses that N(eV) < N < N(eI) (where N = population census number) due to a recent demographic bottleneck. N(eV) was estimated using a stochastic demographic model based on detailed life-history data spanning a 28-year period. Using empirical estimates of age-specific rates of survival and fertility for both sexes, individual-based simulations were used to estimate the variance in lifetime reproductive success. The resultant values translated into an N(eV)/N estimate of 0.329 (SD = 0.116, 95% CI = 0.172-0.537). Historical N(eI), was estimated from 14-locus microsatellite genotypes using a coalescent-based simulation model. Estimates of N(eI) were 2.2 to 7.2 times higher than the contemporary census number of the Amboseli baboon population. In addition to the effects of immigration, the disparity between historical N(eI) and contemporary N is likely attributable to the time lag between the recent drop in census numbers and the rate of increase in the average probability of allelic identity-by-descent. Thus, observed levels of genetic diversity may primarily reflect the population's prebottleneck history rather than its current demography.     

A novel mode of embryonic nutrition in the tiger shark,Galeocerdo cuvier

How are tiger shark embryos nourished to large size without a placental connection? Tiger sharks belong to the family Carcharhinidae, and all carcharhinid sharks are placental with the exception of the tiger shark. The aim of this study was to test the hypothesis that tiger shark embryos are nourished to large size by imbibing a clear uterine fluid found in their egg cases. Based on weights of fertilized eggs and of term embryos, the tiger shark is a matrotrophic species, and its embryos appear to reach gains of 2119% in wet weight and 1092% in dry weight during gestation. By measuring the total energy content of the fluid in the egg case by chemical oxygen demand (COD), the authors demonstrate that clear liquid in the tiger shark egg case is an energy-rich embryotrophe that nourishes the embryos to large size. We suggest that the process be termed ‘embryotrophy'. The process appears to be an adaptation for producing large broods of large embryos in a species lacking a placental connection.     

Extracting DNA from ‘jaws’: high yield and quality from archived tiger shark (Galeocerdo cuvier) skeletal material

Archived specimens are highly valuable sources of DNA for retrospective genetic/genomic analysis. However, often limited effort has been made to evaluate and optimize extraction methods, which may be crucial for downstream applications. Here, we assessed and optimized the usefulness of abundant archived skeletal material from sharks as a source of DNA for temporal genomic studies. Six different methods for DNA extraction, encompassing two different commercial kits and three different protocols, were applied to material, so‐called bio‐swarf, from contemporary and archived jaws and vertebrae of tiger sharks (Galeocerdo cuvier). Protocols were compared for DNA yield and quality using a qPCR approach. For jaw swarf, all methods provided relatively high DNA yield and quality, while large differences in yield between protocols were observed for vertebrae. Similar results were obtained from samples of white shark (Carcharodon carcharias). Application of the optimized methods to 38 museum and private angler trophy specimens dating back to 1912 yielded sufficient DNA for downstream genomic analysis for 68% of the samples. No clear relationships between age of samples, DNA quality and quantity were observed, likely reflecting different preparation and storage methods for the trophies. Trial sequencing of DNA capture genomic libraries using 20 000 baits revealed that a significant proportion of captured sequences were derived from tiger sharks. This study demonstrates that archived shark jaws and vertebrae are potential high‐yield sources of DNA for genomic‐scale analysis. It also highlights that even for similar tissue types, a careful evaluation of extraction protocols can vastly improve DNA yield.     

Teeth penetration force of the tiger shark Galeocerdo cuvier and sandbar shark Carcharhinus plumbeus

This study examined the minimum force required of functional teeth and replacement teeth in the tiger shark Galeocerdo cuvier and the sandbar shark Carcharhinus plumbeus to penetrate the scales and muscle of sheepshead Archosargus probatocephalus and pigfish Orthopristis chrysoptera. Penetration force ranged from 7·7–41·9 and 3·2–26·3 N to penetrate A. probatocephalus and O. chrysoptera, respectively. Replacement teeth required significantly less force to penetrate O. chrysoptera for both shark species, most probably due to microscopic wear of the tooth surfaces supporting the theory shark teeth are replaced regularly to ensure sharp teeth that are efficient for prey capture.     

Genetic evidence does not support an etruscan origin in Anatolia

The debate on the origins of Etruscans, documented in central Italy between the eighth century BC and the first century AD, dates back to antiquity. Herodotus described them as a group of immigrants from Lydia, in Western Anatolia, whereas for Dionysius of Halicarnassus they were an indigenous population. Dionysius' view is shared by most modern archeologists, but the observation of similarities between the (modern) mitochondrial DNAs (mtDNAs) of Turks and Tuscans was interpreted as supporting an Anatolian origin of the Etruscans. However, ancient DNA evidence shows that only some isolates, and not the bulk of the modern Tuscan population, are genetically related to the Etruscans. In this study, we tested alternative models of Etruscan origins by Approximate Bayesian Computation methods, comparing levels of genetic diversity in the mtDNAs of modern and ancient populations with those obtained by millions of computer simulations. The results show that the observed genetic similarities between modern Tuscans and Anatolians cannot be attributed to an immigration wave from the East leading to the onset of the Etruscan culture in Italy. Genetic links between Tuscany and Anatolia do exist, but date back to a remote stage of prehistory, possibly but not necessarily to the spread of farmers during the Neolithic period. Am J Phys Anthropol 152:11–18, 2013. © 2013 Wiley Periodicals, Inc.     

Age and growth of the tiger shark Galeocerdo cuvier off the east coast of Australia

Total lengths (L_T) at age and growth rates for south‐west Pacific Galeocerdo cuvier were estimated from vertebral growth‐band counts of 202 sagitally sectioned centra from 112 females (71–430 cm L_T), 79 males (72–351 cm L_T) and 11 of unknown sex. Captive growth data were also examined to complement vertebral age estimations. The sexes combined modelled growth coefficient (k = 0·08) was smaller than previously reported for G. cuvier populations elsewhere. Split‐band and narrow banding patterns were identified as potential sources of age underestimation in this species.     

Reconstruction of a beech population bottleneck using archival demographic information and Bayesian analysis of genetic data

Range expansion and contraction has occurred in the history of most species and can seriously impact patterns of genetic diversity. Historical data about range change are rare and generally appropriate for studies at large scales, whereas the individual pollen and seed dispersal events that form the basis of geneflow and colonization generally occur at a local scale. In this study, we investigated range change in Fagus sylvatica on Mont Ventoux, France, using historical data from 1838 to the present and approximate Bayesian computation (ABC) analyses of genetic data. From the historical data, we identified a population minimum in 1845 and located remnant populations at least 200 years old. The ABC analysis selected a demographic scenario with three populations, corresponding to two remnant populations and one area of recent expansion. It also identified expansion from a smaller ancestral population but did not find that this expansion followed a population bottleneck, as suggested by the historical data. Despite a strong support to the selected scenario for our data set, the ABC approach showed a low power to discriminate among scenarios on average and a low ability to accurately estimate effective population sizes and divergence dates, probably due to the temporal scale of the study. This study provides an unusual opportunity to test ABC analysis in a system with a well-documented demographic history and identify discrepancies between the results of historical, classical population genetic and ABC analyses. The results also provide valuable insights into genetic processes at work at a fine spatial and temporal scale in range change and colonization.     

Fine-scale population structure in a desert amphibian: landscape genetics of the black toad (Bufo exsul)

Environmental variables can strongly influence a variety of intra- and inter-population processes, including demography, population structure and gene flow. When environmental conditions are particularly harsh for a certain species, investigating these effects is important to understanding how populations persist under difficult conditions. Furthermore, species inhabiting challenging environments present excellent opportunities to examine the effects of complex landscapes on population processes because these effects will often be more pronounced. In this study, I use 16 microsatellite loci to examine population structure, gene flow and demographic history in the black toad, Bufo exsul , which has one of the most restricted natural ranges of any amphibian. Bufo exsul inhabits four springs in the Deep Springs Valley high desert basin and has never been observed more than several meters from any source of water. My results reveal limited gene flow and moderately high levels of population structure ( F _ST = 0.051–0.063) between all but the two closest springs. I found that the geographic distance across the arid scrub habitat between springs is significantly correlated with genetic structure when distance accounts for topography and barriers to dispersal. I also found very low effective population sizes ( N _e = 7–30) and substantial evidence for historical population bottlenecks in all four populations. Together, these results suggest that the desert landscape and B.   exsul 's high habitat specificity contribute significantly to population structure and demography in this species and emphasize the importance of considering behavioural and landscape data in conservation genetic studies of natural systems.     

Ontogenetic dietary shifts and feeding behavior of the tiger shark, Galeocerdo cuvier, in Hawaiian waters

Synopsis Stomach content data from 281 tiger sharks caught during shark control programs in Hawaii between 1967 and 1969, and during 1976 were analyzed to examine feeding habits and ontogenetic shifts in diet. As sharks increased in size, prey diversity and frequency of occurrence of large prey items increased. The percent occurrence of teleosts and cephalopods in stomachs decreased as sharks increased in length, while occurrence of elasmobranchs, turtles, land mammals, crustaceans, and undigestible items increased. Comparisons between the diets of tiger sharks from Hawaii and other locations indicate that ontogenetic shifts are universal in this species and that tiger sharks may be opportunistic feeders that prey heavily on abundant, easy to capture prey. Small tiger sharks may be spatially segregated from medium and large sharks and appear to be primarily nocturnal, bottom feeders. Large tiger sharks feed near the bottom at night, but also feed at the surface during the day. Prey, similar in size to humans, begin to occur in the diet of tiger sharks approximately 230 cm TL, and therefore sharks of this size and larger may pose the greatest threat to humans. Ontogenetic shifts in diet may be attributed to increased size of sharks, expanded range and exploitation of habitats of larger sharks, and/or improved hunting skill of larger sharks.Deceased 1974     

Bighorn Sheep Genetic Structure in Wyoming Reflects Geography and Management

Aligning wildlife management boundaries with accurate biological units promotes effective conservation and management practices that reflect ecological and evolutionary processes. Neutral genetic markers allow for quantitative delineation of population structure without a priori assumptions or biases. In the United States, bighorn sheep (Ovis canadensis) are a charismatic component of Wyoming's biodiversity and a species that provides important viewing and hunting opportunities. Bighorn sheep abundances are relatively stable throughout Wyoming, and the species is managed by administrative units identified using expert knowledge, distribution and movement data, and geographic and administrative boundaries. We used a panel of 38 variable microsatellite loci and 512 base pairs of mitochondrial DNA sequence to identify the genetic structure throughout the state and in translocation source herds, quantify the extent of genetic diversity within each genetic cluster, and estimate the degree of gene flow among herds using blood and tissue samples collected 1989–2017. We identified genetic structure of Rocky Mountain bighorn sheep in the major mountain ranges of Wyoming, with strong support for ≥5 genetic clusters using microsatellite loci. These genetic clusters generally aligned with current management units, whereas mitochondrial data showed a more complex mosaic that was not geographically patterned. Genetic variation estimated from both markers was high within each herd and comparable among herds. The assignment of individuals reflected a combination of geographic isolation and translocation, which has been extensive. Our results provide a state-wide assessment of genetic diversity and structure that will enhance management by understanding the outcomes of translocation, identifying the source of unknown individuals, and parameterizing disease ecology models. © 2020 The Wildlife Society.     

The genetic structure of Nautilus pompilius populations surrounding Australia and the Philippines

Understanding the distribution of genetic diversity in exploited species is fundamental to successful conservation. Genetic structure and the degree of gene flow among populations must be assessed to design appropriate strategies to prevent the loss of distinct populations. The cephalopod Nautilus pompilius is fished unsustainably in the Philippines for the ornamental shell trade and has limited legislative protection, despite the species' recent dramatic decline in the region. Here, we use 14 microsatellite markers to evaluate the population structure of N. pompilius around Australia and the Philippines. Despite their relative geographical proximity, Great Barrier Reef individuals are genetically isolated from Osprey Reef and Shark Reef in the Coral Sea (F_ST = 0.312, 0.229, respectively). Conversely, despite the larger geographical distances between the Philippines and west Australian reefs, samples display a small degree of genetic structure (F_ST = 0.015). Demographic scenarios modelled using approximate Bayesian computation analysis indicate that this limited divergence is not due to contemporary gene flow between the Philippines and west Australia. Instead, present‐day genetic similarity can be explained by very limited genetic drift that has occurred due to large average effective population sizes that persisted at both locations following their separation. The lack of connectivity among populations suggests that immigrants from west Australia would not facilitate natural recolonization if Philippine populations were fished to extinction. These data help to rectify the paucity of information on the species' biology currently inhibiting their conservation classification. Understanding population structure can allow us to facilitate sustainable harvesting, thereby preserving the diversity of genetically distinct stocks.     

In situ population structure and ex situ representation of the endangered Amur tiger

The Amur tiger ( Panthera tigris altaica ) is a critically endangered felid that suffered a severe demographic contraction in the 1940s. In this study, we sampled 95 individuals collected throughout their native range to investigate questions relative to population genetic structure and demographic history. Additionally, we sampled targeted individuals from the North American ex situ population to assess the genetic representation found in captivity. Population genetic and Bayesian structure analyses clearly identified two populations separated by a development corridor in Russia. Despite their well-documented 20th century decline, we failed to find evidence of a recent population bottleneck, although genetic signatures of a historical contraction were detected. This disparity in signal may be due to several reasons, including historical paucity in population genetic variation associated with postglacial colonization and potential gene flow from a now extirpated Chinese population. Despite conflicting signatures of a bottleneck, our estimates of effective population size ( N _e = 27–35) and N _e /N ratio (0.07–0.054) were substantially lower than the only other values reported for a wild tiger population. Lastly, the extent and distribution of genetic variation in captive and wild populations were similar, yet gene variants persisted ex situ that were lost in situ . Overall, our results indicate the need to secure ecological connectivity between the two Russian populations to minimize loss of genetic diversity and overall susceptibility to stochastic events, and support a previous study suggesting that the captive population may be a reservoir of gene variants lost in situ .     

The genetic effective and adult census size of an Australian population of tiger prawns (Penaeus esculentus)

This study compares estimates of the census size of the spawning population with genetic estimates of effective current and long-term population size for an abundant and commercially important marine invertebrate, the brown tiger prawn (Penaeus esculentus). Our aim was to focus on the relationship between genetic effective and census size that may provide a source of information for viability analyses of naturally occurring populations. Samples were taken in 2001, 2002 and 2003 from a population on the east coast of Australia and temporal allelic variation was measured at eight polymorphic microsatellite loci. Moments-based and maximum-likelihood estimates of current genetic effective population size ranged from 797 to 1304. The mean long-term genetic effective population size was 9968. Although small for a large population, the effective population size estimates were above the threshold where genetic diversity is lost at neutral alleles through drift or inbreeding. Simulation studies correctly predicted that under these experimental conditions the genetic estimates would have non-infinite upper confidence limits and revealed they might be overestimates of the true size. We also show that estimates of mortality and variance in family size may be derived from data on average fecundity, current genetic effective and census spawning population size, assuming effective population size is equivalent to the number of breeders. This work confirms that it is feasible to obtain accurate estimates of current genetic effective population size for abundant Type III species using existing genetic marker technology.     

Genetic structure of the Danish red deer (Cervus elaphus)

The red deer (Cervus elaphus) population in Denmark became almost extinct in recent historical times due to over‐hunting. The species has subsequently recovered within remote areas, but non‐Danish individuals have been introduced at several localities. To assess genetic structure, past demographic history, and the possibility of a still existing original stock, we analysed 349 specimens from 11 geographically separate areas and from three enclosed areas, genotyping 11 microsatellite loci. Moreover, an 826‐bp fragment of the control region of the mitochondrial DNA was sequenced for 116 recent specimens and seven museum specimens. There was a significant difference in mean expected heterozygosity (H_E) between the three enclosed areas and the 11 unenclosed areas. Significant departures from Hardy–Weinberg equilibrium were observed in the three enclosed areas and in nine of the unenclosed areas. The overall degree of genetic differentiation among all 14 areas was significant (F_ST = 0.09, P < 0.01), primarily because the mean pairwise F_ST for the three enclosed areas was significantly higher than that for the 11 unenclosed areas. A Bayesian clustering procedure detected three genetically distinct populations and indicated reduced gene flow between the enclosed and unenclosed areas. The individuals in the unenclosed areas show genotypic mixture, presumably as a result of gene flow among them. Markov Chain Monte Carlo simulations, based on the genealogical history of the microsatellite alleles, suggest a drastic decline in the effective population size of the enclosed areas some 188–474 years ago. Mitochondrial DNA analysis of the recent specimens showed seven haplotypes. Individuals from the enclosed Jægersborg Dyrehave contain haplotypes that occur all over Denmark and also are found in Western Europe. A close relationship between Scandinavian and Western European red deer is most likely. Only individuals from the unenclosed Lindenborg Estate and the enclosed Tofte Skov did not group with any other Danish individuals. As six of seven museum specimens had haplotypes also found in modern Danish samples, the current population of red deer in Denmark is genetically close to the original Danish red deer. © 2008 The Linnean Society of London, Biological Journal of the Linnean Society, 2008, 95 , 688–701.     

IUCN classification zones concord with, but underestimate, the population genetic structure of the zebra shark Stegostoma fasciatum in the Indo-West Pacific

The Indo-West Pacific (IWP), from South Africa in the western Indian Ocean to the western Pacific Ocean, contains some of the most biologically diverse marine habitats on earth, including the greatest biodiversity of chondrichthyan fishes. The region encompasses various densities of human habitation leading to contrasts in the levels of exploitation experienced by chondrichthyans, which are targeted for local consumption and export. The demersal chondrichthyan, the zebra shark, Stegostoma fasciatum , is endemic to the IWP and has two current regional International Union for the Conservation of Nature (IUCN) Red List classifications that reflect differing levels of exploitation: 'Least Concern' and 'Vulnerable'. In this study, we employed mitochondrial ND4 sequence data and 13 microsatellite loci to investigate the population genetic structure of 180 zebra sharks from 13 locations throughout the IWP to test the concordance of IUCN zones with demographic units that have conservation value. Mitochondrial and microsatellite data sets from samples collected throughout northern Australia and Southeast Asia concord with the regional IUCN classifications. However, we found evidence of genetic subdivision within these regions, including subdivision between locations connected by habitat suitable for migration. Furthermore, parametric F _ST analyses and Bayesian clustering analyses indicated that the primary genetic break within the IWP is not represented by the IUCN classifications but rather is congruent with the Indonesian throughflow current. Our findings indicate that recruitment to areas of high exploitation from nearby healthy populations in zebra sharks is likely to be minimal, and that severe localized depletions are predicted to occur in zebra shark populations throughout the IWP region.