Philopatry influences the genetic population structure of the blacktip shark (Carcharhinus limbatus) at multiple spatial scales.

Understanding how interactions among microevolutionary forces generate genetic population structure of exploited species is vital to the implementation of management policies that facilitate persistence. Philopatry displayed by many coastal shark species can impact gene flow and facilitate selection, and has direct implications for the spatial scales of management. Here, genetic structure of the blacktip shark (Carcharhinus limbatus) was examined using a mixed-marker approach employing mitochondrial control region sequences and 4339 SNP-containing loci generated using ddRAD-Seq. Genetic variation was assessed among young-of-the-year sampled in 11 sites in waters of the United States in the western North Atlantic Ocean, including the Gulf of Mexico. Spatial and environmental analyses detected 68 nuclear loci putatively under selection, enabling separate assessments of neutral and adaptive genetic structure. Both mitochondrial and neutral SNP data indicated three genetically distinct units-the Atlantic, eastern Gulf, and western Gulf-that align with regional stocks and suggest regional philopatry by males and females. Heterogeneity at loci putatively under selection, associated with temperature and salinity, was observed among sites within Gulf units, suggesting local adaptation. Furthermore, five pairs of siblings were identified in the same site across timescales corresponding with female reproductive cycles. This indicates that females re-used a site for parturition, which has the potential to facilitate the sorting of adaptive variation among neighbouring sites. The results demonstrate differential impacts of microevolutionary forces at varying spatial scales and highlight the importance of conserving essential habitats to maintain sources of adaptive variation that may buffer species against environmental change.

Interrelationships and biogeography of the New World pufferfish genus Sphoeroides (Tetraodontiformes: Tetraodontidae) inferred using ultra-conserved DNA elements.

Colonization of the New World by marine taxa has been hypothesized to have occurred through the Tethys Sea or by crossing the East Pacific Barrier. To better understand patterns and timing of diversification, geological events can be coupled with time calibrated phylogenetic hypotheses to infer major drivers of diversification. Phylogenetic relationships among members of Sphoeroides, a genus of four toothed pufferfishes (Tetraodontiformes: Tetraodontidae) which are found nearly exclusively in the New World (eastern Pacific and western Atlantic), were reconstructed using sequences from ultra-conserved DNA elements, nuclear markers with clear homology among many vertebrate taxa. Hypotheses derived from concatenated maximum-likelihood and species tree summary methods support a paraphyletic Sphoeroides, with Colomesus deeply nested within the genus. Analyses also revealed S. pachygaster, a pelagic species with a cosmopolitan distribution, as the sister taxon to the remainder of Sphoeroides and recovered distinct lineages within S. pachygaster, indicating that this cosmopolitan species may represent a species complex. Ancestral range reconstruction may suggest the genus colonized the New World through the eastern Pacific before diversifying in the western Atlantic, though date estimates for these events are uncertain due to the lack of reliable fossil record for the genus.

Temperature-associated selection linked to putative chromosomal inversions in king scallop (Pecten maximus).

The genomic landscape of divergence-the distribution of differences among populations or species across the genome-is increasingly characterized to understand the role that microevolutionary forces such as natural selection and recombination play in causing and maintaining genetic divergence. This line of inquiry has also revealed chromosome structure variation to be an important factor shaping the landscape of adaptive genetic variation. Owing to a high prevalence of chromosome structure variation and the strong pressure for local adaptation necessitated by their sessile nature, bivalve molluscs are an ideal taxon for exploring the relationship between chromosome structure variation and local adaptation. Here, we report a population genomic survey of king scallop (Pecten maximus) across its natural range in the northeastern Atlantic Ocean, using a recent chromosome-level genome assembly. We report the presence of at least three large (12-22 Mb), putative chromosomal inversions associated with sea surface temperature and whose frequencies are in contrast to neutral population structure. These results highlight a potentially large role for recombination-suppressing chromosomal inversions in local adaptation and suggest a hypothesis to explain the maintenance of differences in reproductive timing found at relatively small spatial scales across king scallop populations.

Hexabothriidae and Monocotylidae (Monogenoidea) from the gills of neonate hammerhead sharks (Sphyrnidae) Sphyrna gilberti, Sphyrna lewini and their hybrids from the western North Atlantic Ocean.

