Thierry Backeljau

Dispersal and gene flow in free-living marine nematodes

Dispersal and gene flow determine connectivity among populations, and can be studied through population genetics and phylogeography. We here review the results of such a framework for free-living marine nematodes. Although field experiments have illustrated substantial dispersal in nematodes at ecological time scales, analysis of the genetic diversity illustrated the importance of priority effects, founder effects and genetic bottlenecks for population structuring between patches <1 km apart. In contrast, only little genetic structuring was observed within an estuary (<50 km), indicating that these small scale fluctuations in genetic differentiation are stabilized over deeper time scales through extensive gene flow. Interestingly, nematode species with contrasting life histories (extreme colonizers vs persisters) or with different habitat preferences (algae vs sediment) show similar, low genetic structuring. Finally, historical events have shaped the genetic pattern of marine nematodes and show that gene flow is restricted at large geographical scales. We also discuss the presence of substantial cryptic diversity in marine nematodes, and end with highlighting future important steps to further unravel nematode evolution and diversity.

Rapid large-scale evolutionary divergence in morphology and performance associated with exploitation of a different dietary resource

Although rapid adaptive changes in morphology on ecological time scales are now well documented in natural populations, the effects of such changes on whole-organism performance capacity and the consequences on ecological dynamics at the population level are often unclear. Here we show how lizards have rapidly evolved differences in head morphology, bite strength, and digestive tract structure after experimental introduction into a novel environment. Despite the short time scale (≈36 years) since this introduction, these changes in morphology and performance parallel those typically documented among species and even families of lizards in both the type and extent of their specialization. Moreover, these changes have occurred side-by-side with dramatic changes in population density and social structure, providing a compelling example of how the invasion of a novel habitat can evolutionarily drive multiple aspects of the phenotype.

Limits to gene flow in a cosmopolitan marine planktonic diatom

The role of geographic isolation in marine microbial speciation is hotly debated because of the high dispersal potential and large population sizes of planktonic microorganisms and the apparent lack of strong dispersal barriers in the open sea. Here, we show that gene flow between distant populations of the globally distributed, bloom-forming diatom species Pseudo-nitzschia pungens (clade I) is limited and follows a strong isolation by distance pattern. Furthermore, phylogenetic analysis implies that under appropriate geographic and environmental circumstances, like the pronounced climatic changes in the Pleistocene, population structuring may lead to speciation and hence may play an important role in diversification of marine planktonic microorganisms. A better understanding of the factors that control population structuring is thus essential to reveal the role of allopatric speciation in marine microorganisms.

Detection of the East and West African kdr mutation in Anopheles gambiae and Anopheles arabiensis from Uganda using a new assay based on FRET/Melt Curve analysis

BackgroundAppropriate monitoring of vector resistance to insecticides is an integral component of planning and evaluation of insecticide use in malaria control programmes. The malaria vectors Anopheles gambiae s.s. and Anopheles arabiensis have developed resistance to pyrethroid insecticides as a result of a mechanism conferring reduced nervous system sensitivity, better known as knockdown resistance (kdr). In An. gambiae s.s. and An. arabiensis, two different substitutions in the para-type sodium channel, a L1014F substitution common in West Africa and a L1014S replacement found in Kenya, are linked with kdr. Two different allele-specific polymerase chain reactions (AS-PCR) are needed to detect these known kdr mutations. However, these AS-PCR assays rely on a single nucleotide polymorphism mismatch, which can result in unreliable results.MethodsHere, a new assay for the detection of knockdown resistance in An. gambiae s.s. and An. arabiensis based on Fluorescence Resonance Energy Transfer/Melt Curve analysis (FRET/MCA) is presented and compared with the existing assays.ResultsThe new FRET/MCA method has the important advantage of detecting both kdr alleles in one assay. Moreover, results show that the FRET/MCA is more reliable and more sensitive than the existing AS-PCR assays and is able to detect new genotypes. By using this technique, the presence of the East African kdr mutation (L1014S) is shown for the first time in An. arabiensis specimens from Uganda. In addition, a new kdr genotype is reported in An. gambiae s.s. from Uganda, where four An. gambiaes.s. mosquitoes possess both, the West (L1014F) and East (L1014S) African kdr allele, simultaneously.ConclusionThe presence of both kdr mutations in the same geographical region shows the necessity of a reliable assay that enables to detect both mutations in one single assay. Hence, this new assay based on FRET/MCA will improve the screening of the kdr frequencies in An. gambiae s.s. and An. arabiensis.

Synchronization of hatching date with budburst of individual host trees (Quercus robur) in the winter moth (Operophtera brumata) and its fitness consequences

