High population connectivity across the Indo-Pacific: Congruent lack of phylogeographic structure in three reef fish congeners

Geographic genetic variation in the Coral Hawkfish, Cirrhitichthys oxycephalus (Cirrhitidae), in relation to biogeographic barriers across the Tropical Indo-Pacific

The Tropical Indo-Pacific (TIP) includes about two thirds of the world's tropical oceans and harbors an enormous number of marine species. The distributions of those species within the region is affected by habitat discontinuities and oceanographic features. As well as many smaller ones, the TIP contains seven large recognized biogeographic barriers that separate the Red Sea and Indian Ocean, the Indian from the Pacific Ocean, the central and eastern Pacific, the Hawaiian archipelago, the Marquesas and Easter Islands. We examined the genetic structuring of populations of Cirrhitichthys oxycephalus, a small cryptic species of reef fish, across its geographic range, which spans the longitudinal limits of the TIP. We assessed geographic variation in the mitochondrial cytb gene and the nuclear RAG1 gene, using 166 samples collected in 46 localities from the western to eastern edges of the TIP. Sequences from cytb show three well-structured groups that are separated by large genetic distances (1.58-2.96%): two in the Tropical Eastern Pacific (TEP), one at Clipperton Atoll another occupying the rest of that region and the third that ranges across the remainder of the TIP, from the central Pacific to the Red Sea and South Africa. These results indicate that the ~4,000 km wide Eastern Pacific Barrier between the central and eastern Pacific is an efficient barrier separating the two main groups. Further, the ~950 km of open ocean that isolates Clipperton Atoll from the rest of the TEP is also an effective barrier. Contrary to many other cases, various major and minor barriers from the Central Indo-Pacific to the Red Sea are not effective against dispersal by C. oxycephalus, although this species has not colonized the Hawiian islands and Easter Island. The nuclear gene partially supports the genetic structure evident in cytb, although all haplotypes are geographically mixed.

Genetic diversity of Flaviviridae and Rhabdoviridae EVEs in Aedes aegypti and Aedes albopictus on Hainan Island and the Leizhou Peninsula, China

BACKGROUND

Hainan Island and the Leizhou Peninsula, the southernmost part of mainland China, are areas where Aedes aegypti and Ae. albopictus are sympatric and are also high-incidence areas of dengue outbreaks in China. Many studies have suggested that Aedes endogenous viral components (EVEs) are enriched in piRNA clusters which can silence incoming viral genomes. Investigation the EVEs present in the piRNA clusters associated with viral infection of Aedes mosquitoes in these regions may provide a theoretical basis for novel transmission-blocking vector control strategies.

METHODS

In this study, specific primers for endogenous Flaviviridae elements (EFVEs) and endogenous Rhabdoviridae elements (ERVEs) were used to detect the distribution of Zika virus infection associated EVEs in the genomes of individuals of the two Aedes mosquitoes. Genetic diversity of EVEs with a high detection rate was also analyzed.

RESULTS

The results showed that many EVEs associated with Zika virus infection were detected in both Aedes species, with the detection rates were 47.68% to 100% in Ae. aegypti and 36.15% to 92.31% in sympatric Ae. albopictus populations. EVEs detection rates in another 17 Ae. albopictus populations ranged from 29.39% to 89.85%. Genetic diversity analyses of the four EVEs (AaFlavi53, AaRha61, AaRha91 and AaRha100) of Ae. aegypti showed that each had high haplotype diversity and low nucleotide diversity. The number of haplotypes in AaFlavi53 was 8, with the dominant haplotype being Hap_1 and the other 7 haplotypes being further mutated from Hap_1 in a lineage direction. In contrast, the haplotype diversity of the other three ERVEs (AaRha61, AaRha91 and AaRha100) was more diverse and richer, with the haplotype numbers were 9, 15 and 19 respectively. In addition, these EVEs all showed inconsistent patterns of both population differentiation and dispersal compared to neutral evolutionary genes such as the Mitochondrial COI gene.

CONCLUSION

The EFVEs and ERVEs tested were present at high frequencies in the field Aedes mosquito populations. The haplotype diversity of the EFVE AaFlavi53 was relatively lower and the three ERVEs (AaRha61, AaRha91, AaRha100) were higher. None of the four EVEs could be indicative of the genetic diversity of the Ae. aegypti population. This study provided theoretical support for the use of EVEs to block arbovirus transmission, but further research is needed into the mechanisms by which these EVEs are antiviral to Aedes mosquitoes.

Gene flow throughout the evolutionary history of a colour polymorphic and generalist clownfish

Even seemingly homogeneous on the surface, the oceans display high environmental heterogeneity across space and time. Indeed, different soft barriers structure the marine environment, which offers an appealing opportunity to study various evolutionary processes such as population differentiation and speciation. Here, we focus on Amphiprion clarkii (Actinopterygii; Perciformes), the most widespread of clownfishes that exhibits the highest colour polymorphism. Clownfishes can only disperse during a short pelagic larval phase before their sedentary adult lifestyle, which might limit connectivity among populations, thus facilitating speciation events. Consequently, the taxonomic status of A. clarkii has been under debate. We used whole-genome resequencing data of 67 A. clarkii specimens spread across the Indian and Pacific Oceans to characterize the species' population structure, demographic history and colour polymorphism. We found that A. clarkii spread from the Indo-Pacific Ocean to the Pacific and Indian Oceans following a stepping-stone dispersal and that gene flow was pervasive throughout its demographic history. Interestingly, colour patterns differed noticeably among the Indonesian populations and the two populations at the extreme of the sampling distribution (i.e. Maldives and New Caledonia), which exhibited more comparable colour patterns despite their geographic and genetic distances. Our study emphasizes how whole-genome studies can uncover the intricate evolutionary past of wide-ranging species with diverse phenotypes, shedding light on the complex nature of the species concept paradigm.

Phylogenetic and Genomic Characterization of Whole Genome Sequences of Ocular Herpes Simplex Virus Type 1 Isolates Identifies Possible Virulence Determinants in Humans

Purpose

There are limited data on the prevalence and genetic diversity of herpes simplex virus type 1 (HSV-1) virulence genes in ocular isolates. Here, we sequenced 36 HSV-1 ocular isolates, collected by the Bascom Palmer Eye Institute, a university-based eye hospital, from three different ocular anatomical sites (conjunctiva, cornea, and eyelid) and carried out a genomic and phylogenetic analyses.

Methods

The PacBio Sequel II long read platform was used for genome sequencing. Phylogenetic analysis and genomic analysis were performed to help better understand genetic variability among common virulence genes in ocular herpetic disease.

Results

A phylogenetic network generated using the genome sequences of the 36 Bascom Palmer ocular isolates, plus 174 additional strains showed that ocular isolates do not group together phylogenetically. Analysis of the thymidine kinase and DNA polymerase protein sequences from the Bascom Palmer isolates showed multiple novel single nucleotide polymorphisms, but only one, BP-K14 encoded a known thymidine kinase acyclovir resistance mutation. An analysis of the multiple sequence alignment comprising the 51 total ocular isolates versus 159 nonocular strains detected several possible single nucleotide polymorphisms in HSV-1 genes that were found significantly more often in the ocular isolates. These genes included UL6, gM, VP19c, VHS, gC, VP11/12, and gG.

Conclusions

