Population subdivision in marine environments: the contributions of biogeography, geographical distance and discontinuous habitat to genetic differentiation in a blennioid fish, Axoclinus nigricaudus

A Comparative Phylogeography of Three Marine Species with Different PLD Modes Reveals Two Genetic Breaks across the Southern Caribbean Sea

The comparative phylogeography of marine species with contrasting dispersal potential across the southern Caribbean Sea was evaluated by the presence of two putative barriers: the Magdalena River plume (MRP) and the combination of the absence of a rocky bottom and the almost permanent upwelling in the La Guajira Peninsula (ARB + PUG). Three species with varying biological and ecological characteristics (i.e., dispersal potentials) that inhabit shallow rocky bottoms were selected: Cittarium pica (PLD < 6 days), Acanthemblemaria rivasi (PLD < 22 days), and Nerita tessellata (PLD > 60 days). We generated a set of SNPs for the three species using the ddRad-seq technique. Samples of each species were collected in five locations from Capurganá to La Guajira. For the first time, evidence of a phylogeographic break caused by the MRP is provided, mainly for A. rivasi (AMOVA: ΦCT = 0.420). The ARB + PUG barrier causes another break for A. rivasi (ΦCT = 0.406) and C. pica (ΦCT = 0.224). Three populations (K = 3) were identified for A. rivasi and C. pica, while N. tessellata presented one population (K = 1). The Mantel correlogram indicated that A. rivasi and C. pica fit the hierarchical population model, and only the A. rivasi and C. pica comparisons showed phylogeographic congruence. Our results demonstrate how the biological traits of these three species and the biogeographic barriers have influenced their phylogeographic structure.

Coastal dunefields maintain pre-Holocene genetic structure in a rocky shore red alga

Most intertidal algae have limited dispersal potential, and areas that lack hard substratum suitable for attachment are thus expected to isolate regional populations from each other. Here, we used nuclear and mitochondrial genetic data to compare genetic structure in two co-distributed intertidal red algae with different dispersal potential along the South African coastline. Gelidium pristoides is divided into a south-eastern and a south-western evolutionary lineage separated by extensive, continuous sandy shoreline habitat adjacent to coastal dunefields. In contrast, Hypnea spicifera is genetically homogeneous throughout its range. In G. pristoides, the genetic breaks are associated with contemporary coastal dunefields. The age of the divergence event suggests that this may reflect the effect of older dispersal barriers, and that genetic structure was subsequently maintained by the formation of contemporary coastal dunefields.

A multi-locus approach to elucidating the evolutionary history of the clingfish Tomicodon petersii (Gobiesocidae) in the Tropical Eastern Pacific

Marine species that are widely distributed in the Tropical Eastern Pacific (TEP) has served as a model for studying biogeographic patterns resulting from the effects of intraregional habitat discontinuities and oceanographic processes on the diversification and evolution of cryptobenthic reef fishes. Tomicodon petersii, a clingfish (Gobiesocidae) endemic to the TEP, is found on very shallow rocky reefs from central Mexico to northern Peru, and in the Cocos and Galapagos islands. We evaluated the effect of likely biogeographic barriers in different parts of the TEP on the diversification process of this species. We used one mitochondrial and three nuclear DNA markers from 112 individuals collected across the distribution range of T. petersii. Our phylogenetic results showed the samples constituted a monophyletic group, with three well-supported, allopatric subgroups: in the Mexican province, the Panamic province (from El Salvador to Ecuador), and the Galapagos Islands. The split between the Mexican and more southerly clades was estimated to occur at the end of the Miocene ca. 5.74 Mya, and the subsequent cladogenetic event separating the Galapagos population from the Panamic population at the junction of the Pliocene and Pleistocene, ca. 2.85 Mya. The species tree, Bayesian species delimitation tests (BPP), STACEY, and substantial genetic distances separating these three populations indicate that these three independent evolutionary units likely include two unnamed species. The cladogenetic events that promoted the formation of those genetically differentiated groups are consistent with disruptive effects on gene flow of habitat discontinuities and oceanographic processes along the mainland shoreline in the TEP and of ocean-island isolation, in conjunction with the species intrinsic life-history characteristics.

Genomic methods reveal independent demographic histories despite strong morphological conservatism in fish species

Human overexploitation of natural resources has placed conservation and management as one of the most pressing challenges in modern societies, especially in regards to highly vulnerable marine ecosystems. In this context, cryptic species are particularly challenging to conserve because they are hard to distinguish based on morphology alone, and thus it is often unclear how many species coexist in sympatry, what are their phylogenetic relationships and their demographic history. We answer these questions using morphologically similar species of the genus Mugil that are sympatric in the largest coastal Marine Protected Area in the Tropical Southwestern Atlantic marine province. Using a sub-representation of the genome, we show that individuals are assigned to five highly differentiated genetic clusters that are coincident with five mitochondrial lineages, but discordant with morphological information, supporting the existence of five species with conserved morphology in this region. A lack of admixed individuals is consistent with strong genetic isolation between sympatric species, but the most likely species tree suggests that in one case speciation has occurred in the presence of interspecific gene flow. Patterns of genetic diversity within species suggest that effective population sizes differ up to two-fold, probably reflecting differences in the magnitude of population expansions since species formation. Together, our results show that strong morphologic conservatism in marine environments can lead to species that are difficult to distinguish morphologically but that are characterized by an independent evolutionary history, and thus that deserve species-specific management strategies.