Neonates of hammerhead sharks (Sphyrnidae), Sphyrna lewini (Griffith and Smith, 1834), the sympatric cryptic species, Sphyrna gilberti Quattro et al., 2013, and their hybrids were captured in the western North Atlantic, along the coast of South Carolina, USA, between 2018 and 2019 and examined for gill monogenoids. Parasites were identified and redescribed from the gills of 79 neonates, and DNA sequences from partial fragments of the nuclear 28S ribosomal RNA (rDNA) and cytochrome c oxidase I mitochondrial DNA (COI) genes were generated to confirm species identifications. Three species of monogenoids from Hexabothriidae Price, 1942 and Monocotylidae Taschenberg, 1879 were determined and redescribed. Two species of Hexabothriidae, Erpocotyle microstoma (Brooks, 1934) and Erpocotyle sphyrnae (MacCallum, 1931), infecting both species of Sphyrna and hybrids; and 1 species of Monocotylidae, Loimosina wilsoni Manter, 1944, infecting only S. lewini and hybrids. Loimosina wilsoni 28S rDNA sequences matched those of Loimosina sp. from the southern coast of Brazil. Based on limited morphological analysis, Loimosina parawilsoni is likely a junior synonym of L. wilsoni. This is the first taxonomic study of monogenoids infecting S. gilberti and hybrids of S. gilberti and S. lewini.

Genetic diversity and structure of circumtropical almaco jack, Seriola rivoliana: tool for conservation and management.

The almaco jack, Seriola rivoliana, is a circumtropical pelagic fish of importance both in commercial fisheries and in aquaculture. To understand levels of genetic diversity within and among populations in the wild, population genetic structure and the relative magnitude of migration were assessed using mtDNA sequence data and single nucleotide polymorphisms (SNPs) from individuals sampled from locations in the Pacific and Atlantic Oceans. A total of 25 variable sites of cytochrome c oxidase subunit 1 and 3678 neutral SNPs were recovered. Three genetic groups were identified, with both marker types distributed in different oceanic regions: Pacific-1 in central Pacific, Pacific-2 in eastern Pacific and Atlantic in western Atlantic. Nonetheless, the analysis of SNP identified a fourth population in the Pacific coast of Baja California Sur, Mexico (Pacific-3), whereas that of mtDNA did not. This mito-nuclear discordance is likely explained by a recently diverged Pacific-3 population. In addition, two mtDNA haplogroups were found within the western Atlantic, likely indicating that the species came into the Atlantic from the Indian Ocean with historical gene flow from the eastern Pacific. Relative gene flow among ocean basins was low with r m < 0.2, whereas in the eastern Pacific it was asymmetric and higher from south to north (r m > 0.79). The results reflect the importance of assessing genetic structure and gene flow of natural populations for the purposes of sustainable management.

Hybridization between sympatric hammerhead sharks in the western North Atlantic Ocean.

Hybridization between closely related species has been documented across a wide range of taxa but has not been well studied in elasmobranchs. Hammerhead sharks have drawn global conservation concern because they experience some of the highest mortality rates among sharks when interacting with fisheries. Here we report on the detection of hybrids between the globally distributed scalloped hammerhead ( Sphyrna lewini) and recently described Carolina hammerhead ( S. gilberti) which are only known from the western Atlantic Ocean. Using a genomics approach, 10 first-generation hybrids and 15-17 backcrosses were detected from 554 individuals. The identification of backcrosses demonstrates hybrids are viable, and all backcrosses but one involved a scalloped hammerhead. All hybrids but one possessed Carolina hammerhead mtDNA, indicating sex-biased gene flow between species. Repeated hybridization and backcrossing with scalloped hammerheads could lead to the loss of endemic Carolina hammerheads.

Genetic mapping and comparative genomics to inform restoration enhancement and culture of southern flounder, Paralichthys lethostigma.

