Next Article in Journal
Searching for “Environmentally-Benign” Antifouling Biocides
Previous Article in Journal
The Application of an Emerging Technique for Protein–Protein Interaction Interface Mapping: The Combination of Photo-Initiated Cross-Linking Protein Nanoprobes with Mass Spectrometry
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Isolation among the Northwestern, Southwestern and Central-Eastern Indian Ocean Populations of the Pronghorn Spiny Lobster Panulirus penicillatus

1
Graduate School of Engineering and Science, University of the Ryukyus, Nishihara, Okinawa 903-0213, Japan
2
Faculty of Fisheries and Marine Science, Bogor Agricultural University, Bogor 16680, Indonesia
3
Ministry of Fisheries and Agriculture, Malé 20125, Maldives
4
Faculty of Marine Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia
5
Faculty of Science, University of the Ryukyus, Nishihara, Okinawa 903-0213, Japan
*
Author to whom correspondence should be addressed.
Present address: Kooddoo Fisheries Maldives Ltd., Block 389 Ujaalaa Higun, Malé, Maldives
Int. J. Mol. Sci. 2014, 15(6), 9242-9254; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms15069242
Submission received: 10 April 2014 / Revised: 14 May 2014 / Accepted: 16 May 2014 / Published: 26 May 2014
(This article belongs to the Section Biochemistry)

Abstract

:
The pronghorn spiny lobster Panulirus penicillatus is a highly valuable species which is widely distributed in Indo-West Pacific and Eastern Pacific regions. Mitochondrial DNA control region sequences (566–571 bp) were determined to investigate the population genetic structure of this species in the Indian Ocean. In total, 236 adult individuals of Panulirus penicillatus were collected from five locations in the Indian Ocean region. Almost all individuals had a unique haplotype. Intrapopulation haplotype (h) and nucleotide (π) diversities were high for each locality, ranging from h = 0.9986–1.0000 and π = 0.031593–0.043441. We observed distinct genetic isolation of population located at the northwestern and southwestern edge of the species range. Gene flow was found within localities in the central and eastern region of the Indian Ocean, probably resulting from an extended planktonic larval stage and prevailing ocean currents.

Graphical Abstract

1. Introduction

The pronghorn spiny lobster Panulirus penicillatus is widely distributed in the Indo-West Pacific and Eastern Pacific regions. It is found in tropical and adjacent regions from South-Eastern Africa, the Red Sea, Southern India, the Southeast Asian Archipelago, Japan, Northern Australia, and the Southern and Western Pacific Islands, to Hawaii, the Galápagos Islands, and other islands of the Eastern Pacific [1,2,3]. P. penicillatus supports considerable fisheries in Indo-Pacific regions. The species is fished wherever it occurs. Much attention has been paid to biological investigation of young and adult stages. Recognition of stocks is complicated by the life cycle of pronghorn spiny lobster, primarily because of the potential for long-range dispersal during its planktonic larvae. The early life history of spiny lobsters consists of a drifting larval period adapted for a relatively long-term stay in the open ocean, extending from several months to more than a year, with many possibilities for dispersal through ocean currents [4,5].
Several molecular techniques, including allozymes, restriction fragment length polymorphism (RFLP), mitochondrial DNA sequencing, and microsatellites, have been used in studying genetic variability in marine organisms at the population level [6,7]. Mitochondrial DNA sequencing, particularly of the most rapidly evolving and highly variable control region, has been a useful tool for population genetic studies of many terrestrial and aquatic organisms [8]. The highly polymorphic mitochondrial control region sequence has been used previously for population genetics studies in spiny lobster species [9,10,11,12].
There recently been studies of P. penicillatus, East-West Pacific differentiation using 16s and COI (cytochrome oxidase I) marker was described by Chow et al. [13], and this study was corroborated with the study of genetic population structure of P. penicillatus in the Pacific region using highly polymorphic marker, mitochondrial control region by Abdullah et al. [3]. In this study, we focused on the Indian Ocean region of P. penicillatus.
There have been no previous population genetic studies of P. penicillatus in the Indian Ocean region. The aim of the present study is to investigate population genetic structuring in a widely distributed P. penicillatus in the Indian Ocean region. We address this question using mitochondrial control region sequences data from Red Sea, Jeddah, Saudi Arabia; Fort Dauphin, Madagascar; Nilandhoo atoll, Maldives; Aceh and Java Sea, Indonesia (Figure 1). Given the economic importance of P. penicillatus, this genetic information may be valuable for long-term fisheries’ management decisions.
Figure 1. Map of the Indian Ocean (left): showing SEC—South Equatorial Current, SECC—South Equatorial Countercurrent, SEMC (NEMC)—Southeast (Northeast) Madagascar Current, NMC (SMC)—Northeast (Southwest) Monsoon Current, JC—South Java Current, ocean current was drawn referring to Schott and McCreary [14]. Panulirus penicillatus haplotype network (right); Sizes of the circles are proportional to the frequency of each haplotype; lengths of the lines are relative to the number of mutations between haplotypes.
Figure 1. Map of the Indian Ocean (left): showing SEC—South Equatorial Current, SECC—South Equatorial Countercurrent, SEMC (NEMC)—Southeast (Northeast) Madagascar Current, NMC (SMC)—Northeast (Southwest) Monsoon Current, JC—South Java Current, ocean current was drawn referring to Schott and McCreary [14]. Panulirus penicillatus haplotype network (right); Sizes of the circles are proportional to the frequency of each haplotype; lengths of the lines are relative to the number of mutations between haplotypes.
Ijms 15 09242 g001

