Bart Hellemans

Late glacial history of the cold‐adapted freshwater fish Cottus gobio, revealed by microsatellites

The distribution of genetic diversity at 10 highly polymorphic microsatellite loci within the European freshwater fish, Cottus gobio, L. was examined. The sampling range comprised a large geographical scale including lineages known to be highly divergent at both mitochondrial DNA (mtDNA) and allozymes. An analysis of genetic variability within populations showed that expected heterozygosity and allelic richness could be explained largely by current effective population sizes. Evidence was found, however, that historical processes predating the last major glaciation affected allelic richness. In addition to confirming the large‐scale patterns from earlier studies, the microsatellite data revealed new insights into recent processes by analysing genetic structure within ancient lineages defined by mtDNA data. Stepwise mutation model (SMM) and nonSMM‐based methods demonstrated a clear genetic structuring within the Northwestern European lineage comprising populations from Britain and Belgium, and within the Central European lineage populations from the rivers Danube, Elbe and Main. Supported by an analysis of genetic variability within populations these results showed that the bullhead populations most probably persisted throughout the last major glaciation within the British Isles and within the drainages of the rivers Elbe and Main. Such observations provide the first genetic evidence for a glacial refugium in such close proximity to the European glacial margins.

Phylogeography of the common goby, Pomatoschistus microps, with particular emphasis on the colonization of the Mediterranean and the North Sea

The phylogeographical patterns of a small marine fish, the common goby, Pomatoschistus microps, were assessed at 12 sites along the northeastern Atlantic coasts and the western Mediterranean Sea. A combination of two genetic markers was employed: cellulose acetate allozyme electrophoresis (CAGE) and sequence analysis of a 289 bp fragment of the mitochondrial locus cytochrome b. Both markers were congruent in revealing significant differences between samples (global FST = 0.247 for the allozymes and ΦST = 0.437 for the mitochondrial DNA data) and a pattern of isolation‐by‐distance. Phylogeographical analyses yielded a shallow branching structure with four groups. Three of those were confined to the Atlantic basin and showed a star‐like pattern. The fourth group contained a central haplotype occurring at the edges of the species’ distribution, accompanied by a few more rare variants, which were restricted to the Mediterranean Sea. A genetic break was observed around the British Isles, with distinct haplotypes dominating at either side of the English Channel. A significantly negative correlation between the degree of genetic diversity and latitude was recorded both for mitochondrial DNA (mtDNA) and allozymes in the Atlantic basin. Gene flow analysis suggested that recolonization of the North Sea and the coasts of western Scotland and Ireland may have taken place from a glacial refugium in the Southern Bight of the North Sea. These results are discussed in the perspective of possible postglacial migration routes of marine fish along the northeastern Atlantic coasts.

Timing of the population dynamics of bullhead Cottus gobio (Teleostei : Cottidae) during the Pleistocene

Abstract Pleistocene genetic structure of the bullhead, Cottus gobio, was evaluated across the western Palearctic using a 771‐bp long fragment of the mitochondrial control region in 123 individuals collected at 35 sites (data set I). In total, 59 haplotypes that differed at 73 positions (9.3%) were detected. Data analysis also included sequences from Englbrecht et al. (2000 ), thus increasing the sampling to a more comprehensive data set of 529 fish and 63 control region sequences of 482 bp (data set II). A minimum spanning and phylogenetic tree identified a seventh clade (Brittany–Loire) in addition to the previously identified six clades. The geographical range of the North Sea and Lower Rhine clades was considerably larger than thought previously. Haplotype diversity was generally low, and the total fixation index high (FST  = 0.49). Among‐group differentiation accounted for 41.7% (data set I) of the variation. Contiguous range expansions and restricted gene flow combined with isolation by distance, interspersed with past fragmentation characterize bullhead across its range. New is the knowledge that dated interglacial periods correlated with population expansions; river captures, proglacial lake systems and sea level played a significant role in the dispersal and expansion either in northern or southern direction. Hence it became possible to identify and date the colonization routes and putative palaeorefugia, most of which were located in Central and North‐west Europe. Glacial periods resulted in distinct fragmentation events and lineage sorting.

Evidence for isolation by time in the European eel (Anguilla anguilla L.)

