Population genomics reveals seahorses (Hippocampus erectus) of the western mid-Atlantic coast to be residents rather than vagrants

Meta-Transcriptomic Analysis Reveals Novel RNA Viruses in Hippocampus erectus

Lined seahorse, Hippocampus erectus, is an important aquatic animal due to its medicinal and ornamental purposes. However, our understanding of the viral spectrum in H. erectus is still limited. Here, we studied the viruses in H. erectus using meta-transcriptomic sequencing. A total of 213,770,166 reads were generated and assembled de novo into 539 virus-associated contigs. Three novel RNA viruses from the Astroviridae, Paramyxoviridae, and Picornaviridae families were finally identified. In addition, we identified a strain of nervous necrosis virus from H. erectus. In particular, the unhealthy group showed a higher viral diversity and abundance than the normal group. These results revealed the diversity and cross-species transmission of viruses in H. erectus and highlighted the threat of viral infections to H. erectus.

Influence of the Atlantic Ocean thermal anomaly on the Longsnout seahorse Hippocampus reidi in a Brazilian estuary

One of the consequences of climate change is an increase in the temperature of the oceans, which is considered to be one of the greatest impacts on biodiversity. Fish may respond to this impact in several ways, including shifts in their patterns of occurrence. The present study investigated the variation in the structure of a H. reidi population between 2015 and 2017 in the northern Guaíba Island area, highlighting a possible relationship to thermal anomaly associated with the El Niño phenomenon. The seahorse population monitoring was performed monthly, recording sex ratio, abundance, juvenile and adult proportion, depth of occurrence, total length and the holdfast which the seahorse were found attached. The influence of the El Niño event on the study population was evaluated by the correlation of the thermal anomaly data reported for the Tropical South Atlantic Index. Seahorse density on northern Guaíba island was positively and significantly correlated with water temperature, but the sex ratio and number of juveniles were not. The diversity of holdfasts used increased over the study period and was inversely proportional to the thermal anomaly. These results suggest that the thermal anomalies caused by the El Niño in the South Atlantic might trigger migration behaviour in the study species, providing a large aggregation during that period in Guaíba island.

Diversity of Seahorse Species (Hippocampus spp.) in the International Aquarium Trade

Seahorses (Hippocampus spp.) are threatened as a result of habitat degradation and overfishing. They have commercial value as traditional medicine, curio objects, and pets in the aquarium industry. There are 48 valid species, 27 of which are represented in the international aquarium trade. Most species in the aquarium industry are relatively large and were described early in the history of seahorse taxonomy. In 2002, seahorses became the first marine fishes for which the international trade became regulated by CITES (Convention for the International Trade in Endangered Species of Wild Fauna and Flora), with implementation in 2004. Since then, aquaculture has been developed to improve the sustainability of the seahorse trade. This review provides analyses of the roles of wild-caught and cultured individuals in the international aquarium trade of various Hippocampus species for the period 1997–2018. For all species, trade numbers declined after 2011. The proportion of cultured seahorses in the aquarium trade increased rapidly after their listing in CITES, although the industry is still struggling to produce large numbers of young in a cost-effective way, and its economic viability is technically challenging in terms of diet and disease. Whether seahorse aquaculture can benefit wild populations will largely depend on its capacity to provide an alternative livelihood for subsistence fishers in the source countries. For most species, CITES trade records of live animals in the aquarium industry started a few years earlier than those of dead bodies in the traditional medicine trade, despite the latter being 15 times higher in number. The use of DNA analysis in the species identification of seahorses has predominantly been applied to animals in the traditional medicine market, but not to the aquarium trade. Genetic tools have already been used in the description of new species and will also help to discover new species and in various other kinds of applications.

The population genomics of repeated freshwater colonizations by Gulf pipefish