2. Results and Discussion

2.1. Genetic Diversity Analysis

A total of 566–571 bp sequence data were obtained from 236 specimens of P. penicillatus. Of the 577 aligned base pairs, 246 variable sites (including indels) were found: 172 were parsimony informative and 74 were singletons. In total, 231 haplotypes were identified, and nearly all sequenced individuals had a unique haplotype (Table 1). The nucleotide sequences of the haplotypes were deposited in the DNA Data Bank of Japan (DDBJ) (accession numbers AB922639–AB922816, AB689451–AB689503).
As expected, this fragment is AT-rich in nucleotide composition with ratios of 37.8% adenine, 31.8% thymine, 19.4% cytosine, and 11% guanine, consistent with previous reports of an AT-rich control region of the mitochondrial genome in many invertebrates, including crustaceans [7,15].
The haplotype network showed a geographic structure with two main groups (Figure 1): (1) the northeastern edge of the population; Red Sea cluster and (2) other localities from the southwestern to eastern Indian Ocean (Madagascar, Maldives, Aceh and Java) that are closely together.
Table 1. Sample abbreviation (ID), sample locality, sample size (N), year of collection, numbers of haplotypes, haplotype diversity (h), nucleotide diversity (π), Tajima’s D and Fu’s F Test for Panulirus penicillatus mitochondrial DNA control region.
Table 1. Sample abbreviation (ID), sample locality, sample size (N), year of collection, numbers of haplotypes, haplotype diversity (h), nucleotide diversity (π), Tajima’s D and Fu’s F Test for Panulirus penicillatus mitochondrial DNA control region.
IDLocalityNYearNo. of Haplotypes (Unique)h ± SDπ ± SDTajima’s D (p-Value)Fu’s F (p-Value)
JDRed Sea, Jeddah, Saudi Arabia46201145(45)0.9990 ± 0.00480.031593 ± 0.015851−0.79131 (0.215)−24.20576 (0.000)
MGFort Dauphin, Madagascar49201449(49)1.0000 ± 0.00410.043441 ± 0.021533−0.97966 (0.168)−24.14011 (0.000)
MDNilandhoo atoll, Maldives38201337(36)0.9986 ± 0.00650.034957 ± 0.017561−1.05781 (0.127)−19.70407 (0.000)
ACAceh, Indonesia48201147(47)0.9991 ± 0.00450.038100 ± 0.018974−1.17981 (0.098)−24.14392 (0.000)
JVJava Sea, Indonesia55200853(52)0.9993 ± 0.00360.038855 ± 0.019284−1.01996 (0.153)−24.12306 (0.000)
Total236------
Genetic diversity, as shown by haplotype diversity, was found to be very high in P. penicillatus. The values generated in this study are comparable with those derived from the mtDNA control region sequences for the crustaceans Farfantepenaeus aztecus, Litopenaeus setiferus and Sergia lucens [16,17]. Haplotype diversity (h) within the geographical populations was high, ranging from 0.9986 in Nilandhoo atoll, Maldives to 1.0000 in Fort Dauphin, Madagascar. Nucleotide diversity (π) was generally high, ranging from 0.031593 in Red Sea, Jeddah, Saudi Arabia to 0.043441 in Fort Dauphin, Madagascar (Table 1). There were few shared haplotypes among locations belonging to the central-eastern Indian Ocean.
In the present study, the mtDNA control region exhibited high genetic variability (h = 0.9986−1.0000 and π = 0.031593−0.043441), which are relatively high values compared with the control region diversities of other Panulirus spp. (P. inflatus: h = 0.957–1.0000 and π = 0.019–0.028 [10]; P. argus: h = 0.667–1.0000 and π = 0.003–0.066 [7], P. homarus: h = 0.985–1.0000 and π = 0.031–0.041 [12]) and similar to Pacific region populations of P. penicillatus (h = 0.9996–1.0000 and π = 0.0255–0.0448) [3]. The control region is a non-coding region that includes signals necessary for replication of a molecule and is the most rapidly evolving region of mtDNA [18].
The high genetic variability of the control region in P. penicillatus may be maintained by its large population size. In fact, because nearly all individuals represented a unique haplotype, an investigation based on sequencing the control region would require a very large sample size to detect common haplotypes and associate them with particular geographic areas [10]. The evolutionary pattern of interhaplotypic divergence will almost always produce a biogeographic signature, and an analysis that recognizes genetic distance among alleles will provide additional resolution, especially for hypervariable markers [19].

