Regal phylogeography: Range-wide survey of the marine angelfish Pygoplites diacanthus reveals evolutionary partitions between the Red Sea, Indian Ocean, and Pacific Ocean

The regal angelfish (Pygoplites diacanthus; family Pomacanthidae) occurs on reefs from the Red Sea to the central Pacific, with an Indian Ocean/Rea Sea color morph distinct from a Pacific Ocean morph. To assess population differentiation and evaluate the possibility of cryptic evolutionary partitions in this monotypic genus, we surveyed mtDNA cytochrome b and two nuclear introns (S7 and RAG2) in 547 individuals from 15 locations. Phylogeographic analyses revealed four mtDNA lineages (d=0.006-0.015) corresponding to the Pacific Ocean, the Red Sea, and two admixed lineages in the Indian Ocean, a pattern consistent with known biogeographic barriers. Christmas Island in the eastern Indian Ocean had both Indian and Pacific lineages. Both S7 and RAG2 showed strong population-level differentiation between the Red Sea, Indian Ocean, and Pacific Ocean (ΦST=0.066-0.512). The only consistent population sub-structure within these three regions was at the Society Islands (French Polynesia), where surrounding oceanographic conditions may reinforce isolation. Coalescence analyses indicate the Pacific (1.7Ma) as the oldest extant lineage followed by the Red Sea lineage (1.4Ma). Results from a median-joining network suggest radiations of two lineages from the Red Sea that currently occupy the Indian Ocean (0.7-0.9Ma). Persistence of a Red Sea lineage through Pleistocene glacial cycles suggests a long-term refuge in this region. The affiliation of Pacific and Red Sea populations, apparent in cytochrome b and S7 (but equivocal in RAG2) raises the hypothesis that the Indian Ocean was recolonized from the Red Sea, possibly more than once. Assessing the genetic architecture of this widespread monotypic genus reveals cryptic evolutionary diversity that merits subspecific recognition. We recommend P.d. diacanthus and P.d. flavescens for the Pacific and Indian Ocean/Red Sea forms.

The Three Domains of Conservation Genetics: Case Histories from Hawaiian Waters

The scientific field of conservation biology is dominated by 3 specialties: phylogenetics, ecology, and evolution. Under this triad, phylogenetics is oriented towards the past history of biodiversity, conserving the divergent branches in the tree of life. The ecological component is rooted in the present, maintaining the contemporary life support systems for biodiversity. Evolutionary conservation (as defined here) is concerned with preserving the raw materials for generating future biodiversity. All 3 domains can be documented with genetic case histories in the waters of the Hawaiian Archipelago, an isolated chain of volcanic islands with 2 types of biodiversity: colonists, and new species that arose from colonists. This review demonstrates that 1) phylogenetic studies have identified previously unknown branches in the tree of life that are endemic to Hawaiian waters; 2) population genetic surveys define isolated marine ecosystems as management units, and 3) phylogeographic analyses illustrate the pathways of colonization that can enhance future biodiversity. Conventional molecular markers have advanced all 3 domains in conservation biology over the last 3 decades, and recent advances in genomics are especially valuable for understanding the foundations of future evolutionary diversity.

Genomic signatures of geographic isolation and natural selection in coral reef fishes

The drivers of speciation remain among the most controversial topics in evolutionary biology. Initially, Darwin emphasized natural selection as a primary mechanism of speciation, but the architects of the modern synthesis largely abandoned that view in favour of divergence by geographic isolation. The balance between selection and isolation is still at the forefront of the evolutionary debate, especially for the world's tropical oceans where biodiversity is high, but isolating barriers are few. Here, we identify the drivers of speciation in Pacific reef fishes of the genus Acanthurus by comparative genome scans of two peripheral populations that split from a large Central-West Pacific lineage at roughly the same time. Mitochondrial sequences indicate that populations in the Hawaiian Archipelago and the Marquesas Islands became isolated approximately 0.5 Ma. The Hawaiian lineage is morphologically indistinguishable from the widespread Pacific form, but the Marquesan form is recognized as a distinct species that occupies an unusual tropical ecosystem characterized by upwelling, turbidity, temperature fluctuations, algal blooms and little coral cover. An analysis of 3737 SNPs reveals a strong signal of selection at the Marquesas, with 59 loci under disruptive selection including an opsin Rh2 locus. While both the Hawaiian and Marquesan populations indicate signals of drift, the former shows a weak signal of selection that is comparable with populations in the Central-West Pacific. This contrast between closely related lineages reveals one population diverging due primarily to geographic isolation and genetic drift, and the other achieving taxonomic species status under the influence of selection.

