Secondary contact and asymmetrical gene flow in a cosmopolitan marine fish across the Benguela upwelling zone

Strong philopatry in an estuarine-dependent fish

Understanding fish movement is critical in determining the spatial scales in which to appropriately manage wild populations. Genetic markers provide a natural tagging approach to assess the degree of gene flow and population connectivity across a species distribution. We investigated the genetic structure of black bream Acanthopagrus butcheri across its entire distribution range in Australia, as well as regional scale gene flow across south-eastern Australia by undertaking a comprehensive analysis of the populations in estuaries across the region. We applied genome-wide sampling of single-nucleotide polymorphism (SNP) markers generated from restriction site-associated DNA sequencing. Genetic structure and potential gene flow was assessed using principal component analyses and admixture analyses (STRUCTURE). Using 33,493 SNPs, we detected broad scale genetic structuring, with limited gene flow among regional clusters (i.e. Western Australia, South Australia and western Victoria; and eastern Victoria, Tasmania and New South Wales). This is likely the result of unsuitable habitats, strong ocean currents (e.g. the Leeuwin Current and the East Australian Current), large water bodies (e.g. Bass Strait) and known biogeographical provinces across the continent. Local-scale genetic structuring was also identified across the south-eastern Australian estuaries sampled, reflecting that the coexistence of both migratory and resident individuals within populations (i.e. partial migration), and the movement of fish into coastal waters, still results in strong philopatry across the region. Instances of movement among estuaries at this spatial scale were primarily found between adjacent estuaries and were likely attributed to lone migrants utilising inshore coastal currents for movement beyond nearby habitats. Targeting SNP markers in A. butcheri at this continental scale highlighted how neither spatial proximity of estuaries nor black bream's ability to move into coastal waters reflects increased gene flow. Overall, our findings highlight the importance of location-specific management.

Local adaptation with gene flow in a highly dispersive shark

Adaptive divergence in response to environmental clines are expected to be common in species occupying heterogeneous environments. Despite numerous advances in techniques appropriate for non-model species, gene-environment association studies in elasmobranchs are still scarce. The bronze whaler or copper shark (Carcharhinus brachyurus) is a large coastal shark with a wide distribution and one of the most exploited elasmobranchs in southern Africa. Here, we assessed the distribution of neutral and adaptive genomic diversity in C. brachyurus across a highly heterogeneous environment in southern Africa based on genome-wide SNPs obtained through a restriction site-associated DNA method (3RAD). A combination of differentiation-based genome-scan (outflank) and genotype-environment analyses (redundancy analysis, latent factor mixed models) identified a total of 234 differentiation-based outlier and candidate SNPs associated with bioclimatic variables. Analysis of 26,299 putatively neutral SNPs revealed moderate and evenly distributed levels of genomic diversity across sites from the east coast of South Africa to Angola. Multivariate and clustering analyses demonstrated a high degree of gene flow with no significant population structuring among or within ocean basins. In contrast, the putatively adaptive SNPs demonstrated the presence of two clusters and deep divergence between Angola and all other individuals from Namibia and South Africa. These results provide evidence for adaptive divergence in response to a heterogeneous seascape in a large, mobile shark despite high levels of gene flow. These results are expected to inform management strategies and policy at the national and regional level for conservation of C. brachyurus populations.

A Comparative Phylogeography of Three Marine Species with Different PLD Modes Reveals Two Genetic Breaks across the Southern Caribbean Sea

The comparative phylogeography of marine species with contrasting dispersal potential across the southern Caribbean Sea was evaluated by the presence of two putative barriers: the Magdalena River plume (MRP) and the combination of the absence of a rocky bottom and the almost permanent upwelling in the La Guajira Peninsula (ARB + PUG). Three species with varying biological and ecological characteristics (i.e., dispersal potentials) that inhabit shallow rocky bottoms were selected: Cittarium pica (PLD < 6 days), Acanthemblemaria rivasi (PLD < 22 days), and Nerita tessellata (PLD > 60 days). We generated a set of SNPs for the three species using the ddRad-seq technique. Samples of each species were collected in five locations from Capurganá to La Guajira. For the first time, evidence of a phylogeographic break caused by the MRP is provided, mainly for A. rivasi (AMOVA: ΦCT = 0.420). The ARB + PUG barrier causes another break for A. rivasi (ΦCT = 0.406) and C. pica (ΦCT = 0.224). Three populations (K = 3) were identified for A. rivasi and C. pica, while N. tessellata presented one population (K = 1). The Mantel correlogram indicated that A. rivasi and C. pica fit the hierarchical population model, and only the A. rivasi and C. pica comparisons showed phylogeographic congruence. Our results demonstrate how the biological traits of these three species and the biogeographic barriers have influenced their phylogeographic structure.