2.2. Population Structure Analysis

Genetic analysis of mitochondrial control region revealed distinct population at the northwestern and southwestern margin of geographic range of P. penicillatus with other localities. This suggests that the geographic isolation may be due to contemporary isolation/gene flow or prior isolation with recent gene flow, although the overall level of genetic exchange must be below that required to homogenize the populations.
Congruent with the pattern observed in the network haplotypes (Figure 1), we found strong regional structuring between Red Sea and other Indian Ocean populations of P. penicillatus, which explained high genetic variations (51%) between the regions. High genetic variation (85%) was also observed in previous study between eastern Pacific and west-central Pacific population of P. penicillatus [3]. To analyze the sub-structure in the other area we attempted three types of regional grouping: (1) Madagascar and Maldives/Aceh/Java; (2) Madagascar/Maldives and Aceh/Java; and (3) Madagascar/Aceh and Maldives/Java. Regarding these groupings, genetic variation was relatively greater between Madagascar and Maldives/Aceh/Java (2.41%) than the variations between other combinations (%var −0.49 and −0.75) (Table 2).
Table 2. AMOVA results showing degree of freedom (df), variance components (Var.), percent variation (% Var.) and Φ-statistics for Panulirus penicillatus (1000 permutations).
Table 2. AMOVA results showing degree of freedom (df), variance components (Var.), percent variation (% Var.) and Φ-statistics for Panulirus penicillatus (1000 permutations).
Source of VariationdfVar.% Var.Φ-Statisticsp-Value
All sites (Red Sea/Madagascar-Maldives-Aceh-Java)
Among groups111.5081951.31ΦCT = 0.51313p = 0.19746
Among populations30.188860.84ΦSC = 0.01730p = 0.00098
Within populations23110.7303047.84ΦST = 0.51313p = 0.00000
Substructure grouping test: (Madagascar-Maldives/Aceh/Java)
Among groups10.275912.41ΦCT = 0.02408p = 0.24047
Among populations20.038500.34ΦSC = 0.00344p = 0.14467
Within populations18611.1441397.26ΦST = 0.02744p = 0.00098
(Madagascar/Maldives-Aceh/Java)
Among groups1−0.05540−0.49ΦCT = −0.00490p = 0.67742
Among populations20.216971.92ΦSC = 0.01910p = 0.00293
Within populations18611.1441398.57ΦST = 0.01429p = 0.00098
(Madagascar/Aceh-Maldives/Java)
Among groups1−0.08441−0.75ΦCT = −0.00747p = 1.00000
Among populations20.236662.09ΦSC = 0.02079p = 0.00098
Within populations18611.1441398.65ΦST = 0.01348p = 0.00098
Central-eastern region (Maldives/Aceh/Java) samples were pooled, suggested by population tree, haplotype network and ocean division. In addition, our preliminary analyses of pairwise ΦST comparisons of each data set found no heterogeneity as observed within samples from the central-eastern region (Maldives/Aceh/Java) (data not shown), although it was a limited ΦCT (Table 2, grouping Madagascar-Maldives/Aceh/Java). This supports the notion of a genetically homogeneous central-eastern region population that includes Maldives, Aceh and Java Sea, justified by combining the samples within the locality.
Populations of P. penicillatus in the Indian Ocean region do not represent a genetically homogenous assemblage. Pairwise ΦST comparisons detected genetic structure between the central-eastern region (Maldives, Aceh, and Java) and other populations (Madagascar and Red Sea) (Table 3). A relatively small structure was found between the southwestern population (Madagascar) and central-eastern populations (Maldives/Aceh/Java) of P. penicillatus (ΦST value = 0.02629, p < 0.005; see Table 3). This pattern is in agreement with was previously observed from this region in reef fish Myripritis berndti [20] and amphidromous prawn Macrobrachium lar [21]. The prevailing ocean currents in the southwest Indian Ocean splits the coast of Madagascar into the southeast Madagascar current and the northeast Madagascar current, which may have limited genetic exchange with other regions [21,22].
Table 3. Pairwise ΦST values (above the diagonal) and pairwise ΦST p-values (below the diagonal) of mitochondrial DNA control region among populations of Panulirus penicillatus.
Table 3. Pairwise ΦST values (above the diagonal) and pairwise ΦST p-values (below the diagonal) of mitochondrial DNA control region among populations of Panulirus penicillatus.
PopulationNorthwestern Region (Red Sea)Southwestern Region (Madagascar)Central-Eastern Region (Maldives-Aceh-Java)
Northwestern Region (Red Sea)-0.510580.53463
Southwestern Region (Madagascar)0.00000 *-0.02629
Central-Eastern Region (Maldives-Aceh-Java)0.00000 *0.00098 *-
* Significant p-value after Bonferroni correction (α = 0.05).
The greatest significance is found in northwestern region (Red Sea) with supported by significant and high pairwise ΦST values (Table 3), UPGMA tree based on ΦST values showed a prominent genetic break between Red Sea and other populations (Figure 2) this result also corroborated with the distinction of the Red Sea cluster in the haplotype network (Figure 1), other studies on invertebrates, mud crab Scylla serrata [23,24], have shown a genetic isolation of Red Sea populations from the Indian Ocean. Oceanographic currents are likely to create barriers and impact routes and directions of larval dispersal [25,26,27]. More importantly, the lack of shared mtDNA haplotypes and large ΦST values suggest that there has been prolonged historical isolation between the Red Sea and others populations. In addition, high degree of endemism and biodiversity reef fishes were reported previously from the Red Sea, despite a geological history characterized by intermittent isolation and multiple salinity crises, mtDNA indicate that some widespread reef fish species originate in the Red Sea with subsequent contributing species to adjacent regions [28].
Further analysis on recruitment of P. penicillatus phyllosoma larvae may be necessary since genetic distinction in Red Sea populations may indicate the potential self-recruitment of phyllosoma larvae in this region. Although phyllosoma larvae are weak swimmers, they have the ability to maintain their vertical position in the water column [29]. Larval dispersal for spiny lobsters has been previously modeled. Griffin et al. [30], discovered that P. cygnus was a self-recruiting species, restricted to Western Australia. A more recent simulation study of P. argus found the average dispersal of P. argus in the Caribbean Sea to be possibly only 200–400 km, suggesting that larval behavior such as vertical migration, together with a retentive oceanographic environment, may increase the potential for self-recruitment [31].
No population structure was found in the localites of P. penicillatus inhabiting the central-eastern region of Indian Ocean: Nilandhoo atoll, Maldives; Aceh and Java Sea, Indonesia. The early life history of spiny lobsters consists of a drifting larval period adapted for a relatively long-term stay in the open ocean, extending from several months to more than a year, with many possibilities for dispersal through ocean currents [4,5]. Long larval life (8.3–9.4 months) of the P. penicillatus phyllosoma [32], coupled with the considerable impact of oceanic current systems or gyres (e.g., Equatorial currents and South Java current) (Figure 1) may explain the observed patterns. Furthermore, transportation of larvae likely occurred indirectly via stepping stones by the currents, and a small amount of larval dispersal over years might result in genetic homogeneity. Indeed, pelagic larval duration was considerably correlated with dispersal distance, but there were many exceptions; larval behavior can play a crucial role in determining dispersal distance [33], while type of eggs and life-history parameters are also important predictors of connectivity in fishes [27,34].
Figure 2. UPGMA tree based on ΦST values, illustrating the most probable geographical structure in the analysis of molecular variance (AMOVA).
Figure 2. UPGMA tree based on ΦST values, illustrating the most probable geographical structure in the analysis of molecular variance (AMOVA).
Ijms 15 09242 g002
The neutrality of mutations in the mtDNA control region was rejected on the basis of Tajima’s D and Fu’ F tests. The significant negative values obtained in these P. penicillatus are often associated which experienced a population expansion during their geographic history [35]. This pattern has also been observed in other spiny lobster [12,36,37].
Although mtDNA has been proven to be a very useful marker, its use is not without complications. Due to the exclusively maternal inheritance mtDNA represents a marker strictly for historical process in females; they can potentially introgress between species [38,39], and symbiont-driven changes in mtDNA variation over space specifically in arthropods can occur [40]. Therefore, further studies incorporating supporting evidence from nuclear markers is required to obtain a more robust indication.