There does not seem to be a specific genetic feature of viruses causing ocular infection. The identification of novel and common recurrent polymorphisms may help to understand the drivers of herpetic pathogenicity and specific factors that may influence the virulence of ocular disease.

The complete mitochondrial genome of sleek unicornfish, Naso hexacanthus (Acanthuridae, Perciformes)

The complete mitochondrial DNA sequence of sleek unicornfish, Naso hexacanthus was first determined in this study. The complete mitogenome is 16,611 bp in length composed of 13 protein-coding genes, 2 ribosomal RNAs, 22 transfer RNAs, and a control region. The nucleotides consist of 33.8% A, 20.6% C, 25.0% G, 20.6% T. The gene order and direction are identical to those of N. lopezi and the species of Acanthuridae. The result would be useful to investigate genetic relationships among the species of Naso.

Mitochondrial marker implies fishery separate management units for spotted sardinella, Amblygaster sirm (Walbaum, 1792) populations in the South China Sea and the Andaman Sea

The spotted sardinella, Amblygaster sirm (Walbaum, 1792), is a commercial sardine commonly caught in Malaysia. Lack of management of these marine species in Malaysian waters could lead to overfishing and potentially declining fish stock populations. Therefore, sustainable management of this species is of paramount importance to ensure its longevity. As such, molecular information is vital in determining the A. sirm population structure and management strategy. In the present study, mitochondrial DNA Cytochrome b was sequenced from 10 A. sirm populations: the Andaman Sea (AS) (two), South China Sea (SCS) (six), Sulu Sea (SS) (one), and Celebes Sea (CS) (one). Accordingly, the intra-population haplotype diversity (Hd) was high (0.91-1.00), and nucleotide diversity (π) was low (0.002-0.009), which suggests a population bottleneck followed by rapid population growth. Based on the phylogenetic trees, minimum spanning network (MSN), population pairwise comparison, and F _ST,and supported by analysis of molecular variance (AMOVA) and spatial analysis of molecular variance (SAMOVA) tests, distinct genetic structures were observed (7.2% to 7.6% genetic divergence) between populations in the SCS and its neighboring waters, versus those in the AS. Furthermore, the results defined A. sirm stock boundaries and evolutionary between the west and east coast (which shares the same waters as western Borneo) of Peninsular Malaysia. In addition, genetic homogeneity was revealed throughout the SCS, SS, and CS based on the non-significant F _STpairwise comparisons. Based on the molecular evidence, separate management strategies may be required for A. sirm of the AS and the SCS, including its neighboring waters.

Molecular phylogeography reveals multiple Pleistocene divergence events in estuarine crabs from the tropical West Pacific

Due to the lack of visible barriers to gene flow, it was a long-standing assumption that marine coastal species are widely distributed, until molecular studies revealed geographically structured intraspecific genetic differentiation in many taxa. Historical events of sea level changes during glacial periods are known to have triggered sequential disjunctions and genetic divergences among populations, especially of coastal organisms. The Parasesarma bidens species complex so far includes three named plus potentially cryptic species of estuarine brachyuran crabs, distributed along East to Southeast Asia. The aim of the present study is to address phylogeography and uncover real and hidden biological diversity within this complex, by revealing the underlying genetic structure of populations and species throughout their distribution ranges from Japan to West Papua, with a comparison of mitochondrial COX1 and 16S rRNA gene sequences. Our results reveal that the P. bidens species complex consists of at least five distinct clades, resulting from four main cladogenesis events during the mid to late Pleistocene. Among those clades, P. cricotum and P. sanguimanus are recovered as monophyletic taxa. Geographically restricted endemic clades are encountered in southeastern Indonesia, Japan and China respectively, whereas the Philippines and Taiwan share two clades. As individuals of the Japanese clade can also be found in Taiwan, we provide evidence of a third lineage and the occurrence of a potential cryptic species on this island. Ocean level retreats during Pleistocene ice ages and present oceanic currents appear to be the main triggers for the divergences of the five clades that are here addressed as the P. bidens complex. Secondary range expansions converted Taiwan into the point of maximal overlap, sharing populations with Japan and the Philippines, but not with mainland China.

Comparative Phylogeography and Phylogeny of Pennah Croakers (Teleostei: Sciaenidae) in Southeast Asian Waters

A broad-scale comparative phylogeographic and phylogenetic study of pennah croakers, mainly Pennahia anea, P. macrocephalus, and P. ovata was conducted to elucidate the mechanisms that may have driven the diversification of marine organisms in Southeast Asian waters. A total of 316 individuals from the three species, and an additional eight and six individuals of P. argentata and P. pawak were employed in this study. Two genetically divergent lineages each of P. argentata and P. anea (lineages L1 and L2) were respectively detected from the analyses based on mitochondrial cytochrome b gene data. Historical biogeography analysis with a multi-gene dataset revealed that Pennahia species most likely originated in the South China Sea and expanded into the eastern Indian Ocean, East China Sea, and northwestern Pacific Ocean through three separate range expansions. The main diversifications of Pennahia species occurred during Miocene and Pliocene periods, and the occurrences of lineage divergences within P. anea and P. argentata were during the Pleistocene, likely as a consequence of cyclical glaciations. The population expansions that occurred after the sea level rise might be the reason for the population homogeneity observed in P. macrocephalus and most P. anea L2 South China Sea populations. The structure observed between the two populations of P. ovata, and the restricted distributions of P. anea lineage L1 and P. ovata in the eastern Indian Ocean, might have been hampered by the northward flowing ocean current at the Malacca Strait and by the distribution of coral reefs or rocky bottoms. While our results support S. Ekman’s center-of-origin hypothesis taking place in the South China Sea, the Malacca Strait serving as the center of overlap is a supplementary postulation for explaining the present-day high diversity of pennah croakers centered in these waters.

A paradigm for the recognition of cryptic trematode species in tropical Indo-west Pacific fishes: the problematic genus Preptetos (Trematoda: Lepocreadiidae)

Molecular data have transformed approaches to trematode taxonomy by providing objective evidence for the delineation of species. However, although the data are objective, the interpretation of these data regarding species boundaries is subjective, especially when different markers conflict. Conserved markers can lead to an underestimation of richness and those used for finer species delineation have the capacity to inflate species recognition, perhaps unrealistically. Here we examine molecular and morphological evidence for species recognition in an especially confusing system, the lepocreadiid genus Preptetos Pritchard, 1960 in acanthuriform fishes of the tropical Indo-west Pacific. We consider species boundaries within this genus based on combined data (ITS2 and 28S rDNA; cox1 mtDNA and morphometrics) for substantial new collections. Delineation of species using only morphological data suggest fewer species than analysis of the sequence data; the latter suggests the presence of potential cryptic species and analysis of different markers suggests the presence of differing numbers of species. We conclude that an integrative interpretation creates the most satisfying taxonomic hypothesis. In the light of the new data, we have chosen and propose a model of trematode species recognition that demands reciprocal monophyly in the most discriminating available molecular marker plus distinction in morphology or host distribution. By invoking these criteria, we distinguish eight species in our new tropical Indo-west Pacific collections. Six of these are new (Preptetos allocaballeroi n. sp., Preptetos paracaballeroi n. sp., Preptetos pearsoni n. sp., Preptetos prudhoei n. sp., Preptetos quandamooka n. sp. and Preptetos zebravaranus n. sp.) and we continue to recognise Preptetos cannoni Barker, Bray & Cribb, 1993 and Preptetos laguncula Bray and Cribb, 1996. Notably; two of the new species, P. allocaballeroi n. sp. and P. paracaballeroi n. sp., are morphologically cryptic relative to each other. Our criteria lead us to recognise, as species, populations with unvarying morphology and similar host relationships but which may have a complex population structure over their range. In our view, this paradigm has the capacity to render tractable the interpretation of the species status of the huge trematode fauna of the tropical Indo-west Pacific.

