Skip to main content
  • Research article
  • Open access
  • Published:

Genome-wide SNP analyses reveal high gene flow and signatures of local adaptation among the scalloped spiny lobster (Panulirus homarus) along the Omani coastline

Abstract

Background

The scalloped spiny lobster (Panulirus homarus) is a popular seafood commodity worldwide and an important export item from Oman. Annual catches in commercial fisheries are in serious decline, which has resulted in calls for the development of an integrated stock management approach. In Oman, the scalloped spiny lobster is currently treated as a single management unit (MU) or stock and there is an absence of information on the genetic population structure of the species that can inform management decisions, particularly at a fine-scale level. This work is the first to identify genome-wide single nucleotide polymorphisms (SNPs) for P. homarus using Diversity Arrays Technology sequencing (DArT-seq) and to elucidate any stock structure in the species.

Results

After stringent filtering, 7988 high utility SNPs were discovered and used to assess the genetic diversity, connectivity and structure of P. homarus populations from Al Ashkharah, Masirah Island, Duqm, Ras Madrakah, Haitam, Ashuwaymiyah, Mirbat and Dhalkut landing sites. Pairwise FST estimates revealed low differentiation among populations (pairwise FST range = − 0.0008 - 0.0021). Analysis of genetic variation using putatively directional FST outliers (504 SNPs) revealed higher and significant pairwise differentiation (p < 0.01) for all locations, with Ashuwaymiyah being the most diverged population (Ashuwaymiyah pairwise FST range = 0.0288–0.0736). Analysis of population structure using Discriminant Analysis of Principal Components (DAPC) revealed a broad admixture among P. homarus, however, Ashuwaymiyah stock appeared to be potentially under local adaptive pressures. Fine scale analysis using Netview R provided further support for the general admixture of P. homarus.

Conclusions

Findings here suggested that stocks of P. homarus along the Omani coastline are admixed. Yet, fishery managers need to treat the lobster stock from Ashuwaymiyah with caution as it might be subject to local adaptive pressures. We emphasize further study with larger number of samples to confirm the genetic status of the Ashuwaymiyah stock. The approach utilised in this study has high transferability in conservation and management of other marine stocks with similar biological and ecological attributes.

Background

Severe decline of many commercial fish stocks highlights the emerging need for sustainable management plans for regulation and conservation of marine biodiversity. Managing marine stocks sustainably is a dynamic process and requires an in-depth understanding of the stock and its spatial boundaries, along with biological, ecological, evolutionary, economic, social or even political factors that influence the fishery [1, 2]. While traditional fishery management plans rely on morphological and demographic aspects of a population such as growth, size, and mortality rates [3, 4], appropriate management should also consider evolutionary criteria, including conservation of genetic diversity and maintenance of sustainable spawning stock biomass [2]. Recent studies have shown the complexity in the population genetic structure of many marine species [5,6,7]. Generally, marine organisms possess high genetic diversity and show weak population differentiation due to highly dispersive larval stages and relative absence of barriers to dispersal in the marine environment [8,9,10]. However, seascape factors (e.g. water currents, seafloor features and bathymetry) and environmental attributes can significantly influence rates of gene flow, connectivity and genetic structure in some species. Further, evolutionary processes like genetic drift and selection [11,12,13] continuously shape the genomes of marine organisms. For these reasons, defining the population structure of such organisms is challenging, but important for their conservation and management [14,15,16]. Current progress in the fields of genomics and computational biology doubtlessly offers a versatile platform for fishery managers to answer questions and issues related to population structure, stock boundaries and the level of divergence of marine organisms [1, 16, 17]. Recent reports support the successful application of genomic approaches to identify conservation or management units (MUs) of marine species [8, 18, 19]. Many of these utilise advanced genomic approaches, using high-throughput genotyping technologies i.e. Next generation sequencing (e.g. Illumina HiSeq and MiSeq platforms) and third generation sequencing (e.g. PacBio and Nanopore technology), to isolate large number of genetic markers suitable for inference of population differentiation and structure [20,21,22,23]. These technologies have enabled the development of panels of SNP markers to investigate interspecific hybrids [24, 25], to assign individuals to populations, or to identify MUs [8, 11, 26]. Harnessing of genomic wide SNPs in the assessment of commercial marine stocks is a successful approach that could address many questions related to the level of genetic diversity and stability of these stocks [27].

The scalloped spiny lobster P. homarus (Linnaeus 1758) is characterised by a relatively long pelagic larval duration (PLD) of about 4.5–6.5 months [28], during which the larvae is exposed to oceanic dispersal as a result of currents and wind-shear, before metamorphosing into the puerulus stage and continuing its life as a benthic organism [28, 29]. The species is distributed throughout the Indo-Pacific [30] and in the region supports valuable fisheries of considerable socio-economic importance. There are major concerns about the future of spiny lobster fisheries owing to a general decline in catch [31, 32], emphasizing the need for serious efforts towards sustainable fishery management and regulation of the species. It is essential to introduce comprehensive fishery management guidelines for the species, considering a wide range of biological aspects e.g. demographic interactions of individuals and genetic structuring. Characterisation of stock boundaries and identification of population divergence will greatly support managers in deciding whether two populations should be managed together, or as separate stocks [27, 33]. Many recent works of P. homarus have studied sub-species resolution, phylogeography throughout its wide range [34,35,36,37] or dispersal capabilities [38]. However, the fine-scale genetic structure of the species in many regions is still remains unrevealed [38].

Commercial spiny lobster fisheries have a long tradition in Oman, with the country currently being one of the major suppliers to the global market [39]. Of concern, however, is the observation that the annual harvest of Omani lobsters has declined dramatically from 2000 tons/year in the 1980’s, to less than 485 tons in 2016 [40]. Presently, the lobster fishery management in Oman is primarily based on data from growth, mortality and catch rates and aims to increase population densities [39, 41]. Despite regulations implemented by the government, they are not regularly reinforced (i.e. high incidence of illegal catch) and no clear legislation system against illegal practices [42]. There is a lack of knowledge surrounding the population genetic structure of P. homarus in Oman, its levels of genetic fitness and relatedness in this region. Hence, the lobster stock along the Omani coastline is currently treated as one single management unit. This study is the first to assess the genetic structure of P. homarus in Oman using high resolution genome-wide SNPs genotyping. The findings provide valuable insight into the connectivity of the Omani P. homarus population and will aid in the identification of management units for the fishery of this commercially important crustacean.

Results

SNPs quality control and filtering

A primary dataset of 48,140 SNPs was filtered to retain 7988 SNPs (Additional file 1) suitable for genomic analysis (Table 1). Significant deviations from Hardy Weinberg Equilibrium (HWE) were observed across all populations (p < 0.000004 after Bonferroni correction). After the removal of these SNP loci deviating from HWE, significant skew in estimation of genetic diversity indices (Fis and Ho) was observed. Thus, indicating the status of those SNPs as putative null alleles or genotypic artifacts. Additionally, nine individuals were excluded from the dataset due to poor genotyping coverage < 80%.

Table 1 Filtering steps and SNPs counts retained after each step

Population genetic diversity