1. Due to varying selection pressures among host individuals, herbivorous insects may show local adaptation at the individual level. In the winter moth, synchrony of egg hatching and host tree (Quercus robur L.) budburst may have important fitness consequences and, therefore, may result in a local adaptation to the hosts phenology. However, relatively high gene flow levels may disrupt such a fine-scale adaption. 2. We determined hatching dates of clutches laid by females collected during copulation, and variation in budburst dates of the trees on which these couples were collected. 3. In our study area, budburst date showed substantial variation within years and areas, as the first and last tree to start leafing were separated by up to 26 days. Relative budburst dates of individual trees were constant between years. Average hatching dates differed significantly among clutches and varied over 30 days, a range comparable to that of budburst. 4. Hatching dates were positively associated with budburst of the respective trees indicating individual synchrony, which may be mediated by two mechanisms. First, early active adults are captured more frequently on early trees and late adults on late trees. As adult activity period is positively correlated with date of egg hatch, early clutches will tend to be laid on earlier trees while late clutches will be laid more frequently on late trees. Secondly, a presumed active choice mechanism may additionally increases individual host synchronization. Yet, the latter requires further study. 5. Our data support that synchronization of larval hatching with host budburst is adaptive as it increases adult size and thus expected fitness.

Phylogeography of the **Rhabditis (Pellioditis) marina** species complex: evidence for long-distance dispersal, and for range expansions and restricted gene flow in the northeast Atlantic

Pinpointing processes that structure the geographical distribution of genetic diversity of marine species and lead to speciation is challenging because of the lack of obvious dispersal barriers and the likelihood of substantial (passive) dispersal in oceans. In addition, cryptic radiations with sympatric distributions abound in marine species, challenging the allopatric speciation mechanism. Here, we present a phylogeographical study of the marine nematode species complex Rhabditis (Pellioditis) marina to investigate processes shaping genetic structure and speciation. Rhabditis (P.) marina lives on decaying macroalgae in the intertidal, and may therefore disperse over considerable distances. Rhabditis (P.) marina consists of several cryptic species sympatrically distributed at a local scale. Genetic variation in the COI gene was screened in 1362 specimens from 45 locations around the world. Two nuclear DNA genes (ITS and D2D3) were sequenced to infer phylogenetic species. We found evidence for ten sympatrically distributed cryptic species, seven of which show a strong genetic structuring. A historical signature showed evidence for restricted gene flow with occasional long‐distance dispersal and range expansions pre‐dating the last glacial maximum. Our data also point to a genetic break around the British Isles and a contact zone in the Southern Bight of the North Sea. We provide evidence for the transoceanic distribution of at least one cryptic species (PmIII) and discuss the dispersal capacity of marine nematodes. The allopatric distribution of some intraspecific phylogroups and of closely related cryptic species points to the potential for allopatric speciation in R. (P.) marina.

Exploring the Use of Cytochrome Oxidase c Subunit 1 (COI) for DNA Barcoding of Free-Living Marine Nematodes

Background The identification of free-living marine nematodes is difficult because of the paucity of easily scorable diagnostic morphological characters. Consequently, molecular identification tools could solve this problem. Unfortunately, hitherto most of these tools relied on 18S rDNA and 28S rDNA sequences, which often lack sufficient resolution at the species level. In contrast, only a few mitochondrial COI data are available for free-living marine nematodes. Therefore, we investigate the amplification and sequencing success of two partitions of the COI gene, the M1-M6 barcoding region and the I3-M11 partition. Methodology Both partitions were analysed in 41 nematode species from a wide phylogenetic range. The taxon specific primers for the I3-M11 partition outperformed the universal M1-M6 primers in terms of amplification success (87.8% vs. 65.8%, respectively) and produced a higher number of bidirectional COI sequences (65.8% vs 39.0%, respectively). A threshold value of 5% K2P genetic divergence marked a clear DNA barcoding gap separating intra- and interspecific distances: 99.3% of all interspecific comparisons were >0.05, while 99.5% of all intraspecific comparisons were <0.05 K2P distance. Conclusion The I3-M11 partition reliably identifies a wide range of marine nematodes, and our data show the need for a strict scrutiny of the obtained sequences, since contamination, nuclear pseudogenes and endosymbionts may confuse nematode species identification by COI sequences.

Rangewide phylogeography of a terrestrial slug in Europe: evidence for Alpine refugia and rapid colonization after the Pleistocene glaciations

