Search

Studying range expansions is central for understanding genetic variation through space and time as well as for identifying refugia and biological invasions. Range expansions are characterized by serial founder events causing clines of decreasing genetic diversity away from the center of origin and asymmetries in the two-dimensional allele frequency spectra. These asymmetries, summarized by the directionality index (ψ), are sensitive to range expansions and persist for longer than clines in genetic diversity. In continuous and finite meta-populations, genetic drift tends to be stronger at the edges of the species distribution in equilibrium populations and populations undergoing range expansions alike. Such boundary effects are expected to affect geographic patterns in genetic diversity and ψ. Here we demonstrate that boundary effects cause high false positive rates in equilibrium meta-populations when testing for range expansions. In the simulations, the absolute value of ψ (| ψ |) in equilibrium data sets was proportional to the fixation index (F ST). By fitting signatures of range expansions as a function of ɛ | ψ |/ F ST and geographic clines in ψ , strong evidence for range expansions could be detected in data from a recent rapid invasion of the cane toad, Rhinella marina , in Australia, but not in 28 previously published empirical data sets from Australian scincid lizards that were significant for the standard range expansion tests. Thus, while clinal variation in ψ is still the most sensitive statistic to range expansions, to detect true signatures of range expansions in natural populations, its magnitude needs to be considered in relation to the overall levels of genetic structuring in the data. [ABSTRACT FROM AUTHOR]

When populations repeatedly adapt to similar environments they can evolve similar phenotypes based on shared genetic mechanisms (parallel evolution). The likelihood of parallel evolution is affected by demographic history, as it depends on the standing genetic variation of the source population. The three-spined stickleback (Gasterosteus aculeatus) repeatedly colonized and adapted to brackish and freshwater. Most parallel evolution studies in G. aculeatus were conducted at high latitudes, where freshwater populations maintain connectivity to the source marine populations. Here, we analysed southern and northern European marine and freshwater populations to test two hypotheses. First, that southern European freshwater populations (which currently lack connection to marine populations) lost genetic diversity due to bottlenecks and inbreeding compared to their northern counterparts. Second, that the degree of genetic parallelism is higher among northern than southern European freshwater populations, as the latter have been subjected to strong drift due to isolation. The results show that southern populations exhibit lower genetic diversity but a higher degree of genetic parallelism than northern populations. Hence, they confirm the hypothesis that southern populations have lost genetic diversity, but this loss probably happened after they had already adapted to freshwater conditions, explaining the high degree of genetic parallelism in the south. [ABSTRACT FROM AUTHOR]

The Earth's polar regions are low rates of inter‐ and intraspecific diversification. An extreme mammalian example is the Arctic ringed seal (Pusa hispida hispida), which is assumed to be panmictic across its circumpolar Arctic range. Yet, local Inuit communities in Greenland and Canada recognize several regional variants; a finding supported by scientific studies of body size variation. It is however unclear whether this phenotypic variation reflects plasticity, morphs or distinct ecotypes. Here, we combine genomic, biologging and survey data, to document the existence of a unique ringed seal ecotype in the Ilulissat Icefjord (locally 'Kangia'), Greenland; a UNESCO World Heritage site, which is home to the most productive marine‐terminating glacier in the Arctic. Genomic analyses reveal a divergence of Kangia ringed seals from other Arctic ringed seals about 240 kya, followed by secondary contact since the Last Glacial Maximum. Despite ongoing gene flow, multiple genomic regions appear under strong selection in Kangia ringed seals, including candidate genes associated with pelage coloration, growth and osmoregulation, potentially explaining the Kangia seal's phenotypic and behavioural uniqueness. The description of 'hidden' diversity and adaptations in yet another Arctic species merits a reassessment of the evolutionary processes that have shaped Arctic diversity and the traditional view of this region as an evolutionary freezer. Our study highlights the value of indigenous knowledge in guiding science and calls for efforts to identify distinct populations or ecotypes to understand how these might respond differently to environmental change. [ABSTRACT FROM AUTHOR]

