Demographic history has shaped the strongly differentiated corkwing wrasse populations in Northern Europe

Understanding the biological processes involved in genetic differentiation and divergence between populations within species is a pivotal aim in evolutionary biology. One particular phenomenon that requires clarification is the maintenance of genetic barriers despite the high potential for gene flow in the marine environment. Such patterns have been attributed to limited dispersal or local adaptation, and to a lesser extent to the demographic history of the species. The corkwing wrasse (Symphodus melops) is an example of a marine fish species where regions of particular strong divergence are observed. One such genetic break occurred at a surprisingly small spatial scale (F_ST ~0.1), over a short coastline (<60 km) in the North Sea‐Skagerrak transition area in southwestern Norway. Here, we investigate the observed divergence and purported reproductive isolation using genome resequencing. Our results suggest that historical events during the post‐glacial recolonization route can explain the present population structure of the corkwing wrasse in the northeast Atlantic. While the divergence across the break is strong, we detected ongoing gene flow between populations over the break suggesting recent contact or negative selection against hybrids. Moreover, we found few outlier loci and no clear genomic regions potentially being under selection. We concluded that neutral processes and random genetic drift e.g., due to founder events during colonization have shaped the population structure in this species in Northern Europe. Our findings underline the need to take into account the demographic process in studies of divergence processes.     

A continuous genome assembly of the corkwing wrasse (Symphodus melops)

The wrasses (Labridae) are one of the most successful and species-rich families of the Perciformes order of teleost fish. Its members display great morphological diversity, and occupy distinct trophic levels in coastal waters and coral reefs. The cleaning behaviour displayed by some wrasses, such as corkwing wrasse (Symphodus melops), is of particular interest for the salmon aquaculture industry to combat and control sea lice infestation as an alternative to chemicals and pharmaceuticals. There are still few genome assemblies available within this fish family for comparative and functional studies, despite the rapid increase in genome resources generated during the past years. Here, we present a highly continuous genome assembly of the corkwing wrasse using PacBio SMRT sequencing (x28.8) followed by error correction with paired-end Illumina data (x132.9). The present genome assembly consists of 5040 contigs (N50?=?461,652?bp) and a total size of 614 Mbp, of which 8.5% of the genome sequence encode known repeated elements. The genome assembly covers 94.21% of highly conserved genes across ray-finned fish species. We find evidence for increased copy numbers specific for corkwing wrasse possibly highlighting diversification and adaptive processes in gene families including N-linked glycosylation (ST8SIA6) and stress response kinases (HIPK1). By comparative analyses, we discover that de novo repeats, often not properly investigated during genome annotation, encode hundreds of immune-related genes. This new genomic resource, together with the ballan wrasse (Labrus bergylta), will allow for in-depth comparative genomics as well as population genetic analyses for the understudied wrasses.     

Impact of harvesting cleaner fish for salmonid aquaculture assessed from replicated coastal marine protected areas

Wrasse (Labridae) fisheries have increased markedly in Norway since 2010. Wrasse are being used as cleaner fish in salmonid aquaculture to control sea-lice infestations. However, fundamental knowledge on the demography and abundance of the targeted wrasse populations in Norwegian waters is lacking, and the consequences of harvesting at the current intensity have not been assessed. Here, we compared catch per unit effort (CPUE), size, age and sex ratio of goldsinny wrasse (Ctenolabrus rupestris) and corkwing wrasse (Symphodus melops) between marine protected areas (MPAs) and control areas open for fishing at four localities on the Skagerrak coast in Southern Norway. The CPUE of goldsinny larger than the minimum size limit was 33–65% higher within MPAs, while for corkwing three of four MPAs had higher CPUE with the relative difference between MPAs and control areas ranging from ?16% to 92%. Moreover, corkwing, but not goldsinny, was significantly older and larger within MPAs than in control areas. Sex ratios did not differ between MPAs and control areas for either species. Our study suggests that despite its short history, the wrasse fisheries have considerable impacts on the target populations and, further, that small MPAs hold promise as a management tool for maintaining natural population sizes and size structure. Goldsinny, being a smaller-sized species, also seems to benefit from the traditional minimum size limit management tool, which applies outside MPAs.     

Cleaner wrasse influence habitat selection of young damselfish