Three mitochondrial lineages and no Atlantic-Mediterranean barrier for the bogue Boops boops across its widespread distribution

Marine species exhibiting wide distributional ranges are frequently subdivided into discrete genetic units over limited spatial scales. This is often due to specific life-history traits or oceanographic barriers that prevent gene flow. Fine-scale sampling studies revealed distinct phylogeographic patterns in the northeastern Atlantic and the Mediterranean, ranging from panmixia to noticeable population genetic structure. Here, we used mitochondrial sequence data to analyse connectivity in the bogue Boops boops throughout most of its widespread distribution. Our results identified the existence of three clades, one comprising specimens from the Azores and eastern Atlantic/Mediterranean, another with individuals from the Canary Islands, Madeira and Cape Verde archipelagos, and the third with samples from Mauritania only. One of the branches of the northern subtropical gyre (Azores Current) that drifts towards the Gulf of Cádiz promotes a closer connection between the Azores, southern Portugal and the Mediterranean B. boops populations. The Almería-Oran Front, widely recognised as an oceanographic barrier for many organisms to cross the Atlantic-Mediterranean divide, does not seem to affect the dispersal of this benthopelagic species. The southward movement of the Cape Verde Frontal Zone during the winter, combined with the relatively short duration of the pelagic larval stage of B. boops, may be potential factors for preventing the connectivity between the Atlantic oceanic archipelagos and Mauritania shaping the genetic signature of this species.

Spatial and temporal stability in the genetic structure of a marine crab despite a biogeographic break

Elucidating the processes responsible for maintaining the population connectivity of marine benthic species mediated by larval dispersal remains a fundamental question in marine ecology and fishery management. Understanding these processes becomes particularly important in areas with a biogeographic break and unidirectional water movement along the sides of the break. Based on variability at 4209 single-nucleotide polymorphisms in 234 individuals, we determine the genetic structure, temporal genetic stability, and gene flow among populations of the commercially important mola rock crab Metacarcinus edwardsii in a system in southern Chile with a biogeographic break at latitude 42°S. Specimens were collected at eight sites within its geographic distribution, with collection at four of these sites was performed twice. Using population genetic approaches, we found no evidence of geographic or temporal population differentiation. Similarly, we found no evidence of an effect on gene flow of the biogeographic break caused by the the West Wind Drift Current. Moreover, migration analyses supported gene flow among all sites but at different rates for different pairs of sites. Overall, our findings indicate that M. edwardsii comprises a single large population with high levels of gene flow among sites separated by over 1700 km and demonstrate temporal stability in its genetic structure.

Combining population genomics with demographic analyses highlights habitat patchiness and larval dispersal as determinants of connectivity in coastal fish species

Gene flow shapes spatial genetic structure and the potential for local adaptation. Among marine animals with nonmigratory adults, the presence or absence of a pelagic larval stage is thought to be a key determinant in shaping gene flow and the genetic structure of populations. In addition, the spatial distribution of suitable habitats is expected to influence the distribution of biological populations and their connectivity patterns. We used whole genome sequencing to study demographic history and reduced representation (double‐digest restriction associated DNA) sequencing data to analyse spatial genetic structure in broadnosed pipefish (Syngnathus typhle). Its main habitat is eelgrass beds, which are patchily distributed along the study area in southern Norway. Demographic connectivity among populations was inferred from long‐term (~30‐year) population counts that uncovered a rapid decline in spatial correlations in abundance with distance as short as ~2 km. These findings were contrasted with data for two other fish species that have a pelagic larval stage (corkwing wrasse, Symphodus melops; black goby, Gobius niger). For these latter species, we found wider spatial scales of connectivity and weaker genetic isolation‐by‐distance patterns, except where both species experienced a strong barrier to gene flow, seemingly due to lack of suitable habitat. Our findings verify expectations that a fragmented habitat and absence of a pelagic larval stage promote genetic structure, while presence of a pelagic larvae stage increases demographic connectivity and gene flow, except perhaps over extensive habitat gaps.

Site fidelity and shallow genetic structure in the common smooth-hound shark Mustelus mustelus confirmed by tag-recapture and genetic data