We used mitochondrial DNA (mtDNA) gene sequences and nuclear microsatellite loci to investigate the extent and outcome of hybridization between the Black‐billed Gull Chroicocephalus bulleri and the Red‐billed Gull Chroicocephalus novaehollandiae scopulinus in New Zealand. Six of 26 sampled Black‐billed Gulls possessed mtDNA typical of Red‐billed Gulls, but allele frequencies at six polymorphic microsatellites provided little evidence of mixed ancestry expected in very recent hybrids. None of the Red‐billed Gulls sampled from different colonies possessed Black‐billed Gull mtDNA expected in the reciprocal cross, suggesting that hybridization in the two species typically occurs between female Red‐billed Gulls and Black‐billed Gull males. The lack of any hybrid signal in the nuclear loci indicates that there has been extensive backcrossing with Black‐billed Gulls, effectively diluting the Red‐billed Gull nuclear DNA contribution. Divergence of Red‐billed Gulls and Black‐billed Gulls occurred approximately 250 000 years ago, indicating that unsorted ancestral polymorphism is an unlikely alternative to hybridization. Comparing demographic models within an approximate Bayesian computation (ABC) framework, we confirm that the observed patterns cannot result from incomplete lineage sorting. Using an ABC random forest approach, we determined that the most likely model explaining the data is a recent introgression scenario, whereby unidirectional gene flow is re‐established following a period of strict isolation. The ability of Black‐billed and Red‐billed Gulls to successfully interbreed shows that despite significant differentiation (FST > 0.3), there has been insufficient time for the two species to develop complete reproductive isolation. The apparent one‐way transfer of Red‐billed Gull mtDNA into Black‐billed Gulls and extensive backcrossing argues against cytoplasmic–nuclear genome incompatibilities between the two species. We hypothesize that the specific mate recognition system cued on colours of soft parts normally functions to prevent hybridization, but that it can break down under demographic conditions where there is a shortage of available mates and a surplus of females in the Red‐billed Gull population. The high incidence of introgression in Black‐billed Gulls conflicts with field observations that interbreeding is extremely rare. [ABSTRACT FROM AUTHOR]

The evolution of sexual dimorphism (SD) is driven by intricate interplays between sexual and natural selection. When it comes to SD variation within populations, however, environmental factors play a major role. Sexually selected traits are expected to be strongly dependent on individual body condition, which is influenced by the local environment that individuals experience. As a consequence, the degree of SD may also depend on resource availability. Here, we investigated the potential drivers of SD expression at two sexually dimorphic morphometric traits, body size (snout vent length) and head shape (head geometric morphometrics), in the Italian wall lizard (Podarcis siculus). We assessed the existence of condition- and context-dependent SD across ten islands of the Aeolian archipelago (southern Italy), at within- and among-population scales. We observed strong geographical variation of SD among islands, and tested three potential SD predictors related to resource availability (individual body condition, ecosystem productivity, temperature). Body condition and ecosystem productivity were the main drivers of body size SD variation, and body condition was also the main driver for head shape SD. Our results highlight that the expression of SD in the Italian wall lizard is both condition- and context-dependent. These results are congruent at within- and among-populations scales highlighting that spatial multi-scale analysis represents a useful approach to understand patterns of SD expression. [ABSTRACT FROM AUTHOR]

The European whitefish (Coregonus lavaretus) species complex is a classic example of recent adaptive radiation. Here, we examine a whitefish population introduced to northern Finnish Lake Tsahkal in the late 1960s, where three divergent morphs (viz. littoral, pelagic, and profundal feeders) were found 10 generations after. Using demographic modeling based on genomic data, we show that whitefish morphs evolved during a phase of strict isolation, refuting a rapid sympatric divergence scenario. The lake is now an artificial hybrid zone between morphs originated in allopatry. Despite their current syntopy, clear genetic differentiation remains between two of the three morphs. Using admixture mapping, we identify five SNPs associated with gonad weight variation, a proxy for sexual maturity and spawning time. We suggest that ecological adaptations in spawning time evolved in allopatry are currently maintaining partial reproductive isolation in the absence of other barriers to gene flow. [ABSTRACT FROM AUTHOR]