How organisms adapt to the novel challenges imposed by the colonization of a new habitat has long been a central question in evolutionary biology. When multiple populations of the same species independently adapt to similar environmental challenges, the question becomes whether the populations have arrived at their adaptations through the same genetic mechanisms. In recent years, genetic techniques have been used to tackle these questions by investigating the genome-level changes underlying local adaptation. Here, we present a genomic analysis of colonization of freshwater habitats by a primarily marine fish, the Gulf pipefish (Syngnathus scovelli). We sample pipefish from four geographically distinct freshwater locations and use double-digest restriction site associated DNA sequencing to compare them to 12 previously studied saltwater populations. The two most geographically distant and isolated freshwater populations are the most genetically distinct, although demographic analysis suggests that these populations are experiencing ongoing migration with their saltwater neighbours. Additionally, outlier regions were found genome-wide, showing parallelism across ecotype pairs. We conclude that these multiple freshwater colonizations involve similar genomic regions, despite the large geographical distances and different underlying mechanisms. These similar patterns are probably facilitated by the interacting effects of intrinsic barriers, gene flow among populations and ecological selection in the Gulf pipefish.

Genomic and morphological evidence of distinct populations in the endemic common (weedy) seadragon Phyllopteryx taeniolatus (Syngnathidae) along the east coast of Australia

The common or weedy seadragon, Phyllopteryx taeniolatus, is an iconic and endemic fish found across temperate reefs of southern Australia. Despite its charismatic nature, few studies have been published, and the extent of population sub-structuring remains poorly resolved. Here we used 7462 single nucleotide polymorphisms (SNPs) to identify the extent of population structure in the weedy seadragon along the temperate southeast coast of Australia. We identified four populations, with strong genetic structure (F_ST = 0.562) between them. Both Discriminant Analysis of Principle Components (DAPC) and Bayesian clustering analyses support four distinct genetic clusters (north to south: central New South Wales, southern NSW, Victoria and Tasmania). In addition to these genetic differences, geographical variation in external morphology was recorded, with individuals from New South Wales shaped differently for a few measurements to those from the Mornington Peninsula (Victoria). We posit that these genetic and morphological differences suggest that the Victorian population of P. taeniolatus was historically isolated by the Bassian Isthmus during the last glacial maximum and should now be considered at least a distinct population. We also recorded high levels of genetic structure among the other locations. Based on the genomic and to a degree morphological evidence presented in this study, we recommend that the Victorian population be managed separately from the eastern populations (New South Wales and Tasmania).

Adaptive divergence, neutral panmixia, and algal symbiont population structure in the temperate coral Astrangia poculata along the Mid-Atlantic United States

Astrangia poculata is a temperate scleractinian coral that exists in facultative symbiosis with the dinoflagellate alga Breviolum psygmophilum across a range spanning the Gulf of Mexico to Cape Cod, Massachusetts. Our previous work on metabolic thermal performance of Virginia (VA) and Rhode Island (RI) populations of A. poculata revealed physiological signatures of cold (RI) and warm (VA) adaptation of these populations to their respective local thermal environments. Here, we used whole-transcriptome sequencing (mRNA-Seq) to evaluate genetic differences and identify potential loci involved in the adaptive signature of VA and RI populations. Sequencing data from 40 A. poculata individuals, including 10 colonies from each population and symbiotic state (VA-white, VA-brown, RI-white, and RI-brown), yielded a total of 1,808 host-associated and 59 algal symbiont-associated single nucleotide polymorphisms (SNPs) post filtration. Fst outlier analysis identified 66 putative high outlier SNPs in the coral host and 4 in the algal symbiont. Differentiation of VA and RI populations in the coral host was driven by putatively adaptive loci, not neutral divergence (Fst = 0.16, p = 0.001 and Fst = 0.002, p = 0.269 for outlier and neutral SNPs respectively). In contrast, we found evidence of neutral population differentiation in B. psygmophilum (Fst = 0.093, p = 0.001). Several putatively adaptive host loci occur on genes previously associated with the coral stress response. In the symbiont, three of four putatively adaptive loci are associated with photosystem proteins. The opposing pattern of neutral differentiation in B. psygmophilum, but not the A. poculata host, reflects the contrasting dynamics of coral host and algal symbiont population connectivity, dispersal, and gene by environment interactions.

Genotyping-by-sequencing and ecological niche modeling illuminate phylogeography, admixture, and Pleistocene range dynamics in quaking aspen (Populus tremuloides)