Vertical and horizontal genetic connectivity in Chromis verater, an endemic damselfish found on shallow and mesophotic reefs in the Hawaiian Archipelago and adjacent Johnston Atoll

Understanding vertical and horizontal connectivity is a major priority in research on mesophotic coral ecosystems (30-150 m). However, horizontal connectivity has been the focus of few studies, and data on vertical connectivity are limited to sessile benthic mesophotic organisms. Here we present patterns of vertical and horizontal connectivity in the Hawaiian Islands-Johnston Atoll endemic threespot damselfish, Chromis verater, based on 319 shallow specimens and 153 deep specimens. The mtDNA markers cytochrome b and control region were sequenced to analyze genetic structure: 1) between shallow (< 30 m) and mesophotic (30-150 m) populations and 2) across the species' geographic range. Additionally, the nuclear markers rhodopsin and internal transcribed spacer 2 of ribosomal DNA were sequenced to assess connectivity between shallow and mesophotic populations. There was no significant genetic differentiation by depth, indicating high levels of vertical connectivity between shallow and deep aggregates of C. verater. Consequently, shallow and deep samples were combined by location for analyses of horizontal connectivity. We detected low but significant population structure across the Hawaiian Archipelago (overall cytochrome b: ΦST = 0.009, P = 0.020; control region: ΦST = 0.012, P = 0.009) and a larger break between the archipelago and Johnston Atoll (cytochrome b: ΦST = 0.068, P < 0.001; control region: ΦST = 0.116, P < 0.001). The population structure within the archipelago was driven by samples from the island of Hawaii at the southeast end of the chain and Lisianski in the middle of the archipelago. The lack of vertical genetic structure supports the refugia hypothesis that deep reefs may constitute a population reservoir for species depleted in shallow reef habitats. These findings represent the first connectivity study on a mobile organism that spans shallow and mesophotic depths and provide a reference point for future connectivity studies on mesophotic fishes.

Phylogeographic analyses of submesophotic snappers Etelis coruscans and Etelis "marshi" (family Lutjanidae) reveal concordant genetic structure across the Hawaiian Archipelago

The Hawaiian Archipelago has become a natural laboratory for understanding genetic connectivity in marine organisms as a result of the large number of population genetics studies that have been conducted across this island chain for a wide taxonomic range of organisms. However, population genetic studies have been conducted for only two species occurring in the mesophotic or submesophotic zones (30+m) in this archipelago. To gain a greater understanding of genetic connectivity in these deepwater habitats, we investigated the genetic structure of two submesophotic fish species (occurring ∼200-360 m) in this archipelago. We surveyed 16 locations across the archipelago for submesophotic snappers Etelis coruscans (N = 787) and E. "marshi" (formerly E. carbunculus; N = 770) with 436-490 bp of mtDNA cytochrome b and 10-11 microsatellite loci. Phylogeographic analyses reveal no geographic structuring of mtDNA lineages and recent coalescence times that are typical of shallow reef fauna. Population genetic analyses reveal no overall structure across most of the archipelago, a pattern also typical of dispersive shallow fishes. However some sites in the mid-archipelago (Raita Bank to French Frigate Shoals) had significant population differentiation. This pattern of no structure between ends of the Hawaiian range, and significant structure in the middle, was previously observed in a submesophotic snapper (Pristipomoides filamentosus) and a submesophotic grouper (Hyporthodus quernus). Three of these four species also have elevated genetic diversity in the mid-archipelago. Biophysical larval dispersal models from previous studies indicate that this elevated diversity may result from larval supplement from Johnston Atoll, ∼800 km southwest of Hawaii. In this case the boundaries of stocks for fishery management cannot be defined simply in terms of geography, and fishery management in Hawaii may need to incorporate external larval supply into management plans.

One size does not fit all: the emerging frontier in large-scale marine conservation