The ongoing decline of large marine vertebrates must be urgently mitigated, particularly under increasing levels of climate change and other anthropogenic pressures. However, characterizing the connectivity among populations remains one of the greatest challenges for the effective conservation of an increasing number of endangered species. Achieving conservation targets requires an understanding of which seascape features influence dispersal and subsequent genetic structure. This is particularly challenging for adult‐disperser species, and when distribution‐wide sampling is difficult. Here, we developed a two‐step modelling framework to investigate how seascape features drive the genetic connectivity of marine species without larval dispersal, to better guide the design of marine protected area networks and corridors. We applied this framework to the endangered grey reef shark, Carcharhinus amblyrhynchos, a reef‐associated shark distributed across the tropical Indo‐Pacific. In the first step, we developed a seascape genomic approach based on isolation‐by‐resistance models involving circuit theory applied to 515 shark samples, genotyped for 4991 nuclear single‐nucleotide polymorphisms. We show that deep oceanic areas act as strong barriers to dispersal, while proximity to habitat facilitates dispersal. In the second step, we predicted the resulting genetic differentiation across the entire distribution range of the species, providing both local and global‐scale conservation units for future management guidance. We found that grey reef shark populations are more fragmented than expected for such a mobile species, raising concerns about the resilience of isolated populations under high anthropogenic pressures. We recommend the use of this framework to identify barriers to gene flow and to help in the delineation of conservation units at different scales, together with its integration across multiple species when considering marine spatial planning. [ABSTRACT FROM AUTHOR]

The three‐spined stickleback (Gasterosteus aculeatus) has repeatedly and independently adapted to freshwater habitats from standing genetic variation (SGV) following colonization from the sea. However, in the Mediterranean Sea G. aculeatus is believed to have gone extinct, and thus the spread of locally adapted alleles between different freshwater populations via the sea since then has been highly unlikely. This is expected to limit parallel evolution, that is the extent to which phylogenetically related alleles can be shared among independently colonized freshwater populations. Using whole genome and 2b‐RAD sequencing data, we compared levels of genetic differentiation and genetic parallelism of 15 Adriatic stickleback populations to 19 Pacific, Atlantic and Caspian populations, where gene flow between freshwater populations across extant marine populations is still possible. Our findings support previous studies suggesting that Adriatic populations are highly differentiated (average FST ≈ 0.45), of low genetic diversity and connectivity, and likely to stem from multiple independent colonizations during the Pleistocene. Linkage disequilibrium network analyses in combination with linear mixed models nevertheless revealed several parallel marine–freshwater differentiated genomic regions, although still not to the extent observed elsewhere in the world. We hypothesize that current levels of genetic parallelism in the Adriatic lineages are a relic of freshwater adaptation from SGV prior to the extinction of marine sticklebacks in the Mediterranean that has persisted despite substantial genetic drift experienced by the Adriatic stickleback isolates. [ABSTRACT FROM AUTHOR]