3. Experimental Section

3.1. Lobster Samples

Adults of Panulirus penicillatus were collected from five localities in the Indian Ocean region including Red Sea, Jeddah, Saudi Arabia; Fort Dauphin, Madagascar; Nilandhoo atoll, Maldives; Aceh and Java Sea, Indonesia (Figure 1; Table 1). Lobsters were purchased from local commercial fishers and fully complied with local fisheries management and marine protected area controls. Tissues samples from walking legs (pereiopods) or abdomen were dissected on site, immediately fixed in 70%–99% ethanol, and transferred to the laboratory.

3.2. DNA Analysis

About 50 mg of finely minced tissue sample was added to a 1.5 mL plastic test tube containing 0.5 mL TNES-8 M urea buffer [41]. After adding 10 μL proteinase K, it was incubated at 38 °C for 3 h and then genomic DNA was extracted using phenol-chloroform procedure and precipitated with absolute ethanol according method modified in Imai et al. [42]. The mtDNA control region was amplified using the polymerase chain reaction (PCR). PCR primers previously designed for mitochondrial control region (Panulirus12s: 5'-TATAGCAAGAATCAAACTATAG; and Penicillatus-R: 5'-CATAGG(T/C)GTG(T/C)GAGGGAACAAAGTC) were used for amplification of the control region [3]. PCR amplifications were performed in 50 μL reaction mixtures containing 1 μL template DNA, 12.5 pmol of each primer and 25 μL EmeraldAmp PCR Master Mix (2× Premix) (Takara Bio Inc., Shiga, Japan); the final volume of the reaction mixture was adjusted to 50 μL with sterile water. Reactions were performed in a thermal cycler (Perkin Elmer GeneAmp PCR System 9700) under conditions of an initial denaturation step at 94 °C for 2 min, followed by 30 cycles of 30 s at 95 °C, 30 s at 50 °C, and 1 min at 72 °C, with a final 7 min extension at 72 °C. PCR products were purified using a PCR Product Pre-sequencing kit (Exosap; USB Co., Cleveland, OH, USA). The cleaned products were sequenced on an ABI 3730xl Genetic Analyzer (Applied Biosystems, Foster City, CA, USA) using a Big Dye Terminator Cycle Sequencing kit (ver. 3.1; Applied Biosystems).