Populus tremuloides is the widest-ranging tree species in North America and an ecologically important component of mesic forest ecosystems displaced by the Pleistocene glaciations. Using phylogeographic analyses of genome-wide SNPs (34,796 SNPs, 183 individuals) and ecological niche modeling, we inferred population structure, ploidy levels, admixture, and Pleistocene range dynamics of P. tremuloides, and tested several historical biogeographical hypotheses. We found three genetic lineages located mainly in coastal-Cascades (cluster 1), east-slope Cascades-Sierra Nevadas-Northern Rockies (cluster 2), and U.S. Rocky Mountains through southern Canadian (cluster 3) regions of the P. tremuloides range, with tree graph relationships of the form ((cluster 1, cluster 2), cluster 3). Populations consisted mainly of diploids (86%) but also small numbers of triploids (12%) and tetraploids (1%), and ploidy did not adversely affect our genetic inferences. The main vector of admixture was from cluster 3 into cluster 2, with the admixture zone trending northwest through the Rocky Mountains along a recognized phenotypic cline (Utah to Idaho). Clusters 1 and 2 provided strong support for the "stable-edge hypothesis" that unglaciated southwestern populations persisted in situ since the last glaciation. By contrast, despite a lack of clinal genetic variation, cluster 3 exhibited "trailing-edge" dynamics from niche suitability predictions signifying complete northward postglacial expansion. Results were also consistent with the "inland dispersal hypothesis" predicting postglacial assembly of Pacific Northwestern forest ecosystems, but rejected the hypothesis that Pacific-coastal populations were colonized during outburst flooding from glacial Lake Missoula. Overall, congruent patterns between our phylogeographic and ecological niche modeling results and fossil pollen data demonstrate complex mixtures of stable-edge, refugial locations, and postglacial expansion within P. tremuloides. These findings confirm and refine previous genetic studies, while strongly supporting a distinct Pacific-coastal genetic lineage of quaking aspen.

Morphological and molecular evidence for first records and range extension of the Japanese seahorse, Hippocampus mohnikei (Bleeker 1853) in a bay-estuarine system of Goa, central west coast of India

Accurate information of taxonomy and geographic range of seahorse species (genus Hippocampus) is the first step in preparing threat assessments and designing effective conservation measures. Here, we report first records and a range extension of the Japanese seahorse, Hippocampus mohnikei (Bleeker, 1853) from the Mandovi estuarine ecosystem of Goa, central west coast of India (CWCI) based on morphological and molecular analyses. The morphometric and meristic traits, particularly short snout (29–35% head length), double cheek spine, low coronet, long tail (51.2–57.9% of standard length), 11 trunk rings, 37–39 tail rings, 15–16 dorsal and 12–14 fin rays observed in four collected specimens matched with the reported key diagnostic morphological criteria of vouchered specimens of H. mohnikei. The seahorse mitochondrial cytochrome oxidase subunit I (COI) and cytochrome b (Cyt b) genes were partially sequenced for conclusive genetic identification of the species under study. Molecular analysis showed that all four individuals clustered together suggesting a monophyletic lineage. Using the maximum similarity with GenBank database, maximum likelihood network and subsequent morphological analysis, the identity of the collected seahorse species was reconfirmed as H. mohnikei. With this new report, the geographic range of H. mohnikei extended significantly to the west from its previously known range. This new sighting of H. mohnikei could indicate a long-distance dispersal facilitated by the prevailing oceanic circulation in the Indo-Pacific region or increased habitat suitability in bay-estuarine systems of Goa, CWCI. Comparison of the pair-wise genetic distances (Kimura 2-parameter) based on COI and Cyt b sequences revealed that the specimens examined in this study are genetically closer to H. mohnikei populations from Vietnam and Thailand than they are to those in Japan and China. To test the hypothesis whether H. mohnikei are vagrants or previously unreported established population, long-term inter-annual sampling and analyses are warranted.

A haplotype-resolved draft genome of the European sardine (Sardina pilchardus)

Background

The European sardine (Sardina pilchardus Walbaum, 1792) is culturally and economically important throughout its distribution. Monitoring studies of sardine populations report an alarming decrease in stocks due to overfishing and environmental change, which has resulted in historically low captures along the Iberian Atlantic coast. Important biological and ecological features such as population diversity, structure, and migratory patterns can be addressed with the development and use of genomics resources.