"A cleaner break": Genetic divergence between geographic groups and sympatric phenotypes revealed in ballan wrasse (Labrus bergylta)

Capture and long-distance translocation of cleaner fish to control lice infestations on marine salmonid farms has the potential to influence wild populations via overexploitation in source regions, and introgression in recipient regions. Knowledge of population genetic structure is therefore required. We studied the genetic structure of ballan wrasse, a phenotypically diverse and extensively used cleaner fish, from 18 locations in Norway and Sweden, and from Galicia, Spain, using 82 SNP markers. We detected two very distinct genetic groups in Scandinavia, northwestern and southeastern. These groups were split by a stretch of sandy beaches in southwest Norway, representing a habitat discontinuity for this rocky shore associated benthic egg-laying species. Wrasse from Galicia were highly differentiated from all Scandinavian locations, but more similar to northwestern than southeastern locations. Distinct genetic differences were observed between sympatric spotty and plain phenotypes in Galicia, but not in Scandinavia. The mechanisms underlying the geographic patterns between phenotypes are discussed, but not identified. We conclude that extensive aquaculture-mediated translocation of ballan wrasse from Sweden and southern Norway to western and middle Norway has the potential to mix genetically distinct populations. These results question the sustainability of the current cleaner fish practice.

Unraveling hierarchical genetic structure in a marine metapopulation: A comparison of three high-throughput genotyping approaches

Marine metapopulations often exhibit subtle population structure that can be difficult to detect. Given recent advances in high-throughput sequencing, an emerging question is whether various genotyping approaches, in concert with improved sampling designs, will substantially improve our understanding of genetic structure in the sea. To address this question, we explored hierarchical patterns of structure in the coral reef fish Elacatinus lori using a high-resolution approach with respect to both genetic and geographic sampling. Previously, we identified three putative E. lori populations within Belize using traditional genetic markers and sparse geographic sampling: barrier reef and Turneffe Atoll; Glover's Atoll; and Lighthouse Atoll. Here, we systematically sampled individuals at ~10 km intervals throughout these reefs (1,129 individuals from 35 sites) and sequenced all individuals at three sets of markers: 2,418 SNPs; 89 microsatellites; and 57 nonrepetitive nuclear loci. At broad spatial scales, the markers were consistent with each other and with previous findings. At finer spatial scales, there was new evidence of genetic substructure, but our three marker sets differed slightly in their ability to detect these patterns. Specifically, we found subtle structure between the barrier reef and Turneffe Atoll, with SNPs resolving this pattern most effectively. We also documented isolation by distance within the barrier reef. Sensitivity analyses revealed that the number of loci (and alleles) had a strong effect on the detection of structure for all three marker sets, particularly at small spatial scales. Taken together, these results illustrate empirically that high-throughput genotyping data can elucidate subtle genetic structure at previously-undetected scales in a dispersive marine fish.

Fluctuating thermal environments of shallow-water rocky reefs in the Gulf of California, Mexico

As part of a broad-scale study of the biogeography of rocky reefs in the Gulf of California, Mexico (GOC), we collected a continuous 1-yr temperature time series at ~5 m water depth at 16 sites spanning 5° of latitude and ~700 km along the western boundary of the basin. Throughout the region, thermal conditions were most variable in summer with fluctuations concentrated at diurnal and semi-diurnal frequencies, likely associated with solar and wind forcing and vertical water column oscillations forced by internal waves. Temperatures in winter were less variable than in summer, and minimum temperatures also differed among sites. Thermal variability integrated across the diurnal and semi-diurnal frequency bands was greatest near the Midriff Islands in the northern GOC and decreased toward the southern sites. Diurnal variability was greater than semi-diurnal variability at 13 of the 16 sites. A statistic-of-extremes analysis indicated shortest return times for cooling events in summer, and reef organisms at many of the sites may experience anomalous 2 to 5 °C cooling events multiple times per month. The significant extent of local temperature variability may play important roles in limiting species occurrences among sites across this biogeographic region.

Spatial differences in genetic diversity and northward migration suggest genetic erosion along the boreal caribou southern range limit and continued range retraction

With increasing human activities and associated landscape changes, distributions of terrestrial mammals become fragmented. These changes in distribution are often associated with reduced population sizes and loss of genetic connectivity and diversity (i.e., genetic erosion) which may further diminish a species' ability to respond to changing environmental conditions and lead to local population extinctions. We studied threatened boreal caribou (Rangifer tarandus caribou) populations across their distribution in Ontario/Manitoba (Canada) to assess changes in genetic diversity and connectivity in areas of high and low anthropogenic activity. Using data from >1,000 caribou and nine microsatellite loci, we assessed population genetic structure, genetic diversity, and recent migration rates using a combination of network and population genetic analyses. We used Bayesian clustering analyses to identify population genetic structure and explored spatial and temporal variation in those patterns by assembling networks based on R ST and F ST as historical and contemporary genetic edge distances, respectively. The Bayesian clustering analyses identified broad-scale patterns of genetic structure and closely aligned with the R ST network. The F ST network revealed substantial contemporary genetic differentiation, particularly in areas presenting contemporary anthropogenic disturbances and habitat fragmentation. In general, relatively lower genetic diversity and greater genetic differentiation were detected along the southern range limit, differing from areas in the northern parts of the distribution. Moreover, estimation of migration rates suggested a northward movement of animals away from the southern range limit. The patterns of genetic erosion revealed in our study suggest ongoing range retraction of boreal caribou in central Canada.