Life history traits of highly vagile marine species, such as adult reproductive success and larval dispersal, are strongly determined by oceanographic and climatic forces. Nevertheless, marine organisms may show restricted dispersal in time and space. Patterns of isolation by distance (IBD) have been repeatedly observed in marine species. If spawning time is a function of geographical location, temporal and spatial isolation, can easily be confounded or misinterpreted. In this study, we aimed at discriminating between various forces shaping the genetic composition of recruiting juveniles of the European eel (Anguilla anguilla L.). By controlling for geographical variation, we assessed temporal variation and tested for possible isolation by time (IBT) between recruitment waves within and between years. Using 12 polymorphic allozyme and six variable microsatellite loci, we show that genetic differentiation was low (FST = 0.01–0.002) and significant between temporal samples. Regression analysis between genetic and temporal distance, was consistent with a subtle interannual pattern of IBT. Our data suggest that the population dynamics of the European eel may be governed by a double pattern of temporal variance in genetic composition: (i) a broad‐scale IBT of spawning cohorts, possibly as a consequence of the large migration loop in anguillids and strong variance in annual adult reproductive contribution; and (ii) a smaller‐scale variance in reproductive success (genetic patchiness) within cohorts among seasonally separated spawning groups, most likely originating from fluctuating oceanic and climatic forces. The consistency of both mechanisms remains to be verified with fine‐scale analyses of both spawning/migrating aged adults and their offspring to confirm the stochastic/deterministic nature of the IBT pattern in eel.

QTL for body weight, morphometric traits and stress response in European sea bass Dicentrarchus labrax

Natural mating and mass spawning in the European sea bass (Dicentrarchus labrax L., Moronidae, Teleostei) complicate genetic studies and the implementation of selective breeding schemes. We utilized a two-step experimental design for detecting QTL in mass-spawning species: 2122 offspring from natural mating between 57 parents (22 males, 34 females and one missing) phenotyped for body weight, eight morphometric traits and cortisol levels, had been previously assigned to parents based on genotypes of 31 DNA microsatellite markers. Five large full-sib families (five sires and two dams) were selected from the offspring (570 animals), which were genotyped with 67 additional markers. A new genetic map was compiled, specific to our population, but based on the previously published map. QTL mapping was performed with two methods: half-sib regression analysis (paternal and maternal) and variance component analysis accounting for all family relationships. Two significant QTL were found for body weight on linkage group 4 and 6, six significant QTL for morphometric traits on linkage groups 1B, 4, 6, 7, 15 and 23 and three suggestive QTL for stress response on linkage groups 3, 14 and 23. The QTL explained between 8% and 38% of phenotypic variance. The results are the first step towards identifying genes involved in economically important traits like body weight and stress response in European sea bass.

Phylogenetic relationships among Palearctic and Nearctic burbot (Lota lota): Pleistocene extinctions and recolonization

The burbot (Lota lota Linnaeus, 1758) is the only freshwater species from the cod family. Various taxonomic hypotheses were tested against molecular data by sequencing the mitochondrial cytochrome b locus of 120 burbot from 41 populations together with the related species Molva molva (ling) and Brosme brosme (tusk), which represented the other Lotinae genera. Within the genus Lota two distinct phylogroups were observed: one in North America south of the Great Slave Lakes (Lota lota maculosa) and one in Eurasia and the remainder of the Nearctic region (Lota lota lota). The burbot lineage separated 10 Myr BP from the other Lotinae, while the genetic variation within burbot appeared to be approximately 1 Myr old. However, fossil evidence suggested that burbot already existed in the Early Pliocene in Europe, from were it probably colonized North America in the Early Pleistocene. While Nearctic burbot survived climatic oscillations and diverged in several refugia, the Eurasian form became extinct or was reduced to a very small population. In the Late Pleistocene the species recolonized the Palearctic region to establish its present distribution range.

Gynogenesis in the African catfish (Clarias gariepinus). I: Induction of meiogynogenesis with thermal and pressure shocks

Abstract Meiogynogenesis of the African catfish (Clarias gariepinus) was optimized by altering the intensity, duration and timing of application of pressure, heat or cold shocks shortly after activation of the egg with sperm which had been genetically inactivated with UV irradiation. The distributions of the survival rates showed a unimodal (heat and hydrostatic pressure shock) or bimodal (cold shock) pattern and were specific to each treatment. The timing of the expulsion of the second polar body invariably occurred 3 to 5 min after fertilization. Of the three methods tested, cold (5°C, during 40 min at 3 min post-activation) and pressure (55.00 MPa during 2 min at 3 min post-activation) shocks turned out to be the most efficient in terms of survival rate. A relative hatching rate of up to 81% of gynogenetic diploid fish was achieved. The final proof of gynogenetic status was made with a paternity test by means of multiple-locus DNA fingerprinting.