Population genetic theory predicts that small effective population sizes (N e) and restricted gene flow limit the potential for local adaptation. In particular, the probability of evolving similar phenotypes based on shared genetic mechanisms (i.e. parallel evolution), is expected to be reduced. We tested these predictions in a comparative genomic study of two ecologically similar and geographically codistributed stickleback species (viz. Gasterosteus aculeatus and Pungitius pungitius). We found that P. pungitius harbors less genetic diversity and exhibits higher levels of genetic differentiation and isolation-by-distance than G. aculeatus. Conversely, G. aculeatus exhibits a stronger degree of genetic parallelism across freshwater populations than P. pungitius : 2,996 versus 379 single nucleotide polymorphisms located within 26 versus 9 genomic regions show evidence of selection in multiple freshwater populations of G. aculeatus and P. pungitius , respectively. Most regions involved in parallel evolution in G. aculeatus showed increased levels of divergence, suggestive of selection on ancient haplotypes. In contrast, haplotypes involved in freshwater adaptation in P. pungitius were younger. In accordance with theory, the results suggest that connectivity and genetic drift play crucial roles in determining the levels and geographic distribution of standing genetic variation, providing evidence that population subdivision limits local adaptation and therefore also the likelihood of parallel evolution. [ABSTRACT FROM AUTHOR]

Testing among competing demographic models of divergence has become an important component of evolutionary research in model and non-model organisms. However, the effect of unaccounted demographic events on model choice and parameter estimation remains largely unexplored. Using extensive simulations, we demonstrate that under realistic divergence scenarios, failure to account for population size (N e) changes in daughter and ancestral populations leads to strong biases in divergence time estimates as well as model choice. We illustrate these issues reconstructing the recent demographic history of North Sea and Baltic Sea turbots (Scophthalmus maximus) by testing 16 isolation with migration (IM) and 16 secondary contact (SC) scenarios, modeling changes in N e as well as the effects of linked selection and barrier loci. Failure to account for changes in N e resulted in selecting SC models with long periods of strict isolation and divergence times preceding the formation of the Baltic Sea. In contrast, models accounting for N e changes suggest recent (<6 kya) divergence with constant gene flow. We further show how interpreting genomic landscapes of differentiation can help discerning among competing models. For example, in the turbot data, islands of differentiation show signatures of recent selective sweeps, rather than old divergence resisting secondary introgression. The results have broad implications for the study of population divergence by highlighting the potential effects of unmodeled changes in N e on demographic inference. Tested models should aim at representing realistic divergence scenarios for the target taxa, and extreme caution should always be exercised when interpreting results of demographic modeling. [ABSTRACT FROM AUTHOR]

Oyster reef habitats are critical to coastal biodiversity and their decline has prompted restoration efforts in Australia. Knowledge gaps exist regarding the population structure and diversity of key species in these habitats. This may be critical information for the design of effective restoration programs. Sydney rock oysters (Saccostrea glomerata) are the dominant reef-forming bivalve in eastern Australia. Wild populations of S. glomerata have declined due to overharvesting, disease outbreaks, coastal development and reduced water quality. Here, we use genetic markers identified by genome-wide sequencing to investigate the genetic structure and diversity of wild Sydney rock oysters throughout their distribution in eastern Australia. We examine evidence for past population bottlenecks and spatial genetic structure associated with the East Australian Current. Analysis of 3, 400 neutral single-nucleotide polymorphisms (SNPs) revealed a single population, and an overlap with two other Saccostrea sp. at the northernmost boundary of the distribution. We detected signals of asymmetric gene flow consistent with the direction of the East Australian Current, and spatial structure patterns of limited genetic isolation by distance and spatial autocorrelation in the northern region (which experiences stronger effects of the East Australian Current) but not in the southern region of the distribution. We found no evidence of significant recent bottlenecks, with high effective population size throughout the species' range. This information will provide a baseline against which to assess the impact of restoration projects, and guide strategies for sourcing stock for the enhancement of wild oyster populations. Our results provide a positive outlook for the resilience and adaptive capacity of Sydney rock oysters, and highlight wild populations as valuable resources for aquaculture and restoration initiatives. [ABSTRACT FROM AUTHOR]