Dissimilarity of individual microsatellite profiles under different mutation models: Empirical approach

Microsatellites (simple sequence repeats, SSRs) still remain popular molecular markers for studying neutral genetic variation. Two alternative models outline how new microsatellite alleles evolve. Infinite alleles model (IAM) assumes that all possible alleles are equally likely to result from a mutation, while stepwise mutation model (SMM) describes microsatellite evolution as stepwise adding or subtracting single repeat units. Genetic relationships between individuals can be analyzed in higher precision when assuming the SMM scenario with allele size differences as a proxy of genetic distance. If population structure is not predetermined in advance, an empirical data analysis usually includes (a) estimating proximity between individual SSR profiles with a selected dissimilarity measure and (b) determining putative genetic structure of a given set of individuals using methods of clustering and/or ordination for the obtained dissimilarity matrix. We developed new dissimilarity indices between SSR profiles of haploid, diploid, or polyploid organisms assuming different mutation models and compared the performance of these indices for determining genetic structure with population data and with simulations. More specifically, we compared SMM with a constant or variable mutation rate at different SSR loci to IAM using data from natural populations of a freshwater bryozoan Cristatella mucedo (diploid), wheat leaf rust Puccinia triticina (dikaryon), and wheat powdery mildew Blumeria graminis (monokaryon). We show that inferences about population genetic structure are sensitive to the assumed mutation model. With simulations, we found that Bruvo's distance performs generally poorly, while the new metrics are capturing the differences in the genetic structure of the 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.

The evolutionary history of the goby Elacatinus puncticulatus in the tropical eastern pacific: Effects of habitat discontinuities and local environmental variability

Habitat discontinuities, temperature gradients, upwelling systems, and ocean currents, gyres and fronts, can affect distributions of species with narrow environmental tolerance or motility and influence the dispersal of pelagic larvae, with effects ranging from the isolation of adjacent populations to connections between them. The coast of the Tropical Eastern Pacific (TEP) is a highly dynamic environment, with various large gyres and upwelling systems, alternating currents and large rocky-habitat discontinuities, which may greatly influence the genetic connectivity of populations in different parts of the coast. Elacatinus puncticulatus is a cryptic, shallow-living goby that is distributed along the continental shore of virtually the entire TEP, which makes it a good model for testing the influence of these environmental characteristics in the molecular evolution of widespread species in this region. A multilocus phylogeny was used to evaluate the influence of habitat gaps, and oceanographic processes in the evolutionary history of E. puncticulatus throughout its geographical range in the TEP. Two well-supported allopatric clades (one with two allopatric subclades) were recovered, the geographic distribution of which does not correspond to any previously proposed major biogeographic provinces. These populations show strong genetic structure and substantial genetic distances between clades and sub-clades (cytb 0.8-7.3%), with divergence times between them ranging from 0.53 to 4.88 Mya, and recent population expansions dated at 170-130 Kya. The ancestral area of all populations appears to be the Gulf of Panama, while several isolation events have formed the phylogeographic patterns evident in this species. Local and regional oceanographic processes as well as habitat discontinuities have shaped the distribution patterns of the genetic lineages along the continental TEP. Large genetic distances, high genetic differentiation, and the results of species-tree and phylogenetic analyses indicate that E. puncticulatus comprises a complex of three allopatric species with an unusual geographic arrangement.

Different patterns of colonization of Oxalis alpina in the Sky Islands of the Sonoran desert via pollen and seed flow

Historical factors such as climatic oscillations during the Pleistocene epoch have dramatically impacted species distributions. Studies of the patterns of genetic structure in angiosperm species using molecular markers with different modes of inheritance contribute to a better understanding of potential differences in colonization and patterns of gene flow via pollen and seeds. These markers may also provide insights into the evolution of reproductive systems in plants. Oxalis alpina is a tetraploid, herbaceous species inhabiting the Sky Island region of the southwestern United States and northern Mexico. Our main objective in this study was to analyze the influence of climatic oscillations on the genetic structure of O. alpina and the impact of these oscillations on the evolutionary transition from tristylous to distylous reproductive systems. We used microsatellite markers and compared our results to a previous study using chloroplast genetic markers. The phylogeographic structure inferred by both markers was different, suggesting that intrinsic characteristics including the pollination system and seed dispersal have influenced patterns of gene flow. Microsatellites exhibited low genetic structure, showed no significant association between geographic and genetic distances, and all individual genotypes were assigned to two main groups. In contrast, chloroplast markers exhibited a strong association between geographic and genetic distance, had higher levels of genetic differentiation, and were assigned to five groups. Both types of DNA markers showed evidence of a northward expansion as a consequence of climate warming occurring in the last 10,000 years. The data from both types of markers support the hypothesis for several independent transitions from tristyly to distyly.

Role of oceanography in shaping the genetic structure in the North Pacific hake Merluccius productus

