Ecological genomics of adaptation to unpredictability in experimental rotifer populations.

Elucidating the genetic basis of phenotypic variation in response to different environments is key to understanding how populations evolve. Facultatively sexual rotifers can develop adaptive responses to fluctuating environments. In a previous evolution experiment, diapause-related traits changed rapidly in response to two selective regimes (predictable vs unpredictable) in laboratory populations of the rotifer Brachionus plicatilis. Here, we investigate the genomic basis of adaptation to environmental unpredictability in these experimental populations. We identified and genotyped genome-wide polymorphisms in 169 clones from both selective regimes after seven cycles of selection using genotyping by sequencing (GBS). Additionally, we used GBS data from the 270 field clones from which the laboratory populations were established. This GBS dataset was used to identify candidate SNPs under selection. A total of 76 SNPs showed divergent selection, three of which are candidates for being under selection in the particular unpredictable fluctuation pattern studied. Most of the remaining SNPs showed strong signals of adaptation to laboratory conditions. Furthermore, a genotype-phenotype association approach revealed five SNPs associated with two key life-history traits in the adaptation to unpredictability. Our results contribute to elucidating the genomic basis for adaptation to unpredictable environments and lay the groundwork for future evolution studies in rotifers.

Founder effects drive the genetic structure of passively dispersed aquatic invertebrates.

Populations of passively dispersed organisms in continental aquatic habitats typically show high levels of neutral genetic differentiation despite their high dispersal capabilities. Several evolutionary factors, including founder events, local adaptation, and life cycle features such as high population growth rates and the presence of propagule banks, have been proposed to be responsible for this paradox. Here, we have modeled the colonization process to assess the impact of migration rate, population growth rate, population size, local adaptation and life-cycle features on the population genetic structure in these organisms. Our simulations show that the strongest effect on population structure are persistent founder effects, resulting from the interaction of a few population founders, high population growth rates, large population sizes and the presence of diapausing egg banks. In contrast, the role of local adaptation, genetic hitchhiking and migration is limited to small populations in these organisms. Our results indicate that local adaptation could have different impact on genetic structure in different groups of zooplankters.

Genomic signatures of local adaptation to the degree of environmental predictability in rotifers.

Environmental fluctuations are ubiquitous and thus essential for the study of adaptation. Despite this, genome evolution in response to environmental fluctuations -and more specifically to the degree of environmental predictability- is still unknown. Saline lakes in the Mediterranean region are remarkably diverse in their ecological conditions, which can lead to divergent local adaptation patterns in the inhabiting aquatic organisms. The facultatively sexual rotifer Brachionus plicatilis shows diverging local adaptation in its life-history traits in relation to estimated environmental predictability in its habitats. Here, we used an integrative approach -combining environmental, phenotypic and genomic data for the same populations- to understand the genomic basis of this diverging adaptation. Firstly, a novel draft genome for B. plicatilis was assembled. Then, genome-wide polymorphisms were studied using genotyping by sequencing on 270 clones from nine populations in eastern Spain. As a result, 4,543 high-quality SNPs were identified and genotyped. More than 90 SNPs were found to be putatively under selection with signatures of diversifying and balancing selection. Over 140 SNPs were correlated with environmental or phenotypic variables revealing signatures of local adaptation, including environmental predictability. Putative functions were associated to most of these SNPs, since they were located within annotated genes. Our results reveal associations between genomic variation and the degree of environmental predictability, providing genomic evidence of adaptation to local conditions in natural rotifer populations.

Barcoding rotifer biodiversity in Mediterranean ponds using diapausing egg banks.