On the 20th anniversary of the Convention on Biological Diversity, a network of very large marine protected areas (the Big Ocean network) has emerged as a key strategy in the move to arrest marine decline and conserve some of the last remaining relatively undisturbed marine areas on the globe. Here we outline the ecological, economic and policy benefits of very large-scale MPAs and show their disproportionate value to global marine conservation targets. In particular we point out that very large-scale MPAs are a critical component of reaching the Aichi targets of protecting 10% of global marine habitats by 2020, because in addition to encompassing entire ecosystems, they will bring forward the expected date of achievement by nearly three decades (2025 as opposed to 2054). While the need for small MPAs remains critical, large MPAs will complement and enhance these conservation efforts. Big Ocean sites currently contain more than 80% of managed area in the sea, and provide our best hope for arresting the global decline in marine biodiversity.

Population structure in the native range predicts the spread of introduced marine species

Forecasting invasion success remains a fundamental challenge in invasion biology. The effort to identify universal characteristics that predict which species become invasive has faltered in part because of the diversity of taxa and systems considered. Here, we use an alternative approach focused on the spread stage of invasions. FST, a measure of alternative fixation of alleles, is a common proxy for realized dispersal among natural populations, summarizing the combined influences of life history, behaviour, habitat requirements, population size, history and ecology. We test the hypothesis that population structure in the native range (FST) is negatively correlated with the geographical extent of spread of marine species in an introduced range. An analysis of the available data (29 species, nine phyla) revealed a significant negative correlation (R(2) = 0.245-0.464) between FST and the extent of spread of non-native species. Mode FST among pairwise comparisons between populations in the native range demonstrated the highest predictive power (R(2) = 0.464, p < 0.001). There was significant improvement when marker type was considered, with mtDNA datasets providing the strongest relationship (n = 21, R(2) = 0.333-0.516). This study shows that FST can be used to make qualitative predictions concerning the geographical extent to which a non-native marine species will spread once established in a new area.

An invasive fish and the time-lagged spread of its parasite across the Hawaiian archipelago

Efforts to limit the impact of invasive species are frustrated by the cryptogenic status of a large proportion of those species. Half a century ago, the state of Hawai'i introduced the Bluestripe Snapper, Lutjanus kasmira, to O'ahu for fisheries enhancement. Today, this species shares an intestinal nematode parasite, Spirocamallanus istiblenni, with native Hawaiian fishes, raising the possibility that the introduced fish carried a parasite that has since spread to naïve local hosts. Here, we employ a multidisciplinary approach, combining molecular, historical, and ecological data to confirm the alien status of S. istiblenni in Hawai'i. Using molecular sequence data we show that S. istiblenni from Hawai'i are genetically affiliated with source populations in French Polynesia, and not parasites at a geographically intermediate location in the Line Islands. S. istiblenni from Hawai'i are a genetic subset of the more diverse source populations, indicating a bottleneck at introduction. Ecological surveys indicate that the parasite has found suitable intermediate hosts in Hawai'i, which are required for the completion of its life cycle, and that the parasite is twice as prevalent in Hawaiian Bluestripe Snappers as in source populations. While the introduced snapper has spread across the entire 2600 km archipelago to Kure Atoll, the introduced parasite has spread only half that distance. However, the parasite faces no apparent impediments to invading the entire archipelago, with unknown implications for naïve indigenous Hawaiian fishes and the protected Papahānaumokuākea Marine National Monument.

Microsatellites for next-generation ecologists: a post-sequencing bioinformatics pipeline

Microsatellites are the markers of choice for a variety of population genetic studies. The recent advent of next-generation pyrosequencing has drastically accelerated microsatellite locus discovery by providing a greater amount of DNA sequencing reads at lower costs compared to other techniques. However, laboratory testing of PCR primers targeting potential microsatellite markers remains time consuming and costly. Here we show how to reduce this workload by screening microsatellite loci via bioinformatic analyses prior to primer design. Our method emphasizes the importance of sequence quality, and we avoid loci associated with repetitive elements by screening with repetitive sequence databases available for a growing number of taxa. Testing with the Yellowstripe Goatfish Mulloidichthys flavolineatus and the marine planktonic copepod Pleuromamma xiphias we show higher success rate of primers selected by our pipeline in comparison to previous in silico microsatellite detection methodologies. Following the same pipeline, we discover and select microsatellite loci in nine additional species including fishes, sea stars, copepods and octopuses.