BACKGROUND: Southern flounder, Paralichthys lethostigma, historically support a substantial fishery along the Atlantic and Gulf coasts of the southern United States. Low year-class strengths over the past few years in the western Gulf of Mexico have raised concern that spawning stocks may be overfished. Current management of the resource includes releasing hatchery-raised juveniles to restock bays and estuaries; additionally, there is a growing interest in the potential for commercial aquaculture of the species. Currently, genomic resources for southern flounder do not exist. Here, we used two hatchery-reared families and double-digest, restriction-site-associated DNA (ddRAD) sequencing to create a reduced-representation genomic library consisting of several thousand single nucleotide polymorphisms (SNPs) located throughout the genome. RESULTS: The relative position of each SNP-containing locus was determined to create a high-density genetic map spanning the 24 linkage groups of the southern flounder genome. The consensus map was used to identify regions of shared synteny between southern flounder and seven other fish species for which genome assemblies are available. Finally, syntenic blocks were used to localize genes identified from transcripts in European flounder as potentially being involved in ecotoxicological and osmoregulatory responses, as well as QTLs associated with growth and disease resistance in Japanese flounder, on the southern flounder linkage map. CONCLUSIONS: The information provided by the linkage map will enrich restoration efforts by providing a foundation for interpreting spatial genetic variation within the species, ultimately furthering an understanding of the adaptive potential and resilience of southern flounder to future changes in local environmental conditions. Further, the map will facilitate the use of genetic markers to enhance restoration and commercial aquaculture.

These aren't the loci you'e looking for: Principles of effective SNP filtering for molecular ecologists.

Sequencing reduced-representation libraries of restriction site-associated DNA (RADseq) to identify single nucleotide polymorphisms (SNPs) is quickly becoming a standard methodology for molecular ecologists. Because of the scale of RADseq data sets, putative loci cannot be assessed individually, making the process of filtering noise and correctly identifying biologically meaningful signal more difficult. Artefacts introduced during library preparation and/or bioinformatic processing of SNP data can create patterns that are incorrectly interpreted as indicative of population structure or natural selection. Therefore, it is crucial to carefully consider types of errors that may be introduced during laboratory work and data processing, and how to minimize, detect and remove these errors. Here, we discuss issues inherent to RADseq methodologies that can result in artefacts during library preparation and locus reconstruction resulting in erroneous SNP calls and, ultimately, genotyping error. Further, we describe steps that can be implemented to create a rigorously filtered data set consisting of markers accurately representing independent loci and compare the effect of different combinations of filters on four RAD data sets. At last, we stress the importance of publishing raw sequence data along with final filtered data sets in addition to detailed documentation of filtering steps and quality control measures.

Molecular phylogenetics of New World searobins (Triglidae; Prionotinae).

Phylogenetic relationships among members of the New World searobin genera Bellator and Prionotus (Family Triglidae, Subfamily Prionotinae) and among other searobins in the families Triglidae and Peristediidae were investigated using both mitochondrial and nuclear DNA sequences. Phylogenetic hypotheses derived from maximum likelihood and Bayesian methodologies supported a monophyletic Prionotinae that included four well resolved clades of uncertain relationship; three contained species in the genus Prionotus and one contained species in the genus Bellator. Bellator was always recovered within the genus Prionotus, a result supported by post hoc model testing. Two nominal species of Prionotus (P. alatus and P. paralatus) were not recovered as exclusive lineages, suggesting the two may comprise a single species. Phylogenetic hypotheses also supported a monophyletic Triglidae but only if armored searobins (Family Peristediidae) were included. A robust morphological assessment is needed to further characterize relationships and suggest classification of clades within Prionotinae; for the time being we recommend that Bellator be considered a synonym of Prionotus. Relationships between armored searobins (Family Peristediidae) and searobins (Family Triglidae) and relationships within Triglidae also warrant further study.

Complex patterns of genetic population structure in the mouthbrooding marine catfish, Bagre marinus, in the Gulf of Mexico and U.S. Atlantic.

Patterns of genetic variation reflect interactions among microevolutionary forces that vary in strength with changing demography. Here, patterns of variation within and among samples of the mouthbrooding gafftopsail catfish (Bagre marinus, Family Ariidae) captured in the U.S. Atlantic and throughout the Gulf of Mexico were analyzed using genomics to generate neutral and non-neutral SNP data sets. Because genomic resources are lacking for ariids, linkage disequilibrium network analysis was used to examine patterns of putatively adaptive variation. Finally, historical demographic parameters were estimated from site frequency spectra. The results show four differentiated groups, corresponding to the (1) U.S. Atlantic, and the (2) northeastern, (3) northwestern, and (4) southern Gulf of Mexico. The non-neutral data presented two contrasting signals of structure, one due to increases in diversity moving west to east and north to south, and another to increased heterozygosity in the Atlantic. Demographic analysis suggested that recently reduced long-term effective population size in the Atlantic is likely an important driver of patterns of genetic variation and is consistent with a known reduction in population size potentially due to an epizootic. Overall, patterns of genetic variation resemble that of other fishes that use the same estuarine habitats as nurseries, regardless of the presence/absence of a larval phase, supporting the idea that adult/juvenile behavior and habitat are important predictors of contemporary patterns of genetic structure.