The biodiversity of Mediterranean freshwater bodies is among the most threatened worldwide; therefore, its accurate estimation is an urgent issue. However, traditional methods are likely to underestimate freshwater zooplankton biodiversity due to its high species seasonality and cryptic diversity. We test the value of applying DNA barcoding to diapausing egg banks, in combination with the creation of a reference collection of DNA barcodes using adult individual samples, to characterize rotifer communities. We use monogonont rotifers from two lakes in Doñana National Park and one from Ruidera Natural Park in Spain as models to create a reference collection of DNA barcodes for taxonomically diagnosed adult individuals sampled from the water column, to compare with the sequences obtained from individual eggs from the diapausing egg banks. We apply two different approaches to carry out DNA taxonomy analyses, the generalized mixed Yule coalescent method (GMYC) and the Automatic Barcode Gap Discovery (ABGD), to the obtained sequences and to publicly available rotifer sequences. We obtained a total of 210 new rotifer COI sequences from all three locations (151 diapausing eggs and 59 adults). Both GMYC and ABGD generated the same 35 operational taxonomic units (OTUs), revealing four potential cryptic species. Most sequences obtained from diapausing eggs (85%) clustered with sequences obtained from morphologically diagnosed adults. Our approach, based on a single sediment sample, retrieved estimates of rotifer biodiversity higher than or similar to those of previous studies based on a number of seasonal samples. This study shows that DNA barcoding of diapausing egg banks is an effective aid to characterize rotifer diversity in Mediterranean freshwater bodies.

A new molecular diagnostic tool for surveying and monitoring Triops cancriformis populations.

The tadpole shrimp, Triops cancriformis, is a freshwater crustacean listed as endangered in the UK and Europe living in ephemeral pools. Populations are threatened by habitat destruction due to land development for agriculture and increased urbanisation. Despite this, there is a lack of efficient methods for discovering and monitoring populations. Established macroinvertebrate monitoring methods, such as net sampling, are unsuitable given the organism's life history, that include long lived diapausing eggs, benthic habits and ephemerally active populations. Conventional hatching methods, such as sediment incubation, are both time consuming and potentially confounded by bet-hedging hatching strategies of diapausing eggs. Here we develop a new molecular diagnostic method to detect viable egg banks of T. cancriformis, and compare its performance to two conventional monitoring methods involving diapausing egg hatching. We apply this method to a collection of pond sediments from the Wildfowl & Wetlands Trust Caerlaverock National Nature Reserve, which holds one of the two remaining British populations of T. cancriformis. DNA barcoding of isolated eggs, using newly designed species-specific primers for a large region of mtDNA, was used to estimate egg viability. These estimates were compared to those obtained by the conventional methods of sediment and isolation hatching. Our method outperformed the conventional methods, revealing six ponds holding viable T. cancriformis diapausing egg banks in Caerlaverock. Additionally, designed species-specific primers for a short region of mtDNA identified degraded, inviable eggs and were used to ascertain the levels of recent mortality within an egg bank. Together with efficient sugar flotation techniques to extract eggs from sediment samples, our molecular method proved to be a faster and more powerful alternative for assessing the viability and condition of T. cancriformis diapausing egg banks.

Isolation mediates persistent founder effects on zooplankton colonisation in new temporary ponds.

Understanding the colonisation process in zooplankton is crucial for successful restoration of aquatic ecosystems. Here, we analyzed the clonal and genetic structure of the cyclical parthenogenetic rotifer Brachionus plicatilis by following populations established in new temporary ponds during the first three hydroperiods. Rotifer populations established rapidly after first flooding, although colonisation was ongoing throughout the study. Multilocus genotypes from 7 microsatellite loci suggested that most populations (10 of 14) were founded by few clones. The exception was one of the four populations that persisted throughout the studied hydroperiods, where high genetic diversity in the first hydroperiod suggested colonisation from a historical egg bank, and no increase in allelic diversity was detected with time. In contrast, in another of these four populations, we observed a progressive increase of allelic diversity. This population became less differentiated from the other populations suggesting effective gene flow soon after its foundation. Allelic diversity and richness remained low in the remaining two, more isolated, populations, suggesting little gene flow. Our results highlight the complexity of colonisation dynamics, with evidence for persistent founder effects in some ponds, but not in others, and with early immigration both from external source populations, and from residual, historical diapausing egg banks.

Laboratory generation of new parthenogenetic lineages supports contagious parthenogenesis in Artemia.