3.3. Genetic Data and Analyses

Sequence data were aligned using ClustalX [43], with default alignment parameters and were checked manually for misalignments. The nucleotide compositions and numbers of variable sites were assessed with MEGA6 [44]. Haplotype and nucleotide diversity for each location were estimated using Arlequin (ver. 3.5) software [45]. AMOVA were performed to test the geographic divisions among population experimented with various groupings. This approach is a hierarchical approach that computes the proportion of variations among groups (ΦCT), the proportion of variation among population within groups (ΦSc), and the proportion of variation within populations (ΦST). We ran preliminary analysis of pairwise ΦST comparisons on each data sheet. The statistical significance of ΦST value was tested by 1000 permutations in Arlequin (ver. 3.5) software [45]. The Bonferroni test [46] was used to correct for multiple tests of the hypothesis that pairwise ΦST statistics did not differ from zero. An UPGMA tree based on ΦST values was constructed using Neighbor in Phylip (ver. 3.6) [47,48]. The median-joining network [49], for the haplotypes was estimated using Network (ver. 4.611) [50]. Past demographic patterns were inferred by Tajima’s D [51], and Fu’s F [52], estimated using Arlequin (ver. 3.5). Both tests are commonly used to test neutrality; however, they can also be used to examine population growth because population expansion may result in rejection of the null hypothesis of neutrality (significant negative value).

4. Conclusions

The present study detected genetic structure in P. penicillatus populations located at the northwestern and southwestern edge of the Indian Ocean region with other populations. Relatively small structure was found between southwestern population and central-eastern, while vast majority structures were observed between the northwestern edge population and the other populations. Despite a lengthy larval period, oceanographic currents are likely to create barriers and impact routes and directions of larval dispersal in the northwestern edge population that may have limited genetic exchange with other regions. These results have implications for fisheries management in the region; this should essentially provide information for ascertaining stock boundaries to better evaluate the long-term needs of spiny lobster management for this valuable species.
At the scale of the central-eastern Indian Ocean, results suggested that P. penicillatus was belonging to a single unique panmictic population or at least several breeding grounds with significant exchange of genetic material. Long larval life (8.3–9.4 months) of the P. penicillatus coupled with the considerable impact of oceanic current systems or gyres (e.g., Equatorial currents and South Java current) may explain the observed patterns. Furthermore, transportation of larvae likely occurred indirectly via stepping stones by the currents, and a small amount of larval dispersal over some years might result in genetic homogeneity.

Acknowledgments

We thank the Merci Co., Ltd. (Shizuoka, Japan) for the help in sampling of spiny lobster from Madagascar. This study was supported in part by grants from the University of the Ryukyus grant for “International Research Hub Project for Climate Change and Coral Reef/Island Dynamics” and “the 2012 Nissui Research Fund” from Nippon Suisan Kaisha., Ltd. (Tokyo, Japan).

Author Contributions