Population Structure, Genetic Connectivity, and Signatures of Local Adaptation of the Giant Black Tiger Shrimp (Penaeus monodon) throughout the Indo-Pacific Region

The giant black tiger shrimp (Penaeus monodon) is native to the Indo-Pacific and is the second most farmed penaeid shrimp species globally. Understanding genetic structure, connectivity, and local adaptation among Indo-Pacific black tiger shrimp populations is important for informing sustainable fisheries management and aquaculture breeding programs. Population genetic and outlier detection analyses were undertaken using 10,593 genome-wide single nucleotide polymorphisms (SNPs) from 16 geographically disparate Indo-Pacific P. monodon populations. Levels of genetic diversity were highest for Southeast Asian populations and were lowest for Western Indian Ocean (WIO) populations. Both neutral (n = 9,930) and outlier (n = 663) loci datasets revealed a pattern of strong genetic structure of P. monodon corresponding with broad geographical regions and clear genetic breaks among samples within regions. Neutral loci revealed seven genetic clusters and the separation of Fiji and WIO clusters from all other clusters, whereas outlier loci revealed six genetic clusters and high genetic differentiation among populations. The neutral loci dataset estimated five migration events that indicated migration to Southeast Asia from the WIO, with partial connectivity to populations in both oceans. We also identified 26 putatively adaptive SNPs that exhibited significant Pearson correlation (P < 0.05) between minor allele frequency and maximum or minimum sea surface temperature. Matched transcriptome contig annotations suggest putatively adaptive SNPs involvement in cellular and metabolic processes, pigmentation, immune response, and currently unknown functions. This study provides novel genome-level insights that have direct implications for P. monodon aquaculture and fishery management practices.

Conspicuous genetic similarity within a widely distributed and newly described species of Parasesarma De Man, 1895 from Western Pacific oceanic islands, with notes on the allied P. calypso group (Crustacea : Brachyura : Sesarmidae)

Some brachyuran crab species of the Western Pacific appear to be widespread throughout the region and distributed across a large geographic area, without obvious phylogeographic structuring. In the present study, we describe a new species of Parasesarma that appears to be restricted to Western Pacific islands (so far Guam, Palau, Vanuatu, Fiji, Wallis and New Caledonia). Comparisons of partial sequences of the COX1 gene show that individuals of this species, though from relatively isolated and widely separated islands, are monophyletic and, surprisingly, genetically uniform. These results give credence to the hypothesis that these oceanic islands serve as ‘stepping stones’ for the current-mediated dispersal and genetic homogenisation of coastal–littoral marine species. Morphologically, the new species differs most significantly from similar congeners in the tuberculation pattern of the chelar dactyli, whereas genetically it is markedly divergent from other morphologically similar species of Parasesarma, with a minimum COX1 p-distance of 6.9%. With such evidence, the new species is here formally described as Parasesarma daviei sp. nov. It is the fifth species of Parasesarma reported from oceanic islands of the Western Pacific. Compared to other congeners, P. daviei sp. nov. shows a close relationship with a clade including P. calypso. Therefore, P. calypso (De Man, 1895), and three of its former subspecies or varieties were subjected to a closer examination and are here rediagnosed and illustrated. In consequence, we suggest full species status for P. kuekenthali (De Man, 1902), P. lanchesteri (Tweedie, 1936), and P. ellenae (Pretzmann, 1968).

Population Genetic Structure of A Marine Pelagic Egg Producer and Popular Marine Aquarium Species, the Mandarinfish Synchiropus splendidus

The mandarinfish Synchiropus splendidus is extensively collected in Southeast Asia (mainly in the Philippines) and highly favoured for the marine aquarium trade. Males are more popular than females for their large first dorsal fins and the fishery is not managed. To examine possible population replenishment dynamics arising as a result of selective fishing, the effects of sex-selective fishing on sex ratios and population connectivity were considered. This study determined the sex ratios and analyzed the population genetic structure from mandarinfish collected at six locations: one from Palau, where the species is not exploited, and five from Bohol in the Philippines, where the species has long been heavily fished. The findings reported very low male to female ratios (0.12 to 0.30) from four of the five locations in Bohol, with relatively more males to females in the specimens collected from Palau (2.3). The analyses from allozymes (43 alleles from 10 loci) and microsatellites (118 alleles from 5 loci) revealed that genetic connectivity was high among the five locations in the Philippines as well as with the specimens collected from the more-distant Palau. The genetic homogeneity observed across the geographical range considered is inconsistent with the hypothesized limited dispersal ability of the species and could be explained by recent species range expansion associated with sea level rise in the region. The results suggest that the present genetic structure, at least in the geographic region considered, may not be determined by current patterns of gene flow, but may, instead, be driven by recent sea-level changes associated with periods of glaciation. Caution is suggested to ensure that heavily localized fishing does not produce excessively biased adult sex ratios.

New microsatellite DNA markers to resolve population structure of the convict surgeonfish, Acanthurus triostegus, and cross-species amplifications on thirteen other Acanthuridae

Microsatellites are widely used to investigate connectivity and parentage in marine organisms. Despite surgeonfish (Acanthuridae) being dominant members of most reef fish assemblages and having an ecological key role in coral reef ecosystems, there is limited information describing the scale at which populations are connected and very few microsatellite markers have been screened. Here, we developed fourteen microsatellite markers for the convict surgeonfish Acanthurus triostegus with the aim to infer its genetic connectivity throughout its distribution range. Genetic diversity and variability was tested over 152 fishes sampled from four locations across the Indo-Pacific: Mayotte (Western Indian Ocean), Papua New Guinea and New Caledonia (Southwestern Pacific Ocean), and Moorea (French Polynesia). Over all locations, the number of alleles per locus varied from 5 to 24 per locus, and expected heterozygosities ranged from 0.468 to 0.941. Significant deviations from Hardy-Weinberg equilibrium were detected for two loci in two to three locations and were attributed to the presence of null alleles. These markers revealed for the first time a strong and significant distinctiveness between Indian Ocean and Pacific Ocean A. triostegus populations. We further conducted cross-species amplification tests in 13 Pacific congener species to investigate the possible use of these microsatellites in other Acanthuridae species. The phylogenetic placement of A. triostegus branching off from the clade containing nearly all Acanthurus + Ctenochaetus species likely explain the rather good transferability of these microsatellite markers towards other Acanthuridae species. This suggests that this fourteen new microsatellite loci will be helpful tools not only for inferring population structure of various surgeonfish but also to clarify systematic relationships among Acanthuridae.

Coping with Pleistocene climatic fluctuations: Demographic responses in remote endemic reef fishes