Contagious parthenogenesis-a process involving rare functional males produced by a parthenogenetic lineage which mate with coexisting sexual females resulting in fertile parthenogenetic offspring-is one of the most striking mechanisms responsible for the generation of new parthenogenetic lineages. Populations of the parthenogenetic diploid brine shrimp Artemia produce fully functional males in low proportions. The evolutionary role of these so-called Artemia rare males is, however, unknown. Here we investigate whether new parthenogenetic clones could be obtained in the laboratory by mating these rare males with sexual females. We assessed the survival and sex ratio of the hybrid ovoviviparous offspring from previous crosses between rare males and females from all Asiatic sexual species, carried out cross-mating experiments between F1 hybrid individuals to assess their fertility, and estimated the viability and the reproductive mode of the resulting F2 offspring. Molecular analysis confirmed the parentage of hybrid parthenogenetic F2. Our study documents the first laboratory synthesis of new parthenogenetic lineages in Artemia and supports a model for the contagious spread of parthenogenesis. Our results suggest recessive inheritance but further experiments are required to confirm the likelihood of the contagious parthenogenesis model.

Origin and genetic diversity of diploid parthenogenetic Artemia in Eurasia.

There is wide interest in understanding how genetic diversity is generated and maintained in parthenogenetic lineages, as it will help clarify the debate of the evolution and maintenance of sexual reproduction. There are three mechanisms that can be responsible for the generation of genetic diversity of parthenogenetic lineages: contagious parthenogenesis, repeated hybridization and microorganism infections (e.g. Wolbachia). Brine shrimps of the genus Artemia (Crustacea, Branchiopoda, Anostraca) are a good model system to investigate evolutionary transitions between reproductive systems as they include sexual species and lineages of obligate parthenogenetic populations of different ploidy level, which often co-occur. Diploid parthenogenetic lineages produce occasional fully functional rare males, interspecific hybridization is known to occur, but the mechanisms of origin of asexual lineages are not completely understood. Here we sequenced and analysed fragments of one mitochondrial and two nuclear genes from an extensive set of populations of diploid parthenogenetic Artemia and sexual species from Central and East Asia to investigate the evolutionary origin of diploid parthenogenetic Artemia, and geographic origin of the parental taxa. Our results indicate that there are at least two, possibly three independent and recent maternal origins of parthenogenetic lineages, related to A. urmiana and Artemia sp. from Kazakhstan, but that the nuclear genes are very closely related in all the sexual species and parthenogegetic lineages except for A. sinica, who presumable took no part on the origin of diploid parthenogenetic strains. Our data cannot rule out either hybridization between any of the very closely related Asiatic sexual species or rare events of contagious parthenogenesis via rare males as the contributing mechanisms to the generation of genetic diversity in diploid parthenogenetic Artemia lineages.

Bird migratory flyways influence the phylogeography of the invasive brine shrimp Artemia franciscana in its native American range.

Since Darwin's time, waterbirds have been considered an important vector for the dispersal of continental aquatic invertebrates. Bird movements have facilitated the worldwide invasion of the American brine shrimp Artemia franciscana, transporting cysts (diapausing eggs), and favouring rapid range expansions from introduction sites. Here we address the impact of bird migratory flyways on the population genetic structure and phylogeography of A. franciscana in its native range in the Americas. We examined sequence variation for two mitochondrial gene fragments (COI and 16S for a subset of the data) in a large set of population samples representing the entire native range of A. franciscana. Furthermore, we performed Mantel tests and redundancy analyses (RDA) to test the role of flyways, geography and human introductions on the phylogeography and population genetic structure at a continental scale. A. franciscana mitochondrial DNA was very diverse, with two main clades, largely corresponding to Pacific and Atlantic populations, mirroring American bird flyways. There was a high degree of regional endemism, with populations subdivided into at least 12 divergent, geographically restricted and largely allopatric mitochondrial lineages, and high levels of population structure (Φ ST of 0.92), indicating low ongoing gene flow. We found evidence of human-mediated introductions in nine out of 39 populations analysed. Once these populations were removed, Mantel tests revealed a strong association between genetic variation and geographic distance (i.e., isolation-by-distance pattern). RDA showed that shared bird flyways explained around 20% of the variance in genetic distance between populations and this was highly significant, once geographic distance was controlled for. The variance explained increased to 30% when the factor human introduction was included in the model. Our findings suggest that bird-mediated transport of brine shrimp propagules does not result in substantial ongoing gene flow; instead, it had a significant historical role on the current species phylogeography, facilitating the colonisation of new aquatic environments as they become available along their main migratory flyways.

Multiple global radiations in tadpole shrimps challenge the concept of 'living fossils'.