Findings

The genome of a single female individual was sequenced using Illumina HiSeq X Ten 10x Genomics linked reads, generating 113.8 gigabase pairs of data. Three draft genomes were assembled: 2 haploid genomes with a total size of 935 megabase pairs (N50 103 kilobase pairs) each, and a consensus genome of total size 950 megabase pairs (N50 97 kilobase pairs). The genome completeness assessment captured 84% of Actinopterygii Benchmarking Universal Single-Copy Orthologs. To obtain a more complete analysis, the transcriptomes of 11 tissues were sequenced to aid the functional annotation of the genome, resulting in 40,777 genes predicted. Variant calling on nearly half of the haplotype genome resulted in the identification of >2.3 million phased single-nucleotide polymorphisms with heterozygous loci.

Conclusions

A draft genome was obtained, despite a high level of sequence repeats and heterozygosity, which are expected genome characteristics of a wild sardine. The reference sardine genome and respective variant data will be a cornerstone resource of ongoing population genomics studies to be integrated into future sardine stock assessment modelling to better manage this valuable resource.

Range-wide population genomics of the Mexican fruit fly: Toward development of pathway analysis tools

Recurrently invading pests provide unique challenges for pest management, but also present opportunities to utilize genomics to understand invasion dynamics and inform regulatory management through pathway analysis. In the southern United States, the Mexican fruit fly Anastrepha ludens is such a pest, and its incursions into Texas and California represent major threats to the agricultural systems of those regions. We developed a draft genome assembly for A. ludens, conducted range-wide population genomics using restriction site-associated DNA sequencing, and then developed and demonstrated a panel of highly differentiated diagnostic SNPs for source determination of intercepted flies in this system. Using 2,081 genomewide SNPs, we identified four populations across the range of A. ludens, corresponding to western Mexico, eastern Mexico/Texas, Guatemala/Belize/Honduras, and Costa Rica/Panama, with some intergradation present between clusters, particularly in Central America. From this population genomics framework, we developed a diagnostic panel of 28 highly differentiated SNPs that were able to recreate the genomewide population structure in this species. We demonstrated this panel on a set of test specimens, including specimens intercepted as part of regular trapping surveillance in Texas and California, and we were able to predict populations of origin for these specimens. This methodology presents a highly applied use of genomic techniques and can be implemented in any group of recurrently invading pests.

Temperature-induced physiological stress and reproductive characteristics of the migratory seahorse Hippocampus erectus during a thermal stress simulation

Inshore-offshore migration occurs frequently in seahorse species, either because of prey opportunities or because they are driven by reproduction, and variations in water temperature may dramatically change migratory seahorse behavior and physiology. The present study investigated the behavioral and physiological responses of the lined seahorse Hippocampus erectus under thermal stress and evaluated the potential effects of different temperatures on its reproduction. The results showed that the thermal tolerance of the seahorses was time dependent. Acute thermal stress (30°C, 2-10 h) increased the basal metabolic rate (breathing rate) and the expression of stress response genes (Hsp genes) significantly and further stimulated seahorse appetite. Chronic thermal treatment (30°C, 4 weeks) led to a persistently higher basal metabolic rate, higher stress response gene expression and higher mortality rates, indicating that the seahorses could not acclimate to chronic thermal stress and might experience massive mortality rates due to excessively high basal metabolic rates and stress damage. Additionally, no significant negative effects on gonad development or reproductive endocrine regulation genes were observed in response to chronic thermal stress, suggesting that seahorse reproductive behavior could adapt to higher-temperature conditions during migration and within seahorse breeding grounds. In conclusion, this simulation experiment indicates that temperature variations during inshore-offshore migration have no effect on reproduction, but promote significantly high basal metabolic rates and stress responses. Therefore, we suggest that the observed high tolerance of seahorse reproduction is in line with the inshore-offshore reproductive migration pattern of lined seahorses.This article has an associated First Person interview with the first author of the paper.

Evidence for panmixia despite barriers to gene flow in the hooked mussel, Ischadium recurvum (Mytilidae; Brachidontinae) along the North American coastline