Determining the relative influence of biotic and abiotic factors on genetic connectivity among populations remains a major challenge in evolutionary biology and in the management and conservation of species. North Pacific hake (Merluccius productus) inhabits upwelling regions in the California Current ecosystem from the Gulf of California to the Gulf of Alaska. In this study, we examined mitochondrial DNA (mtDNA) and microsatellite variation to estimate levels of genetic differentiation of M. productus in relation to the role of oceanographic features as potential barriers to gene flow. Samples were obtained from nine sites spanning a large part of the geographic range of the species, from Puget Sound, Washington to Costa Rica. The microsatellite results revealed three genetically discrete populations: one spanning the eastern Pacific coast, and two apparently resident populations circumscribed to the Puget Sound and the northern Gulf of California (F_ST = 0.032, p = 0.036). Cytochrome b sequence data indicated that isolation between the Puget Sound and northern Gulf of California populations from the coastal Pacific were recent phenomena (18.5 kyr for Puget Sound and 40 kyr for the northern Gulf of California). Oceanographic data obtained from the Gulf of California support the hypothesis that permanent fronts within the region, and strong gradients at the entrance to the Gulf of California act as barriers to gene flow. A seascape genetics approach found significant genetic–environment associations, where the daytime sea surface temperature and chlorophyll concentrations were the best predictive variables for the observed genetic differentiation. Considering the potential causes of genetic isolation among the three populations, e.g. spawning areas in different latitudes associated with upwelling processes, oceanographic barriers, asymmetric migration and specialized diet, oceanographic barriers appear to be a likely mechanism restricting gene flow.

Eco-Evolutionary Processes Generating Diversity Among Bottlenose Dolphin, Tursiops truncatus, Populations off Baja California, Mexico

For highly mobile species that nevertheless show fine-scale patterns of population genetic structure, the relevant evolutionary mechanisms determining structure remain poorly understood. The bottlenose dolphin (Tursiops truncatus) is one such species, exhibiting complex patterns of genetic structure associated with local habitat dependence in various geographic regions. Here we studied bottlenose dolphin populations in the Gulf of California and Pacific Ocean off Baja California where habitat is highly structured to test associations between ecology, habitat dependence and genetic differentiation. We investigated population structure at a fine geographic scale using both stable isotope analysis (to assess feeding ecology) and molecular genetic markers (to assess population structure). Our results show that there are at least two factors affecting population structure for both genetics and feeding ecology (as indicated by stable isotope profiles). On the one hand there is a signal for the differentiation of individuals by ecotype, one foraging more offshore than the other. At the same time, there is differentiation between the Gulf of California and the west coast of Baja California, meaning that for example, nearshore ecotypes were both genetically and isotopically differentiated either side of the peninsula. We discuss these data in the context of similar studies showing fine-scale population structure for delphinid species in coastal waters, and consider possible evolutionary mechanisms.

Population differentiation or species formation across the Indian and the Pacific Oceans? An example from the brooding marine hydrozoan Macrorhynchia phoenicea

Assessing population connectivity is necessary to construct effective marine protected areas. This connectivity depends, among other parameters, inherently on species dispersal capacities. Isolation by distance (IBD) is one of the main modes of differentiation in marine species, above all in species presenting low dispersal abilities. This study reports the genetic structuring in the tropical hydrozoan Macrorhynchia phoenicea α (sensu Postaire et al., 2016a), a brooding species, from 30 sampling sites in the Western Indian Ocean and the Tropical Southwestern Pacific, using 15 microsatellite loci. At the local scale, genet dispersal relied on asexual propagation at short distance, which was not found at larger scales. Considering one representative per clone, significant positive F_IS values (from −0.327*** to 0.411***) were found within almost all sites. Gene flow was extremely low at all spatial scales, among sites within islands (<10 km distance) and among islands (100 to >11,000 km distance), with significant pairwise F_ST values (from 0.035*** to 0.645***). A general pattern of IBD was found at the Indo‐Pacific scale, but also within ecoregions in the Western Indian Ocean province. Clustering and network analyses identified each island as a potential independent population, while analysis of molecular variance indicated that population genetic differentiation was significant at small (within island) and intermediate (among islands within province) spatial scales. As shown by this species, a brooding life cycle might be corollary of the high population differentiation found in some coastal marine species, thwarting regular dispersal at distances more than a few kilometers and probably leading to high cryptic diversity, each island housing independent evolutionary lineages.

Phenotypic variation in dorsal fin morphology of coastal bottlenose dolphins (Tursiops truncatus) off Mexico

Geographic variation in external morphology is thought to reflect an interplay between genotype and the environment. Morphological variation has been well-described for a number of cetacean species, including the bottlenose dolphin (Tursiops truncatus). In this study we analyzed dorsal fin morphometric variation in coastal bottlenose dolphins to search for geographic patterns at different spatial scales. A total of 533 dorsal fin images from 19 available photo-identification catalogs across the three Mexican oceanic regions (Pacific Ocean n = 6, Gulf of California n = 6 and, Gulf of Mexico n = 7) were used in the analysis. Eleven fin shape measurements were analyzed to evaluate fin polymorphism through multivariate tests. Principal Component Analysis on log-transformed standardized ratios explained 94% of the variance. Canonical Discriminant Function Analysis on factor scores showed separation among most study areas (p < 0.05) with exception of the Gulf of Mexico where a strong morphometric cline was found. Possible explanations for the observed differences are related to environmental, biological and evolutionary processes. Shape distinction between dorsal fins from the Pacific and those from the Gulf of California were consistent with previously reported differences in skull morphometrics and genetics. Although the functional advantages of dorsal fin shape remains to be assessed, it is not unlikely that over a wide range of environments, fin shape may represent a trade-off among thermoregulatory capacity, hydrodynamic performance and the swimming/hunting behavior of the species.

Large-scale gene flow in the barnacle Jehlius cirratus and contrasts with other broadly-distributed taxa along the Chilean coast