The presence of bluestreak cleaner wrasse, Labroides dimidiatus, on coral reefs increases total abundance and biodiversity of reef fishes. The mechanism(s) that cause such shifts in population structure are unclear, but it is possible that young fish preferentially settle into microhabitats where cleaner wrasse are present. As a first step to investigate this possibility, we conducted aquarium experiments to examine whether settlement-stage and young juveniles of ambon damselfish, Pomacentrus amboinensis, selected a microhabitat near a cleaner wrasse (adult or juvenile). Both settlement-stage (0 d post-settlement) and juvenile (~5 weeks post-settlement) fish spent a greater proportion of time in a microhabitat adjacent to L. dimidiatus than in one next to a control fish (a non-cleaner wrasse, Halichoeres melanurus) or one where no fish was present. This suggests that cleaner wrasse may serve as a positive cue during microhabitat selection. We also conducted focal observations of cleaner wrasse and counts of nearby damselfishes (1 m radius) to examine whether newly settled fish obtained direct benefits, in the form of cleaning services, from being near a cleaner wrasse. Although abundant, newly settled recruits (<20 mm total length) were rarely (2 %) observed being cleaned in 20 min observations compared with larger damselfishes (58 %). Individual damselfish that were cleaned were significantly larger than the median size of the surrounding nearby non-cleaned conspecifics; this was consistent across four species. The selection by settlement-stage fish of a microhabitat adjacent to cleaner wrasse in the laboratory, despite only being rarely cleaned in the natural environment, suggests that even rare cleaning events and/or indirect benefits may drive their settlement choices. This behaviour may also explain the decreased abundance of young fishes on reefs from which cleaner wrasse had been experimentally removed. This study reinforces the potentially important role of mutualism during the processes of settlement and recruitment of young reef fishes.     

Interbreeding between local and translocated populations of a cleaner fish in an experimental mesocosm predicts risk of disrupted local adaptation

Translocation of organisms within or outside its native range carries the risk of modifying the community of the recipient ecosystems and induces gene flow between locally adapted populations or closely related species. In this study, we evaluated the genetic consequences of large‐scale translocation of cleaner wrasses that has become a common practice within the salmon aquaculture industry in northern Europe to combat sea lice infestation. A major concern with this practice is the potential for hybridization of escaped organisms with the local, recipient wrasse population, and thus potentially introduce exogenous alleles and breaking down coadapted gene complexes in local populations. We investigated the potential threat for such genetic introgressions in a large seminatural mesocosm basin. The experimental setting represented a simulated translocation of corkwing wrasse (Symphodus melops) that occurs on a large scale in the Norwegian salmon industry. Parentage assignment analysis of mesocosm's offspring revealed 30% (195 out of 651 offspring) interbreeding between the two populations, despite their being genetically (F_ST = 0.094, p < 0.05) and phenotypically differentiated. Moreover, our results suggest that reproductive fitness of the translocated western population doubled that of the local southern population. Our results confirm that human translocations may overcome the impediments imposed by natural habitat discontinuities and urge for immediate action to manage the genetic resources of these small benthic wrasses.     

Geographic patterns in morphometric and genetic variation for coyote populations with emphasis on southeastern coyotes

Prior to 1900, coyotes (Canis latrans) were restricted to the western and central regions of North America, but by the early 2000s, coyotes became ubiquitous throughout the eastern United States. Information regarding morphological and genetic structure of coyote populations in the southeastern United States is limited, and where data exist, they are rarely compared to those from other regions of North America. We assessed geographic patterns in morphology and genetics of coyotes with special consideration of coyotes in the southeastern United States. Mean body mass of coyote populations increased along a west‐to‐east gradient, with southeastern coyotes being intermediate to western and northeastern coyotes. Similarly, principal component analysis of body mass and linear body measurements suggested that southeastern coyotes were intermediate to western and northeastern coyotes in body size but exhibited shorter tails and ears from other populations. Genetic analyses indicated that southeastern coyotes represented a distinct genetic cluster that differentiated strongly from western and northeastern coyotes. We postulate that southeastern coyotes experienced lower immigration from western populations than did northeastern coyotes, and over time, genetically diverged from both western and northeastern populations. Coyotes colonizing eastern North America experienced different selective pressures than did stable populations in the core range, and we offer that the larger body size of eastern coyotes reflects an adaptation that improved dispersal capabilities of individuals in the expanding range.     