The discovery of phylogeographic patterns within broadly distributed marine species can be particularly challenging because absolute physical barriers to dispersal can be inconspicuous. Genetic boundaries often lie where ocean currents meet, forming sharp physical and ecological gradients, which may act as barriers to successful migrants. In eastern North America, coastal species often show phylogeographic differentiation associated with two recognized genetic barriers: the Gulf/Atlantic and the Virginia/Carolina discontinuities. We examined 185 specimens of the intertidal hooked mussel Ischadium recurvum collected from 15 locations along the eastern coastline of North America to examine phylogeographic, migration and historical demographic patterns associated climate change associated with Pleistocene glacial patterns. Hypothesis testing using Bayes factors in Migrate-n rejected the presence of phylogeographic breaks consistent with either maritime discontinuity and favoured a panmictic population model. The migration rate from the Gulf to the Atlantic was approximately three times higher than the migration from the Atlantic to the Gulf whereas the Carolina-Virginia migration rates were nearly equal. The summary statistics (Tajima's D, Fu's Fs) were significant and the demographic analyses (mismatch distributions, Bayesian skyline plot) were consistent with patterns of population expansion following glacial retreat during the Pleistocene epoch.

Exploring the phylogeography of a hexaploid freshwater fish by RAD sequencing

The KwaZulu‐Natal yellowfish (Labeobarbus natalensis) is an abundant cyprinid, endemic to KwaZulu‐Natal Province, South Africa. In this study, we developed a single‐nucleotide polymorphism (SNP) dataset from double‐digest restriction site‐associated DNA (ddRAD) sequencing of samples across the distribution. We addressed several hidden challenges, primarily focusing on proper filtering of RAD data and selecting optimal parameters for data processing in polyploid lineages. We used the resulting high‐quality SNP dataset to investigate the population genetic structure of L. natalensis. A small number of mitochondrial markers present in these data had disproportionate influence on the recovered genetic structure. The presence of singleton SNPs also confounded genetic structure. We found a well‐supported division into northern and southern lineages, with further subdivision into five populations, one of which reflects north–south admixture. Approximate Bayesian Computation scenario testing supported a scenario where an ancestral population diverged into northern and southern lineages, which then diverged to yield the current five populations. All river systems showed similar levels of genetic diversity, which appears unrelated to drainage system size. Nucleotide diversity was highest in the smallest river system, the Mbokodweni, which, together with adjacent small coastal systems, should be considered as a key catchment for conservation.

Landscape genetic approaches to guide native plant restoration in the Mojave Desert

Restoring dryland ecosystems is a global challenge due to synergistic drivers of disturbance coupled with unpredictable environmental conditions. Dryland plant species have evolved complex life-history strategies to cope with fluctuating resources and climatic extremes. Although rarely quantified, local adaptation is likely widespread among these species and potentially influences restoration outcomes. The common practice of reintroducing propagules to restore dryland ecosystems, often across large spatial scales, compels evaluation of adaptive divergence within these species. Such evaluations are critical to understanding the consequences of large-scale manipulation of gene flow and to predicting success of restoration efforts. However, genetic information for species of interest can be difficult and expensive to obtain through traditional common garden experiments. Recent advances in landscape genetics offer marker-based approaches for identifying environmental drivers of adaptive genetic variability in non-model species, but tools are still needed to link these approaches with practical aspects of ecological restoration. Here, we combine spatially explicit landscape genetics models with flexible visualization tools to demonstrate how cost-effective evaluations of adaptive genetic divergence can facilitate implementation of different seed sourcing strategies in ecological restoration. We apply these methods to Amplified Fragment Length Polymorphism (AFLP) markers genotyped in two Mojave Desert shrub species of high restoration importance: the long-lived, wind-pollinated gymnosperm Ephedra nevadensis, and the short-lived, insect-pollinated angiosperm Sphaeralcea ambigua. Mean annual temperature was identified as an important driver of adaptive genetic divergence for both species. Ephedra showed stronger adaptive divergence with respect to precipitation variability, while temperature variability and precipitation averages explained a larger fraction of adaptive divergence in Sphaeralcea. We describe multivariate statistical approaches for interpolating spatial patterns of adaptive divergence while accounting for potential bias due to neutral genetic structure. Through a spatial bootstrapping procedure, we also visualize patterns in the magnitude of model uncertainty. Finally, we introduce an interactive, distance-based mapping approach that explicitly links marker-based models of adaptive divergence with local or admixture seed sourcing strategies, promoting effective native plant restoration.