We evaluate the population genetic structure of the intertidal barnacle Jehlius cirratus across a broad portion of its geographic distribution using data from the mitochondrial cytochrome oxidase I (COI) gene region. Despite sampling diversity from over 3,000 km of the linear range of this species, there is only slight regional structure indicated, with overall Φ CT of 0.036 (p < 0.001) yet no support for isolation by distance. While these results suggest greater structure than previous studies of J. cirratus had indicated, the pattern of diversity is still far more subtle than in other similarly-distributed species with similar larval and life history traits. We compare these data and results with recent findings in four other intertidal species that have planktotrophic larvae. There are no clear patterns among these taxa that can be associated with intertidal depth or other known life history traits.

Spatial and temporal genetic homogeneity of the Monterey Spanish mackerel, Scomberomorus concolor, in the Gulf of California

The genetic homogeneity of the Monterey Spanish mackerel Scomberomorus concolor population in the Gulf of California was confirmed using nine nuclear microsatellite loci in combination with mitochondrial cytochrome b gene sequences. Samples were collected from the upper and central Gulf areas, representing the two main biogeographical regions of the Gulf. The analyses support the existence of a single panmictic population of S. concolor inhabiting the Gulf of California which in terms of fishery management represents a single genetic stock. Additionally, the contemporary effective population size estimated for the S. concolor population (Ne = 3056.9) was high and similar to another pelagic species. The gene flow seems to be bidirectional between the upper and central Gulf, which coincides with the seasonal movements between both regions related to spawning and feeding activities. A population expansion event was detected, which agrees with a colonization-expansion hypothesis of the S. concolor population in the Gulf.

Habitat Discontinuities Separate Genetically Divergent Populations of a Rocky Shore Marine Fish

Habitat fragmentation has been suggested to be responsible for major genetic differentiations in a range of marine organisms. In this study, we combined genetic data and environmental information to unravel the relative role of geography and habitat heterogeneity on patterns of genetic population structure of corkwing wrasse (Symphodus melops), a rocky shore species at the northern limit of its distribution range in Scandinavia. Our results revealed a major genetic break separating populations inhabiting the western and southern coasts of Norway. This genetic break coincides with the longest stretch of sand in the whole study area, suggesting habitat fragmentation as a major driver of genetic differentiation of this obligate rocky shore benthic fish in Scandinavia. The complex fjords systems extending along the western coast of Norway appeared responsible for further regional genetic structuring. Our findings indicate that habitat discontinuities may lead to significant genetic fragmentation over short geographical distances, even for marine species with a pelagic larval phase, as for this rocky shore fish.

Panmixia in a Critically Endangered Fish: The Totoaba (Totoaba macdonaldi) in the Gulf of California

Conservation of the evolutionary legacy of endangered species is a key component for long-term persistence. Totoaba is a long-lived fish endemic to the Gulf of California and is considered critically endangered. There is currently a debate concerning its conservation status and whether it can be used as a fishery resource. Unfortunately, basic information on biological and genetic population structure of the species is lacking. We sampled 313 individuals and employed 16 microsatellite loci and 3 mitochondrial DNA markers (16S, 547 pb; COI, 619 pb; control region, 650 pb) to assess population structure and demography of totoaba in the Gulf of California, with samples from locations that encompass nearly all of its recognized geographic distribution. We could not reject a hypothesis of panmixia for totoaba, using nuclear or mitochondrial markers. Demographic analysis of mtDNA suggests a sudden population expansion model. The results have important implications for totoaba conservation because poaching is a significant conservation challenge and could have additive negative effects over the single population of totoaba in the Gulf of California.

Post-glacial population expansion of the Monterey Spanish mackerel Scomberomorus concolor in the Gulf of California

The level of genetic homogeneity and demographic history of the Monterey Spanish mackerel Scomberomorus concolor was assessed by analyses using sequences of the mitochondrial (mt)DNA-control region of samples from the two oceanographic regions of the Gulf of California in order to define the stock structure for this exploited vulnerable species. The data were consistent with a single homogeneous population and revealed the hallmark of fluctuations in population size; these fluctuations appear to have occurred during glacial events of the middle to late Pleistocene, which may in turn be related to the colonization and expansion of S. concolor populations in the gulf.

Population structure of Hirundichthys oxycephalus in the northwestern Pacific inferred from mitochondrial cytochrome oxidase I gene

BACKGROUND

Hirundichthys oxycephalus is an important flyingfish resource in eastern Taiwan and northwestern Japan. A substantial catch decline in Taiwan has caused serious concerns on stock status of the fish, prompting the government to impose a set of regulations on flyingfish egg fishery since 2008. However, the regulations were set in a precautionary manner, without considering the fundamental understanding of the population genetic structure. This study aims to investigate the population genetic structure of H. oxycephalus in the region based on mtDNA cytochrome oxidase I (COI) gene and to thus provide scientific information for sustainable management of theresource.

RESULTS

Tissue samples (156) from six localities of eastern Taiwan and western Japan were collected, and 616 bpof mtDNA COI gene were sequenced. Seventy haplotypes were determined, and the haplotype diversity and nucleotide diversity were estimated as 0.93% and 0.57%, respectively. Results of various statistical analyses suggested that the genetic differentiations among the six localities were small and most variation occurred within populations, indicating a high gene flow in the region with undergoing population expansion. Although the study showed that the fishes were genetically divided into two groups, the support was low and the separation was not geologically evident.

CONCLUSIONS