Presence of cleaner wrasse increases the recruitment of damselfishes to coral reefs

Mutualisms affect the biodiversity, distribution and abundance of biological communities. However, ecological processes that drive mutualism-related shifts in population structure are often unclear and must be examined to elucidate how complex, multi-species mutualistic networks are formed and structured. In this study, we investigated how the presence of key marine mutualistic partners can drive the organisation of local communities on coral reefs. The cleaner wrasse, Labroides dimidiatus, removes ectoparasites and reduces stress hormones for multiple reef fish species, and their presence on coral reefs increases fish abundance and diversity. Such changes in population structure could be driven by increased recruitment of larval fish at settlement, or by post-settlement processes such as modified levels of migration or predation. We conducted a controlled field experiment to examine the effect of cleaners on recruitment processes of a common group of reef fishes, and showed that small patch reefs (61–285 m²) with cleaner wrasse had higher abundances of damselfish recruits than reefs from which cleaner wrasse had been removed over a 12-year period. However, the presence of cleaner wrasse did not affect species diversity of damselfish recruits. Our study provides evidence of the ecological processes that underpin changes in local population structure in the presence of a key mutualistic partner.     

Cleaning by the wrasse Thalassoma noronhanum, with two records of predation by its grouper client Cephalopholis fulva

The Noronha wrasse Thalassoma noronhanum was recorded cleaning 19 client fish species at Fernando de Noronha Archipelago, o. north–eastern Brazil. The preferred clients were non–dangerous, mostly planktivorous species, whereas the potentially dangerous, predatory species were rarely cleaned. T. noronhanum acts as a cleaner in two distinct ecological situations, at and outside the cleaning stations, and attends different client species in each of them. Potentially dangerous clients were mostly attended outside the cleaning stations. Many attacks and two instances of predation on the cleaner wrasse by the grouper client Cephalopholis fulva were recorded. The attacks occurred on individual wrasses foraging near the bottom outside the cleaning stations.     

Primula farinosa in Denmark; genetic diversity and population management

During the past centuries Danish populations of Primula farinosa have seriously declined in number. We investigated the genetic structure and genetic diversity of plants of seven populations from two different regions, Zealand and Bornholm in Denmark, using three AFLP markers. Two populations from nearby Scania, Sweden were included as reference. We found 54 unambiguously polymorphic loci. The genetic structure analysis suggested division of the 268 plants into three distinct groups, to a large extent matching the geographical distribution of the populations. Analysis of molecular variance (AMOVA) indicated significant genetic differentiation of 67% within populations and 33% among the populations. Our results suggest that genetic differentiation among regions and unique local genetic diversity should carefully be considered in future conservation attempts if we are to maintain as much genetic variation as possible. We present a historical overview of the decline in Danish populations and discuss conservation management and restoration strategies.     

Molecular phylogeny and geography of Korean medaka fish (Oryzias latipes)

The phylogeny and geography of the medaka (Oryzias latipes) populations of Korea were investigated by analyzing sequence data for the mitochondrial control region. From the 41 haplotypes including 25 Korean haplotypes detected in 64 Korean specimens and data for the Japanese and Chinese populations, phylogenetic and nested clade analyses were executed to examine the phylogeny of haplogroups and the relation of the genetic architecture of the haplotypes to the historical geography of the Korean medaka fish. The analyses suggest that there are two very distinct lineages of Korean medaka, and that these result from reproductive isolation mechanisms due to geographic barriers. The southeastern lineage has experienced recent range expansion to the western region. The northwestern lineage, sister to Chinese populations, showed evidence of internal range expansion with shared haplotypes.     

The viability of a plant “on the edge”: Glaucium flavum (Papaveraceae) in Norway

The prospects for persistence of Glaucium flavum in Norway were evaluated by integrating data concerning spatio‐temporal distribution, plant growth and development, reproduction, dispersal, seed bank and genetic variability. The distribution of Glaucium flavum in the Oslo Fjord region is much more restricted now than before. Out of the formerly recorded 35 populations, only eight populations were present in 1998, all in the outer parts of the fjord. The species is a short‐lived iteroparous perennial in Norway. This growth habit gives high reproductive output from fertile individuals. G. flavum is self‐compatible and can easily reproduce even in small populations. The ability to disperse long distances on water is limited. The observed seed bank will therefore enhance survival. No variation in isoenzymes was observed neither within nor among the G. flavum populations in the Oslo Fjord region. The species should be considered Endangered in Norway, according to the criteria given by IUCN.     