Elucidating demographic history during the settlement of ecological communities is crucial for properly inferring the mechanisms that shape patterns of species diversity and their persistence through time. Here, we used genomic data and coalescent-based approaches to elucidate for the first time the demographic dynamics associated with the settlement by endemic reef fish fauna of one of the most remote peripheral islands of the Pacific Ocean, Rapa Nui (Easter Island). We compared the demographic history of nine endemic species in order to explore their demographic responses to Pleistocene climatic fluctuations. We found that species endemic to Rapa Nui share a common demographic history, as signatures of population expansions were retrieved for almost all of the species studied here, and synchronous demographic expansions initiated during the last glacial period were recovered for more than half of the studied species. These results suggest that eustatic fluctuations associated with Milankovitch cycles have played a central role in species demographic histories and in the final stage of the community assembly of many Rapa Nui reef fishes. Specifically, sea level lowstands resulted in the maximum reef habitat extension for Rapa Nui endemic species; we discuss the potential role of seamounts in allowing endemic species to cope with Pleistocene climatic fluctuations, and we highlight the importance of local historical processes over regional ones. Overall, our results shed light on the mechanisms by which endemism arises and is maintained in peripheral reef fish fauna.

Phylogenetic relationships and cryptic species in the genus Sthenoteuthis (Cephalopoda: Ommastrephidae) in the South China Sea

Large squids of the genus Sthenoteuthis are commercially relevant species that include two truly oceanic squids. They are large nektonic predators being widely distributed throughout tropical and subtropical waters of the Atlantic and Indo-Pacific Ocean. The present study investigates different morphs varying in size at maturity, and assesses the genetic divergence in Sthenotheutis in relation to geographic patterns in the South China Sea. We obtained sequences using a mitochondrial (cytochrome c oxidase subunit I) and a nuclear (Histone H3) gene marker from 111 individuals in 23 locations of the South China Sea. In combination with sequences available in public databases, we performed tests on DNA taxonomy, mostly based on molecular phylogenies. Our results suggest that the genus Sthenoteuthis includes at least three species. The Indo-Pacific purpleback squid Sthenoteuthis oualaniensis contains at least two genetically distinct lineages that can be considered separate species, a dwarf species and a medium-sized species, separated by both the mitochondrial marker and the more conserved nuclear marker. We also assessed whether the few cases of mitonuclear discordance could be the result of genetic introgression and past hybridization or incongruence lineage sorting. The medium-sized species is more widely distributed and dominant in the South China Sea than the dwarf species. The medium-sized species inhabits the whole South China Sea, whereas the dwarf species is restricted to the equatorial waters of the South China Sea. The medium-sized species has two further genetic clades, one distributed in the East Pacific Ocean and the other in the South China Sea. This high level of genetic differentiation is in agreement with the discriminant analysis on the morphological measurements, clearly separating the dwarf and medium-sized species, indicating the presence of a complex of pseudo-cryptic species in S. oualaniensis, clearly identifiable by differences in DNA sequences and in body size, and statistically differentiated in their body measurements.

Population structure and genetic diversity of Aedes aegypti and Aedes albopictus in Penang as revealed by mitochondrial DNA cytochrome oxidase I

The population genetics study is crucial as it helps in understanding the epidemiological aspects of dengue and help improving a vector control measures. This research aims to investigate the population genetics structure of two common species of Aedes mosquitoes in Penang; Aedes aegypti and Aedes albopictus using Cytochrome Oxidase I (COI) mitochondrial DNA (mtDNA) marker. Molecular investigations were derived from 440 bp and 418 bp mtDNA COI on 125 and 334 larvae of Aedes aegypti and Aedes albopictus respectively, from 32 locations in Penang. All samples were employed in the BLASTn for species identification. The haplotype diversity, nucleotide diversity, neutrality test and mismatch distribution analysis were conducted in DnaSP version 5.10.1. AMOVA analysis was conducted in ARLEQUIN version 3.5 and the phylogenetic reconstructions based on maximum likelihood (ML) and neighbor-joining (NJ) methods were implemented in MEGA X. The relationships among haplotypes were further tested by creating a minimum spanning tree using Network version 4.6.1. All samples were genetically identified and clustered into six distinct species. Among the species, Ae. albopictus was the most abundant (67.2%), followed by Ae. aegypti (25.2%) and the rest were counted for Culex sp. and Toxorhynchites sp. Both Ae. aegypti and Ae. albopictus show low nucleotide diversity (π) and high haplotype diversity (h), while the neutrality test shows a negative value in most of the population for both species. There are a total of 39 and 64 haplotypes recorded for Ae. aegypti and Ae. albopictus respectively. AMOVA analysis revealed that most of the variation occurred within population for both species. Mismatch distribution analysis showed bimodal characteristic of population differentiation for Ae. aegypti but Ae. albopictus showed unimodal characteristics of population differentiation. Genetic distance based on Tamura-Nei parameter showed low genetic divergent within population and high genetic divergent among population for both species. The maximum likelihood tree showed no obvious pattern of population genetic structure for both Ae. aegypti and Ae. albopictus from Penang and a moderate to high bootstrap values has supported this conclusion. The minimum spanning network for Ae. aegypti and Ae. albopictus showed five and three dominant haplotypes respectively, which indicates a mixture of haplotypes from the regions analysed. This study revealed that there is no population genetic structure exhibited by both Ae. aegypti and Ae. albopictus in Penang. Mutation has occurred rapidly in both species and this will be challenging in controlling the populations. However, further analysis needed to confirm this statement.

RADseq analyses reveal concordant Indian Ocean biogeographic and phylogeographic boundaries in the reef fish Dascyllus trimaculatus

Population genetic analysis is an important tool for estimating the degree of evolutionary connectivity in marine organisms. Here, we investigate the population structure of the three-spot damselfish Dascyllus trimaculatus in the Red Sea, Arabian Sea and Western Indian Ocean, using 1174 single nucleotide polymorphisms (SNPs). Neutral loci revealed a signature of weak genetic differentiation between the Northwestern (Red Sea and Arabian Sea) and Western Indian Ocean biogeographic provinces. Loci potentially under selection (outlier loci) revealed a similar pattern but with a much stronger signal of genetic structure between regions. The Oman population appears to be genetically distinct from all other populations included in the analysis. While we could not clearly identify the mechanisms driving these patterns (isolation, adaptation or both), the datasets indicate that population-level divergences are largely concordant with biogeographic boundaries based on species composition. Our data can be used along with genetic connectivity of other species to identify the common genetic breaks that need to be considered for the conservation of biodiversity and evolutionary processes in the poorly studied Western Indian Ocean region.

Recombination contributes to population diversification in the polyploid intestinal symbiont Epulopiscium sp. type B

Epulopiscium sp. type B (Lachnospiraceae) is an exceptionally large, highly polyploid, intestinal symbiont of the coral reef dwelling surgeonfish Naso tonganus. These obligate anaerobes do not form mature endospores and reproduce solely through the production of multiple intracellular offspring. This likely makes them dependent on immediate transfer to a receptive host for dispersal. During reproduction, only a small proportion of Epulopiscium mother-cell DNA is inherited. To explore the impact of this unusual viviparous lifestyle on symbiont population dynamics, we investigated Epulopiscium sp. type B and their fish hosts collected over the course of two decades, at island and reef habitats near Lizard Island, Australia. Using multi-locus sequence analysis, we found that recombination plays an important role in maintaining diversity of these symbionts and yet populations exhibit linkage disequilibrium (LD). Symbiont populations showed spatial but not temporal partitioning. Surgeonfish are long-lived and capable of traveling long distances, yet the population structures of Epulopiscium suggest that adult fish tend to not roam beyond a limited locale. Codiversification analyses and traits of this partnership suggest that while symbionts are obligately dependent on their host, the host has a facultative association with Epulopiscium. We suggest that congression of unlinked markers contributes to LD estimates in this and other recombinant populations of bacteria. The findings here inform our understanding of evolutionary processes within intestinal Lachnospiraceae populations.