The role of geological events and Pleistocene climatic fluctuations as drivers of current patterns of genetic variation in extant species has been a topic of continued interest among evolutionary biologists. Nevertheless, comprehensive studies of widely distributed species are still rare, especially from Asia. Using geographically extensive sampling of many individuals and a large number of nuclear single nucleotide polymorphisms (SNPs), we studied the phylogeography and historical demography of Hyla annectans populations in southern China. Thirty‐five sampled populations were grouped into seven clearly defined genetic clusters that closely match phenotype‐based subspecies classification. These lineages diverged 2.32–5.23 million years ago (Ma), a timing that closely aligns with the rapid and drastic uplifting of the Qinghai‐Tibet Plateau and adjacent southwest China. Demographic analyses and species distribution models indicate that different populations of this species have responded differently to past climatic changes. In the Hengduan Mountains, most populations experienced a bottleneck, whereas the populations located outside of the Hengduan Mountains have gradually declined in size since the end of the last glaciation. In addition, the levels of phenotypic and genetic divergence were strongly correlated across major clades. These results highlight the combined effects of geological events and past climatic fluctuations, as well as natural selection, as drivers of contemporary patterns of genetic and phenotypic variation in a widely distributed anuran in Asia. [ABSTRACT FROM AUTHOR]

Recent years have brought the realization that evolutionary changes driven by selection can occur in ecological time scales. However, recent evolutionary events can be hard to detect and may easily go unnoticed. For harvested species, such cryptic diversity may lead to suboptimal management. These points are illustrated by the two flounder species in the Baltic Sea. Although early ecological studies identified two ecotypes of European flounder (Platichthys flesus) based on spawning differences, genomic studies only recently demonstrated that they were reproductively isolated species, P. flesus and P. solemdali , separated through rapid ecologically driven speciation. These morphologically indistinguishable species are harvested within a mixed-stock fishery. In the northern Baltic Sea flounder landings have declined since the mid-1980s, with a drop in the Gulf of Finland (GoF) being particularly dramatic (∼90%). Genetic analyses of historical otolith samples from GoF catches have revealed that back in 1983 the fishery unknowingly targeted primarily P. flesus , whereas thereafter almost exclusively P. solemdali. Hence, the case of two flounder species illustrates (i) how ecological studies stimulated genetic investigations leading to discovery of ecological speciation and (ii) how cryptic species turnover discovered with genetic tools in turn improved ecological understanding with benefits to management and conservation. [ABSTRACT FROM AUTHOR]

Hybridization and convergent evolution are phenomena of broad interest in evolutionary biology, but their occurrence poses challenges for reconstructing evolutionary affinities among affected taxa. Sticklebacks in the genus Pungitius are a case in point: evolutionary relationships and taxonomic validity of different species and populations in this circumpolarly distributed species complex remain contentious due to convergent evolution of traits regarded as diagnostic in their taxonomy, and possibly also due to frequent hybridization among taxa. To clarify the evolutionary relationships among different Pungitius species and populations globally, as well as to study the prevalence and extent of introgression among recognized species, genomic data sets of both reference genome‐anchored single nucleotide polymorphisms and de novo assembled RAD‐tag loci were constructed with RAD‐seq data. Both data sets yielded topologically identical and well‐supported species trees. Incongruence between nuclear and mitochondrial DNA‐based trees was found and suggested possibly frequent hybridization and mitogenome capture during the evolution of Pungitius sticklebacks. Further analyses revealed evidence for frequent nuclear genetic introgression among Pungitius species, although the estimated proportions of autosomal introgression were low. Apart from providing evidence for frequent hybridization, the results challenge earlier mitochondrial and morphology‐based hypotheses regarding the number of species and their affinities in this genus: at least seven extant species can be recognized on the basis of genetic data. The results also shed new light on the biogeographical history of the Pungitius‐complex, including suggestion of several trans‐Arctic invasions of Europe from the Northern Pacific. The well‐resolved phylogeny should facilitate the utility of this genus as a model system for future comparative evolutionary studies. [ABSTRACT FROM AUTHOR]