Observed heterozygosity (Ho) across populations ranged from 0.1660 to 0.1840 and were generally lower than the expected heterozygosity values (Hn.b) (0.2260–0.2333, Table 2). Average individual multilocus heterozygosity (Av.MLH) revealed similar values and distribution to Ho, ranging from 0.1683 to 0.1858 (for Ashuwaymiyah and Haitam, respectively) (Table 2). Average standardized MLH (sMLH) values slightly varied across populations and ranged from 0.9302 to 1.029. Inbreeding coefficient (Fis) was significantly high across all populations, (0.2094–0.2861, Table 2). The estimated parameters of identity disequilibrium (g2) slightly differentiated from zero (0.0003–0.0038). However, this differentiation was statistically significant (i.e. 95% C. I. does not overlap zero) in two locations (Haitam and Dhalkut, Table 2). Estimated effective population size (NeLD) varied from 5507.5 for Haitam and 10,305.8 for Al Ashkharah to an infinite value for other populations.

Table 2 Genetic diversity indices for P. homarus for each sampling site using 7988 SNPs

Population differentiation and genetic structure

In general, pairwise genetic differentiation estimates (FST) using 7988 SNPs indicated very low levels of genetic differentiation, with average FST = 0.0004 (±SD = 0.1843), with only seven out of 28 pairwise comparisons being statistically significant (Table 3). AMOVA indicated an absence of hierarchical genetic structure between populations (variation of 80.47% within individuals; 19.51% among individuals; 0.01% among populations and 0% among groups (AS, MA; group1; DU, RM, group2; HA, SH, group3; MI, DA; group4)). While visualization of population structure using DAPC with all 7988 SNPs revealed two admixed genetic clusters (Additional file 2), population Network analysis with Netview R displayed only one genetic cluster (Fig. 1a, b). Similarly, high admixture was observed through the NJ tree, indicating high genetic relatedness among individuals from different geographical locations (similar branch lengths among all individuals, Fig. 1c).

Table 3 Pairwise FST values for 7988 SNPs using Genetix v.4.05.2 with permuted p-values inside brackets
Fig. 1
figure 1

Population structure of 164 individuals of P. homarus samples using 7988 SNPs. Population network was constructed using NetView R v.1.0. The network is visualised at nearest neighbor (k-NN) threshold of 25 (a) and 50 (b). Un-rooted neighbor-joining tree (c) was drawn in MEGA6 using 1-psa genetic distances, scale bar indicates genetic distance

Putatively selective SNPs

The genomic scan using Bayescan v.2.1 identified only one SNP as a directional outlier, which was characterized by Jeffery’s scale as decisive. The frequentist approach with Arlequin v.3.5.2.2 indicated out of 7988 SNPs, 504 SNPs (6.31%) potentially under divergent selection, 168 SNPs (2.10%) under balancing selection and 7316 SNPs (91.59%) likely to be neutral (Genotypes of 504 SNP loci are available as Additional file 3). Global FST for the 504 SNPs candidate outliers was 0.0423 (Table 4), more than 100 times greater than for all 7988 SNPs (0.0004). All of the pairwise estimates with these directional outliers were significant (p < 0.001), with highest differentiation between Ashuwaymiyah and Ras Madrakah (0.0736) and the lowest between Masirah and Dhalkut (0.0288) (Fig. 2; Table 4). Visualization of DAPC revealed that Ashuwaymiyah samples represented a distinct genetic cluster, while all other samples were grouped into a second admixed genetic cluster (Fig. 3). The same pattern was also observed by population Netview at a range from k-NN =10 to 30 and visualized at k-NN = 15 (Fig. 4).

Table 4 Pairwise FST values for 504 directional SNPs using Genetix v.4.05.2
Fig. 2
figure 2