QTL mapping of freezing tolerance: links to fitness and adaptive trade‐offs

Local adaptation, defined as higher fitness of local vs. nonlocal genotypes, is commonly identified in reciprocal transplant experiments. Reciprocally adapted populations display fitness trade‐offs across environments, but little is known about the traits and genes underlying fitness trade‐offs in reciprocally adapted populations. We investigated the genetic basis and adaptive significance of freezing tolerance using locally adapted populations of Arabidopsis thaliana from Italy and Sweden. Previous reciprocal transplant studies of these populations indicated that subfreezing temperature is a major selective agent in Sweden. We used quantitative trait locus (QTL) mapping to identify the contribution of freezing tolerance to previously demonstrated local adaptation and genetic trade‐offs. First, we compared the genomic locations of freezing tolerance QTL to those for previously published QTL for survival in Sweden, and overall fitness in the field. Then, we estimated the contributions to survival and fitness across both field sites of genotypes at locally adaptive freezing tolerance QTL. In growth chamber studies, we found seven QTL for freezing tolerance, and the Swedish genotype increased freezing tolerance for five of these QTL. Three of these colocalized with locally adaptive survival QTL in Sweden and with trade‐off QTL for overall fitness. Two freezing tolerance QTL contribute to genetic trade‐offs across environments for both survival and overall fitness. A major regulator of freezing tolerance, CBF2, is implicated as a candidate gene for one of the trade‐off freezing tolerance QTL. Our study provides some of the first evidence of a trait and gene that mediate a fitness trade‐off in nature.     

Caribbean Cleaning Gobies Prefer Client Ectoparasites Over Mucus

If cooperation often involves investment, then what specific conditions prevent selection from acting on cheaters that do not invest? The mutualism between the Indo‐Pacific cleaner wrasse Labroides dimidiatus and its reef fish clients has been a model system to study conflicts of interest and their resolution. These cleaners prefer client mucus over ectoparasites – that is, they prefer to cheat – but punishment and partner switching by clients enforce cooperative behaviour by cleaners. By contrast, clients of Caribbean cleaning gobies (Elacatinus spp.) do not to use punishment or partner switching. Here, we test the hypothesis that the behavioural differences between these two cleaner fish systems are caused by differences in cleaner foraging preferences. In foraging choice experiments, we offered broadstripe cleaning gobies Elacatinus prochilos client‐derived parasitic isopods, client mucus and a control food item. The cleaning gobies significantly preferred ectoparasites over mucus or the control item, which contrasts with cleaner wrasses. We propose that the low level of cleaner–client conflict arising from cleaning goby foraging preferences explains the observed lack of strategic partner control behaviour in the clients of cleaning gobies.     

Embryonic development in ballan wrasse Labrus bergylta

Eight primary embryonic developmental stages were assigned to eggs of ballan wrasse Labrus bergylta using key morphological features following standardized nomenclature: Ia, Ib, II, III, IV, V, VI and VI+, reared from single family clutches under comparable environmental conditions in Ireland and Norway. Development in L. bergylta is typical of demersal marine finfish species with a short egg stage. Hatching occurred c. 123 h post-fertilization (hpf) equivalent to 62·5 degree days at 12·2 ± 1·10° C (mean ±s.d.), after which the larvae swam intermittently near the surface of the water column.     

Cleaning behaviour of the goldsinny (Pisces,Labridae) in Swedish waters

The goldsinny, $\text{Ctenolabrus}$$\text{rupestris}$, on the Swedish west-coast has been shown to be a facultative cleaner with the ballan wrasse, $\text{Labrus}$$\text{berggylta}$, as the cleanee. Twenty-four cleanings were observed during 1975 to 1981. The cleaning pattern of the symbiosis is described and the low rate of cleaners within the species is discussed.     

Y chromosome haplotype distribution of brown bears (Ursus arctos) in Northern Europe provides insight into population history and recovery