Genomic structure and expression pattern of MHC IIα and IIβ genes reveal an unusual immune trait in lined seahorse Hippocampus erectus

The major histocompatibility complex (MHC) genes are crucial in the adaptive immune system, and the gene duplication of MHC in animals can generally result in immune flexibility. In this study, we found that the lined seahorse (Hippocampus erectus) has only one gene copy number (GCN) of MHC IIα and IIβ, which is different from that in other teleosts. Together with the lack of spleen and gut-associated lymphatic tissue (GALT), the seahorse may be referred to as having a partial but natural "immunodeficiency". Highly variable amino acid residues were found in the IIα and IIβ domains, especially in the α1 and β1 domains with 9.62% and 8.43% allelic variation, respectively. Site models revealed seven and ten positively selected positions in the α1 and β1 domains, respectively. Real-time PCR experiments showed high expression levels of the MHC II genes in intestine (In), gill (Gi) and trunk kidney (TK) and medium in muscle (Mu) and brood pouch (BP), and the expression levels were significantly up-regulated after bacterial infection. Specially, relative higher expression level of both MHC IIα and IIβ was found in Mu and BP when compared with other fish species, in which MHC II is expressed negligibly in Mu. These results indicate that apart from TK, Gi and In, MU and BP play an important role in the immune response against pathogens in the seahorse. In conclusion, high allelic variation and strong positive selection in PBR and relative higher expression in MU and BP are speculated to partly compensate for the immunodeficiency.

Large-Scale Genotyping-by-Sequencing Indicates High Levels of Gene Flow in the Deep-Sea Octocoral Swiftia simplex (Nutting 1909) on the West Coast of the United States

Deep-sea corals are a critical component of habitat in the deep-sea, existing as regional hotspots for biodiversity, and are associated with increased assemblages of fish, including commercially important species. Because sampling these species is so difficult, little is known about the connectivity and life history of deep-sea octocoral populations. This study evaluates the genetic connectivity among 23 individuals of the deep-sea octocoral Swiftia simplex collected from Eastern Pacific waters along the west coast of the United States. We utilized high-throughput restriction-site associated DNA (RAD)-tag sequencing to develop the first molecular genetic resource for the deep-sea octocoral, Swiftia simplex. Using this technique we discovered thousands of putative genome-wide SNPs in this species, and after quality control, successfully genotyped 1,145 SNPs across individuals sampled from California to Washington. These SNPs were used to assess putative population structure across the region. A STRUCTURE analysis as well as a principal coordinates analysis both failed to detect any population differentiation across all geographic areas in these collections. Additionally, after assigning individuals to putative population groups geographically, no significant F_ST values could be detected (F_ST for the full data set 0.0056), and no significant isolation by distance could be detected (p = 0.999). Taken together, these results indicate a high degree of connectivity and potential panmixia in S. simplex along this portion of the continental shelf.

Dealing with paralogy in RADseq data: in silico detection and single nucleotide polymorphism validation in Robinia pseudoacacia L

The RADseq technology allows researchers to efficiently develop thousands of polymorphic loci across multiple individuals with little or no prior information on the genome. However, many questions remain about the biases inherent to this technology. Notably, sequence misalignments arising from paralogy may affect the development of single nucleotide polymorphism (SNP) markers and the estimation of genetic diversity. We evaluated the impact of putative paralog loci on genetic diversity estimation during the development of SNPs from a RADseq dataset for the nonmodel tree species Robinia pseudoacacia L. We sequenced nine genotypes and analyzed the frequency of putative paralogous RAD loci as a function of both the depth of coverage and the mismatch threshold allowed between loci. Putative paralogy was detected in a very variable number of loci, from 1% to more than 20%, with the depth of coverage having a major influence on the result. Putative paralogy artificially increased the observed degree of polymorphism and resulting estimates of diversity. The choice of the depth of coverage also affected diversity estimation and SNP validation: A low threshold decreased the chances of detecting minor alleles while a high threshold increased allelic dropout. SNP validation was better for the low threshold (4×) than for the high threshold (18×) we tested. Using the strategy developed here, we were able to validate more than 80% of the SNPs tested by means of individual genotyping, resulting in a readily usable set of 330 SNPs, suitable for use in population genetics applications.