'Living fossils', a phrase first coined by Darwin, are defined as species with limited recent diversification and high morphological stasis over long periods of evolutionary time. Morphological stasis, however, can potentially lead to diversification rates being underestimated. Notostraca, or tadpole shrimps, is an ancient, globally distributed order of branchiopod crustaceans regarded as 'living fossils' because their rich fossil record dates back to the early Devonian and their morphology is highly conserved. Recent phylogenetic reconstructions have shown a strong biogeographic signal, suggesting diversification due to continental breakup, and widespread cryptic speciation. However, morphological conservatism makes it difficult to place fossil taxa in a phylogenetic context. Here we reveal for the first time the timing and tempo of tadpole shrimp diversification by inferring a robust multilocus phylogeny of Branchiopoda and applying Bayesian divergence dating techniques using reliable fossil calibrations external to Notostraca. Our results suggest at least two bouts of global radiation in Notostraca, one of them recent, so questioning the validity of the 'living fossils' concept in groups where cryptic speciation is widespread.

High lability of sexual system over 250 million years of evolution in morphologically conservative tadpole shrimps.

BACKGROUND: Sexual system is a key factor affecting the genetic diversity, population structure, genome structure and the evolutionary potential of species. The sexual system androdioecy - where males and hermaphrodites coexist in populations - is extremely rare, yet is found in three crustacean groups, barnacles, a genus of clam shrimps Eulimnadia, and in the order Notostraca, the tadpole shrimps. In the ancient crustacean order Notostraca, high morphological conservatism contrasts with a wide diversity of sexual systems, including androdioecy. An understanding of the evolution of sexual systems in this group has been hampered by poor phylogenetic resolution and confounded by the widespread occurrence of cryptic species. Here we use a multigene supermatrix for 30 taxa to produce a comprehensive phylogenetic reconstruction of Notostraca. Based on this phylogenetic reconstruction we use character mapping techniques to investigate the evolution of sexual systems. We also tested the hypothesis that reproductive assurance has driven the evolution of androdioecy in Notostraca. RESULTS: Character mapping analysis showed that sexual system is an extremely flexible trait within Notostraca, with repeated shifts between gonochorism and androdioecy, the latter having evolved a minimum of five times. In agreement with the reproductive assurance hypothesis androdioecious notostracans are found at significantly higher latitudes than gonochoric ones indicating that post glacial re-colonisation may have selected for the higher colonisation ability conferred by androdioecy. CONCLUSIONS: In contrast to their conserved morphology, sexual system in Notostraca is highly labile and the rare reproductive mode androdioecy has evolved repeatedly within the order. Furthermore, we conclude that this lability of sexual system has been maintained for at least 250 million years and may have contributed to the long term evolutionary persistence of Notostraca. Our results further our understanding of the evolution of androdioecy and indicate that reproductive assurance is a recurrent theme involved in the evolution of this sexual system.

Long-term coexistence of rotifer cryptic species.

Despite their high morphological similarity, cryptic species often coexist in aquatic habitats presenting a challenge in the framework of niche differentiation theory and coexistence mechanisms. Here we use a rotifer species complex inhabiting highly unpredictable and fluctuating salt lakes to gain insights into the mechanisms involved in stable coexistence in cryptic species. We combined molecular barcoding surveys of planktonic populations and paleogenetic analysis of diapausing eggs to reconstruct the current and historical coexistence dynamics of two highly morphologically similar rotifer species, B. plicatilis and B. manjavacas. In addition, we carried out laboratory experiments using clones isolated from eight lakes where both species coexist to explore their clonal growth responses to salinity, a challenging, highly variable and unpredictable condition in Mediterranean salt lakes. We show that both species have co-occurred in a stable way in one lake, with population fluctuations in which no species was permanently excluded. The seasonal occurrence patterns of the plankton in two lakes agree with laboratory experiments showing that both species differ in their optimal salinity. These results suggest that stable species coexistence is mediated by differential responses to salinity and its fluctuating regime. We discuss the role of fluctuating salinity and a persistent diapausing egg banks as a mechanism for species coexistence in accordance with the 'storage effect'.

Widespread secondary contact and new glacial refugia in the halophilic rotifer Brachionus plicatilis in the Iberian Peninsula.