Intraspecific phylogeographical patterns largely depend on the life history traits of a species. Especially species with a high degree of cold tolerance, limited requirements towards habitat preferences, and relatively low active dispersal capacities may have responded in a different way to the Pleistocene climatological fluctuations than the majority of taxa studied so far. To evaluate this possibility, we studied Arion fuscus (Müller, 1774), a common and widespread European terrestrial slug, from 88 locations (N = 964). Sequence variation was assessed for fragments of the mitochondrial 16S rDNA and COI genes by means of single‐strand conformation polymorphisms (SSCP) and subsequent DNA sequencing. Additionally, eight allozyme loci were scored in 843 individuals. Phylogenetic analysis revealed the presence of two major evolutionary lineages, one in the Balkan region and another in the Alps and the rest of Europe. The sequence divergence between the two lineages was limited (3.3%), but gene flow between the regions was absent, suggesting that the two regions have been isolated since the late Pliocene or early Pleistocene. Allozyme differentiation among geographical regions and mitochondrial DNA (mtDNA) lineages was low. The geographical patterns observed in our data showed that (i) haplotype and nucleotide diversities are very low in northern Europe, suggesting that single haplotypes rapidly colonized large areas; (ii) recently expanded haplotype clades have restricted distribution ranges, suggesting that current gene flow is low; and (iii) genetic diversity in the Alps is much higher than in other regions and estimated past gene flow from the Eastern Alps to other regions was high, suggesting that this was a refugial zone during the Pleistocene. This full‐range phylogeography suggests the existence of an alternative refugial zone, situated north of the refugial areas currently recognized in most other taxa.

Extreme mtDNA divergences in a terrestrial slug (Gastropoda, Pulmonata, Arionidae): accelerated evolution, allopatric divergence and secondary contact

Extremely high levels of intraspecific mtDNA differences in pulmonate gastropods have been reported repeatedly and several hypotheses to explain them have been postulated. We studied the phylogeny and phylogeography of 51 populations (n = 843) of the highly polymorphic terrestrial slug Arion subfuscus ( Draparnaud, 1805 ) across its native distribution range in Western Europe. By combining the analysis of single stranded conformation polymorphisms (SSCP) and nucleotide sequencing, we obtained individual sequence data for a fragment of the mitochondrial 16S rDNA and a fragment of the nuclear ITS1. Additionally, five polymorphic allozyme loci were scored. Based on the 16S rDNA phylogeny, five monophyletic haplotype groups with sequence divergences of 9–21% were found. Despite this deep mitochondrial divergence, the haplotype groups were not monophyletic for the nuclear ITS1 fragment and haplotype group‐specific allozyme alleles were not found. Although there is evidence for an accelerated mtDNA clock, the divergence among the haplotype groups is older than the Pleistocene and their current allopatric ranges probably reflect allopatric divergence and glacial survival in separate refugia from which different post‐glacial colonization routes were established. A range‐overlap of two mtDNA groups (S1 and S2, 21% sequence divergence) stretched from Central France and Belgium up to the North of the British Isles. The nuclear data suggest that this secondary contact resulted in hybridization between the allopatrically diverged groups. Therefore, it seems that, at least for two of the groups, the deep mtDNA divergence was only partially accompanied by the formation of reproductive isolation.

Identifying Insects with Incomplete DNA Barcode Libraries, African Fruit Flies (Diptera: Tephritidae) as a Test Case

We propose a general working strategy to deal with incomplete reference libraries in the DNA barcoding identification of species. Considering that (1) queries with a large genetic distance with their best DNA barcode match are more likely to be misidentified and (2) imposing a distance threshold profitably reduces identification errors, we modelled relationships between identification performances and distance thresholds in four DNA barcode libraries of Diptera (n = 4270), Lepidoptera (n = 7577), Hymenoptera (n = 2067) and Tephritidae (n = 602 DNA barcodes). In all cases, more restrictive distance thresholds produced a gradual increase in the proportion of true negatives, a gradual decrease of false positives and more abrupt variations in the proportions of true positives and false negatives. More restrictive distance thresholds improved precision, yet negatively affected accuracy due to the higher proportions of queries discarded (viz. having a distance query-best match above the threshold). Using a simple linear regression we calculated an ad hoc distance threshold for the tephritid library producing an estimated relative identification error <0.05. According to the expectations, when we used this threshold for the identification of 188 independently collected tephritids, less than 5% of queries with a distance query-best match below the threshold were misidentified. Ad hoc thresholds can be calculated for each particular reference library of DNA barcodes and should be used as cut-off mark defining whether we can proceed identifying the query with a known estimated error probability (e.g. 5%) or whether we should discard the query and consider alternative/complementary identification methods.