A combined AFLP and microsatellite linkage map and pilot comparative genomic analysis of European sea bass Dicentrarchus labrax L.

European sea bass (Dicentrarchus labrax L., Moronidae, Teleostei) sustains a regional fishery and is commonly farmed in the Mediterranean basin, but has not undergone much long-term genetic improvement. An updated genetic linkage map of the European sea bass was constructed using 190 microsatellites, 176 amplified fragment length polymorphisms and two single nucleotide polymorphisms. From the 45 new microsatellite markers (including 31 type I markers) reported in this study, 28 were mapped. A total of 368 markers were assembled into 35 linkage groups. Among these markers, 28 represented type I (coding) markers, including those located within the peptide Y, SOX10, PXN1, ERA and TCRB genes (linkage groups 1, 7, 16, 17 and 27 respectively). The sex-averaged map spanned 1373.1 centimorgans (cM) of the genome. The female map measured 1380.0 cM, whereas the male map measured 1046.9 cM, leading to a female-to-male (F:M) recombination rate ratio of 1.32:1. The intermarker spacing of the second-generation linkage map of the European sea bass was 3.67 cM, which is smaller than that of the first-generation linkage map (5.03 cM). Comparative mapping of microsatellite flanking regions was performed with five model teleosts and this revealed a high percentage (33.6%) of evolutionarily conserved regions with the three-spined stickleback.

Evidence for fine-scale genetic structure and estuarine colonisation in a potential high gene flow marine goby ( Pomatoschistus minutus )

Marine fish seem to experience evolutionary processes that are expected to produce genetically homogeneous populations. We have assessed genetic diversity and differentiation in 15 samples of the sand goby Pomatoschistus minutus (Pallas, 1770) (Gobiidae, Teleostei) from four major habitats within the Southern Bight of the North Sea, using seven microsatellite and 13 allozyme loci. Despite its high dispersal potential, microsatellite loci revealed a moderate level of differentiation (overall FST=0.026; overall RST=0.058). Both hierarchical analysis of molecular variance and multivariate analysis revealed significant differentiation (P<0.01) between estuarine, coastal and marine samples with microsatellites, but not with allozymes. Comparison among the different estimators of differentiation (FST and RST) pointed to possible historical events and contemporary habitat fragmentation. Samples were assigned to two breeding units in the estuary and coastal region. Despite this classification, there were indications of a complex and dynamic spatiotemporal structure, which is, most likely, determined by historical events and local oceanic currents.

Mito-nuclear discordance in the degree of population differentiation in a marine goby

An increasing number of phylogeographic studies on marine species shows discordant patterns in the degree of population differentiation between nuclear and mitochondrial markers. To understand better which factors have the potential to cause these patterns of discordance in marine organisms, a population genetic study was realized on the sand goby Pomatoschistus minutus (Pallas 1770; Gobiidae, Teleostei). Sand gobies from eight European locations were genotyped at eight microsatellite markers. Microsatellites confirmed the global phylogeographical pattern of P. minutus observed with mitochondrial DNA (mtDNA) markers and nuclear allozyme markers. Three groups consistent with the mitochondrial lineages were defined (the Mediterranean, Iberian and North Atlantic groups) and indications of a recent founder event in the northern Baltic Sea were found. Nevertheless, differences in the degree of population differentiation between the nuclear and mitochondrial markers were large (global FST-values for microsatellites=0.0121; for allozymes=0.00831; for mtDNA=0.4293). Selection, sex-biased dispersal, homoplasy and a high effective population size are generally accepted as explanations for this mitonuclear discrepancy in the degree of population differentiation. In this study, selection on mtDNA and microsatellites, male-biased dispersal and homoplasy on microsatellite markers are unlikely to be a main cause for this discrepancy. The most likely reason for the discordant pattern is a recent demographical expansion of the sand goby, resulting in high effective population sizes slowing down the differentiation of nuclear DNA.