Phylogeography of two closely related Indo-Pacific butterflyfishes reveals divergent evolutionary histories and discordant results from mtDNA and microsatellites

Marine biogeographic barriers can have unpredictable consequences, even among closely related species. To resolve phylogeographic patterns for Indo-Pacific reef fauna, we conducted range-wide surveys of sister species, the scrawled butterflyfish (Chaetodon meyeri; N = 134) and the ornate butterflyfish (Chaetodon ornatissimus; N = 296), using mitochondrial DNA cytochrome b sequences and 10 microsatellite loci. The former is distributed primarily in the Indian Ocean but also extends to the Line Islands in the Central Pacific, whereas the latter is distributed primarily in the Central-West Pacific (including Hawaii and French Polynesia) but extends to the eastern margin of the Indian Ocean. Analyses of molecular variance and Bayesian STRUCTURE results revealed 1 range-wide group for C. meyeri and 3 groups for C. ornatissimus: 1) eastern Indian Ocean and western Pacific, 2) Central Pacific, and 3) Hawaii. Estimates of the last population expansion were much more recent for C. meyeri (61 500 to 95 000 years) versus C. ornatissimus (184 700 to 286 300 years). Despite similarities in ecology, morphology, life history, and a broadly overlapping distribution, these sister species have divergent patterns of dispersal and corresponding evolutionary history. The mtDNA and microsatellite markers did not provide concordant results within 1 of our study species (C. meyeri), or in 7 out of 12 other cases of marine fishes in the published literature. This discordance renews caution in relying on one or a few markers for reconstructing historical demography.

Coming out of the starting blocks: extended lag time rearranges genetic diversity in introduced marine fishes of Hawai'i

Biological invasions with known histories are rare, especially in the sea, and empirical studies of the genetic consequences are even rarer. Fifty-five years ago, the state of Hawai'i began a remarkable, if unintentional, 'experiment' with the introduction of three reef fishes, Lutjanus fulvus, Cephalopholis argus and Lutjanus kasmira. All have since expanded from the initial introduction of 2204 to 3163 individuals; however, historical records show that initially L. fulvus remained scarce, C. argus had modest population expansion and L. kasmira experienced rapid population growth. The consequences of differential population growth rates are apparent in F-statistics: Hawaiian L. fulvus demonstrate strong and significant haplotype frequency shifts from the founder location (F(ST) = 0.449), C. argus shows low but significant differentiation (F(ST) = 0.066) and L. kasmira is nearly identical to the founder location (F(ST) = 0.008). All three species had higher mtDNA diversity in the introduced range, which can be explained by multiple sources for L. fulvus and L. kasmira, but not for C. argus. We conclude that lag time before population expansion, in conjunction with genetic drift, has defined the genetic architecture of these three species in the introduced range.

Living in the past: phylogeography and population histories of Indo-Pacific wrasses (genus Halichoeres) in shallow lagoons versus outer reef slopes

Sea level fluctuations during glacial cycles affect the distribution of shallow marine biota, exposing the continental shelf on a global scale, and displacing coral reef habitat to steep slopes on oceanic islands. In these circumstances we expect that species inhabiting lagoons should show shallow genetic architecture relative to species inhabiting more stable outer reefs. Here we test this expectation on an ocean-basin scale with four wrasses (genus Halichoeres): H. claudia (N = 194, with ocean-wide distribution) and H. ornatissimus (N = 346, a Hawaiian endemic) inhabit seaward reef slopes, whereas H. trimaculatus (N = 239) and H. margaritaceus (N = 118) inhabit lagoons and shallow habitats throughout the Pacific. Two mitochondrial markers (cytochrome oxidase I and control region) were sequenced to resolve population structure and history of each species. Haplotype and nucleotide diversity were similar among all four species. The outer reef species showed significantly less population structure, consistent with longer pelagic larval durations. Mismatch distributions and significant negative Fu's F values indicate Pleistocene population expansion for all species, and (contrary to expectations) shallower histories in the outer slope species. We conclude that lagoonal wrasses may persist through glacial habitat disruptions, but are restricted to refugia during lower sea level stands. In contrast, outer reef slope species have homogeneous and well-connected populations through their entire ranges regardless of sea level fluctuations. These findings contradict the hypothesis that shallow species are less genetically diverse as a consequence of glacial cycles.