Small aquatic organisms harbour deep phylogeographic patterns and highly structured populations even at local scales. These patterns indicate restricted gene flow, despite these organisms' high dispersal abilities, and have been explained by a combination of (1) strong founder effects due to rapidly growing populations and very large population sizes, and (2) the development of diapausing egg banks and local adaptation, resulting in low effective gene flow, what is known as the Monopolization hypothesis. In this study, we build up on our understanding of the mitochondrial phylogeography of the halophilic rotifer Brachionus plicatilis in the Iberian Peninsula by both increasing the number of sampled ponds in areas where secondary contact is likely and doubling sample sizes. We analyzed partial mitochondrial sequences of 252 individuals. We found two deep mitochondrial DNA lineages differing in both their genetic diversity and the complexity of their phylogeographic structure. Our analyses suggest that several events of secondary contact between clades occurred after their expansion from glacial refugia. We found a pattern of isolation-by-distance, which we interpret as being the result of historical colonization events. We propose the existence of at least one glacial refugium in the SE of the Iberian Peninsula. Our findings challenge predictions of the Monopolization hypothesis, since coexistence (i.e., secondary contact) of divergent lineages in some ponds in the Iberian Peninsula is common. Our results indicate that phylogeographic structures in small organisms can be very complex and that gene flow between diverse lineages after population establishment can indeed occur.

Evolutionary origin and phylogeography of the diploid obligate parthenogen Artemia parthenogenetica (Branchiopoda: Anostraca).

BACKGROUND: Understanding the evolutionary origin and the phylogeographic patterns of asexual taxa can shed light on the origin and maintenance of sexual reproduction. We assessed the geographic origin, genetic diversity, and phylogeographic history of obligate parthenogen diploid Artemia parthenogenetica populations, a widespread halophilic crustacean. METHODOLOGY/PRINCIPAL FINDINGS: We analysed a partial sequence of the Cytochrome c Oxidase Subunit I mitochondrial gene from an extensive set of localities (including Eurasia, Africa, and Australia), and examined their phylogeographic patterns and the phylogenetic relationships of diploid A. parthenogenetica and its closest sexual relatives. Populations displayed an extremely low level of mitochondrial genetic diversity, with one widespread haplotype shared by over 79% of individuals analysed. Phylogenetic and phylogeographic analyses indicated a multiple and recent evolutionary origin of diploid A. parthenogenetica, and strongly suggested that the geographic origin of parthenogenesis in Artemia was in Central Asia. Our results indicate that the maternal sexual ancestors of diploid A. parthenogenetica were an undescribed species from Kazakhstan and A. urmiana. CONCLUSIONS/SIGNIFICANCE: We found evidence for multiple origin of parthenogenesis in Central Asia. Our results indicated that, shortly after its origin, diploid A. parthenogenetica populations underwent a rapid range expansion from Central Asia towards the Mediterranean region, and probably to the rest of its current geographic distribution. This contrasts with the restricted geographic distribution, strong genetic structure, and regional endemism of sexual Artemia lineages and other passively dispersed sexual continental aquatic invertebrates. We hypothesize that diploid parthenogens might have reached their current distribution in historical times, with a range expansion possibly facilitated by an increased availability of suitable habitat provided by anthropogenic activities, such as the spread of solar saltworks, aided by their natural dispersal vectors (i.e., waterbirds).

Phylogeography and local endemism of the native Mediterranean brine shrimp Artemia salina (Branchiopoda: Anostraca).

There has been a recent appreciation of the ecological impacts of zooplanktonic species invasions. The North American brine shrimp Artemia franciscana is one such alien invader in hyper-saline water ecosystems at a global scale. It has been shown to outcompete native Artemia species, leading to their local extinction. We used partial sequences of the mitochondrial Cytochrome c Oxidase Subunit 1 (COI or cox1) gene to investigate the genetic diversity and phylogeography of A. salina, an extreme halophilic sexual brine shrimp, over its known distribution range (Mediterranean Basin and South Africa) and to assess the extent of local endemism, the degree of population structure and the potential impact of traditional human saltpan management on this species. We also examined the phylogenetic relationships in the genus Artemia using COI sequences. Our results show extensive regional endemism and indicate an early Pleistocene expansion of A. salina in the Mediterranean Basin. Subsequent population isolation in a mosaic of Pleistocene refugia is suggested, with two or three refugia located in the Iberian Peninsula. Two instances of long-distance colonization were also observed. Surprisingly, given its strong phylogeographical structure, A. salina showed a signature of correlation between geographical and genetic distance. Owing to strong 'priority effects', extensive population differentiation is retained, despite dispersal via migrant birds and human management of saltpans. The foreseeable expansion of A. franciscana is likely to be followed by substantial loss of genetic diversity in Mediterranean A. salina. Large genetic divergences between Mediterranean and South African A. salina suggest that the latter deserves species status.