Unobserved diversity, such as undetected genetic structure or the presence of cryptic species, is of concern for the conservation and management of global biodiversity in the face of threatening anthropogenic processes. For instance, unobserved diversity can lead to overestimation of maximum sustainable yields and therefore to overharvesting of the more vulnerable stock components within unrecognized mixed‐stock fisheries. We used DNA from archival (otolith) samples to reconstruct the temporal (1976–2011) genetic makeup of two mixed‐stock flounder fisheries in the Åland Sea (AS) and the Gulf of Finland (GoF). Both fisheries have hitherto been managed as a single stock of European flounders (Platichthys flesus), but were recently revealed to target two closely related species: the pelagic‐spawning P. flesus and the newly described, demersal‐spawning P. solemdali. While the AS and GoF fisheries were assumed to consist exclusively of P. solemdali, P. flesus dominated the GoF flounder assemblage (87% of total) in 1983, had disappeared (0%) by 1993, and remained in low proportions (10%–11%) thereafter. In the AS, P. solemdali dominated throughout the sampling period (>70%), and P. flesus remained in very low proportions after 1983. The disappearance of P. flesus from the GoF coincides in time with a dramatic (~60%) decline in commercial landings and worsening environmental conditions in P. flesus' northernmost spawning ground, the Eastern Gotland Basin, in the preceding 4–6 years. These results are compatible with the hypothesis that P. flesus in the GoF is a sink population relying on larval subsidies from southern spawning grounds and the cause of their disappearance is a cessation of larval supply. Our results highlight the importance of uncovering unobserved genetic diversity and studying spatiotemporal changes in the relative contribution of different stock components, as well as the underlying environmental causes, to manage marine resources in the age of rapid anthropogenic change. [ABSTRACT FROM AUTHOR]

Divergent selection may initiate ecological speciation extremely rapidly. How often and at what pace ecological speciation proceeds to yield strong reproductive isolation is more uncertain. Here, we document a case of extraordinarily rapid speciation associated with ecological selection in the postglacial Baltic Sea. European flounders (Platichthys flesus) in the Baltic exhibit two contrasting reproductive behaviors: pelagic and demersal spawning. Demersal spawning enables flounders to thrive in the low salinity of the Northern Baltic, where eggs cannot achieve neutral buoyancy. We show that demersal and pelagic flounders are a species pair arising from a recent event of speciation. Despite having a parapatric distribution with extensive overlap, the two species are reciprocally monophyletic and show strongly bimodal genotypic clustering and no evidence of contemporary migration, suggesting strong reproductive isolation. Divergence across the genome is weak but shows strong signatures of selection, a pattern suggestive of a recent ecological speciation event. We propose that spawning behavior in Baltic flounders is the trait under ecologically based selection causing reproductive isolation, directly implicating a process of ecological speciation. We evaluated different possible evolutionary scenarios under the approximate Bayesian computation framework and estimate that the speciation process started in allopatry ~2,400 generations ago, following the colonization of the Baltic by the demersal lineage. This is faster than most known cases of ecological speciation and represents the most rapid event of speciation ever reported for any marine vertebrate. [ABSTRACT FROM AUTHOR]

The grey nurse shark (Carcharius taurus) is listed as threatened throughout much of its global distribution, and as critically endangered in eastern Australia. Captive breeding programs have thus far been largely unsuccessful and little is known of its mating system in this context. Here we carry out a paternity analysis to determine if the mating system in captivity is characterised by multiple mating, and whether poor offspring survival is associated with a particular male. Tissue samples from grey nurse sharks were collected from three potential sires, the two dams and nine pups housed at Manly SEA LIFE Sanctuary in eastern Australia. Each individual was genotyped at seven microsatellite markers and three cases of multiple paternity were inferred. No paternal link to stillborn (5), or scoliotic (2) pups was indicated. For the first time, we show the natural wild phenomenon of multiple paternity occurring in a captive environment. [ABSTRACT FROM AUTHOR]