Thestudy revealed two groups of H.oxycephalus inthe northwestern Pacific Ocean. However, due to high gene flow, an association of either group to a spatial distribution was not observed, and so the two groups may be considered as one population. Thus, the results favored the conclusion that H. oxycephalus from eastern Taiwan and western Japan belong to the same population and, consequently, that the management unit of the current regulations only covering eastern Taiwan does not match the spatial structure of the population. Rather, the results suggest that joint efforts from countries within the population boundary are necessary to maintain a sustainable exploitation.

Seascape continuity plays an important role in determining patterns of spatial genetic structure in a coral reef fish

Detecting patterns of spatial genetic structure (SGS) can help identify intrinsic and extrinsic barriers to gene flow within metapopulations. For marine organisms such as coral reef fishes, identifying these barriers is critical to predicting evolutionary dynamics and demarcating evolutionarily significant units for conservation. In this study, we adopted an alternative hypothesis-testing framework to identify the patterns and predictors of SGS in the Caribbean reef fish Elacatinus lori. First, genetic structure was estimated using nuclear microsatellites and mitochondrial cytochrome b sequences. Next, clustering and network analyses were applied to visualize patterns of SGS. Finally, logistic regressions and linear mixed models were used to identify the predictors of SGS. Both sets of markers revealed low global structure: mitochondrial ΦST=0.12, microsatellite FST=0.0056. However, there was high variability among pairwise estimates, ranging from no differentiation between sites on contiguous reef (ΦST=0) to strong differentiation between sites separated by ocean expanses≥20 km (maximum ΦST=0.65). Genetic clustering and statistical analyses provided additional support for the hypothesis that seascape discontinuity, represented by oceanic breaks between patches of reef habitat, is a key predictor of SGS in E. lori. Notably, the estimated patterns and predictors of SGS were consistent between both sets of markers. Combined with previous studies of dispersal in E. lori, these results suggest that the interaction between seascape continuity and the dispersal kernel plays an important role in determining genetic connectivity within metapopulations.

Geographic distance affects dispersal of the patchy distributed greater long-tailed hamster (Tscherskia triton)

Dispersal is a fundamental process in ecology influencing the genetic structure and the viability of populations. Understanding how variable factors influence the dispersal of the population is becoming an important question in animal ecology. To date, geographic distance and geographic barriers are often considered as main factors impacting dispersal, but their effects are variable depending on different conditions. In general, geographic barriers affect more significantly than geographic distance on dispersal. In rapidly expanding populations, however, geographic barriers have less effect on dispersal than geographic distance. The effects of both geographic distance and geographic barriers in low-density populations with patchy distributions are poorly understood. By using a panel of 10 microsatellite loci we investigated the genetic structure of three patchy-distributed populations of the Greater long-tailed hamster (Tscherskia triton) from Raoyang, Guan and Shunyi counties of the North China Plain. The results showed that (i) high genetic diversity and differentiation exist in three geographic populations with patchy distributions; (ii) gene flow occurs among these three populations with physical barriers of Beijing city and Hutuo River, which potentially restricted the dispersal of the animal; (iii) the gene flow is negatively correlated with the geographic distance, while the genetic distance shows the positive correlation. Our results suggest that the effect of the physical barriers is conditional-dependent, including barrier capacity or individual potentially dispersal ability. Geographic distance also acts as an important factor affecting dispersal for the patchy distributed geographic populations. So, gene flow is effective, even at relatively long distances, in balancing the effect of geographic barrier in this study.

Contrasting genetic diversity patterns in two sister kelp species co-distributed along the coast of Brittany, France

We investigated patterns of genetic structure in two sister kelp species to explore how distribution width along the shore, zonation, latitudinal distribution and historical factors contribute to contrasting patterns of genetic diversity. We implemented a hierarchical sampling scheme to compare patterns of genetic diversity and structure in these two kelp species co-distributed along the coasts of Brittany (France) using a total of 12 microsatellites, nine for Laminaria hyperborea and 11 for Laminaria digitata, of which eight amplified in both species. The genetic diversity and connectivity of L. hyperborea populations were greater than those of L. digitata populations in accordance with the larger cross-shore distribution width along the coast and the greater depth occupied by L. hyperborea populations in contrast to L. digitata populations. In addition, marginal populations showed reduced genetic diversity and connectivity, which erased isolation-by-distance patterns in both species. As L. digitata encounters its southern range limit in southern Brittany (SBr) while L. hyperborea extends down to mid-Portugal, it was possible to distinguish the effect of habitat continuity from range edge effects. We found that L. digitata did not harbour high regional diversity at its southern edge, as expected in a typical rear edge, suggesting that refuges from the last glacial maximum for L. digitata were probably not located in SBr, but most likely further north. For both species, the highest levels of genetic diversity were found in the Iroise Sea and Morlaix Bay, the two regions in which they are being currently harvested. Preserving genetic diversity of these two foundation species in these areas should, thus, be a priority for the management of this resource in Brittany.

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.

Testing the genetic predictions of a biogeographical model in a dominant endemic Eastern Pacific coral (Porites panamensis) using a genetic seascape approach