Conceived and designed the experiments: M.F.A., H.I. Performed the experiments: M.F.A., H.I. Analyzed the data: M.F.A., H.I. Contributed to sampling materials/reagents: M.F.A., A., M.M., A.J.S., H.I. Wrote the paper: M.F.A., H.I.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Holthuis, L.B. Marine Lobsters of the World: An Annotated and Illustrated Catalogue of Species of Interest to Fisheries Known to Date; Food and Agriculture Organization of the United Nations: Rome, Italy, 1991; Volume 13. [Google Scholar]
  2. Pitcher, C.R. Spiny lobster. In Nearshore Marine Resources of the South Pacific. Information for Fisheries Development and Management; Wright, A., Hill, L., Eds.; IPS: Suva; Fiji/FFA: Honiara; Solomon/ICOD: Halifax, NS, Canada, 1993; pp. 539–607. [Google Scholar]
  3. Abdullah, M.F.; Chow, S.; Sakai, M.; Cheng, J.; Imai, H. Genetic diversity and population structure of pronghorn spiny lobster Panulirus penicillatus in the Pacific region. Pac. Sci. 2014, 68, 197–211. [Google Scholar] [CrossRef]
  4. Booth, J.D.; Phillips, B.F. Early life history of spiny lobster. Crustaceana 1994, 66, 271–294. [Google Scholar] [CrossRef]
  5. Tolley, K.A.; Groeneveld, J.C.; Gopal, K.; Matthee, C.A. Mitochondrial DNA panmixia in spiny lobster Palinurus gilchristi suggests a population expansion. Mar. Ecol. Prog. Ser. 2005, 297, 225–231. [Google Scholar] [CrossRef]
  6. Shaklee, J.B.; Bentzen, P. Genetic identification of stocks of marine fish and shellfish. Bull. Mar. Sci. 1998, 62, 589–621. [Google Scholar]
  7. Diniz, F.M.; Maclean, N.; Ogawa, M.; Cintra, I.H.A.; Bentzen, P. The hypervariable domain of the mitochondrial control region in Atlantic spiny lobsters and its potential as a marker for investigating phylogeographic structuring. Mar. Biotechnol. 2005, 7, 462–473. [Google Scholar] [CrossRef]
  8. Avise, J.C. Molecular Markers: Natural History and Evolution; Chapman and Hall: London, UK, 1994. [Google Scholar]
  9. Inoue, N.; Watanabe, H.; Kojima, S.; Sekiguchi, H. Population structure of Japanese spiny lobster Panulirus japonicus inferred by nucleotide sequence analysis of mitochondrial COI gene. Fish. Sci. 2007, 73, 550–556. [Google Scholar] [CrossRef]
  10. García-Rodríguez, F.J.; Perez-Enriquez, R. Lack of genetic differentiation of blue spiny lobster Panulirus inflatus along the Pacific coast of Mexico inferred from the mtDNA sequences. Mar. Ecol. Prog. Ser. 2008, 361, 203–212. [Google Scholar] [CrossRef]
  11. Sarver, S.; Silberman, J.D.; Walsh, P.J. Mitochondrial DNA sequence evidence supporting the recognition of two subspecies or species of the Florida spiny lobster Panulirus argus. J. Crust. Biol. 1998, 18, 177–186. [Google Scholar] [CrossRef]
  12. Farhadi, A.; Farhamand, H.; Nematollahi, M.A.; Jeffs, A.; Lavery, S.D. Mitochondrial DNA population structure of the scalloped lobster Panulirus homarus (Linnaeus 1758) from the West Indian Ocean. ICES J. Mar. Sci. 2013, 70, 1491–1498. [Google Scholar] [CrossRef]
  13. Chow, S.; Jeff, A.; Miyaeke, Y.; Konishi, K.; Okazaki, M.; Suzuki, N.; Kimura, S.; Abdullah, M.F.; Imai, H.; Sasaki, M.; et al. Genetic isolation between the western and eastern Pacific populations of pronghorn spiny lobster Panulirus penicillatus. PLoS One 2011, 6, e29280. [Google Scholar]
  14. Schott, F.; McCreary, J.P. The monsoon circulation of the Indian Ocean. Prog. Oceanogr. 2001, 51, 1–123. [Google Scholar] [CrossRef]
  15. McMillen-Jackson, A.L.; Bert, T.M. Genetic diversity in the mtDNA control region and population structure in the pink shrimp Farfantepenaeus duorarum. J. Crust. Biol. 2004, 24, 101–109. [Google Scholar] [CrossRef]
  16. McMillen-Jackson, A.L.; Bert, T.M. Disparate patterns of population genetic structure and population history in two sympatric penaeid shrimp species (Farfantepenaeus aztecus and Litopenaeus setiferus) in the eastern United States. Mol. Ecol. 2003, 12, 2895–2905. [Google Scholar] [CrossRef]
  17. Imai, H.; Hanamura, Y.; Cheng, J.H. Genetic and morphological differentiation in the Sakura shrimp (Sergia lucens) between Japanese and Taiwanese populations. Contrib. Zool. 2013, 82, 123–130. [Google Scholar]
  18. Heyer, E.; Zietkiewicz, E.; Rochowski, A.; Yotova, V.; Puymirat, J.; Labuda, D. Phylogenetic and familial estimates of mitochondrial substitution rates: study of control region mutations in deep-rooting pedigrees. Am. J. Hum. Genet. 2001, 69, 1113–1126. [Google Scholar] [CrossRef]
  19. Bird, C.E.; Smouse, P.E.; Karl, S.A.; Toonen, R.J. Detecting and measuring genetic differentiation. In Crustacean Issues: Phylogeography and Population Genetics in Crustacea; Koenemann, S., Held, C., Schubart, C., Eds.; CRC Press: Boca Raton, FL, USA, 2011; pp. 31–55. [Google Scholar]
  20. Muths, D.; Tessier, E.; Gouws, G.; Craig, M.; Mwale, M.; Mwaluma, J.; Mwandya, A.; Bourjea, J. Restricted dispersal of the reef fish Myripristis berndti at the scale of the SW Indian Ocean. Mar. Ecol. Prog. Ser. 2011, 443, 167–180. [Google Scholar] [CrossRef]
  21. Castelin, M.; Feutry, P.; Hautecoeur, M.; Marquet, G.; Wowor, D.; Zimmermann, G.; Keith, P. New insight on population genetic connectivity of widespread amphidromous prawn Macrobrachium lar (Fabricius, 1798) (Crustacea: Decapoda: Palaemonidae). Mar. Biol. 2013, 160, 1395–1406. [Google Scholar] [CrossRef]
  22. Chapman, P.; di Marco, S.F.; Davis, R.E.; Coward, A.C. Flow at intermediate depths around Madagascar based on ALACE float trajectories. Deep Sea Res. Part II Top. Stud. Oceanogr. 2003, 50, 1957–1986. [Google Scholar] [CrossRef]
  23. Gopurenko, D.; Hughes, J.M.; Keenan, C.P. Mitochondrial DNA evidence for rapid colonisation of the Indo-West Pacific by the mudcrab Scylla serrata. Mar. Biol. 1999, 134, 227–233. [Google Scholar] [CrossRef]
  24. Fratini, S.; Vannini, M. Genetic differentiation in the mud crab Scylla serrata (Decapoda: Portunidae) within the Indian Ocean. J. Exp. Mar. Biol. Ecol. 2002, 272, 103–116. [Google Scholar] [CrossRef]
  25. Baums, I.B.; Paris, C.B.; Cherubin, L.M. A biooceanographic filter to larval dispersal in a reef-building coral. Limnol. Oceanogr. 2006, 51, 1969–1981. [Google Scholar] [CrossRef]
  26. Treml, E.A.; Halpin, P.N.; Urban, D.L.; Pratson, L.F. Modeling population connectivity by ocean currents, a graph-theoretic approach for marine conservation. Landsc. Ecol. 2008, 23, 19–36. [Google Scholar] [CrossRef]
  27. Riginos, C.; Douglas, K.E.; Jin, Y.; Shanahan, D.F.; Treml, E.A. Effects of geography and life history traits on genetic differentiation in benthic marine fishes. Ecography 2011, 34, 566–575. [Google Scholar] [CrossRef]
  28. DiBattista, J.D.; Berumen, M.L.; Gaither, M.R.; Rocha, L.A.; Eble, J.A.; Choat, J.H.; Craig, M.T.; Skillings, D.J.; Bowen, B.W. After continents divide: comparative phylogeography of reef fishes from the Red Sea and Indian Ocean. J. Biogeogr. 2013, 40, 1170–1181. [Google Scholar] [CrossRef]
  29. Johnson, M.W. On the dispersal of lobster larvae into the East Pacific Barrier (Decapoda, Palinuridea). Fish. Bull. 1974, 72, 639–647. [Google Scholar]
  30. Griffin, D.A.; Wilkin, J.L.; Chubb, C.F.; Pearce, A.F.; Caputi, N. Ocean currents and the larval phase of Australian western rock lobster, Panulirus cygnus. Mar. Freshw. Res. 2001, 52, 1187–1199. [Google Scholar] [CrossRef]
  31. Butler, M.J., IV; Paris, C.B.; Goldstein, J.S.; Matsuda, H.; Cowen, R.K. Behavior constrains the dispersal of long-lived spiny lobster larvae. Mar. Ecol. Prog. Ser. 2011, 422, 223–237. [Google Scholar] [CrossRef]
  32. Matsuda, H.; Takenouchi, T.; Goldstein, J.S. The complete larval development of the pronghorn spiny lobster Panulirus penicillatus (Decapoda: Palinuridae) in culture. J. Crust. Biol. 2006, 26, 579–600. [Google Scholar] [CrossRef]
  33. Shanks, A.L. Pelagic larval duration and dispersal distance revisited. Biol. Bull. 2009, 216, 373–385. [Google Scholar]
  34. Bradbury, I.R.; Laurel, B.; Snelgrove, P.V.R.; Bentzen, P.; Campana, S.E. Global patterns in marine dispersal estimates: The influence of geography, taxonomic category and life history. Proc. R. Soc. B 2008, 275, 1803–1809. [Google Scholar] [CrossRef]
  35. Rogers, A.R.; Harpending, H. Population growth waves in the distribution of pairwise genetic differences. Mol. Biol. Evol. 1992, 9, 552–569. [Google Scholar]
  36. Palero, F.; Abelló, P.; Macpherson, E.; Gristina, M.; Pascual, M. Phylogeography of the European spiny lobster (Palinurus elephas): Influence of current oceanographical features and historical processes. Mol. Phylogenet. Evol. 2008, 48, 708–717. [Google Scholar] [CrossRef]
  37. Babbucci, M.; Buccoli, S.; Cau, A.; Cannas, R.; Goñi, R.; Díaz, D.; Marcato, S.; Zane, L.; Patarnello, T. Population structure, demographic history, and selective processes: Contrasting evidences from mitochondrial and nuclear markers in the European spiny lobster Palinurus elephas (Fabricius, 1787). Mol. Phylogenet. Evol. 2010, 56, 1040–1050. [Google Scholar] [CrossRef]
  38. Ballard, J.W.; Whitlock, M.C. The incomplete natural history of mitochondria. Mol. Ecol. 2004, 13, 729–744. [Google Scholar] [CrossRef]
  39. Jiggins, F.M. Male-killing Wolbachia and mitochondrial DNA: Selective sweeps, hybrid introgression and parasite population dynamics. Genetics 2003, 164, 5–12. [Google Scholar]
  40. Hurst, G.D.D.; Jiggins, F.M. Problems with mitochondrial DNA as a marker in population, phylogeographic and phylogenetic studies: The effect of inherited symbionts. Proc. R. Soc. B 2005, 272, 1525–1534. [Google Scholar] [CrossRef]
  41. Asahida, T.; Kobayashi, T.; Saitoh, K.; Nakayama, I. Tissue preservation and total DNA extraction from fish stored at ambient temperature using buffers containing high concentration of urea. Fish. Sci. 1996, 62, 727–730. [Google Scholar] [CrossRef]
  42. Imai, H.; Cheng, J.-H.; Hamasaki, K.; Numachi, K.-I. Identification of four mud crab species (genus Scylla) using ITS-1 and 16S rDNA markers. Aquat. Living Resour. 2004, 17, 31–34. [Google Scholar] [CrossRef]
  43. Thompson, J.D.; Gibson, T.J.; Plewniak, F.; Jeanmougin, F.; Higgins, D.G. The ClustalX windows interface: Flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids. Res. 1997, 25, 4876–4882. [Google Scholar] [CrossRef]
  44. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol. Biol. Evol. 2013, 30, 2725–2729. [Google Scholar] [CrossRef]
  45. Excoffier, L.; Laval, G.; Schneider, S. Arlequin ver. 3.0: An integrated software package for population genetics data analysis. Evol. Bioinform. Online 2005, 1, 47–50. [Google Scholar]
  46. Rice, W.R. Analyzing tables of statistical tests. Evolution 1989, 43, 223–225. [Google Scholar] [CrossRef]
  47. Felsenstein, J. PHYLIP-Phylogeny Inference Package. Cladistics 1989, 5, 164–166. [Google Scholar]
  48. Felsenstein, J. PHYLIP Home Page. Available online: http://evolution.genetics.washington.edu/phylip.html (accessed on 21 May 2014).
  49. Bandelt, H.J.; Forster, P.; Rohl, A. Median-joining networks for inferring intraspecific phylogenies. Mol. Biol. Evol. 1999, 16, 37–48. [Google Scholar] [CrossRef]
  50. Fluxus-engineering.com. Available online: http://www.fluxus-engineering.com (accessed on 21 February 2014).
  51. Tajima, F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 1989, 123, 585–595. [Google Scholar]
  52. Fu, Y.-X. Statistical tests of neutrality of mutations against population growth, hitchhiking and backgroud selection. Genetics 1997, 147, 915–925. [Google Scholar]

Share and Cite

MDPI and ACS Style

Abdullah, M.F.; Alimuddin; Muththalib, M.; Salama, A.J.; 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, 9242-9254. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms15069242

AMA Style

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. International Journal of Molecular Sciences. 2014; 15(6):9242-9254. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms15069242

Chicago/Turabian Style

Abdullah, Muhamad Fadry, Alimuddin, Mohamed Muththalib, Adnan Jameel Salama, and Hideyuki Imai. 2014. "Genetic Isolation among the Northwestern, Southwestern and Central-Eastern Indian Ocean Populations of the Pronghorn Spiny Lobster Panulirus penicillatus" International Journal of Molecular Sciences 15, no. 6: 9242-9254. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms15069242

Article Metrics

Back to TopTop