Sampling sites of P. homarus specimens along the coastline of Oman. AS Al Ashkharah, MA Masirah Island, DU Duqm, RM Ras Madrakah, HA Haitam, SH Ashuwaymiyah, MI Mirbat and DA Dhalkut. Numbers in circles represent genotyped samples in the final analyses. The black circle represents Ras Al-Had, the northern border of the commercial fishery sites of P. homarus. The map was obtained from: (https://www.d-maps.com/carte.php?num_car=516&lang=en) and edited to highlight the sampling sites

Fig. 3
figure 3

DAPC scatter plot created using eight optimum PCs out of 504 outlier SNPs across 164 P. homarus individuals in the R package adegenet. Dots represent individuals. The plot showing Ashuwaymiyah represents a distinct genetic cluster and other sampling locations represent a second admixed genetic cluster

Fig. 4
figure 4

Population network of P. homarus individuals using NetView R v.1.0 based on identity by similarity (IBS) distance matrix calculated in PLINK after Steinig et al. [104]. The network was drawn at nearest neighbor (k-NN) threshold of 15, using 504 potentially directional SNPs and 164 individuals

Discussion

Population genetic diversity

Assessment of genetic diversity and population structure of P. homarus along the Omani coastline is vital for its stock management and conservation. This study is the first to identify a genome-wide SNP dataset for P. homarus utilising DArTseq technology and to determine the genetic stock structure of the species in Oman. The obtained SNP dataset revealed unique insights into genetic diversity and stocks boundaries of the spiny lobster in Oman. Significant deviations from Hardy Weinberg Equilibrium (HWE) were observed across all populations, apparently caused by a heterozygote deficit as represented by significantly high Fis values (all Fis values > 0, p < 0.001). Despite repetitive culling of SNPs deviating from HWE, notably lower, but still significantly high Fis values (0.2094–0.2861) were observed. Heterozygote deficits have been observed in many marine invertebrates [43] and are believed to be caused by genotypic artifacts, null alleles, population stratification caused through the Wahlund effect, biological and behavioural traits and selection [44,45,46]. In other cases, genotype calling errors have been reported to cause deviations from HWE [47,48,49]. To overcome these in our dataset, SNP loci which had call rate (< 0.7) and those significantly deviated from HWE or had low minor allele frequencies (MAF) were removed [48]. In this study, estimations of identity disequilibrium (g2) slightly differentiated from zero (0.0003–0.0038, Table 2) with no significance except for Haitam and Dhalkut, thus, the obtained high Fis values for P. homarus samples in this study are less likely due to inbreeding.

Additionally, average multi locus heterozygosity estimates (Av.MLH) showed a moderate to high level of genomic heterozygosity (0.1683–0.1858) compared to some marine invertebrates e.g. Hapalochlaena maculosa (0.0800–0.1720) [50], Pinctada margaritifera (0.0520–0.1030) [11] and P. maxima (0.3030–0.3110) [51]. This may support that the significantly high Fis values are unlikely due to reduction in genomic heterozygosity. Multi locus heterozygosity may reflect genome wide heterozygosity especially when thousands of genomic markers are used [52, 53].

It is also unlikely that the high Fis values are caused by the Wahlund effect, as Netview analyses indicated high levels of admixture (Fig. 1a, b). Possibly, the observed heterozygote deficits are caused by null alleles or selection. In our study, high missingness rate (> 0.2) of SNPs was observed for some samples from Ashuwaymiyah and were excluded, as they did not match our quality criteria. In the remaining samples (with genotyping coverage ≥0.8), monomorphic loci were observed for Ashuwaymiyah only, but not other sites. The monomorphic loci observed in Ashuwaymiyah samples could be attributed to mutations in restriction enzyme sites during DArTseq genotyping causing null alleles in populations [54]. It is also possible, that these observations were due to small sample size effect from Ashuwaymiyah. A study in our group with microsatellite markers revealed divergence of Ashuawymiyah stock in a larger sample size of 20 individuals (Delghandi et al., under prep).

Generally, null alleles may indicate occurrence of genetic variation in the form of point mutations or structural mutations (i.e. insertions/deletions) and contribute to the organism’s fitness and adaptation [55]. It has been reported that the frequency of missing data is correlated with the level of genetic divergence between populations [54, 56]. Other studies, however, have reported that null alleles are frequently encountered in SNP datasets and generate biases in estimation of diversity indices [57].

Population differentiation and genetic structure

In our study, pairwise differentiation estimates (FST) analysis using 7988 SNP loci, showed low but statistically significant (p < 0.05) differentiation between Ashuwaymiyah and other locations, except for Al Ashkharah and Mirbat (Table 3). Earlier reports showed that P. homarus possesses a genetic structure at the large Indo-Pacific scale [36] and relatively small scales, within both the north west Indian Ocean and south west Indian Ocean [38, 58].

Other studies demonstrated different patterns and levels of genetic structure in spiny lobsters e.g. P. ornatus, exhibit low levels of genetic structure across an expansive distribution of the Indo-Pacific [59], but not at smaller scales i.e. South East Asia [60]. Similarly, no genetic structure was observed in Hawaiian P. penicillatus [61], while at the scale of the Indian Ocean and across its Indo-Pacific range, significant structure was revealed [62, 63]. These contrasting patterns are referred to, mainly, environmental (e.g. pattern of water circulation) [38], bioecological (e.g. larval retention) [64], behavioral (e.g. spawning migrations) [60], or geographical factors (e.g. habitats patchiness) [59]. In this study, the observed genetic divergence among Omani P. homarus could be caused by distinct environmental and geographical factors in the region. The Omani coast in the Arabian sea is known to be influenced by complex water circulation which varies seasonally with the Monsoon and results in a series of eddies along the coast of Oman [65, 66]. Generally, eddies might act as a larval retention system [67], limiting larval dispersal and maintaining divergence in marine populations including spiny lobsters, e.g. Jasus edwardsii [68] and P. h. rubellus [38]. Additionally, fragmentation of marine habitats across the Omani coastline by sandy stretches and absence of corals and rock reefs [44] could have contributed to this observation. Similar observations have been reported for other marine organisms with larval life stage such as Corkwing wrasse across Norwegian coastline [69] and the Omani clownfish [70]. In fact, Ashuwaymiyah is a shallow bay characterized by large rocky reefs, and its extended sea shelf (about 50 km) is less affected by the Monsoon currents, even during its extremes from June to September (National Survey Authority, Ministry of Defense, Oman, unpublished local data). This unique geography of Ashuwaymiyah could have limited the gene flow among lobsters from Ashuwaymiyah and other sites. This could be a possible explanation for the putative genetic divergence of P. homarus samples from Ashuwaymiyah. Moreover, it is not surprising to observe a genetically heterogeneous stock among other admixed stocks over a relatively small geographical scale. Such observation of fine scale genetic differentiation within relatively high gene flow environment has been widely described in a variety of marine species with planktonic early life stage and in a phenomenon known as chaotic genetic patchiness [71, 72]. Examples include crown-of-thorns starfish, Acanthaster planci [73], clam, Spisula ovalis, [74], the sea urchin, Strongylocentrotus purpuratus [75], marine goby, Coryphopterus personatus [76], marine goby, pulmonate limpets, Siphonaria sp. [71], bicolour damselfish, Stegastes partitus [77] and spiny lobster, P. interruptus [78]. Therefore, it is possible that the observed slight divergence of Ashuwaymiyah is due to chaotic genetic patchiness. A typical feature of chaotic genetic patchiness is being temporal therefore, repetitive sampling from Ashuwaymiyah would be useful to clarify the current status of the genetic heterogeneity in Ashuwaymiyah.

Detecting potential selective SNPs

Genome-wide scan for FST outliers using Bayescan v2.1 identified only one directional outlier, hence further analysis of population structure was not possible using this approach. In contrast, applying Arlequin resulted in more candidate outlier loci being identified, allowing further investigation for adaptive variation. It is common to detect less candidate outliers with Bayescan as it is a conservative approach and may fail to detect relatively low signals of selection [79,80,81]. Therefore other studies have used only the frequentist approach to perform FST outliers analysis [82].

Assessment of population structure with putatively directional SNPs using both population Network and DAPC, revealed Ashuwaymiyah to be potentially under local adaptive pressures (Figs. 3 and 4). The utilization of the whole dataset of SNPs could not capture the low levels of genetic structure, while analyses based on FST outliers allowed detection of selective divergence and identification of possible discrete stock. Similarly, many other studies showed that the use of FST outliers could detect adaptive variation in the absence of broader analyses based on neutral markers [83,84,85]. A possible explanation for the observed divergence is the heterogeneous environmental attributes i.e. massive rocky reefs and shallow waters in Ashuwaymiyah, which might be the driver of this differentiation. In addition to the morphological/biological differences among P. homarus stocks [42], genetic studies with microsatellite markers revealed a significant divergence of lobsters from the Dhofar governorate, including Ashuwamiyah from Al Sharqiyah and Al Wusta governorates (Delghandi et al., under prep).

Implications for fishery management

The population of P. homarus along Oman is currently considered as a single homogenous stock with a single management regulation. This study shows that the samples from Ashuwaymiyah are genetically distinct from other broadly admixed samples, albeit at low levels. An earlier study reported that stock of P. homarus in Ashuwaymiyah was significantly differentiated in body size and that lobsters reach maturity at significantly lower sizes when compared to two geographically close locations [86]. The same study suggested a need to investigate the current fishery management further and indicated that Ashuwaymiyah site might require separate management. Additionally, a recent biological study of P. homarus in Oman revealed that the stock in Dhofar governorate differs from other stocks in size, time of spawning and number of spawning peaks/year, suggesting a need to consider spatial management of P. homarus along the Omani coastline [42]. Our study delivers for the first time genetic support for possible differentiation of the Ashuwamiyah stock from other locations across the coastline of Oman, using genome-wide markers, and that this stock might need to be considered for regional management. We recommend the conduct of further studies with larger number of samples, coupled with environmental and ecological data, to aid integrated assessment studies and potential discovery of unique management units.

Conclusions

Utilisation of genome wide SNPs to study the genetic status of P. homarus stocks in Oman provided valuable insights into the genetic status of the stock. This genomic resource is the first of its kind in P. homarus and the SNP dataset obtained in this study has allowed deep characterization of the lobster population genetic diversity, connectivity and structure in Oman. This study has revealed general admixture and high connectivity of P. homarus across the Omani coastline. Additionally, the study highlighted the potential prevalence of local adaptive pressures in Ashuwaymiyah. These findings indicate the importance of considering spatially customized management strategies for P. homarus across the coastline of Oman. Further studies of the genetic status of Ashuwaymiyah and stocks from other locations in Dhofar with adequate sampling based on different temporal periods together with ecological and environmental data about Omani coastline is required before any conclusive decision on the stock structure can be inferred.

Methods

Sampling and genomic DNA extraction

The commercial P. homarus fishery in Oman is situated between Ras Al-Hadd and Dhalkut (a distance of approximately 1100 km) (Fig. 2). Samples were obtained from eight commercial landing-sites, covering most of the distribution range of the lobster in Oman (Fig. 2). All samples were euthanized and purchased from local fishermen in March 2015 during the legal fishing season. A single walking leg was excised from wild caught P. homarus (n = 172) and preserved immediately in 95% ethanol until DNA extraction. All samples were purchased from local fishermen in March 2015 during the legal fishing season. Genomic DNA was extracted from tissue samples using a modified cetyl trimethyl ammonium bromide (CTAB)/ Chloroform-Isoamyl method [87]. DNA extracts were further purified through a Sephadex G50 (GE, 2007) column prior to quantification with a Nanodrop 1000 Spectrophotometer (Thermo Scientific).

Library preparation and genotyping

Genomic DNA extracts were standardised to 50 ng/μl, and sent for sequencing and genotyping using DArTseq™ technology, with Diversity Arrays Technology, Canberra, Australia [88, 89]. Library preparation was completed as described by Kilian et al. [89] and Sansaloni et al. [90] with all P. homarus DNA samples being digested using a combination of PtsI and HpaII restriction enzymes. Multiplexed reduced representation libraries were then sequenced on the Illumina HiSeq2500 platform for 77 cycles.

To call SNPs and genotype each individual, raw Illumina HiSeq2500 data was first de-multiplexed into individual samples, based on sample-specific barcode sequences. De-multiplexed samples were then assessed for overall sequence quality, with any fragments with an average Q-score of < 25 being removed from the dataset. Sequences were also compared to public databases for identification of contaminant sequences, and any non-target sequences (including bacterial and viral fragments) were removed. SNP calling was conducted using the DArTsoft14 algorithm within the KDCompute framework developed by Diversity Arrays Technology (http://www.kddart.org/kdcompute.html), with initial calling parameters and filtering methods as described in Morse et al. [50] and Lind et al. [91].

SNPs quality control and filtering

To eliminate potentially aberrant SNPs, stringent quality controls were applied using custom python scripts within the DArTQC pipeline (https://github.com/esteinig/dartQC) [92]. Initially, all duplicated sequences with > 95% similarity were identified using CD-HIT and collapsed into a single cluster, or removed [93]. Further, SNPs with a call rate < 70% and those where technical replicates did not return a repeatability value of > 95% were also removed. Additionally, individuals and SNPs with > 20% missing data and SNPs with a Minor Allele Frequency (MAF) < 0.02 were excluded using Plink v1.07 [94].

To investigate the effect of sequencing depth, Fis and Ho were calculated for each population at different reads depth (Average SNP Counts) thresholds (3, 5, 7 and 10) to discover the degree of potential bias caused by lower call depths. Accordingly, four subsets of SNPs were generated at these sequencing depths. To detect potential genotyping artifacts, SNPs were tested for significant deviation from Hardy-Weinberg equilibrium (HWE) using Arlequin v.3.5.2.2 [95]. Any SNP loci which significantly deviated from HWE were excluded following Bonferroni correction (p < 0.000004). To assess the impact of deviation from HWE, Fis and Ho were calculated before and after removal of significantly deviated SNPs.

Population genetic diversity

To estimate the genetic diversity within populations, standard allelic diversity indices including average observed heterozygosity (Ho), average expected heterozygosity corrected for population sample size (Hn.b.) and inbreeding coefficient (Fis) were calculated using Genetix v.4.05.2 [96]. Effective population size, using a linkage disequilibrium method (NeLD) was computed with NeEstimator [97]. To examine individual genome wide diversity and individual inbreeding, multi-locus heterozygosity (MLH) and identity disequilibrium parameter (g2) were calculated for all individuals using the R package inbreedR [52].

Population differentiation and genetic structure

To assess population differentiation and genetic structure, a number of different statistical approaches were conducted. The extent of pairwise population differentiation was evaluated using Weir and Cockerham’s unbiased F-statistics [98] through Genetix v.4.05.2 [96]. To assess hierarchical levels of population structuring, an analysis of molecular variance (AMOVA) using Arlequin v.3.5.2.2 [99] was calculated between sampling locations based on grouping samples in four groups (AS, MA; group1; DU,RM, group2; HA, SH, group3; MI, DA; group4). The grouping criterion was based on habitat similarity i.e. abundance of corals and rock substrates. Obviously, Al Wusta coast is dominated by sandy habitat and sparse coral colonies [100, 101]. Region surrounding DU has been exposed for the last seven years to major construction and industrial influence, most probably having an impact on fragmentation of marine populations. Hence, HA being not affected of these factors, was grouped with SH, which is just below the border line between Al Wusta and Dhofar governates (Fig. 2). In addition, the function find.clusters in the R package adegenet [102] was used to determine the optimal number of clusters with the Bayesian Information Criterion (BIC) method. To assess levels of differentiation between the obtained genetic clusters, Discriminant Analysis of Principal Components (DAPC) was used. DAPC was performed using the optimum number of principal components (PCs) calculated using the α-score function in adegenet [103].

Finally, a network analysis with no prior population assumptions was performed to assess both broad and fine scale population structure using NetView R [104]. NetView was run through the R implementation of NetView P [104, 105] at a k-NN range from 25 to 50 as determined by a k-NN selection plot. Similarly, to visualise the extent of relatedness between individuals within each population and divergence among populations, a Neighbour-Joining (NJ) tree was constructed in MEGA6 [106]. The NJ tree was constructed using 1-proportion of shared alleles (1-psa) genetic distance matrix calculated in the R package adegenet using propShared function [102].

Identifying potential selective SNPs

To detect possible signatures of directional and balancing selection, detection of putatively selective outlier SNPs was performed using an FST approach. To minimize false positive rates in identifying SNPs under selection, two independent statistical approaches were used. A Bayesian approach was implemented in Bayescan v.2.1 [107] and a frequentist approach [108] was implemented in Arlequin v.3.5.2.2 [99]. Bayescan estimates population-specific FST coefficients by the Bayesian method described in [109] and uses a statistical cut-off based on the mode of the posterior distribution to detect SNPs under selection [107]. Bayescan v2.1 was used with 1:10 prior odds for a neutral model and all other parameters were kept as default (20 pilot runs of 5000 iterations followed by 100,000 iterations with an additional burn-in of 50,000) [110]. Once probabilities had been calculated for each locus, they were ranked from largest to smallest. SNPs with posterior probabilities ≥0.91–1, which are categorised as strong to decisive according to the Jeffery’s scale [110], were retained. In addition, the Bayescan v2.1 function, plot R.r in the R v.3.3.1 was used to control the false discovery rate (FDR) of the selective markers at FDR of 0.05. SNPs were considered as outliers if their probability was > 0.9 at FDR of 0.05.

The frequentist approach in Arlequin v.3.5.2.2 was executed under a finite island model with 200,000 simulations and 100 demes simulated [99]. SNPs were considered as outliers based on their FST and p values. SNPs were considered as directional loci if their FST values fell within the upper 5% quantile and p < 0.05. They were considered as balancing SNPs if their FST values fell in the lower 5% quantile and p < 0.05.

To assess the population structure based on directional outliers, a dataset of the putative outlier SNPs was generated. Broad scale population differentiation based on this SNP dataset was examined by calculating magnitude and significance of pairwise FST comparisons using Genetix v.4.05.2. Population structure and network were examined based on the putative outlier SNPs using DAPC and Netview respectively. DAPC was visualised after retaining the optimum number of PCs and NetView was run at k-NN range = 10 to 30.

Abbreviations

1-psa:

1-proportion of shared alleles

AMOVA:

Analysis of molecular variance

AS:

Al Ashkharah

BIC:

Bayesian information criterion

CTAB:

Cetyl trimethyl ammonium bromide

DA:

Dhalkut

DAPC:

Discriminant analysis of principal components

DartQC:

Diversity arrays technology SNPs quality control

DArT-seq:

Diversity arrays technology sequencing

Du:

Duqm

FDR:

False discovery rate

gDNA:

Genomic deoxyribonucleic acid

HA:

Haitam

HWE:

Hardy-Weinberg equilibrium

IBS:

Identity by similarity

k-NN:

Number of nearest neighbours k-threshold

Ma:

Masirah Island

MAF:

Minor allele frequency

MI:

Mirbat

MLH:

Multi-locus heterozygosity

NJ:

Neighbour-joining

PCR:

Polymerase chain reaction

PCs:

Principal components

PLD:

Pelagic larval duration

RM:

Ras Madrakah

SH:

Ashuwaymiyah

sMLH:

Standardized multi-locus heterozygosity

SNP:

Single nucleotide polymorphism

References

  1. Reiss H, Hoarau G. Ckey-Collas M, Wolff WJ. Genetic population structure of marine fish: mismatch between biological and fisheries management units. Fish Fish. 2009;10(4):361–95.

    Article  Google Scholar 

  2. Lal MM, Southgate PC, Jerry DR, Bosserelle C, Zenger KR. A parallel population genomic and hydrodynamic approach to fishery Management of Highly-Dispersive Marine Invertebrates: the case of the Fijian black-lip pearl oyster Pinctada margaritifera. PLoS One. 2016;11(8):e0161390.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  3. Carvalho GR, Hauser L. Molecular-genetics and the stock concept in fisheries. Rev Fish Biol Fish. 1994;4(3):326–50.

    Article  Google Scholar 

  4. McKeown NJ, Arkhipkin AI, Shaw PW. Regional genetic population structure and fine scale genetic cohesion in the southern blue whiting Micromesistius australis. Fish Res. 2017;185:176–84.

    Article  Google Scholar 

  5. Selkoe KA, D'Aloia CC, Crandall ED, Iacchei M, Liggins L, Puritz JB, et al. A decade of seascape genetics: contributions to basic and applied marine connectivity. Mar Ecol Prog Ser. 2016;554:1–19.

    Article  Google Scholar 

  6. Riginos C, Liggins L. Seascape genetics: populations, individuals, and genes marooned and adrift. Geogr Compass. 2013;7(3):197–216.

    Article  Google Scholar 

  7. Roy D, Hurlbut TR, Ruzzante DE. Biocomplexity in a demersal exploited fish, white hake (Urophycis tenuis): depth-related structure and inadequacy of current management approaches. Can J Fish Aquat Sci. 2012;69(3):415–29.

    Article  Google Scholar 

  8. Lal MM, Southgate PC, Jerry DR, Zenger KR. Fishing for divergence in a sea of connectivity: the utility of ddRADseq genotyping in a marine invertebrate, the black-lip pearl oyster Pinctada margaritifera. Mar Genomics. 2016;25:57–68.

    Article  PubMed  Google Scholar 

  9. Waples RS. Separating the wheat from the chaff: patterns of genetic differentiation in high gene flow species. J Hered. 1998;89(5):438–50.

    Article  Google Scholar 

  10. Waples RS, Gaggiotti O. What is a population? An empirical evaluation of some genetic methods for identifying the number of gene pools and their degree of connectivity. Mol Ecol. 2006;15(6):1419–39.

    Article  CAS  PubMed  Google Scholar 

  11. Lal MM, Southgate PC, Jerry DR, Bosserelle C, Zenger KR. Swept away: ocean currents and seascape features influence genetic structure across the 18,000 km indo-Pacific distribution of a marine invertebrate, the black-lip pearl oyster Pinctada margaritifera. BMC Genomics. 2017;18(1):66.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Limborg MT, Helyar SJ, de Bruyn M, Taylor MI, Nielsen EE, Ogden R, et al. Environmental selection on transcriptome-derived SNPs in a high gene flow marine fish, the Atlantic herring (Clupea harengus). Mol Ecol. 2012;21(15):3686–703.

    Article  CAS  PubMed  Google Scholar 

  13. Nosil P, Funk DJ, Ortiz-Barrientos D. Divergent selection and heterogeneous genomic divergence. Mol Ecol. 2009;18(3):375–402.

    Article  PubMed  Google Scholar 

  14. Andre C, Larsson LC, Laikre L, Bekkevold D, Brigham J, Carvalho GR, et al. Detecting population structure in a high gene-flow species, Atlantic herring (Clupea harengus): direct, simultaneous evaluation of neutral vs putatively selected loci. Heredity. 2011;106(2):270–80.

    Article  CAS  PubMed  Google Scholar 

  15. Benestan L, Gosselin T, Perrier C, Sainte-Marie B, Rochette R, Bernatchez L. RAD genotyping reveals fine-scale genetic structuring and provides powerful population assignment in a widely distributed marine species, the American lobster (Homarus americanus). Mol Ecol. 2015;24(13):3299–315.

    Article  PubMed  Google Scholar 

  16. Gaggiotti OE, Bekkevold D, Jorgensen HBH, Foll M, Carvalho GR, Andre C, et al. Disentangling the effects of evolutionary, demographic, and environmental factors influencing genetic structure of natural populations: Atlantic herring as a case study. Evolution. 2009;63(11):2939–51.

    Article  PubMed  Google Scholar 

  17. Garner BA, Hand BK, Amish SJ, Bernatchez L, Foster JT, Miller KM, et al. Genomics in conservation: case studies and bridging the gap between data and application. Trends Ecol Evol. 2016;31(2):81–3.

    Article  PubMed  Google Scholar 

  18. Larson WA, Seeb LW, Everett MV, Waples RK, Templin WD, Seeb JE. Genotyping by sequencing resolves shallow population structure to inform conservation of Chinook salmon ( Oncorhynchus tshawytscha). Evol Appl. 2014;7(3):355–69.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Benestan L, Quinn BK, Maaroufi H, Laporte M, Clark FK, Greenwood SJ, et al. Seascape genomics provides evidence for thermal adaptation and current-mediated population structure in American lobster (Homarus americanus). Mol Ecol. 2016;25(20):5073–92.

    Article  PubMed  Google Scholar 

  20. Aykanat T, Lindqvist M, Pritchard VL, Primmer CR. From population genomics to conservation and management: a workflow for targeted analysis of markers identified using genome-wide approaches in Atlantic salmon Salmo salar. J Fish Biol. 2016;89(6):2658–79.

    Article  CAS  PubMed  Google Scholar 

  21. Loman NJ, Misra RV, Dallman TJ, Constantinidou C, Gharbia SE, Wain J, et al. Performance comparison of benchtop high-throughput sequencing platforms. Nat Biotechnol. 2012;30(5):434–9.

    Article  CAS  PubMed  Google Scholar 

  22. Rusk N. Cheap third-generation sequencing. Nat Methods. 2009;6(4):244–5.

    Article  CAS  Google Scholar 

  23. Rhoads A, Au KF. PacBio sequencing and its applications. Geno Proteo Bioinform. 2015;13(5):278–89.

    Article  Google Scholar 

  24. Campbell NR, Amish SJ, Pritchard VL, McKelvey KS, Young MK, Schwartz MK, et al. Development and evaluation of 200 novel SNP assays for population genetic studies of westslope cutthroat trout and genetic identification of related taxa. Mol Ecol Resour. 2012;12(5):942–9.

    Article  CAS  PubMed  Google Scholar 

  25. Pujolar JM, Jacobsen MW, Als TD, Frydenberg J, Munch K, Jonsson B, et al. Genome-wide single-generation signatures of local selection in the panmictic European eel. Mol Ecol. 2014;23(10):2514–28.

    Article  CAS  PubMed  Google Scholar 

  26. Candy JR, Campbell NR, Grinnell MH, Beacham TD, Larson WA, Narum SR. Population differentiation determined from putative neutral and divergent adaptive genetic markers in eulachon (Thaleichthys pacificus, Osmeridae), an anadromous Pacific smelt. Mol Ecol Resour. 2015;15(6):1421–34.

    Article  PubMed  Google Scholar 

  27. Saha A, Hauser L, Kent M, Planque B, Neat F, Kirubakaran TG, et al. Seascape genetics of saithe (Pollachius virens) across the North Atlantic using single nucleotide polymorphisms. ICES J Mar Sci. 2015;72(9):2732–41.

    Article  Google Scholar 

  28. Phillips BAGR of SLA, Matsuda H. A Global Review of Spiny Lobster Aquaculture. In: Fotedar K, Phillips BF, editors. Recent Advances and New Species in Aquaculture. UK: Blackwell Publishing Ltd. p. 22–84.

    Chapter  Google Scholar 

  29. Booth JD, Phillips BF. Early-life history of spiny lobster. Crustaceana. 1994;66:271–94.

    Article  Google Scholar 

  30. Berry PF. A revision of the Panulirus Homarus-Group of Spiny Lobsters (Decapoda, Palinuridae). Crustaceana. 1974;27(1):31–42.

    Article  Google Scholar 

  31. Delghandi M, Goddard S, Jerry DR, Dao HT, Al Hinai MSN, Al-Amry W, et al. Novel genomic microsatellite markers for genetic population and diversity studies of tropical scalloped spiny lobster (Panulirus homarus) and their potential application in related Panulirus species. Genet Mol Res. 2016;15(2):gmr7846. https://doi.org/10.4238/gmr.15027846.

  32. Rogers PP, Barnard RM, Johnston MD. Lobster aquaculture a commercial reality: a review. J Mar Biol Ass Ind. 2010;52(2):327–35.

    Google Scholar 

  33. Fu CH, Fanning LP. Spatial considerations in the management of Atlantic cod off Nova Scotia. Canada North Amer J Fish Man. 2004;24(3):775–84.

    Article  Google Scholar 

  34. Lavery SD, Farhadi A, Farahmand H, Chan TY, Azhdehakoshpour A, Thakur V, et al. Evolutionary divergence of geographic subspecies within the scalloped spiny lobster Panulirus homarus (Linnaeus 1758). PLoS One. 2014;9(6):e97247.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  35. Singh SP, Groeneveld JC, Al-Marzouqi A, Willows-Munro S. A molecular phylogeny of the spiny lobster Panulirus homarus highlights a separately evolving lineage from the Southwest Indian Ocean. Peer J. 2017;5:e3356.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Farhadi A, Jeffs AG, Farahmand H, Rejiniemon TS, Smith G, Lavery SD. Mechanisms of peripheral phylogeographic divergence in the indo-Pacific: lessons from the spiny lobster Panulirus homarus. BMC Evol Biol. 2017;17(1):195.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Al-Marzouqi A, Al-Amry W, Al-Hadabi A, Al-Senaidi R. DNA barcoding of Panulirus homarus from Oman and Yemen. J Aquac Mar Biol. 2015;2(5):0043.

    Google Scholar 

  38. Reddy MM, Macdonald AHH, Groeneveld JC, Schleyer MH. Phylogeography of the scalloped spiny-lobster Panulirus Homarus Rubellus in the Southwest Indian Ocean. J Crust Biol. 2014;34(6):773–81.

    Article  Google Scholar 

  39. Mehanna S, Al-Shijibi S, Al-Jafary J, Al-Senaidi R. Population dynamics and management of scalloped spiny lobster Panulirus homarus in Oman coastal waters. J Biol Agr Heal. 2012;2(10):184–94.

    Google Scholar 

  40. Book Fishery Statistics. Year Book Issued by Department of Fisheries Statistics. Ministry of Agriculture and Fisheries Wealth. Sultanate of Oman; 2016.

  41. Al-Marzouqi A, Al-Nahdi A, Jayabalan N, Groeneveld J. An assessment of the spiny lobster Panulirus homarus fishery in Oman - another decline in the western Indian Ocean? Western Indian Ocean J Mar Sci. 2007;6(2):159–74.

    Google Scholar 

  42. Al-Marzouqi A, Chesalin M, Al-Shajibi S, Al-Hadabi A, Al SR. Changes in the scalloped spiny lobster, Panulirus Homarus biological structure after a shift of the fishing season. J Aquac Mar Biol. 2015;3(5):1–14.

    Google Scholar 

  43. Raymond M, Vaanto RL, Thomas F. Rousset F, deMeeus T. Renaud F Heterozygote deficiency in the mussel species complex revisited Mar Ecol Prog Ser. 1997;156:225–37.

    Google Scholar 

  44. Addison JA, Hart MW. Spawning, copulation and inbreeding coefficients in marine invertebrates. Biol Lett. 2005;1(4):450–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Zouros E. Heterozygosity and growth in marine bivalves response. Evolution. 1990;44(1):216–8.

    Article  CAS  PubMed  Google Scholar 

  46. Valles-Jimenez R, Cruz P, Perez-Enriquez R. Population genetic structure of pacific white shrimp (Litopenaeus vannamei) from Mexico to Panama: microsatellite DNA variation. Mar Biotechnol. 2004;6(5):475–84.

    Article  CAS  Google Scholar 

  47. Brumfield RT, Beerli P, Nickerson DA, Edwards SV. The utility of single nucleotide polymorphisms in inferences of population history. Trends Ecol Evol. 2003;18(5):249–56.

    Article  Google Scholar 

  48. Laurie CC, Doheny KF, Mirel DB, Pugh EW, Bierut LJ, Bhangale T, et al. Quality control and quality Assurance in Genotypic Data for genome-wide association studies. Genet Epidemiol. 2010;34(6):591–602.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Li H. Toward better understanding of artifacts in variant calling from high-coverage samples. Bioinformatics. 2014;30(20):2843–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Morse P, Kjeldsen SR, Meekan MG, McCormick MI, Finn JK, Huffard CL, et al. Genome-wide comparisons reveal a clinal species pattern within a holobenthic octopodthe Australian southern blue-ringed octopus, Hapalochlaena maculosa (Cephalopoda: Octopodidae). Ecol Evol. 2018;8(4):2253–67.

  51. Nayfa MG, Zenger KR. Unravelling the effects of gene flow and selection in highly connected populations of the silver-lip pearl oyster (Pinctada maxima). Mar Genomics. 2016;28:99–106.

    Article  PubMed  Google Scholar 

  52. Stoffel MA, Esser M, Kardos M, Humble E, Nichols H, David P, et al. inbreedR: an R package for the analysis of inbreeding based on genetic markers. Methods Ecol Evol. 2016;7(11):1331–9.

    Article  Google Scholar 

  53. Herdegen M, Dudka K, Radwan J. Heterozygosity and orange coloration are associated in the guppy (Poecilia reticulata). J Evol Biol. 2014;27(1):220–5.

    Article  CAS  PubMed  Google Scholar 

  54. Andrews KR, Luikart G. Recent novel approaches for population genomics data analysis. Mol Ecol. 2014;23(7):1661–7.

    Article  PubMed  Google Scholar 

  55. Crooks L, Carlborg O, Marklund S, Johansson AM. Identification of Null Alleles and Deletions from SNP Genotypes for an Intercross Between Domestic and Wild Chickens. G3-Gen Geno Genet. 2013;3(8):1253–60.

    Google Scholar 

  56. Kjeldsen SR, Zenger KR, Leigh K, Ellis W, Tobey J, Phalen D, et al. Genome-wide SNP loci reveal novel insights into koala (Phascolarctos cinereus) population variability across its range. Conserv Genet. 2016;17(2):337–53.

    Article  Google Scholar 

  57. Carlson CS, Smith JD, Stanaway IB, Rieder MJ, Nickerson DA. Direct detection of null alleles in SNP genotyping data. Hum Mol Genet. 2006;15(12):1931–7.

    Article  CAS  PubMed  Google Scholar 

  58. Farhadi A, Farhamand H, Nematollahi MA, Jeffs A, Lavery SD. Mitochondrial DNA population structure of the scalloped lobster Panulirus homarus (Linnaeus 1758) from the West Indian Ocean. ICES J Mar Sci. 2013;70(7):1491–8.

    Article  Google Scholar 

  59. Yellapu B, Jeffs A, Battaglene S, Lavery SD. Population subdivision in the tropical spiny lobster Panulirus ornatus throughout its indo-West Pacific distribution. ICES J Mar Sci. 2017;74(3):759–68.

    Google Scholar 

  60. Dao HT, Smith-Keune C, Wolanski E, Jones CM, Jerry DR. Oceanographic Currents and Local Ecological Knowledge Indicate, and Genetics Does Not Refute, a Contemporary Pattern of Larval Dispersal for The Ornate Spiny Lobster, Panulirus ornatus in the South-East Asian Archipelago. PLoS One. 2015;10(5):e0124568.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  61. Iacchei M, O'Malley JM, Toonen RJ. After the gold rush: population structure of spiny lobsters in Hawaii following a fishery closure and the implications for contemporary spatial management. Bull Mar Sci. 2014;90(1):331–57.

    Article  Google Scholar 

  62. Abdullah MF. Alimuddin, Muththalib M, Salama AJ, Imai H. Genetic isolation among the northwestern, southwestern and central-eastern Indian Ocean populations of the pronghorn spiny lobster Panulirus penicillatus. Int J Mol Sci. 2014;15(6):9242–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Iacchei M, Gaither MR, Bowen BW, Toonen RJ. Testing dispersal limits in the sea: range-wide phylogeography of the pronghorn spiny lobster Panulirus penicillatus. J Biogeogr. 2016;43(5):1032–44.

    Article  Google Scholar 

  64. George RW. Evolution of life cycles, including migration, in spiny lobsters (Palinuridae). N Z J Mar Freshw Res. 2005;39(3):503–14.

    Article  Google Scholar 

  65. Ivanochko, T. S. 2005. Sub-orbital scale variations in the intensity of the Arabian Sea monsoon. University of Edinburgh, PhD Thesis.

  66. Schott FA, McCreary JP. The monsoon circulation of the Indian Ocean. Prog Oceano. 2001;51(1):1–123.

    Article  Google Scholar 

  67. Condie S, Condie R. Retention of plankton within ocean eddies. Glob Ecol Biogeog. 2016;25(10):1264–77.

    Article  Google Scholar 

  68. Chiswell SM, Roemmich D. The east cape current and two eddies: a mechanism for larval retention? N Z J Mar Freshw Res. 1998;32(3):385–97.

    Article  Google Scholar 

  69. Gonzalez EB, Knutsen H, Jorde PE. Habitat discontinuities separate genetically divergent populations of a rocky shore marine fish. PLoS One. 2016;11(10):e0163052.

    Article  Google Scholar 

  70. Simpson SD, Harrison HB, Claereboudt MR, Planes S. Long-distance dispersal via ocean currents connects omani clownfish populations throughout entire species range. PLoS One. 2014;9(9):e107610.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  71. Johnson MS, Black R. Chaotic genetic patchiness in an intertidal limpet. Siphonaria sp Mar Biol. 1982;70(2):157–64.

    Article  Google Scholar 

  72. Larson RJ, Julian RM. Spatial and temporal genetic patchiness in marine populations and their implications for fisheries management. Calif Coop Ocean Fish Investig Reports. 1999;40:94–9.

    Google Scholar 

  73. Nash WJ, Goddard M, Lucas JS. Population genetic studies of the crown-of-thorns starfish, Acanthaster planci (L.), in the great barrier reef region. Coral Reefs. 1988;7910:11–8.

    Article  Google Scholar 

  74. David P, Perdieu M. A, Pernot a F, Jarne P. Fine-grained spatial and temporal population genetic structure in the marine bivalve Spisula ovalis. Evolution. 1997;51(4):1318–22.

    Article  PubMed  Google Scholar 

  75. Flowers JM, Schroeter SC, Burton RS. The recruitment sweepstakes has many winners: genetic evidence from the sea urchin Strongylocentrotus purpuratus. Evolution. 2002;56(7):1445–53.

    Article  CAS  PubMed  Google Scholar 

  76. Selwyn JD, Hogan JD, Downey-Wall AM, Gurski LM, Portnoy DS, Heath DD. Kin-aggregations explain chaotic genetic patchiness, a commonly observed genetic pattern, in a marine fish. PLoS One. 2016;11(4):1–11.

    Article  CAS  Google Scholar 

  77. Christie MR, Johnson DW, Stallings CD, Hixon MA. Self-recruitment and sweepstakes reproduction amid extensive gene flow in a coral-reef fish. Mol Ecol. 2010;19(5):1042–57.

    Article  PubMed  Google Scholar 

  78. Iacchei M, Ben-Horin T, Selkoe KA, Bird CE, García-Rodríguez FJ, Toonen RJ. Combined analyses of kinship and FSTsuggest potential drivers of chaotic genetic patchiness in high gene-flow populations. Mol Ecol. 2013;22(13):3476–94.

    Article  PubMed  PubMed Central  Google Scholar 

  79. Zhao S, Zheng P, Dong S, Zhan X, Wu Q, Guo X, et al. Whole-genome sequencing of giant pandas provides insights into demographic history and local adaptation. Nat Genet. 2013;45(1):67–71.

    Article  CAS  PubMed  Google Scholar 

  80. Dillon S, McEvoy R, Baldwin DS, Rees GN, Parsons Y, Southerton S. Characterisation of adaptive genetic diversity in environmentally contrasted populations of Eucalyptus camaldulensis Dehnh. (river red gum). Plos one. 2014; 9(8):e103515.

    Article  PubMed  PubMed Central  Google Scholar 

  81. Zhan X, Dixon A, Batbayar N, Bragin E, Ayas Z, Deutschova L, et al. Exonic versus intronic SNPs: contrasting roles in revealing the population genetic differentiation of a widespread bird species. Heredity. 2015;114(1):1–9.

    Article  CAS  PubMed  Google Scholar 

  82. Drury C, Dale KE, Panlilio JM, Miller SV, Lirman D, Larson EA, et al. Genomic variation among populations of threatened coral: Acropora cervicornis. BMC Genomics. 2016;17(1):286.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  83. Nielsen EE, Hemmer-Hansen J, Larsen PF, Bekkevold D. Population genomics of marine fishes: identifying adaptive variation in space and time. Mol Ecol. 2009;18(15):3128–50.

    Article  PubMed  Google Scholar 

  84. Tsumura Y, Uchiyama K, Moriguchi Y, Ueno S, Ihara-Ujino T. Genome scanning for detecting adaptive genes along environmental gradients in the Japanese conifer, Cryptomeria japonica. Heredity. 2012;109(6):349–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  85. Orsini L, Mergeay J, Vanoverbeke J, De Meester L. The role of selection in driving landscape genomic structure of the waterflea Daphnia magna. Mol Ecol. 2013;22(3):583–601.

    Article  PubMed  Google Scholar 

  86. Mohan R. Size structure and reproductive variation of the spiny lobster Panulirus homarus over a relatively small geographic range along the Dhofar coast in the Sultanate of Oman. R Mar Freshw Res. 1997;48(8):1085–91.

    Article  Google Scholar 

  87. Adamkewicz SL, Harasewych MG. Systematics and biogeography of the genus Donax (Bivalvia: Donacidae) in eastern North America. Amer Malacol Bulletin. 1996;13(1–2):97–103.

    Google Scholar 

  88. Jaccoud D, Peng K, Feinstein D, Kilian A. Diversity arrays: a solid state technology for sequence information independent genotyping. Nuc Acid Res. 2001;29:e25.

    Article  CAS  Google Scholar 

  89. Kilian A, Wenzl P, Huttner E, Carling J, Xia L, Blois H, et al. Diversity arrays technology: a generic genome profiling technology on open platforms. Methods Mol Biol. 2012;888:67–89.

    Article  PubMed  Google Scholar 

  90. Sansaloni C, Petroli C, Jaccoud D, Carling J, Detering F, Grattapaglia D, Kilian A. Diversity arrays technology (DArT) and next-generation sequencing combined: genome-wide, high throughput, highly informative genotyping for molecular breeding of Eucalyptus. BMC Proc. 2011;5(Suppl 7):54–P54.

    Article  Google Scholar 

  91. Lind C, Kilian A, Benzie J. Development of diversity arrays technology markers as a tool for rapid genomic assessment in Nile tilapia. Anim Genet. 2017;48(3):362–4.

    Article  CAS  PubMed  Google Scholar 

  92. Steinig E, Guppy J, Jones D, Zenger K. 2016. DartQC pipeline. https://github.com/esteinig/dartQC (2018). Accessed 5 Apr 2018.

  93. Li W, Godzik A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006;22(13):1658–9.

    Article  CAS  PubMed  Google Scholar 

  94. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. J Hum Genet. 2007;81(3):559–75.

    Article  CAS  PubMed  Google Scholar 

  95. Excoffier L, Laval G, Schneider S. Arlequin ( version 3 . 0 ): An integrated software package for population genetics data analysis. 2005;23(1):47–50.

  96. Belkhir K, Borsa P, Chikhi L, Raufaste N, Bonhomme F. 1996. GENETIX 4.05, logiciel sous Windows TM pour la génétique des populations. Universite´ de Montpellier II. http://www.genetix.univ-montp2.fr/. Accessed 5 Apr 2018.

  97. Do C, Waples RS, Peel D, Macbeth GM, Tillett BJ, Ovenden JR. NEESTIMATOR v2: re-implementation of software for the estimation of contemporary effective population size (N e) from genetic data. Mol Ecol Resour. 2014;14(1):209–14.

    Article  PubMed  Google Scholar 

  98. Weir B, Cockerham C. Estimating F statistics for the analysis of population structure. Evolution. 1984;38(6):1358–70.

    CAS  PubMed  Google Scholar 

  99. Excoffier L, Lischer HEL. Arlequin suite ver 3.5. 5 a new Ser. Programs to perform Popul. Genet. Anal. Under {Linux} {windows}. {molecular} Ecol. Resour. 2010;10:564–7.

    Google Scholar 

  100. Sheppard CRC, Salm RV. Reef and coral communities of Oman, with a description of a new coral species (order scleractinia, genus acanthastrea). J Nat Hist. 1988;22(1):263–79.

    Article  Google Scholar 

  101. Schils T, Wilson SC. Temperature threshold as a biogeographic barrier in northern Indian ocean macroalgae. J Phycol. 2006;42(4):749–56.

    Article  Google Scholar 

  102. Jombart T. Adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics. 2008;24(11):1403–5.

    Article  CAS  PubMed  Google Scholar 

  103. Jombart T, Devillard S, Balloux F. Discriminant analysis of principal components: a new method for the analysis of genetically structured populations. BMC Genet. 2010;11(1):94.

    Article  PubMed  PubMed Central  Google Scholar 

  104. Steinig EJ, Neuditschko M, Khatkar MS, Raadsma HW, Zenger KR. NETVIEW P: a network visualization tool to unravel complex population structure using genome-wide SNPs. Mol Ecol Resour. 2016;16(1):216–27.

    Article  CAS  PubMed  Google Scholar 

  105. Neuditschko M, Khatkar MS, Raadsma HW. NETVIEW: a high-definition network-visualization approach to detect fine-scale population structures from genome-wide patterns of variation. PLoS One. 2012;7(10):e48375.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  106. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30(12):2725–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  107. Foll M, Gaggiotti O. A genome-scan method to identify selected loci appropriate for both dominant and Codominant markers: a Bayesian perspective. Genetics. 2008;180(2):977–93.

    Article  PubMed  PubMed Central  Google Scholar 

  108. Beaumont MA, Nichols RA. Evaluating loci for use in the genetic analysis of population structure. Proc R Soc B Biol Sci. 1996;263(1377):1619–26.

    Article  Google Scholar 

  109. Beaumont MA, Balding DJ. Identifying adaptive genetic divergence among populations from genome scans. Mol Ecol. 2004;13(4):969–80.

    Article  CAS  PubMed  Google Scholar 

  110. Foll M. BayeScan v2.1 user manual. Ecology. 2012;20:1450–62.

Download references

Acknowledgements

We would like to thank Monal Lal and Eike Steinig for their advice and support on different statistical analysis and software including Netview R and DArTQC pipeline. Sincere thanks are also extended for Abdulaziz Al Marzouqi, Ministry of Agriculture and Fisheries Wealth for facilitating samples collection. We thank also Mohammed Al Wahaibi for his assistance in samples collection.

Funding

The present study was supported by the Research Council of Oman (TRC) (ORG/SQU/EBR/13/027).

Availability of data and materials

Genotypes data are provided as Additional files.

Author information

Authors and Affiliations

Authors

Contributions

RDB participated in the investigation design of the study, optimization of DNA extraction, analyzed and interpreted the data, drafted, revised and edited the manuscript. SRK carried out the DNA purification and optimized samples prior to library preparation, provided support in the data analysis and contributed to the writing of the manuscript. HA and MSH, extracted the DNA and provided support in the data analysis. KRZ provided distinct advices in the data quality control and analysis. DRJ participated in the investigation design, provided support in the data analysis and writing the manuscript. MAA provided support in the data analysis and reviewing the manuscript. MD developed the concept of the project and the investigation design, collected samples, provided supervisory support, contributed to interpretation of data, writing and revision of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Madjid Delghandi.

Ethics declarations

Ethics approval

All animal samples in this study complied with the Omani law on Sea Fishing and Protection of the Marine Biological Wealth RD 53–81 and fully complied with local fisheries management and marine protected area controls. All samples were purchased from local fishermen, hence no specific permits were required for the described filed sampling as the fishermen were required to comply with local laws regarding capture. The species sampled are not endangered or protected.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Genotypic data of P. homarus for 164 individuals from Oman using 7988 SNPs. Genotypes are in genetix format. (TXT 576 kb)

Additional file 2:

A plot of the Discriminant Analysis of Principal Components (DAPC) against the discriminant function retained. The plot is indicating presence of two genetic clusters of P. homarus in Oman. The plot was generated using the most informative 13 PCs identified from all 7988 SNPs dataset across 164 P. homarus individuals in the R package adegenet. (TXT 9105 kb)

Additional file 3:

Genotypic data of P. homarus for 164 individuals from Oman using 504 putatively directional SNPs. Genotypes are in genetix format. (PDF 389 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Al-Breiki, R.D., Kjeldsen, S.R., Afzal, H. et al. Genome-wide SNP analyses reveal high gene flow and signatures of local adaptation among the scalloped spiny lobster (Panulirus homarus) along the Omani coastline. BMC Genomics 19, 690 (2018). https://doi.org/10.1186/s12864-018-5044-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-018-5044-8

Keywords