The common smooth-hound shark, Mustelus mustelus, is a widely distributed demersal shark under heavy exploitation from various fisheries throughout its distribution range. To assist in the development of appropriate management strategies, the authors evaluate stock structure, site fidelity and movement patterns along the species' distribution in southern Africa based on a combination of molecular and long-term tag-recapture data. Eight species-specific microsatellite markers (N = 73) and two mitochondrial genes, nicotinamide adenine dehydrogenase subunit 4 and control region (N = 45), did not reveal any significant genetic structure among neighbouring sites. Nonetheless, tagging data demonstrate a remarkable degree of site fidelity with 76% of sharks recaptured within 50 km of the original tagging location. On a larger geographic scale, dispersal is governed by oceanographic features as demonstrated by the lack of movements across the Benguela-Agulhas transition zone separating the South-East Atlantic Ocean (SEAO) and South-West Indian Ocean (SWIO) populations. Microsatellite data supported very shallow ocean-based structure (SEAO and SWIO) and historical southward gene flow following the Agulhas Current, corroborating the influence of this dynamic oceanographic system on gene flow. Moreover, no movements between Namibia and South Africa were observed, indicating that the Lüderitz upwelling formation off the Namibian coast acts as another barrier to dispersal and gene flow. Overall, these results show that dispersal and stock structure of M. mustelus are governed by a combination of behavioural traits and oceanographic features such as steep temperature gradients, currents and upwelling systems.

Range-Wide Population Structure of 3 Deepwater Eteline Snappers Across the Indo-Pacific Basin

Deep-sea habitats may drive unique dispersal and demographic patterns for fishes, but population genetic analyses to address these questions have rarely been conducted for fishes in these environments. This study investigates the population structure of 3 tropical deepwater snappers of the genus Etelis that reside at 100-400 m depth, with broad and overlapping distributions in the Indo-Pacific. Previous studies showed little population structure within the Hawaiian Archipelago for 2 of these species: Etelis coruscans and E. carbunculus. Here we extend sampling to the entire geographic range of each species to resolve the population genetic architecture for these 2 species, as well as a recently exposed cryptic species (Etelis sp.). One goal was to determine whether deepwater snappers are more dispersive than shallow-water fishes. A second goal was to determine whether submesophotic fishes have older, more stable populations than shallow reef denizens that are subject to glacial sea-level fluctuations. Both goals are pertinent to the management of these valuable food fishes. A total of 1153 specimens of E. coruscans from 15 geographic regions were analyzed, along with 1064 specimens of E. carbunculus from 11 regions, and 590 specimens of E. sp. from 16 regions. The first 2 species were analyzed with mtDNA and 9-11 microsatellite loci, while E. sp. was analyzed with mtDNA only. Etelis coruscans had a non-significant microsatellite global FST, but significant global mtDNA Ф ST = 0.010 (P = 0.0007), with the isolation of Seychelles in the western Indian Ocean, and intermittent signals of isolation for the Hawaiian Archipelago. Etelis carbunculus had a non-significant microsatellite global FST, and significant global mtDNA Ф ST = 0.021 (P = 0.0001), with low but significant levels of isolation for Hawai'i, and divergence between Tonga and Fiji. Etelis sp. had mtDNA Ф ST = 0.018 (P = 0.0005), with a strong pattern of isolation for both Seychelles and Tonga. Overall, we observed low population structure, shallow mtDNA coalescence (similar to near-shore species), and isolation at the fringes of the Indo-Pacific basin in Hawai'i and the western Indian Ocean. While most shallow-water species have population structure on the scale of biogeographic provinces, deepwater snapper populations are structured on the wider scale of ocean basins, more similar to pelagic fishes than to shallow-water species. This population structure indicates the capacity for widespread dispersal throughout the Indo-Pacific region.

Phylogeography, Population Structure, and Species Delimitation in Rockhopper Penguins (Eudyptes chrysocome and Eudyptes moseleyi)

Rockhopper penguins are delimited as 2 species, the northern rockhopper (Eudyptes moseleyi) and the southern rockhopper (Eudyptes chrysocome), with the latter comprising 2 subspecies, the western rockhopper (Eudyptes chrysocome chrysocome) and the eastern rockhopper (Eudyptes chrysocome filholi). We conducted a phylogeographic study using multilocus data from 114 individuals sampled across 12 colonies from the entire range of the northern/southern rockhopper complex to assess potential population structure, gene flow, and species limits. Bayesian and likelihood methods with nuclear and mitochondrial DNA, including model testing and heuristic approaches, support E. moseleyi and E. chrysocome as distinct species lineages with a divergence time of 0.97 Ma. However, these analyses also indicated the presence of gene flow between these species. Among southern rockhopper subspecies, we found evidence of significant gene flow and heuristic approaches to species delimitation based on the genealogical diversity index failed to delimit them as species. The best-supported population models for the southern rockhoppers were those where E. c. chrysocome and E. c. filholi were combined into a single lineage or 2 lineages with bidirectional gene flow. Additionally, we found that E. c. filholi has the highest effective population size while E. c. chrysocome showed similar effective population size to that of the endangered E. moseleyi. We suggest that the current taxonomic definitions within rockhopper penguins be upheld and that E. chrysocome populations, all found south of the subtropical front, should be treated as a single taxon with distinct management units for E. c. chrysocome and E. c. filholi.