Global isolation by distance despite strong regional phylogeography in a small metazoan.

BACKGROUND: Small vagile eukaryotic organisms, which comprise a large proportion of the Earth's biodiversity, have traditionally been thought to lack the extent of population structuring and geographic speciation observed in larger taxa. Here we investigate the patterns of genetic diversity, amongst populations of the salt lake microscopic metazoan Brachionus plicatilis s. s. (sensu stricto) (Rotifera: Monogononta) on a global scale. We examine the phylogenetic relationships of geographic isolates from four continents using a 603 bp fragment of the mitochondrial COI gene to investigate patterns of phylogeographic subdivision in this species. In addition we investigate the relationship between genetic and geographic distances on a global scale to try and reconcile the paradox between the high vagility of this species and the previously reported patterns of restricted gene flow, even over local spatial scales. RESULTS: Analysis of global sequence diversity of B. plicatilis s. s. reveals the presence of four allopatric genetic lineages: North American-Far East Asian, Western Mediterranean, Australian, and an Eastern Mediterranean lineage represented by a single isolate. Geographically orientated substructure is also apparent within the three best sampled lineages. Surprisingly, given this strong phylogeographic structure, B. plicatilis s. s. shows a significant correlation between geographic and genetic distance on a global scale ('isolation by distance' - IBD). CONCLUSION: Despite its cosmopolitan distribution and potential for high gene flow, B. plicatilis s. s. is strongly structured at a global scale. IBD patterns have traditionally been interpreted to indicate migration-drift equilibrium, although in this system equilibrium conditions are incompatible with the observed genetic structure. Instead, we suggest the pattern may have arisen through persistent founder effects, acting in a similar fashion to geographic barriers for larger organisms. Our data indicates that geographic speciation, contrary to historical views, is likely to be very important in microorganisms. By presenting compelling evidence for geographic speciation in a small eukaryote we add to the growing body of evidence that is forcing us to rethink our views of global biodiversity.

Recent evolution of alternative reproductive modes in the 'living fossil' Triops cancriformis.

BACKGROUND: The Notostraca is a small but ancient crustacean order with a contrasting combination of a conservative morphology and a wide range of reproductive modes. The tadpole shrimp Triops cancriformis, includes bisexual - the putatively ancestral state -, androdioecious and hermaphrodite populations. As hermaphroditism and androdioecy confer a colonisation advantage, we expect the postglacial colonisation of northern Europe to have been effected by lineages with such reproductive modes. Therefore, N European populations should be composed of closely related lineages reflecting a recent range expansion. In contrast, glacial refugia in the south should contain bisexual populations with high haplotype diversity and more population structuring. To test these hypotheses, we analysed the geographic distribution of reproductive modes based on new and published sex ratio data. In addition, we investigated the European phylogeography of T. cancriformis by sequencing over a 1000 bp of mitochondrial DNA (mtDNA) in individuals from a large sample of populations of the three recognised subspecies. RESULTS: Bisexual populations were only found in the Iberian Peninsula, with the rest of European populations showing low male proportions or no males. Androdioecious populations were found in Central and Eastern Europe. Regarding mtDNA diversity, Spanish and Moroccan populations of T. c. mauritanicus were highly divergent, and showed strong population structure. In contrast, Triops c. cancriformis and T. c. simplex formed a single mtDNA lineage with low haplotype diversity. This diversity was structured into two phylogenetic clades (A, B), coexisting in E Germany. Basal haplotypes of both lineages were found in the Iberian Peninsula. Most of the populations in clade A and B are either hermaphroditic or androdioecious, with the only bisexual population in these clades found in the Iberian Peninsula. The genetic divergence between these two clades suggests a split in the Late Pleistocene and their geographic distribution reflects a complex evolutionary history of European Triops populations, with possibly two episodes of range expansions - one of them by clade A - involving androdioecious and hermaphroditic populations. CONCLUSION: As we predicted, N European populations of T. cancriformis are closely related, with few widely distributed haplotypes and indications of a recent range expansion involving hermaphroditic/androdioecious lineages. A possible second range expansion or long distance colonisation may have created the secondary contact zone between T. c. cancriformis/simplex clades A and B. The large haplotype diversity and strong genetic subdivision in the Iberian Peninsula, which is known to contain only bisexual populations, strongly suggest that this area was a Pleistocene refugium for T. cancriformis, although the occurrence of additional eastern refugia cannot be ruled out. Our data support the status of T. c. mauritanicus as a separate species and the colonisation of N Africa from the Iberian Peninsula. We suggest that hermaphroditism/androdioecy has evolved recently in T. cancriformis and has facilitated the postglacial colonisation of northern Europe.