Contrasting population genetic structure in three aggregating groupers (Percoidei: Epinephelidae) in the Indo-West Pacific: the importance of reproductive mode

Background

Understanding the factors shaping population genetic structure is important for evolutionary considerations as well as for management and conservation. While studies have revealed the importance of palaeogeographic changes in shaping phylogeographic patterns in multiple marine fauna, the role of reproductive behaviour is rarely considered in reef fishes. We investigated the population genetics of three commercially important aggregating grouper species in the Indo-West Pacific, namely the camouflage grouper Epinephelus polyphekadion, the squaretail coral grouper Plectropomus areolatus, and the common coral trout P. leopardus, with similar life histories but distinct spatio-temporal characteristics in their patterns of forming spawning aggregations.

Results

By examining their mitochondrial control region and 9–11 microsatellite markers, we found an overarching influence of palaeogeographic events in the population structure of all species, with genetic breaks largely coinciding with major biogeographic barriers. The divergence time of major lineages in these species coincide with the Pleistocene glaciations. Higher connectivity is evident in E. polyphekadion and P. areolatus that assemble in larger numbers at fewer spawning aggregations and in distinctive offshore locations than in P. leopardus which has multiple small, shelf platform aggregations.

Conclusions

While palaeogeographic events played an important role in shaping the population structure of the target species, the disparity in population connectivity detected may be partly attributable to differences in their reproductive behaviour, highlighting the need for more investigations on this characteristic and the need to consider reproductive mode in studies of connectivity and population genetics.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1284-0) contains supplementary material, which is available to authorized users.

Genetic population structure of the convict surgeonfish Acanthurus triostegus: a phylogeographic reassessment across its range

This study investigates the genetic population structure and connectivity of Acanthurus triostegus in five Indo-Pacific biogeographic regions (western and eastern Indian Ocean, western, central and eastern Pacific Ocean), using a mitochondrial DNA marker spanning the ATPase8 and ATPase6 gene regions. In order to assess the phylogeography and genetic population structure of A. triostegus across its range, 35 individuals were sampled from five localities in the western Indian Ocean and complemented with 227 sequences from two previous studies. Results from the overall analysis of molecular variance (AMOVA) without a priori grouping showed evidence of significant differentiation in the Indo-Pacific, with 25 (8.3%) out of 300 pairwise Φ_ST comparisons being significant. However, the hierarchical AMOVA grouping of Indian and Pacific Ocean populations failed to support the vicariance hypothesis, showing a lack of a genetic break between the two ocean basins. Instead, the correlation between pairwise Φ_ST values and geographic distance showed that dispersal of A. triostegus in the Indo-Pacific Ocean follows an isolation-by-distance model. Three haplogroups could be deduced from the haplotype network and phylogenetic tree, with haplogroup 1 and 2 dominating the Indian and the Pacific Ocean, respectively, while haplogroup 3 exclusively occurring in the Hawaiian Archipelago of the central Pacific Ocean.

Geography and island geomorphology shape fish assemblage structure on isolated coral reef systems

We quantify the relative importance of multi-scale drivers of reef fish assemblage structure on isolated coral reefs at the intersection of the Indian and Indo-Pacific biogeographical provinces. Large (>30 cm), functionally-important and commonly targeted species of fish, were surveyed on the outer reef crest/front at 38 coral reef sites spread across three oceanic coral reef systems (i.e. Christmas Island, Cocos (Keeling) Islands and the Rowley Shoals), in the tropical Indian Ocean (c. 1.126 x 106 km2). The effects of coral cover, exposure, fishing pressure, lagoon size and geographical context, on observed patterns of fish assemblage structure were modelled using Multivariate Regression Trees. Reef fish assemblages were clearly separated in space with geographical location explaining ~53 % of the observed variation. Lagoon size, within each isolated reef system was an equally effective proxy for explaining fish assemblage structure. Among local-scale variables, 'distance from port', a proxy for the influence of fishing, explained 5.2% of total variation and separated the four most isolated reefs from Cocos (Keeling) Island, from reefs with closer boating access. Other factors were not significant. Major divisions in assemblage structure were driven by sister taxa that displayed little geographical overlap between reef systems and low abundances of several species on Christmas Island corresponding to small lagoon habitats. Exclusion of geographical context from the analysis resulted in local processes explaining 47.3% of the variation, highlighting the importance of controlling for spatial correlation to understand the drivers of fish assemblage structure. Our results suggest reef fish assemblage structure on remote coral reef systems in the tropical eastern Indian Ocean reflects a biogeographical legacy of isolation between Indian and Pacific fish faunas and geomorphological variation within the region, more than local fishing pressure or reef condition. Our findings re-emphasise the importance that historical processes play in structuring contemporary biotic communities.

Sea surface currents and geographic isolation shape the genetic population structure of a coral reef fish in the Indian Ocean

In this contribution, we determine the genetic population structure in the Skunk Clownfish (Amphiprion akallopsisos) across the Indian Ocean, and on a smaller geographic scale in the Western Indian Ocean (WIO). Highly restricted gene flow was discovered between populations on either side of the Indian Ocean using the control region as a mitochondrial marker (mtDNA). We verify this conclusion using 13 microsatellite markers and infer fine scale genetic structuring within the WIO. In total 387 samples from 21 sites were analysed using mtDNA and 13 microsatellite loci. Analysis included estimation of genetic diversity and population differentiation. A haplotype network was inferred using mtDNA. Nuclear markers were used in Bayesian clustering and a principal component analysis. Both markers confirmed strong genetic differentiation between WIO and Eastern Indian Ocean (EIO) populations, and a shallower population structure among Malagasy and East African mainland populations. Limited gene flow across the Mozambique Channel may be explained by its complex oceanography, which could cause local retention of larvae, limiting dispersal between Madagascar and the East African coast. Two other potential current-mediated barriers to larval dispersal suggested in the WIO, the split of the SEC at approximately 10° S and the convergence of the Somali Current with the East African Coast Current at approximately 3° S, were not found to form a barrier to gene flow in this species.

Evidence of shallow mitochondrial divergence in the slender armorhead, Pentaceros wheeleri (Pisces, Pentacerotidae) from the Emperor Seamount Chain

Competitive overexploitation of the slender armorhead, Pentaceros wheeleri, a deep-sea fish inhabiting the Emperor Seamount Chain caused a serious population decline. Therefore, it is urgently necessary to clarify its genetic diversity and connectivity among populations of P. wheeleri for appropriate stock management. For this, we compared 677 base pairs (bp) of mitochondrial (mt) DNA control region (CR) sequences of 80 individuals from three seamounts (the Milwaukee, Kinmei, and Koko Seamounts) in the southern part of the Emperor Seamount Chain. Contrary to our expectation, the three seamount populations showed high genetic diversity, not yet reflecting effects from the recent population decline or due to mixed two clades. Analysis of molecular variance indicated no significant genetic differentiation between seamount populations, however, the neighbour-joining tree and minimum spanning network showed significant separation into two clades (K2P distance= 1.2-3.2%, ϕ_st = 0.5739, p < .05) regardless of seamount. The divergence time between the two clades was estimated to be 0.3-0.8 Mya, during the period of Pleistocene glacial cycles, suggesting that associated environmental changes and the unique life history traits of Pentaceros spp. might have resulted in the initiation of divergence between these clades.