Deep phylogeographic structure may indicate cryptic species within the Sparid genus Spondyliosoma

Two geographically nonoverlapping species are currently described within the sparid genus Spondyliosoma: Spondyliosoma cantharus (Black Seabream) occurring across Mediterranean and eastern Atlantic waters from NW Europe to Angola and S. emarginatum (Steentjie) considered endemic to southern Africa. To address prominent knowledge gaps this study investigated range-wide phylogeographic structure across both species. Mitochondrial DNA sequences revealed deep phylogeographic structuring with four regionally partitioned reciprocally monophyletic clades, a Mediterranean clade and three more closely related Atlantic clades [NE Atlantic, Angola and South Africa (corresponding to S. emarginatum)]. Divergence and distribution of the lineages reflects survival in, and expansion from, disjunct glacial refuge areas. Cytonuclear differentiation of S. emarginatum supports its validity as a distinct species endemic to South African waters. However, the results also indicate that S. cantharus may be a cryptic species complex wherein the various regional lineages represent established/incipient species. A robust multilocus genetic assessment combining morphological data and detailing interactions among lineages is needed to determine the full diversity within Spondyliosoma and the most adequate biological and taxonomic status.

Highly regional population structure of Spondyliosoma cantharus depicted by nuclear and mitochondrial DNA data

Resolution of population structure represents an effective way to define biological stocks and inform efficient fisheries management. In the present study, the phylogeography of the protogynous sparid Spondyliosoma cantharus, in the East Atlantic and Mediterranean Sea, was investigated with nuclear (S7) and mitochondrial (cytochrome b) DNA markers. Significant divergence of four regional genetic groups was observed: North Eastern Atlantic, Mediterranean Sea, Western African Transition (Cape Verde) and Gulf of Guinea (Angola). The two southern populations (Cape Verde and Angola) each comprised reciprocally monophyletic mtDNA lineages, revealed low levels of diversity in Cape Verde and high diversity for Angola despite being represented by only 14 individuals. A complete divergence between North Atlantic and Mediterranean populations was depicted by the mitochondrial marker, but a highly shared nuclear haplotype revealed an incomplete lineage sorting between these regions. Bayesian skyline plots and associated statistics revealed different dynamics among the four regions. Cape Verde showed no expansion and the expansion time estimated for Angola was much older than for the other regions. Mediterranean region seems to have experienced an early population growth but has remained with a stable population size for the last 30000 years while the North Atlantic population has been steadily growing. The lack of genetic structuring within these regions should not be taken as evidence of demographic panmixia in light of potential resolution thresholds and previous evidence of intra-regional phenotypic heterogeneity.

Genetic analysis provides insights into species distribution and population structure in East Atlantic horse mackerel (Trachurus trachurus and T. capensis)

Two sister species of horse mackerel (Trachurus trachurus and T. capensis) are described that are intensively harvested in East Atlantic waters. To address long-standing uncertainties as to their respective geographical ranges, overlap and intraspecific population structure this study combined genetic (mitochondrial DNA and microsatellite) analysis and targeted sampling of the hitherto understudied West African coast. mtDNA revealed two reciprocally monophyletic clades corresponding to each species with interspecies nuclear differentiation supported by FST values. The T. trachurus clade was found across the north-east Atlantic down to Ghana but was absent from Angolan and South African samples. The T. capensis clade was found only in South Africa, Angola and a single Ghanaian individual. This pattern suggests that both species may overlap in the waters around Ghana. The potential for cryptic hybridization and/or indiscriminate harvesting of both species in the region is discussed. For T. capensis mtDNA supports high gene flow across the Benguela upwelling system, which fits with the species' ecology. The data add to evidence of a lack of significant genetic structure throughout the range of T. trachurus though the assumption of demographic panmixia is cautioned against. For both species, resolution of stock recruitment heterogeneity relevant to fishery management, as well as potential hybridization, will require more powerful genomic analyses.