Haploid gynogens facilitate disomic marker development in paleotetraploid sturgeons.

Acipenseriformes (sturgeons and paddlefishes) are of substantial conservation concern, and development of genomic resources for these species is difficult due to past whole genome duplication. Development of disomic markers for polyploid organisms can be challenging due to difficulty in resolving alleles at a single locus from those among duplicated loci. In this study, we detail the development of disomic markers for the endangered pallid sturgeon (Scaphirhynchus albus) found in North America. One of the strategies for pallid sturgeon conservation is to stock U.S. rivers with offspring of pure pallid sturgeon, but introgression with the sympatric shovelnose sturgeon (S. platorynchus) threatens pallid sturgeon genetic integrity. Currently, 19 microsatellite loci are used to differentiate between both species and their hybrids, but the markers are insufficient to robustly identify backcrosses. We performed double digest restriction site-associated DNA sequencing (ddRADseq) on shovelnose sturgeon haploid gynogens to produce a reduced-representation genomic reference. Contiguous sequences that were heterozygous within a haploid individual were flagged as potentially encompassing multiple loci. Approximately 60 individuals of each species from two management units were sequenced, and reads were mapped to the haploid reference to identify single nucleotide polymorphisms (SNPs) at individual loci. The final data set contained 11,082 microhaplotyped loci which offer at least an order of magnitude greater resolution for species discrimination than the current panel of 19 microsatellites. These markers will be used to examine a larger sample of Scaphirhynchus individuals throughout their ranges to determine the extent and trajectory of hybridization.

Disentangling complex genomic signals to understand population structure of an exploited, estuarine-dependent flatfish.

Interpreting contemporary patterns of population structure requires an understanding of the interactions among microevolutionary forces and past demographic events. Here, 4,122 SNP-containing loci were used to assess structure in southern flounder (Paralichthys lethostigma) sampled across its range in the US Atlantic Ocean (Atlantic) and Gulf of Mexico (Gulf) and relationships among components of genomic variation and spatial and environmental variables were assessed across estuarine population samples in the Gulf. While hierarchical amova revealed significant heterogeneity within and between the Atlantic and Gulf, pairwise comparisons between samples within ocean basins demonstrated that all significant heterogeneity occurred within the Gulf. The distribution of Tajima's D estimated at a genome-wide scale differed significantly from equilibrium in all estuaries, with more negative values occurring in the Gulf. Components of genomic variation were significantly associated with environmental variables describing individual estuaries, and environment explained a larger component of variation than spatial proximity. Overall, results suggest that there is genetic spatial autocorrelation caused by shared larval sources for proximal nurseries (migration/drift), but that it is modified by environmentally driven differentiation (selection). This leads to conflicting signals in different parts of the genome and creates patterns of divergence that do not correspond to paradigms of strong local directional selection.

Map-like use of Earth's magnetic field in sharks.

Migration is common in marine animals,1-5 and use of the map-like information of Earth's magnetic field appears to play an important role.2,6-9 While sharks are iconic migrants10-12 and well known for their sensitivity to electromagnetic fields,13-20 whether this ability is used for navigation is unresolved.14,17,21,22 We conducted magnetic displacement experiments on wild-caught bonnetheads (Sphyrna tiburo) and show that magnetic map cues can elicit homeward orientation. We further show that use of a magnetic map to derive positional information may help explain aspects of the genetic structure of bonnethead populations in the northwest Atlantic.23-26 These results offer a compelling explanation for the puzzle of how migratory routes and population structure are maintained in marine environments, where few physical barriers limit movements of vagile species. VIDEO ABSTRACT.

World phylogeography and male-mediated gene flow in the sandbar shark, Carcharhinus plumbeus.