Origins of the Greenland shark (Somniosus microcephalus): Impacts of ice-olation and introgression

Herein, we use genetic data from 277 sleeper sharks to perform coalescent‐based modeling to test the hypothesis of early Quaternary emergence of the Greenland shark (Somniosus microcephalus) from ancestral sleeper sharks in the Canadian Arctic‐Subarctic region. Our results show that morphologically cryptic somniosids S. microcephalus and Somniosus pacificus can be genetically distinguished using combined mitochondrial and nuclear DNA markers. Our data confirm the presence of genetically admixed individuals in the Canadian Arctic and sub‐Arctic, and temperate Eastern Atlantic regions, suggesting introgressive hybridization upon secondary contact following the initial species divergence. Conservative substitution rates fitted to an Isolation with Migration (IM) model indicate a likely species divergence time of 2.34 Ma, using the mitochondrial sequence DNA, which in conjunction with the geographic distribution of admixtures and Pacific signatures likely indicates speciation associated with processes other than the closing of the Isthmus of Panama. This time span coincides with further planetary cooling in the early Quaternary period followed by the onset of oscillating glacial‐interglacial cycles. We propose that the initial S. microcephalus–S. pacificus split, and subsequent hybridization events, were likely associated with the onset of Pleistocene glacial oscillations, whereby fluctuating sea levels constrained connectivity among Arctic oceanic basins, Arctic marginal seas, and the North Atlantic Ocean. Our data demonstrates support for the evolutionary consequences of oscillatory vicariance via transient oceanic isolation with subsequent secondary contact associated with fluctuating sea levels throughout the Quaternary period—which may serve as a model for the origins of Arctic marine fauna on a broad taxonomic scale.

Multispecies genetic objectives in spatial conservation planning

Growing threats to biodiversity and global alteration of habitats and species distributions make it increasingly necessary to consider evolutionary patterns in conservation decision making. Yet, there is no clear-cut guidance on how genetic features can be incorporated into conservation-planning processes, despite multiple molecular markers and several genetic metrics for each marker type to choose from. Genetic patterns differ between species, but the potential tradeoffs among genetic objectives for multiple species in conservation planning are currently understudied. We compared spatial conservation prioritizations derived from 2 metrics of genetic diversity (nucleotide and haplotype diversity) and 2 metrics of genetic isolation (private haplotypes and local genetic differentiation) in mitochondrial DNA of 5 marine species. We compared outcomes of conservation plans based only on habitat representation with plans based on genetic data and habitat representation. Fewer priority areas were selected for conservation plans based solely on habitat representation than on plans that included habitat and genetic data. All 4 genetic metrics selected approximately similar conservation-priority areas, which is likely a result of prioritizing genetic patterns across a genetically diverse array of species. Largely, our results suggest that multispecies genetic conservation objectives are vital to creating protected-area networks that appropriately preserve community-level evolutionary patterns.

Wandering behaviour prevents inter and intra oceanic speciation in a coastal pelagic fish

Small pelagic fishes have the ability to disperse over long distances and may present complex evolutionary histories. Here, Old World Anchovies (OWA) were used as a model system to understand genetic patterns and connectivity of fish between the Atlantic and Pacific basins. We surveyed 16 locations worldwide using mtDNA and 8 microsatellite loci for genetic parameters, and mtDNA (cyt b; 16S) and nuclear (RAG1; RAG2) regions for dating major lineage-splitting events within Engraulidae family. The OWA genetic divergences (0-0.4%) are compatible with intra-specific divergence, showing evidence of both ancient and contemporary admixture between the Pacific and Atlantic populations, enhanced by high asymmetrical migration from the Pacific to the Atlantic. The estimated divergence between Atlantic and Pacific anchovies (0.67 [0.53-0.80] Ma) matches a severe drop of sea temperature during the Günz glacial stage of the Pleistocene. Our results support an alternative evolutionary scenario for the OWA, suggesting a coastal migration along south Asia, Middle East and eastern Africa continental platforms, followed by the colonization of the Atlantic via the Cape of the Good Hope.

Phylogeography of the sergeants Abudefduf sexfasciatus and A. vaigiensis reveals complex introgression patterns between two widespread and sympatric Indo-West Pacific reef fishes

On evolutionary timescales, sea level oscillations lead to recurrent spatio-temporal variation in species distribution and population connectivity. In this situation, applying classical concepts of biogeography is challenging yet necessary to understand the mechanisms underlying biodiversity in highly diverse marine ecosystems such as coral reefs. We aimed at studying the outcomes of such complex biogeographic dynamics on reproductive isolation by sampling populations across a wide spatial range of a species-rich fish genus: the sergeants (Pomacentridae: Abudefduf). We generated a mutlilocus data set that included ten morpho-species from 32 Indo-West Pacific localities. We observed a pattern of mito-nuclear discordance in two common and widely distributed species: Abudefduf sexfasciatus and Abudefduf vaigiensis. The results showed three regional sublineages (Indian Ocean, Coral Triangle region, western Pacific) in A. sexfasciatus (0.6-1.5% divergence at cytb). The other species, A. vaigiensis, is polyphyletic and consists of three distinct genetic lineages (A, B and C) (9% divergence at cytb) whose geographic ranges overlap. Although A. vaigiensis A and A. sexfasciatus were found to be distinct based on nuclear information, A. vaigiensis A was found to be nested within A. sexfasciatus in the mitochondrial gene tree. A. sexfasciatus from the Coral Triangle region and A. vaigiensis A were not differentiated from each other at the mitochondrial locus. We then used coalescent-based simulation to characterize a spatially widespread but weak gene flow between the two species. We showed that these fishes are good candidates to investigate the evolutionary complexity of the discrepancies between phenotypic and genetic similarity in closely related species.

Population expansions dominate demographic histories of endemic and widespread Pacific reef fishes

Despite the unique nature of endemic species, their origin and population history remain poorly studied. We investigated the population history of 28 coral reef fish species, close related, from the Gambier and Marquesas Islands, from five families, with range size varying from widespread to small-range endemic. We analyzed both mitochondrial and nuclear sequence data using neutrality test and Bayesian analysis (EBSP and ABC). We found evidence for demographic expansions for most species (24 of 28), irrespective of range size, reproduction strategy or archipelago. The timing of the expansions varied greatly among species, from 8,000 to 2,000,000 years ago. The typical hypothesis for reef fish that links population expansions to the Last Glacial Maximum fit for 14 of the 24 demographic expansions. We propose two evolutionary processes that could lead to expansions older than the LGM: (a) we are retrieving the signature of an old colonization process for widespread, large-range endemic and paleoendemic species or (b) speciation; the expansion reflects the birth of the species for neoendemic species. We show for the first time that the demographic histories of endemic and widespread reef fish are not distinctly different and suggest that a number of processes drive endemism.

Population Genetic Structure of the Tropical Two-Wing Flyingfish (Exocoetus volitans)