High‐resolution, male‐inherited Y‐chromosomal markers are a useful tool for population genetic analyses of wildlife species, but to date have only been applied in this context to relatively few species besides humans. Using nine Y‐chromosomal STRs and three Y‐chromosomal single nucleotide polymorphism markers (Y‐SNPs), we studied whether male gene flow was important for the recent recovery of the brown bear (Ursus arctos) in Northern Europe, where the species declined dramatically in numbers and geographical distribution during the last centuries but is expanding now. We found 36 haplotypes in 443 male extant brown bears from Sweden, Norway, Finland and northwestern Russia. In 14 individuals from southern Norway from 1780 to 1920, we found two Y chromosome haplotypes present in the extant population as well as four Y chromosome haplotypes not present among the modern samples. Our results suggested major differences in genetic connectivity, diversity and structure between the eastern and the western populations in Northern Europe. In the west, our results indicated that the recovered population originated from only four male lineages, displaying pronounced spatial structuring suggestive of large‐scale population size increase under limited male gene flow within the western subpopulation. In the east, we found a contrasting pattern, with high haplotype diversity and admixture. This first population genetic analysis of male brown bears shows conclusively that male gene flow was not the main force of population recovery.     

Population genomics reveals multiple drivers of population differentiation in a sex‐role‐reversed pipefish

A major goal of molecular ecology is to identify the causes of genetic and phenotypic differentiation among populations. Population genomics is suitably poised to tackle these key questions by diagnosing the evolutionary mechanisms driving divergence in nature. Here, we set out to investigate the evolutionary processes underlying population differentiation in the Gulf pipefish, Syngnathus scovelli. We sampled approximately 50 fish from each of 12 populations distributed from the Gulf coast of Texas to the Atlantic coast of Florida and performed restriction‐site‐associated DNA sequencing to identify SNPs throughout the genome. After imposing quality and stringency filters, we selected a panel of 6348 SNPs present in all 12 populations, 1753 of which were not physically linked. We identified a genome‐wide pattern of isolation by distance, in addition to a more substantial genetic break separating populations in the Gulf of Mexico from those in the Atlantic. We also used several divergence outlier approaches and tests for genotype–environment correlations to identify 400 SNPs putatively involved in local adaptation. Patterns of phenotypic differentiation and variation diverged from the overall genomic pattern, suggesting that selection, phenotypic plasticity or demographic factors may be shaping phenotypes in distinct populations. Overall, our results suggest that population divergence is driven by a variety of factors in S. scovelli, including neutral processes and selection on multiple traits.     

Nuclear microsatellites reveal contrasting patterns of genetic structure between western and southeastern European populations of the common ash (Fraxinus excelsior L.)

To determine extant patterns of population genetic structure in common ash and gain insight into postglacial recolonization processes, we applied multilocus-based Bayesian approaches to data from 36 European populations genotyped at five nuclear microsatellite loci. We identified two contrasting patterns in terms of population genetic structure: (1) a large area from the British Isles to Lithuania throughout central Europe constituted effectively a single deme, whereas (2) strong genetic differentiation occurred over short distances in Sweden and southeastern Europe. Concomitant geographical variation was observed in estimates of allelic richness and genetic diversity, which were lowest in populations from southeastern Europe, that is, in regions close to putative ice age refuges, but high in western and central Europe, that is, in more recently recolonized areas. We suggest that in southeastern Europe, restricted postglacial gene flow caused by a rapid expansion of refuge populations in a mountainous topography is responsible for the observed strong genetic structure. In contrast, admixture of previously differentiated gene pools and high gene flow at the onset of postglacial recolonization of western and central Europe would have homogenized the genetic structure and raised the levels of genetic diversity above values in the refuges.     

Implications of life-history strategies for a new wrasse fishery

The variety of life-history patterns exhibited by the five species of wrasse common in Northern Europe are reviewed. The two larger wrasse species, the ballan, Labrus bergylta Ascanius, and cuckoo, Labrus mixtus (L.), are exploited through sport angling. The three smaller species, the corkwing, Symphodus (Crenilabrus) melops (L.), rock cook, Centrolabrus exoletus (L.) and goldsinny, Ctenolabrus rupestris (L.), are being exploited by a new fishery for use as parasite cleaners of farmed salmon.
The nature of salmon farming limits the wrasse fishery to a minimum size, restricted areas and the warmer months of the year. The fishery may be expected to alter population structure through selective removal of larger fish. Removal of dominant territorial males may affect social structures and removal of nest-guarding males would reduce egg survival. Quantitative models incorporating stock size and fishery requirements are now required.