The sandbar shark, Carcharhinus plumbeus, is a large, cosmopolitan, coastal species. Females are thought to show philopatry to nursery grounds while males potentially migrate long distances, creating an opportunity for male-mediated gene flow that may lead to discordance in patterns revealed by mitochondrial DNA (mtDNA) and nuclear markers. While this dynamic has been investigated in elasmobranchs over small spatial scales, it has not been examined at a global level. We examined patterns of historical phylogeography and contemporary gene flow by genotyping 329 individuals from nine locations throughout the species' range at eight nuclear microsatellite markers and sequencing the complete mtDNA control region. Pairwise comparisons often resulted in fixation indices and divergence estimates of greater magnitude using mtDNA sequence data than microsatellite data. In addition, multiple methods of estimation suggested fewer populations based on microsatellite loci than on mtDNA sequence data. Coalescent analyses suggest divergence and restricted migration among Hawaii, Taiwan, eastern and western Australia using mtDNA sequence data and no divergence and high migration rates, between Taiwan and both Australian sites using microsatellite data. Evidence of secondary contact was detected between several localities and appears to be discreet in time rather than continuous. Collectively, these data suggest complex spatial/temporal relationships between shark populations that may feature pulses of female dispersal and more continuous male-mediated gene flow.

Evolution of population structure in an estuarine-dependent marine fish.

Restriction site-associated DNA (RAD) sequencing was used to characterize neutral and adaptive genetic variation among geographic samples of red drum, Sciaenops ocellatus, an estuarine-dependent fish found in coastal waters along the southeastern coast of the United States (Atlantic) and the northern Gulf of Mexico (Gulf). Analyses of neutral and outlier loci revealed three genetically distinct regional clusters: one in the Atlantic and two in the northern Gulf. Divergence in neutral loci indicated gradual genetic change and followed a linear pattern of isolation by distance. Divergence in outlier loci was at least an order of magnitude greater than divergence in neutral loci, and divergence between the regions in the Gulf was twice that of divergence between other regions. Discordance in patterns of genetic divergence between outlier and neutral loci is consistent with the hypothesis that the former reflects adaptive responses to environmental factors that vary on regional scales, while the latter largely reflects drift processes. Differences in basic habitat, initiated by glacial retreat and perpetuated by contemporary oceanic and atmospheric forces interacting with the geomorphology of the northern Gulf, followed by selection, appear to have led to reduced gene flow among red drum across the northern Gulf, reinforcing differences accrued during isolation and resulting in continued divergence across the genome. This same dynamic also may pertain to other coastal or nearshore fishes (18 species in 14 families) where genetically or morphologically defined sister taxa occur in the three regions.

Linkage Mapping and Comparative Genomics of Red Drum (Sciaenops ocellatus) Using Next-Generation Sequencing.

Developments in next-generation sequencing allow genotyping of thousands of genetic markers across hundreds of individuals in a cost-effective manner. Because of this, it is now possible to rapidly produce dense genetic linkage maps for nonmodel species. Here, we report a dense genetic linkage map for red drum, a marine fish species of considerable economic importance in the southeastern United States and elsewhere. We used a prior microsatellite-based linkage map as a framework and incorporated 1794 haplotyped contigs derived from high-throughput, reduced representation DNA sequencing to produce a linkage map containing 1794 haplotyped restriction-site associated DNA (RAD) contigs, 437 anonymous microsatellites, and 44 expressed sequence-tag-linked microsatellites (EST-SSRs). A total of 274 candidate genes, identified from transcripts from a preliminary hydrocarbon exposure study, were localized to specific chromosomes, using a shared synteny approach. The linkage map will be a useful resource for red drum commercial and restoration aquaculture, and for better understanding and managing populations of red drum in the wild.

Haplotyping RAD loci: an efficient method to filter paralogs and account for physical linkage.