Persistent genetic signatures of colonization in Brachionus manjavacas rotifers in the Iberian Peninsula.

Recent phylogeographical assessments have consistently shown that continental zooplankton display high levels of population subdivision, despite the high dispersal capacity of their diapausing propagules. As such, there is an apparent paradox between observed cosmopolitanism in the zooplankton that is associated with long-distance dispersal, and strong phylogeographical structures at a regional scale. Such population dynamics, far from migration-drift equilibrium, have been shown in the rotifer species complex Brachionus plicatilis, a group of over a dozen species inhabiting salt lakes and coastal lagoons worldwide. Here we present the mitochondrial DNA phylogeography of one of these species, Brachionus manjavacas, in the Iberian Peninsula, where it often co-occurs with the morphologically similar species B. plicatilis sensu stricto. We obtained sequences from 233 individuals from diapausing eggs and clonal cultures from 16 lakes in the Iberian Peninsula, and a Tunisian lake. Two strongly supported deep mitochondrial DNA clades were found (A and B). Phylogenetic and nested clade analysis showed that clade A has a strong phylogeographical structure, with a strong similarity of phylogeographical patterns between B. manjavacas clade A and B. plicatilis s.s. These include (i) signatures of allopatric fragmentation between central and southern populations, and (ii) range expansions in the Iberian Peninsula, both likely to have occurred during the Pleistocene. We find evidence for a glacial refugium in the Guadiana basin. Clades A and B co-occurred in several of these lakes because of range expansion and secondary contact between both clades. The co-occurrence between B. plicatilis s.s. and B. manjavacas is not recent, and both species might have experienced similar environmental challenges during the Pleistocene. The strong correlation of genetic and geographical distance found suggests that historical events can lead to such correlation, mirroring the effects of 'isolation by distance' in equilibrium populations.

Mitochondrial DNA phylogeography and mating compatibility reveal marked genetic structuring and speciation in the NE Atlantic bryozoan Celleporella hyalina.

The marine bryozoan Celleporella hyalina is a species complex composed of many highly divergent and mostly allopatric genetic lineages that are reproductively isolated but share a remarkably similar morphology. One such lineage commonly encrusts macroalgae throughout the NE Atlantic coast. To explore the processes leading to geographical diversification, reproductive isolation and speciation in this taxon, we (i) investigated NE Atlantic C. hyalina mitochondrial DNA phylogeography, and (ii) used breeding trials between geographical isolates to ascertain reproductive isolation. We find that haplotype diversity is geographically variable and there is a strong population structure, with significant isolation by distance. NE Atlantic C. hyalina is structured into two main parapatric lineages that appear to have had independent Pleistocene histories. Range expansions have resulted in two contact zones in Spain and W Ireland. Lineage 1 is found from Ireland to Spain and has low haplotype diversity, with closely related haplotypes, suggesting a recent population expansion into the Irish Sea, S Ireland, S England and Spain. Lineage 2 is found from Iceland to Spain and has high haplotype diversity. Complete reproductive isolation was found between some geographical isolates representing both lineages, whereas it was incomplete or asymmetric between others, suggesting these latter phylogeographical groups probably represent incipient species. The phylogeographical distribution of NE Atlantic C. hyalina does not fall easily into a pattern of southern refugia, and we discuss likely differences between terrestrial and marine system responses to Pleistocene glacial cycles.