Common misconceptions in molecular ecology: echoes of the modern synthesis

The field of molecular ecology has burgeoned into a large discipline spurred on by technical innovations that facilitate the rapid acquisition of large amounts of genotypic data, by the continuing development of theory to interpret results, and by the availability of computer programs to analyse data sets. As the discipline grows, however, misconceptions have become enshrined in the literature and are perpetuated by routine citations to other articles in molecular ecology. These misconceptions hamper a better understanding of the processes that influence genetic variation in natural populations and sometimes lead to erroneous conclusions. Here, we consider eight misconceptions commonly appearing in the literature: (i) some molecular markers are inherently better than other markers; (ii) mtDNA produces higher F(ST) values than nDNA; (iii) estimated population coalescences are real; (iv) more data are always better; (v) one needs to do a Bayesian analysis; (vi) selective sweeps influence mtDNA data; (vii) equilibrium conditions are critical for estimating population parameters; and (viii) having better technology makes us smarter than our predecessors. This is clearly not an exhaustive list and many others can be added. It is, however, sufficient to illustrate why we all need to be more critical of our own understanding of molecular ecology and to be suspicious of self-evident truths.

Is the Great Barracuda (Sphyraena barracuda) a reef fish or a pelagic fish? The phylogeographic perspective

Current taxonomy indicates a single global species of the Great Barracuda (Sphyraena barracuda) despite differences in color and behavior between Atlantic and Pacific forms. To investigate these differences and qualify the dispersal characteristics of this unique coastal- pelagic teleost (bony fish), we conducted a global phylogeographic survey of 246 specimens from thirteen sampling locations using a 629-base pair fragment of mtDNA cytochrome b. Data indicate high overall gene flow in the Indo-Pacific over large distances (>16,500 km) bridging several biogeographic barriers. The West Atlantic population contains an mtDNA lineage that is divergent from the Indo-Pacific (d = 1.9%), while the East Atlantic (N = 23) has two mutations (d = 0.6%) apart from the Indo-Pacific. While we cannot rule out distinct evolutionary partitions among ocean basins based on behavior, coloration, and near-monophyly between Atlantic and Indo-Pacific subpopulations, more investigation is required before taxonomic status is revised. Overall, the pattern of high global dispersal and connectivity in S. barracuda more closely resembles those reported for large oceanic predators than reef-associated teleosts.

Global phylogeography with mixed-marker analysis reveals male-mediated dispersal in the endangered scalloped hammerhead shark (Sphyrna lewini)

Background

The scalloped hammerhead shark, Sphyrna lewini, is a large endangered predator with a circumglobal distribution, observed in the open ocean but linked ontogenetically to coastal embayments for parturition and juvenile development. A previous survey of maternal (mtDNA) markers demonstrated strong genetic partitioning overall (global Φ_ST = 0.749) and significant population separations across oceans and between discontinuous continental coastlines.

Methodology/Principal Findings

We surveyed the same global range with increased sample coverage (N = 403) and 13 microsatellite loci to assess the male contribution to dispersal and population structure. Biparentally inherited microsatellites reveal low or absent genetic structure across ocean basins and global genetic differentiation (F_ST = 0.035) over an order of magnitude lower than the corresponding measures for maternal mtDNA lineages (Φ_ST = 0.749). Nuclear allelic richness and heterozygosity are high throughout the Indo-Pacific, while genetic structure is low. In contrast, allelic diversity is low while population structure is higher for populations at the ends of the range in the West Atlantic and East Pacific.

Conclusions/Significance

These data are consistent with the proposed Indo-Pacific center of origin for S. lewini, and indicate that females are philopatric or adhere to coastal habitats while males facilitate gene flow across oceanic expanses. This study includes the largest sampling effort and the most molecular loci ever used to survey the complete range of a large oceanic predator, and findings emphasize the importance of incorporating mixed-marker analysis into stock assessments of threatened and endangered shark species.