Next-generation sequencing of reduced-representation genomic libraries provides a powerful methodology for genotyping thousands of single-nucleotide polymorphisms (SNPs) among individuals of nonmodel species. Utilizing genotype data in the absence of a reference genome, however, presents a number of challenges. One major challenge is the trade-off between splitting alleles at a single locus into separate clusters (loci), creating inflated homozygosity, and lumping multiple loci into a single contig (locus), creating artefacts and inflated heterozygosity. This issue has been addressed primarily through the use of similarity cut-offs in sequence clustering. Here, two commonly employed, postclustering filtering methods (read depth and excess heterozygosity) used to identify incorrectly assembled loci are compared with haplotyping, another postclustering filtering approach. Simulated and empirical data sets were used to demonstrate that each of the three methods separately identified incorrectly assembled loci; more optimal results were achieved when the three methods were applied in combination. The results confirmed that including incorrectly assembled loci in population-genetic data sets inflates estimates of heterozygosity and deflates estimates of population divergence. Additionally, at low levels of population divergence, physical linkage between SNPs within a locus created artificial clustering in analyses that assume markers are independent. Haplotyping SNPs within a locus effectively neutralized the physical linkage issue without having to thin data to a single SNP per locus. We introduce a Perl script that haplotypes polymorphisms, using data from single or paired-end reads, and identifies potentially problematic loci.

Fine-scale partitioning of genomic variation among recruits in an exploited fishery: causes and consequences.

Conservation and management of exploited species depends on accurate knowledge of how genetic variation is partitioned across a fishery, especially as it relates to recruitment. Using double-digest restriction-site associated DNA sequencing, we surveyed variation in 7,382 single nucleotide polymorphisms (SNPs) in red snapper (Lutjanus campechanus) young-of-the-year (YOY) sampled at six localities and in adults sampled at two localities in the northern Gulf of Mexico. Significant genetic heterogeneity was detected between the two adult samples, separated by ~600 km, and at spatial scales less than five kilometers among samples of YOY. Genetic differences between YOY samples and between YOY samples and adult samples were not associated with geographic distance, and a genome scan revealed no evidence of loci under selection. Estimates of the effective number of breeders, allelic richness, and relatedness within YOY samples were not consistent with sweepstakes recruitment. Instead, the data demonstrate, at least within one recruitment season, that multiple pulses of recruits originate from distinct groups of spawning adults, even at small spatial scales. For exploited species with this type of recruitment pattern, protection of spawning adults over wide geographic areas may be critical for ensuring productivity and stability of the fishery by maintaining larval supply and connectivity.

Kin-Aggregations Explain Chaotic Genetic Patchiness, a Commonly Observed Genetic Pattern, in a Marine Fish.

The phenomenon of chaotic genetic patchiness is a pattern commonly seen in marine organisms, particularly those with demersal adults and pelagic larvae. This pattern is usually associated with sweepstakes recruitment and variable reproductive success. Here we investigate the biological underpinnings of this pattern in a species of marine goby Coryphopterus personatus. We find that populations of this species show tell-tale signs of chaotic genetic patchiness including: small, but significant, differences in genetic structure over short distances; a non-equilibrium or "chaotic" pattern of differentiation among locations in space; and within locus, within population deviations from the expectations of Hardy-Weinberg equilibrium (HWE). We show that despite having a pelagic larval stage, and a wide distribution across Caribbean coral reefs, this species forms groups of highly related individuals at small spatial scales (<10 metres). These spatially clustered family groups cause the observed deviations from HWE and local population differentiation, a finding that is rarely demonstrated, but could be more common than previously thought.

A comparison of neutral and immune genetic variation in Atlantic salmon, Salmo salar L. in Chilean aquaculture facilities.

Genetic diversity was assessed in samples of cultured Atlantic salmon, Salmo salar L., obtained from facilities in Chile between 2005 and 2010, a period of time during which the infectious pathogens Infectious Salmon Anemia (ISA) virus, Caligus rogercresseyi (sea lice), and Piscirickettsia salmonis (salmon rickettsial syndrome) were common. Two panels of microsatellite markers were utilized: one with microsatellites with no known gene associations (neutral) and one featuring microsatellites linked to putative immune-related genes (immune-related). Allelic richness and gene diversity across samples were significantly greater in neutral loci as compared to immune-related loci. Both diversity measures were homogeneous among samples for immune-related loci and heterogeneous among samples for neutral loci. Immune-related loci were identified as F(ST) outliers in pairwise comparisons of samples at a 10-fold higher frequency than neutral loci. These results indicate that neutral and immune-related portions of the Atlantic salmon genome may have differed in response to the gauntlet of pathogens and that monitoring of specific, well characterized immune-related loci as well as neutral loci in cultured species could be useful when disease control and prevention is a goal.