Delineating populations of pantropical marine fish is a difficult process, due to widespread geographic ranges and complex life history traits in most species. Exocoetus volitans, a species of two-winged flyingfish, is a good model for understanding large-scale patterns of epipelagic fish population structure because it has a circumtropical geographic range and completes its entire life cycle in the epipelagic zone. Buoyant pelagic eggs should dictate high local dispersal capacity in this species, although a brief larval phase, small body size, and short lifespan may limit the dispersal of individuals over large spatial scales. Based on these biological features, we hypothesized that E. volitans would exhibit statistically and biologically significant population structure defined by recognized oceanographic barriers. We tested this hypothesis by analyzing cytochrome b mtDNA sequence data (1106 bps) from specimens collected in the Pacific, Atlantic and Indian oceans (n = 266). AMOVA, Bayesian, and coalescent analytical approaches were used to assess and interpret population-level genetic variability. A parsimony-based haplotype network did not reveal population subdivision among ocean basins, but AMOVA revealed limited, statistically significant population structure between the Pacific and Atlantic Oceans (ΦST = 0.035, p<0.001). A spatially-unbiased Bayesian approach identified two circumtropical population clusters north and south of the Equator (ΦST = 0.026, p<0.001), a previously unknown dispersal barrier for an epipelagic fish. Bayesian demographic modeling suggested the effective population size of this species increased by at least an order of magnitude ~150,000 years ago, to more than 1 billion individuals currently. Thus, high levels of genetic similarity observed in E. volitans can be explained by high rates of gene flow, a dramatic and recent population expansion, as well as extensive and consistent dispersal throughout the geographic range of the species.

The changing dynamics of coral reef science in Arabia

Six percent of the world's coral reefs occur around the Arabian Peninsula, providing a valuable ecological, economic and scientific resource for the nations bordering its shores. We provide the first region-wide assessment of the current status and historical trends in coral reef research, focusing on research in the Red Sea, Arabian Sea, and Arabian Gulf. In total, 633 regional reef publications have been produced since the 1930s, covering a wide variety of themes and taxa. Our results show a great deal of commonality in regional reef research, but also highlight important differences in research among the various seas as well as knowledge gaps that represent opportunities for future research. A regionally-integrated approach to future research is essential. There is a growing need for large-scale research to guide management of reefs and their stressors, as these operate at much larger scales than the national borders within which most research currently occurs.

Population Connectivity Measures of Fishery-Targeted Coral Reef Species to Inform Marine Reserve Network Design in Fiji

Coral reef fish serve as food sources to coastal communities worldwide, yet are vulnerable to mounting anthropogenic pressures like overfishing and climate change. Marine reserve networks have become important tools for mitigating these pressures, and one of the most critical factors in determining their spatial design is the degree of connectivity among different populations of species prioritized for protection. To help inform the spatial design of an expanded reserve network in Fiji, we used rapidly evolving mitochondrial genes to investigate connectivity patterns of three coral reef species targeted by fisheries in Fiji: Epinephelus merra (Serranidae), Halichoeres trimaculatus (Labridae), and Holothuria atra (Holothuriidae). The two fish species, E. merra and Ha. trimaculatus, exhibited low genetic structuring and high amounts of gene flow, whereas the sea cucumber Ho. atra displayed high genetic partitioning and predominantly westward gene flow. The idiosyncratic patterns observed among these species indicate that patterns of connectivity in Fiji are likely determined by a combination of oceanographic and ecological characteristics. Our data indicate that in the cases of species with high connectivity, other factors such as representation or political availability may dictate where reserves are placed. In low connectivity species, ensuring upstream and downstream connections is critical.

Northern Bobwhite (Colinus virginianus) Mitochondrial Population Genomics Reveals Structure, Divergence, and Evidence for Heteroplasmy

Herein, we evaluated the concordance of population inferences and conclusions resulting from the analysis of short mitochondrial fragments (i.e., partial or complete D-Loop nucleotide sequences) versus complete mitogenome sequences for 53 bobwhites representing six ecoregions across TX and OK (USA). Median joining (MJ) haplotype networks demonstrated that analyses performed using small mitochondrial fragments were insufficient for estimating the true (i.e., complete) mitogenome haplotype structure, corresponding levels of divergence, and maternal population history of our samples. Notably, discordant demographic inferences were observed when mismatch distributions of partial (i.e., partial D-Loop) versus complete mitogenome sequences were compared, with the reduction in mitochondrial genomic information content observed to encourage spurious inferences in our samples. A probabilistic approach to variant prediction for the complete bobwhite mitogenomes revealed 344 segregating sites corresponding to 347 total mutations, including 49 putative nonsynonymous single nucleotide variants (SNVs) distributed across 12 protein coding genes. Evidence of gross heteroplasmy was observed for 13 bobwhites, with 10 of the 13 heteroplasmies involving one moderate to high frequency SNV. Haplotype network and phylogenetic analyses for the complete bobwhite mitogenome sequences revealed two divergent maternal lineages (d_XY = 0.00731; F_ST = 0.849; P < 0.05), thereby supporting the potential for two putative subspecies. However, the diverged lineage (n = 103 variants) almost exclusively involved bobwhites geographically classified as Colinus virginianus texanus, which is discordant with the expectations of previous geographic subspecies designations. Tests of adaptive evolution for functional divergence (MKT), frequency distribution tests (D, F_S) and phylogenetic analyses (RAxML) provide no evidence for positive selection or hybridization with the sympatric scaled quail (Callipepla squamata) as being explanatory factors for the two bobwhite maternal lineages observed. Instead, our analyses support the supposition that two diverged maternal lineages have survived from pre-expansion to post-expansion population(s), with the segregation of some slightly deleterious nonsynonymous mutations.

Extensive gene flow characterizes the phylogeography of a North American migrant bird: Black-headed Grosbeak (Pheucticus melanocephalus)

We describe range-wide phylogeographic variation in the Black-headed Grosbeak (Pheucticus melanocephalus), a songbird that is widely distributed across North American scrublands and forests. Phylogenetic analysis of mitochondrial DNA (mtDNA, n=424) revealed three geographically structured clades. One widespread clade occurs throughout the Rocky Mountains, Great Basin, and Mexican Plateau, a second clade is found on the Pacific coast and in coastal ranges; and, a third in the Sierra Madre del Sur of Oaxaca and Guerrero. Some geographical structuring occurs in Mexican Plateau and Sierra Madre Oriental mtDNA clade, presumably because these populations have been more stable over time than northern populations. Multiple mitochondrial groups are found sympatrically in the Okanogan River Valley in Washington, the eastern Sierra Nevada, and the Transvolcanic Belt across central Mexico, indicating that there is a potential for introgression. Analyses of 12 nuclear loci did not recover the same geographically structured clades. Population analyses show high levels of gene flow in nucDNA from the Interior into the Sierra Madre del Sur and Pacific population groups, possibly indicating expansion of the Interior population at the expense of peripheral populations.

Contrasting genetic structure among populations of two amphidromous fish species (Sicydiinae) in the Central West Pacific