The coral fauna of the Eastern Tropical Pacific (ETP) is depauperate and peripheral; hence, it has drawn attention to the factors allowing its survival. Here, we use a genetic seascape approach and ecological niche modeling to unravel the environmental factors correlating with the genetic variation of Porites panamensis, a hermatypic coral endemic to the ETP. Specifically, we test if levels of diversity and connectivity are higher among abundant than among depauperate populations, as expected by a geographically relaxed version of the Abundant Center Hypothesis (rel-ACH). Unlike the original ACH, referring to a geographical center of distribution of maximal abundance, the rel-ACH refers only to a center of maximum abundance, irrespective of its geographic position. The patterns of relative abundance of P. panamensis in the Mexican Pacific revealed that northern populations from Baja California represent its center of abundance; and southern depauperate populations along the continental margin are peripheral relative to it. Genetic patterns of diversity and structure of nuclear DNA sequences (ribosomal DNA and a single copy open reading frame) and five alloenzymatic loci partially agreed with rel-ACH predictions. We found higher diversity levels in peninsular populations and significant differentiation between peninsular and continental colonies. In addition, continental populations showed higher levels of differentiation and lower connectivity than peninsular populations in the absence of isolation by distance in each region. Some discrepancies with model expectations may relate to the influence of significant habitat discontinuities in the face of limited dispersal potential. Environmental data analyses and niche modeling allowed us to identify temperature, water clarity, and substrate availability as the main factors correlating with patterns of abundance, genetic diversity, and structure, which may hold the key to the survival of P. panamensis in the face of widespread environmental degradation.

Climate Change and Genetic Structure of Leading Edge and Rear End Populations in a Northwards Shifting Marine Fish Species, the Corkwing Wrasse (Symphodus melops)

One mechanism by which marine organisms may respond to climate shifts is range shifts. The corkwing wrasse (Symphodus melops) is a temperate fish species, inhabiting the coasts of Europe, that show strong indications of current as well as historical (ice-age) range shifts towards the north. Nine neutral microsatellite DNA markers were screened to study genetic signatures and spatial population structure over the entire geographic and thermal gradient of the species from Portugal to Norway. A major genetic break (F ST  = 0.159 average among pairs) was identified between Scandinavian and more southern populations, with a marked reduction (30% or more) in levels of genetic variability in Scandinavia. The break is probably related to bottleneck(s) associated with post-glacial colonization of the Scandinavian coasts, and indicates a lack of present gene flow across the North Sea. The lack of gene flow can most likely be attributed to the species' need for rocky substrate for nesting and a relatively short pelagic larval phase, limiting dispersal by ocean currents. These findings demonstrate that long-distance dispersal may be severely limited in the corkwing wrasse, and that successful range-shifts following present climate change may be problematic for this and other species with limited dispersal abilities, even in the seemingly continuous marine environment.

Genetic diversity and population structure of two lancelets along the coast of China

The western Pacific lancelet, once recognized as a monospecies, Branchiostoma belcheri, is a frequently used model in evolutionary and developmental studies, and researchers usually collect samples from the field without consideration of species identification and genetic divergence. However, recent studies found divergence of the lancelets from different localities and divided this monospecies into two separate species (S. belcheri and B. japonicum). To further estimate the genetic diversity of lancelet populations and the cause of their formation, we sampled 70 individuals from four major distribution areas along the coast of China, using both mitochondrial DNA and microsatellite markers in this investigation. Our results demonstrate that the two species possess extremely high genetic diversity at both mtDNA sequence level (h approaches 1.0) and microsatellite loci (He is above 0.8). Further demographic analysis reveals that the lancelets B. japonicum and B. belcheri underwent a recent historical population expansion at approximately 117,000 and 73,000 years ago respectively. Analyses on the population genetic structure revealed weak differentiation among different local populations. No evident differentiation was found among different local populations of the same species using mtDNA sequence data, but certain divergences among them were identified based on the microsatellite data. We suggest that discontinuous habitats may be responsible for the phylogeographic structure of the lancelets along China coasts.

Connectivity of the habitat-forming kelp, Ecklonia radiata within and among estuaries and open coast

With marine protected areas being established worldwide there is a pressing need to understand how the physical setting in which these areas are placed influences patterns of dispersal and connectivity of important marine organisms. This is particularly critical for dynamic and complex nearshore marine environments where patterns of genetic structure of organisms are often chaotic and uncoupled from broad scale physical processes. This study determines the influence of habitat heterogeneity (presence of estuaries) on patterns of genetic structure and connectivity of the common kelp, Ecklonia radiata. There was no genetic differentiation of kelp between estuaries and the open coast and the presence of estuaries did not increase genetic differentiation among open coast populations. Similarly, there were no differences in level of inbreeding or genetic diversity between estuarine and open coast populations. The presence of large estuaries along rocky coastlines does not appear to influence genetic structure of this kelp and factors other than physical heterogeneity of habitat are likely more important determinants of regional connectivity. Marine reserves are currently lacking in this bioregion and may be designated in the future. Knowledge of the factors that influence important habitat forming organisms such as kelp contribute to informed and effective marine protected area design and conservation initiatives to maintain resilience of important marine habitats.

Mitochondrial DNA variability of the pipefish Syngnathus abaster

This study provides data on the genetic structuring of the pipefish Syngnathus abaster in the western Mediterranean and Adriatic Seas. A total of 109 specimens were collected in brackish-water biotopes. The control region and three other regions of the mitochondrial genome were analysed. The most relevant result was the high genetic structuring found by Bayesian inference (BI), maximum likelihood (ML) and network analyses, which were consistent in showing three well-separated clusters of S. abaster populations. Furthermore, BI and ML did not support the monophyly of the taxon S. abaster. These results suggest the occurrence of a species complex in the study area, whose differentiation may have occurred since the Pleistocene. The results also show a very high genetic variability at the inter-population level, with no shared haplotypes among sites. Evolutionary forces due to the fragmented nature of the brackish-water habitats may account for the high genetic divergence found among the groups and populations. Finally, although dispersal by rafting over long distances may occasionally occur, this study suggests linear stepping-stone model of colonization to be most likely. The complexity of the results obtained suggests that further studies are needed to elucidate the phylogeny of S. abaster.