Population panmixia and demographic expansion of a highly piscivorous marine fish Scomberomorus niphonius

Population structure and demographic history of the Japanese Spanish mackerel Scomberomorus niphonius a highly piscivorous and migratory marine fish, were assessed using mitochondrial DNA control region sequences (n = 720) and microsatellite genotypes at five loci (n = 1331) for samples collected on Japanese coasts from 2001 to 2010. The population structure was panmictic and the haplotype and allele frequencies were temporally stable even during the recent recovery process. Demographic expansion was strongly supported throughout the Pleistocene, suggesting that the oscillating glacial and interglacial climate conditions in the Pleistocene had no substantial impact on the demographic history of S. niphonius.

Mechanisms of peripheral phylogeographic divergence in the indo-Pacific: lessons from the spiny lobster Panulirus homarus

Background

There is increasing recognition of the concordance between marine biogeographic and phylogeographic boundaries. However, it is still unclear how population-level divergence translates into species-level divergence, and what are the principal factors that first initiate that divergence, and then maintain reproductive isolation. This study examines the likely forces driving population and lineage divergences in the broadly-distributed Indo-Pacific spiny lobster Panulirus homarus, which has peripheral divergent lineages in the west and east. The study focuses particularly on the West Indian Ocean, which is emerging as a region of unexpected diversity. Mitochondrial control region (mtCR) and COI sequences as well as genotypes of 9 microsatellite loci were examined in 410 individuals from 17 locations grouped into 7 regions from South Africa in the west, and eastward across to Taiwan and the Marquesas Islands. Phylogenetic and population-level analyses were used to test the significance and timing of divergences and describe the genetic relationships among populations.

Results

Analyses of the mtCR revealed high levels of divergence among the seven regions (Ф_ST = 0.594, P < 0.001). Microsatellite analyses also revealed significant divergence among regions, but at a much lower level (F_ST = 0.066, P < 0.001). The results reveal different patterns of mtCR v. nDNA divergence between the two distinct peripheral lineages: a subspecies in South Africa and Madagascar, and a phylogeographically diverged population in the Marquesas. The results also expose a number of other more fine-scale population divergences, particularly in the Indian Ocean.

Conclusions

The divergence of peripheral lineages in the west and east of the species’ range appear to have been initiated and maintained by very different processes. The pattern of mitochondrial and nuclear divergence of the western lineage, implicates processes of parapatric isolation, secondary contact and introgression, and suggests possible maintenance through adaptation and behavioural reproductive isolation. In contrast, the eastern lineage appears to have diverged through a rare colonisation event, maintained through long-term isolation, and matches expectations of the core-periphery hypothesis. The process of active peripheral speciation may be a common force in the Indo-Pacific that helps drive some of the regions’ recognized biogeographic boundaries.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-1050-8) contains supplementary material, which is available to authorized users.

Morphological and genetic divergence between lake and river populations of Triplophysa in Ngangtse Co, Tibet

Ngangtse Co (4535 m a.s.l., 840 km²) is a typical endorheic lake in Tibet. In 2015, we investigated the fish resource in this lake and its rivers and collected two Triplophysa species, including T. brevicauda and T. stewarti. To understand the evolutionary dynamics of Triplophysa in Ngangtse Co, the adaptive and genetic divergence between river and lake populations were investigated by analysing their morphological characteristics and mitochondrial COI and Cytb sequences. Among all of 277 samples, the ratios of the number of T. brevicauda to T. stewarti were 27:8, 34:11, 14:24, 37:0, 16:23 and 0:83 in river populations R_A, R_B, R_C, R_D, lakeside population LS_E and lake population L, respectively. For T. brevicauda, a small but significant genetic divergence (mean F_ST = 0.0890) was detected between population R_D and other three river populations, and significant morphological differences of body length/body height and body length/caudal peduncle length were observed between population R_D and R_A, R_B, suggesting an association between differentiation and geographical distance. For T. stewarti, very high levels of genetic differentiation was observed between population L and other populations with the F_ST values ranging from 0.4737 (L-R_A) to 0.7074 (L-R_C). In addition, significant differences of most of morphological characteristics were observed between the lake (L) and river (R_C) populations, indicating that lake fish of T. stewarti was a population that showed a long-term adaptation to the saltwater lake. These results provide insights on the evolutionary dynamics of Triplophysa in Ngangtse Co and help us to understand the fish speciation in endorheic lakes on the Qinghai-Tibetan Plateau.