Both present-day and past processes can shape connectivity of populations. Pleistocene vicariant events and dispersal have shaped the present distribution and connectivity patterns of aquatic species in the Indo-Pacific region. In particular, the processes that have shaped distribution of amphidromous goby species still remain unknown. Previous studies show that phylogeographic breaks are observed between populations in the Indian and Pacific Oceans where the shallow Sunda shelf constituted a geographical barrier to dispersal, or that the large spans of open ocean that isolate the Hawaiian or Polynesian Islands are also barriers for amphidromous species even though they have great dispersal capacity. Here we assess past and present genetic structure of populations of two amphidromous fish (gobies of the Sicydiinae) that are widely distributed in the Central West Pacific and which have similar pelagic larval durations. We analysed sections of mitochondrial COI, Cytb and nuclear Rhodospine genes in individuals sampled from different locations across their entire known range. Similar to other Sicydiinae fish, intraspecific mtDNA genetic diversity was high for all species (haplotype diversity between 0.9–0.96). Spatial analyses of genetic variation in Sicyopus zosterophorum demonstrated strong isolation across the Torres Strait, which was a geologically intermittent land barrier linking Australia to Papua New Guinea. There was a clear genetic break between the northwestern and the southwestern clusters in Si. zosterophorum (φ_ST = 0.67502 for COI) and coalescent analyses revealed that the two populations split at 306 Kyr BP (95% HPD 79–625 Kyr BP), which is consistent with a Pleistocene separation caused by the Torres Strait barrier. However, this geographical barrier did not seem to affect Sm. fehlmanni. Historical and demographic hypotheses are raised to explain the different patterns of population structure and distribution between these species. Strategies aiming to conserve amphidromous fish should consider the presence of cryptic evolutionary lineages to prevent stock depletion.

Vicariance across major marine biogeographic barriers: temporal concordance and the relative intensity of hard versus soft barriers

The marine tropics contain five major biogeographic regions (East Pacific, Atlantic, Indian Ocean, Indo-Australian Archipelago (IAA) and Central Pacific). These regions are separated by both hard and soft barriers. Reconstructing ancestral vicariance, we evaluate the extent of temporal concordance in vicariance events across three major barriers (Terminal Tethyan Event (TTE), Isthmus of Panama (IOP), East Pacific Barrier, EPB) and two incomplete barriers (either side of the IAA) for the Labridae, Pomacentridae and Chaetodontidae. We found a marked lack of temporal congruence within and among the three fish families in vicariance events associated with the EPB, TTE and IOP. Vicariance across hard barriers separating the Atlantic and Indo-Pacific (TTE, IOP) is temporally diffuse, with many vicariance events preceding barrier formation. In marked contrast, soft barriers either side of the IAA hotspot support tightly concordant vicariance events (2.5 Myr on Indian Ocean side; 6 Myr on Central Pacific side). Temporal concordance in vicariance points to large-scale temporally restricted gene flow during the Late Miocene and Pliocene. Despite different and often complex histories, both hard and soft barriers have comparably strong effects on the evolution of coral reef taxa.

High gene flow in reef fishes and its implications for ad-hoc no-take marine reserves

Integration of genetic connectivity information in effective marine reserve (MR) design is important in sustaining marine biodiversity. Here, genetic connectivity based on mitochondrial DNA (mtDNA) of three reef fish species, namely Epinephelus merra (n = 67; 32 from Bolinao, 14 from Alaminos, and 21 from Masinloc), Parupeneus multifasciatus (n = 23; 12 from Bolinao and 11 from Masinloc), and Odonus niger (n = 35; 21 from Mabini and 14 from Tingloy), sampled across western Luzon, Philippines, was inferred by assessing their genetic diversity, population genetic structure, and historical demography. The results show high haplotype and nucleotide diversity in the three species. Tests for population structure indicate high gene flow and no spatial genetic structuring for the three species. Mismatch analyses suggest unimodal distribution for E. merra and P. multifasciatus, but bimodal distribution for O. niger. Even with differences in mismatch distributions, all the three species exhibit low raggedness index indicating demographic population expansion. The bimodal distribution of O. niger could be attributed to the mixing of two isolated populations. High gene flow between sampling locations implies genetic exchanges and connectivity between many small MRs and fishing grounds in western Luzon, Philippines, at a scale similar to our study. This research is among the first few to elucidate the high genetic connectivity of reef fish communities across the Philippines (here western Luzon), but it also calls for more support (i.e. government and academia) for genetic research that aims to (1) understand the maintenance of megadiversity of the country and (2) search for effective biodiversity conservation options for the coral reefs.

Mitochondrial DNA diversity of mud crab Scylla olivacea (Portunidae) in Peninsular Malaysia: a preliminary assessment

A primary factor in population management and wildlife conservation is the delineation of population units derived from descriptions of population genetic structure. Yet, predicting factors that influence the patterns of gene flow in a population particularly at landscape scales remains a major challenge in evolutionary biology. Here we report a population genetic study of the mud crab Scylla olivacea examined based on a 542 bp segment of the mitochondrial DNA cytochrome c oxidase I gene among 91 individuals from six localities in the west and east coast of Peninsular Malaysia. In total 55 unique haplotypes were distinguished with 45 private haplotypes and a single common haplotype shared among all populations studied. The other ten haplotypes were shared among various populations. The sharing of this haplotype reflects the connection of the mangrove areas between east and west coast of Peninsular Malaysia. High haplotype diversity (h = 0.968 ± 0.021; mean ± SD) and low nucleotide diversity (π = 0.120 ± 0.015; mean ± SD) were displayed, which may be indicative of genetic bottleneck events. No significant phylogenetic lineages were recognized using neighbour-joining and maximum parsimony methods. Hierarchical AMOVA analysis indicated that 99.33 % of the genetic variation was contained within populations and 0.67 % occurred among populations, suggesting no geographical patterning among populations studied, supported by F st test. Mismatch distribution analysis showed that the observed distribution of the pairwise mutation differences among haplotypes was multimodal, which is not concordant with a sudden range expansion scenario. However, neutrality tests showed non-significant negative values suggesting that the populations studied may have experienced past population growth, but the expansion may have been restricted to separate local areas that resulted in the non-significant negative Fu's Fs and Tajima's D value. Overall, this present preliminary study was able to be a reference on the phylogenetic relationships and assessment of genetic structure of Scylla sp. in Malaysia.

Combined analyses of kinship and FST suggest potential drivers of chaotic genetic patchiness in high gene-flow populations

We combine kinship estimates with traditional F-statistics to explain contemporary drivers of population genetic differentiation despite high gene flow. We investigate range-wide population genetic structure of the California spiny (or red rock) lobster (Panulirus interruptus) and find slight, but significant global population differentiation in mtDNA (ΦST = 0.006, P = 0.001; D(est_Chao) = 0.025) and seven nuclear microsatellites (F(ST) = 0.004, P < 0.001; D(est_Chao) = 0.03), despite the species' 240- to 330-day pelagic larval duration. Significant population structure does not correlate with distance between sampling locations, and pairwise FST between adjacent sites often exceeds that among geographically distant locations. This result would typically be interpreted as unexplainable, chaotic genetic patchiness. However, kinship levels differ significantly among sites (pseudo-F(16,988) = 1.39, P = 0.001), and ten of 17 sample sites have significantly greater numbers of kin than expected by chance (P < 0.05). Moreover, a higher proportion of kin within sites strongly correlates with greater genetic differentiation among sites (D(est_Chao), R(2) = 0.66, P < 0.005). Sites with elevated mean kinship were geographically proximate to regions of high upwelling intensity (R(2) = 0.41, P = 0.0009). These results indicate that P. interruptus does not maintain a single homogenous population, despite extreme dispersal potential. Instead, these lobsters appear to either have substantial localized recruitment or maintain planktonic larval cohesiveness whereby siblings more likely settle together than disperse across sites. More broadly, our results contribute to a growing number of studies showing that low F(ST) and high family structure across populations can coexist, illuminating the foundations of cryptic genetic patterns and the nature of marine dispersal.