A DNA-Based Assessment of the Phylogenetic Position of a Morphologically Distinct, Anchialine-Lake-Restricted Seahorse

Isolated populations provide special opportunities to study local adaptation and incipient speciation. In some cases, however, morphological evolution can obscure the taxonomic status of recently founded populations. Here, we use molecular markers to show that an anchialine-lake-restricted population of seahorses, originally identified as Hippocampus reidi, appears on the basis of DNA data to be Hippocampus erectus We collected seahorses from Sweetings Pond, on Eleuthera Island, Bahamas, during the summer of 2014. We measured morphological traits and sequenced 2 genes, cytochrome b and ribosomal protein S7, from 19 seahorses in our sample. On the basis of morphology, Sweetings Pond seahorses could not be assigned definitively to either of the 2 species of seahorse, H. reidi and H. erectus, that occur in marine waters surrounding the Bahamas. However, our DNA-based phylogenetic analysis showed that the Sweetings Pond fish were firmly nested within the H. erectus clade with a Bayesian posterior probability greater than 0.99. Thus, Sweetings Pond seahorses most recently shared a common ancestor with H. erectus populations from the Western Atlantic. Interestingly, the seahorses from Sweetings Pond differ morphologically from other marine populations of H. erectus in having a more even torso to tail length ratio. The substantial habitat differences between Sweetings Pond and the surrounding coastal habitat make Sweetings Pond seahorses particularly interesting from the perspectives of conservation, local adaptation, and incipient speciation.

Contrasting patterns of population connectivity between regions in a commercially important mollusc Haliotis rubra: integrating population genetics, genomics and marine LiDAR data

Estimating contemporary genetic structure and population connectivity in marine species is challenging, often compromised by genetic markers that lack adequate sensitivity, and unstructured sampling regimes. We show how these limitations can be overcome via the integration of modern genotyping methods and sampling designs guided by LiDAR and SONAR data sets. Here we explore patterns of gene flow and local genetic structure in a commercially harvested abalone species (Haliotis rubra) from southeastern Australia, where the viability of fishing stocks is believed to be dictated by recruitment from local sources. Using a panel of microsatellite and genomewide SNP markers, we compare allele frequencies across a replicated hierarchical sampling area guided by bathymetric LiDAR imagery. Results indicate high levels of gene flow and no significant genetic structure within or between benthic reef habitats across 1400 km of coastline. These findings differ to those reported for other regions of the fishery indicating that larval supply is likely to be spatially variable, with implications for management and long-term recovery from stock depletion. The study highlights the utility of suitably designed genetic markers and spatially informed sampling strategies for gaining insights into recruitment patterns in benthic marine species, assisting in conservation planning and sustainable management of fisheries.

Genome-wide SNP analysis explains coral diversity and recovery in the Ryukyu Archipelago

Following a global coral bleaching event in 1998, Acropora corals surrounding most of Okinawa island (OI) were devastated, although they are now gradually recovering. In contrast, the Kerama Islands (KIs) only 30 km west of OI, have continuously hosted a great variety of healthy corals. Taking advantage of the decoded Acropora digitifera genome and using genome-wide SNP analyses, we clarified Acropora population structure in the southern Ryukyu Archipelago (sRA). Despite small genetic distances, we identified distinct clusters corresponding to specific island groups, suggesting infrequent long-distance dispersal within the sRA. Although the KIs were believed to supply coral larvae to OI, admixture analyses showed that such dispersal is much more limited than previously realized, indicating independent recovery of OI coral populations and the necessity of local conservation efforts for each region. We detected strong historical migration from the Yaeyama Islands (YIs) to OI, and suggest that the YIs are the original source of OI corals. In addition, migration edges to the KIs suggest that they are a historical sink population in the sRA, resulting in high diversity. This population genomics study provides the highest resolution data to date regarding coral population structure and history.