The diversity of geographic scales at which marine organisms display genetic variation mirrors the biophysical and ecological complexity of dispersal by pelagic larvae. Yet little is known about the effect of larval ecology on genetic population patterns, partly because detailed data of larval ecology do not yet exist for most taxa. One species for which this data is available is Eleutheronema tetradactylum, a tropical Indo-West Pacific shorefish. Here, we use a partial sequence mitochondrial cytochrome oxidase subunit 1 (COI) marker and five microsatellite loci to survey the genetic structure of E. tetradactylum across northern Australia. Structure was found throughout the range and isolation by distance was strong, explaining approximately 87 and 64% of the genetic variation in microsatellites and mtDNA, respectively. Populations separated by as little as 15 km also showed significant genetic structure, implying that local populations are mainly insular and self-seeding on an ecological time frame. Because the larvae of E. tetradactylum have lower swimming performance and poor orientation compared with other tropical fishes, even modest larval abilities may permit self-recruitment rather than passive dispersal. [ABSTRACT FROM AUTHOR]

Grey reef sharks ( Carcharhinus amblyrhynchos) are important apex predators on coral reefs, and their numbers have declined dramatically as a result of overfishing. Knowledge of environmental factors that shape gene flow is essential for developing appropriate management strategies, but the lack of suitable genetic markers has hindered research on this species. Here, we characterised 15 polymorphic microsatellite loci for grey reef sharks. None of the loci deviated significantly from Hardy-Weinberg equilibrium and there was no evidence of Linkage Disequilibrium. Several loci cross-amplified in other carcharhinid species, and will be useful in future studies of this family. [ABSTRACT FROM AUTHOR]

Grey reef sharks (Carcharhinus amblyrhynchos) can be one of the numerically dominant high order predators on pristine coral reefs, yet their numbers have declined even in the highly regulated Australian Great Barrier Reef (GBR) Marine Park. Knowledge of both large scale and fine scale genetic connectivity of grey reef sharks is essential for their effective management, but no genetic data are yet available. We investigated grey reef shark genetic structure in the GBR across a 1200 km latitudinal gradient, comparing empirical data with models simulating different levels of migration. The empirical data did not reveal any genetic structuring along the entire latitudinal gradient sampled, suggesting regular widespread dispersal and gene flow of the species throughout most of the GBR. Our simulated datasets indicate that even with substantial migrations (up to 25% of individuals migrating between neighboring reefs) both large scale genetic structure and genotypic spatial autocorrelation at the reef scale were maintained. We suggest that present migration rates therefore exceed this level. These findings have important implications regarding the effectiveness of networks of spatially discontinuous Marine Protected Areas to protect reef sharks. [ABSTRACT FROM AUTHOR]

Background:Dinoflagellates are important primary producers, crucial in marine food webs. Toxic strains, however, are the main causative agents of non-bacterial seafood poisoning, a major concern for public health worldwide. Despite their importance, taxonomic uncertainty within many genera of dinoflagellates is still high. The genus Coolia includes potentially harmful species and the diversity within the genus is just starting to become apparent. Methodology/Principal Findings:In the current study, cultures were established from strains of Coolia spp. isolated from the central Great Barrier Reef (GBR). Cultures were identified based on thecal plate morphology and analyses of sequences (18S, ITS and 28S) from the nuclear rRNA operon. We report that the central GBR harbors a high diversity of Coolia species, including two species known to be capable of toxin production (C. tropicalis and C. malayensis), as well as the non-toxic C. canariensis. The strain of C. canariensis isolated from the GBR may in fact be a cryptic species, closely related but nevertheless phylogenetically distinct from the strain on which the holotype of C. canariensis was based. We also found evidence of the occurrence of a cryptic species morphologically very similar to both C. malayensis and C. monotis. The consequences of taxonomic confusion within the genus are discussed. Conclusion/Significance:The central GBR region harbors a previously unreported high diversity of Coolia spp., including two species known to potentially produce toxins. The presence of a cryptic species of unknown toxicity highlights the importance of cryptic diversity within dinoflagellates. [ABSTRACT FROM AUTHOR]