Triton's trident: cryptic Neogene divergences in a marine clam (Lasaea australis) correspond to Australia's three temperate biogeographic provinces

The southern coast of Australia is composed of three distinct biogeographic provinces distinguished primarily by intertidal community composition. Several ecological mechanisms have been proposed to explain their formation and persistence, but no consensus has been reached. The marine clam Lasaea australis is arguably the most common bivalve on southern Australian rocky shores and occurs in all three provinces. Here, we tested if this species exhibits cryptic genetic structuring corresponding to the provinces and if so, what mechanisms potentially drove its divergence. Variation in two mitochondrial genes (16S and COIII) and one nuclear gene (ITS2) was assayed to test for genetic structuring and to reconstruct the clam's phylogenetic history. Our results showed that L. australis is comprised of three cryptic mitochondrial clades, each corresponding almost perfectly to one of the three biogeographic provinces. Divergence time estimates place their cladogenesis in the Neogene. The trident-like topology and Neogene time frame of L. australis cladogenesis are incongruent with Quaternary vicariance predictions: a two-clade topology produced by Pleistocene Bass Strait land bridge formation. We hypothesize that the interaction of the Middle Miocene Climate Transition with the specific geography of the southern coastline of Australia was the primary cladogenic driver in this clam lineage. Additional in-depth studies of the endemic southern Australian marine biota across all three provinces are needed to establish the generality of this proposed older framework for regional cladogenesis.

What maintains the central North Pacific genetic discontinuity in Pacific herring?

Pacific herring show an abrupt genetic discontinuity in the central North Pacific that represents secondary contact between refuge populations previously isolated during Pleistocene glaciations. Paradoxically, high levels of gene flow produce genetic homogeneity among ocean-type populations within each group. Here, we surveyed variability in mtDNA control-region sequences (463 bp) and nine microsatellite loci in Pacific herring from sites across the North Pacific to further explore the nature of the genetic discontinuity around the Alaska Peninsula. Consistent with previous studies, little divergence (Φ(ST)  = 0.011) was detected between ocean-type populations of Pacific herring in the North West Pacific, except for a population in the Yellow Sea (Φ(ST)  = 0.065). A moderate reduction in genetic diversity for both mtDNA and microsatellites in the Yellow Sea likely reflects founder effects during the last colonization of this sea. Reciprocal monophyly between divergent mtDNA lineages (Φ(ST)  = 0.391) across the Alaska Peninsula defines the discontinuity across the North Pacific. However, microsatellites did not show a strong break, as eastern Bering Sea (EBS) herring were more closely related to NE Pacific than to NW Pacific herring. This discordance between mtDNA and microsatellites may be due to microsatellite allelic convergence or to sex-biased dispersal across the secondary contact zone. The sharp discontinuity between Pacific herring populations may be maintained by high-density blocking, competitive exclusion or hybrid inferiority.

Barriers to gene flow in the marine environment: insights from two common intertidal limpet species of the Atlantic and Mediterranean

Knowledge of the scale of dispersal and the mechanisms governing gene flow in marine environments remains fragmentary despite being essential for understanding evolution of marine biota and to design management plans. We use the limpets Patella ulyssiponensis and Patella rustica as models for identifying factors affecting gene flow in marine organisms across the North-East Atlantic and the Mediterranean Sea. A set of allozyme loci and a fragment of the mitochondrial gene cytochrome C oxidase subunit I were screened for genetic variation through starch gel electrophoresis and DNA sequencing, respectively. An approach combining clustering algorithms with clinal analyses was used to test for the existence of barriers to gene flow and estimate their geographic location and abruptness. Sharp breaks in the genetic composition of individuals were observed in the transitions between the Atlantic and the Mediterranean and across southern Italian shores. An additional break within the Atlantic cluster separates samples from the Alboran Sea and Atlantic African shores from those of the Iberian Atlantic shores. The geographic congruence of the genetic breaks detected in these two limpet species strongly supports the existence of transpecific barriers to gene flow in the Mediterranean Sea and Northeastern Atlantic. This leads to testable hypotheses regarding factors restricting gene flow across the study area.

Evolutionary history of almond tree domestication in the Mediterranean basin

Genetic diversity of contemporary domesticated species is shaped by both natural and human-driven processes. However, until now, little is known about how domestication has imprinted the variation of fruit tree species. In this study, we reconstruct the recent evolutionary history of the domesticated almond tree, Prunus dulcis, around the Mediterranean basin, using a combination of nuclear and chloroplast microsatellites [i.e. simple sequence repeat (SSRs)] to investigate patterns of genetic diversity. Whereas conservative chloroplast SSRs show a widespread haplotype and rare locally distributed variants, nuclear SSRs show a pattern of isolation by distance with clines of diversity from the East to the West of the Mediterranean basin, while Bayesian genetic clustering reveals a substantial longitudinal genetic structure. Both kinds of markers thus support a single domestication event, in the eastern side of the Mediterranean basin. In addition, model-based estimation of the timing of genetic divergence among those clusters is estimated sometime during the Holocene, a result that is compatible with human-mediated dispersal of almond tree out of its centre of origin. Still, the detection of region-specific alleles suggests that gene flow from relictual wild preglacial populations (in North Africa) or from wild counterparts (in the Near East) could account for a fraction of the diversity observed.