High connectivity in the deepwater snapper Pristipomoides filamentosus (Lutjanidae) across the Indo-Pacific with isolation of the Hawaiian archipelago

In the tropical Indo-Pacific, most phylogeographic studies have focused on the shallow-water taxa that inhabit reefs to approximately 30 m depth. Little is known about the large predatory fishes, primarily snappers (subfamily Etelinae) and groupers (subfamily Epinephelinae) that occur at 100-400 m. These long-lived, slow-growing species support fisheries across the Indo-Pacific, yet no comprehensive genetic surveys within this group have been conducted. Here we contribute the first range-wide survey of a deepwater Indo-Pacific snapper, Pristipomoides filamentosus, with special focus on Hawai'i. We applied mtDNA cytochrome b and 11 microsatellite loci to 26 samples (N=1,222) collected across 17,000 km from Hawai'i to the western Indian Ocean. Results indicate that P. filamentosus is a highly dispersive species with low but significant population structure (mtDNA Φ(ST)=0.029, microsatellite F(ST)=0.029) due entirely to the isolation of Hawai'i. No population structure was detected across 14,000 km of the Indo-Pacific from Tonga in the Central Pacific to the Seychelles in the western Indian Ocean, a pattern rarely observed in reef species. Despite a long pelagic phase (60-180 days), interisland dispersal as adults, and extensive gene flow across the Indo-Pacific, P. filamentosus is unable to maintain population connectivity with Hawai'i. Coalescent analyses indicate that P. filamentosus may have colonized Hawai'i 26 K-52 K y ago against prevailing currents, with dispersal away from Hawai'i dominating migration estimates. P. filamentosus harbors low genetic diversity in Hawai'i, a common pattern in marine fishes, and our data indicate a single archipelago-wide stock. However, like the Hawaiian Grouper, Hyporthodus quernus, this snapper had several significant pairwise comparisons (F(ST)) clustered around the middle of the archipelago (St. Rogatien, Brooks Banks, Gardner) indicating that this region may be isolated or (more likely) receives input from Johnston Atoll to the south.

Phylogeography of the reef fish Cephalopholis argus (Epinephelidae) indicates Pleistocene isolation across the Indo-Pacific Barrier with contemporary overlap in The Coral Triangle

BACKGROUND

The Coral Triangle (CT), bounded by the Philippines, the Malay Peninsula, and New Guinea, is the epicenter of marine biodiversity. Hypotheses that explain the source of this rich biodiversity include 1) the center of origin, 2) the center of accumulation, and 3) the region of overlap. Here we contribute to the debate with a phylogeographic survey of a widely distributed reef fish, the Peacock Grouper (Cephalopholis argus; Epinephelidae) at 21 locations (N = 550) using DNA sequence data from mtDNA cytochrome b and two nuclear introns (gonadotropin-releasing hormone and S7 ribosomal protein).

RESULTS

Population structure was significant (ΦST = 0.297, P < 0.001; FST = 0.078, P < 0.001; FST = 0.099, P < 0.001 for the three loci, respectively) among five regions: French Polynesia, the central-west Pacific (Line Islands to northeastern Australia), Indo-Pacific boundary (Bali and Rowley Shoals), eastern Indian Ocean (Cocos/Keeling and Christmas Island), and western Indian Ocean (Diego Garcia, Oman, and Seychelles). A strong signal of isolation by distance was detected in both mtDNA (r = 0.749, P = 0.001) and the combined nuclear loci (r = 0.715, P < 0.001). We detected evidence of population expansion with migration toward the CT. Two clusters of haplotypes were detected in the mtDNA data (d = 0.008), corresponding to the Pacific and Indian Oceans, with a low level of introgression observed outside a mixing zone at the Pacific-Indian boundary.

CONCLUSIONS

We conclude that the Indo-Pacific Barrier, operating during low sea level associated with glaciation, defines the primary phylogeographic pattern in this species. These data support a scenario of isolation on the scale of 105 year glacial cycles, followed by population expansion toward the CT, and overlap of divergent lineages at the Pacific-Indian boundary. This pattern of isolation, divergence, and subsequent overlap likely contributes to species richness at the adjacent CT and is consistent with the region of overlap hypothesis.