Babies and bathwater: a comment on the premature obituary for nested clade phylogeographical analysis

Genetic and genomic monitoring with minimally invasive sampling methods

The decreasing cost and increasing scope and power of emerging genomic technologies are reshaping the field of molecular ecology. However, many modern genomic approaches (e.g., RAD-seq) require large amounts of high-quality template DNA. This poses a problem for an active branch of conservation biology: genetic monitoring using minimally invasive sampling (MIS) methods. Without handling or even observing an animal, MIS methods (e.g., collection of hair, skin, faeces) can provide genetic information on individuals or populations. Such samples typically yield low-quality and/or quantities of DNA, restricting the type of molecular methods that can be used. Despite this limitation, genetic monitoring using MIS is an effective tool for estimating population demographic parameters and monitoring genetic diversity in natural populations. Genetic monitoring is likely to become more important in the future as many natural populations are undergoing anthropogenically driven declines, which are unlikely to abate without intensive adaptive management efforts that often include MIS approaches. Here, we profile the expanding suite of genomic methods and platforms compatible with producing genotypes from MIS, considering factors such as development costs and error rates. We evaluate how powerful new approaches will enhance our ability to investigate questions typically answered using genetic monitoring, such as estimating abundance, genetic structure and relatedness. As the field is in a period of unusually rapid transition, we also highlight the importance of legacy data sets and recommend how to address the challenges of moving between traditional and next-generation genetic monitoring platforms. Finally, we consider how genetic monitoring could move beyond genotypes in the future. For example, assessing microbiomes or epigenetic markers could provide a greater understanding of the relationship between individuals and their environment.

Evaluating sample size to estimate genetic management metrics in the genomics era

Inbreeding and relationship metrics among and within populations are useful measures for genetic management of wild populations, but accuracy and precision of estimates can be influenced by the number of individual genotypes analysed. Biologists are confronted with varied advice regarding the sample size necessary for reliable estimates when using genomic tools. We developed a simulation framework to identify the optimal sample size for three widely used metrics to enable quantification of expected variance and relative bias of estimates and a comparison of results among populations. We applied this approach to analyse empirical genomic data for 30 individuals from each of four different free-ranging Rocky Mountain bighorn sheep (Ovis canadensis canadensis) populations in Montana and Wyoming, USA, through cross-species application of an Ovine array and analysis of approximately 14,000 single nucleotide polymorphisms (SNPs) after filtering. We examined intra- and interpopulation relationships using kinship and identity by state metrics, as well as FST between populations. By evaluating our simulation results, we concluded that a sample size of 25 was adequate for assessing these metrics using the Ovine array to genotype Rocky Mountain bighorn sheep herds. However, we conclude that a universal sample size rule may not be able to sufficiently address the complexities that impact genomic kinship and inbreeding estimates. Thus, we recommend that a pilot study and sample size simulation using R code we developed that includes empirical genotypes from a subset of populations of interest would be an effective approach to ensure rigour in estimating genomic kinship and population differentiation.

Quaternary range and demographic expansion of Liolaemus darwinii (Squamata: Liolaemidae) in the Monte Desert of Central Argentina using Bayesian phylogeography and ecological niche modelling

Until recently, most phylogeographic approaches have been unable to distinguish between demographic and range expansion processes, making it difficult to test for the possibility of range expansion without population growth and vice versa. In this study, we applied a Bayesian phylogeographic approach to reconstruct both demographic and range expansion in the lizard Liolaemus darwinii of the Monte Desert in Central Argentina, during the Late Quaternary. Based on analysis of 14 anonymous nuclear loci and the cytochrome b mitochondrial DNA gene, we detected signals of demographic expansion starting at ~55 ka based on Bayesian Skyline and Skyride Plots. In contrast, Bayesian relaxed models of spatial diffusion suggested that range expansion occurred only between ~95 and 55 ka, and more recently, diffusion rates were very low during demographic expansion. The possibility of population growth without substantial range expansion could account for the shared patterns of demographic expansion during the Last Glacial Maxima (OIS 2 and 4) in fish, small mammals and other lizards of the Monte Desert. We found substantial variation in diffusion rates over time, and very high rates during the range expansion phase, consistent with a rapidly advancing expansion front towards the southeast shown by palaeo-distribution models. Furthermore, the estimated diffusion rates are congruent with observed dispersal rates of lizards in field conditions and therefore provide additional confidence to the temporal scale of inferred phylogeographic patterns. Our study highlights how the integration of phylogeography with palaeo-distribution models can shed light on both demographic and range expansion processes and their potential causes.

Conflicting evolutionary patterns due to mitochondrial introgression and multilocus phylogeography of the Patagonian freshwater crab Aegla neuquensis

BACKGROUND

Multiple loci and population genetic methods were employed to study the phylogeographic history of the Patagonian freshwater crab Aegla neuquensis (Aeglidae: Decopoda). This taxon occurs in two large river systems in the Patagonian Steppe, from the foothills of the Andes Mountains east to the Atlantic Ocean.

METHODOLOGY/PRINCIPAL FINDINGS

A nuclear phylogeny and multilocus nested clade phylogeographic analysis detected a fragmentation event between the Negro and Chico-Chubut river systems. This event occurred approximately 137 thousand years ago. An isolation-with-migration analysis and maximum-likelihood estimates of gene flow showed asymmetrical exchange of genetic material between these two river systems exclusively in their headwaters. We used information theory to determine the best-fit demographic history between these two river systems under an isolation-with-migration model. The best-fit model suggests that the Negro and the ancestral populations have the same effective population sizes; whereas the Chico-Chubut population is smaller and shows that gene flow from the Chico-Chubut into the Negro is four times higher than in the reverse direction. Much of the Chico-Chubut system appears to have only been recently colonized while the Negro populations appear to have been in place for most of the evolutionary history of this taxon.

CONCLUSIONS/SIGNIFICANCE

Due to mitochondrial introgression, three nuclear loci provided different phylogeographic resolution than the three mitochondrial genes for an ancient fragmentation event observed in the nuclear phylogeny. However, the mitochondrial locus provided greater resolution on more recent evolutionary events. Our study, therefore, demonstrates the need to include both nuclear and mitochondrial loci for a more complete understanding of evolutionary histories and associated phylogeographic events. Our results suggest that gene flow between these systems, before and after fragmentation was through periodic paleolakes that formed in the headwaters region. Fragmentation between the Negro and Chico-Chubut systems was driven by the disappearance of these paleolakes during the Patagonian Glaciation.

Divergence and dispersal of the nematode-trapping fungus Arthrobotrys oligospora from China

Nematode-trapping fungi are of significant agricultural, forestry and ecological importance. However, relatively little is known about the patterns of genetic variation for any nematode-trapping fungus through its broad geographic and ecological contexts. Here, we analysed DNA sequence variation among strains of the nematode-trapping fungus Arthrobotrys oligospora from China. Our analyses revealed that the Chinese A. oligospora is a species complex with at least three divergent lineages (cryptic species). In addition, there was significant geographic structuring with unambiguous evidence for localized recombination within two of the three lineages in nature. However, evidence for clonal reproduction was also found. We discuss the implications of our results to the conservation and biocontrol application of A. oligospora in agriculture and forestry.

The limits of cryptic diversity in groundwater: phylogeography of the cave shrimp Troglocaris anophthalmus (Crustacea: Decapoda: Atyidae)

Recent studies have revealed high local diversity and endemism in groundwaters, and showed that species with large ranges are extremely rare. One of such species is the cave shrimp Troglocaris anophthalmus from the Dinaric Karst on the western Balkan Peninsula, apparently uniform across a range of more than 500 kilometers. As such it contradicts the paradigm that subterranean organisms form localized, long-term stable populations that cannot disperse over long distances. We tested it for possible cryptic diversity and/or unexpected evolutionary processes, analysing mitochondrial (COI, 16S rRNA) and nuclear (ITS2) genes of 232 specimens from the entire range. The results of an array of phylogeographical procedures congruently suggested that the picture of a widespread, continuously distributed and homogenous T. anophthalmus was wrong. The taxon is composed of four or possibly five monophyletic, geographically defined phylogroups that meet several species delimitation criteria, two of them showing evidence of biological reproductive isolation in sympatry. COI genetic distances between phylogroups turned out to be a poor predictor, as they were much lower than the sometimes suggested crustacean threshold value of 0.16 substitutions per site. Most results confirmed the nondispersal hypothesis of subterranean fauna, but the southern Adriatic phylogroup displayed a paradoxical pattern of recent dispersal across 300 kilometers of hydrographically fragmented karst terrain. We suggest a model of migration under extreme water-level conditions, when flooded poljes could act as stepping-stones. In the north of the range (Slovenia), the results confirmed the existence of a zone of unique biogeographical conflict, where surface fauna is concordant with the current watershed, and subterranean fauna is not.

DIM SUM: demography and individual migration simulated using a Markov chain

An increasing number of studies seek to infer demographic history, often jointly with genetic relationships. Despite numerous analytical methods for such data, few simulations have investigated the methods' power and robustness, especially when underlying assumptions have been violated. DIM SUM (Demography and Individual Migration Simulated Using a Markov chain) is a stand-alone Java program for the simulation of population demography and individual migration while recording ancestor-descendant relationships. It does not employ coalescent assumptions or discrete population boundaries. It is extremely flexible, allowing the user to specify border positions, reactions of organisms to borders, local and global carrying capacities, individual dispersal kernels, rates of reproduction and strategies for sampling individuals. Spatial variables may be specified using image files (e.g., as exported from gis software) and may vary through time. In combination with software for genetic marker simulation, DIM SUM will be useful for testing phylogeographic (e.g., nested clade phylogeographic analysis, coalescent-based tests and continuous-landscape frameworks) and landscape-genetic methods, specifically regarding violations of coalescent assumptions. It can also be used to explore the qualitative features of proposed demographic scenarios (e.g. regarding biological invasions) and as a pedagogical tool. DIM SUM (with user's manual) can be downloaded from http://code.google.com/p/bio-dimsum.

Ginkgo: spatially-explicit simulator of complex phylogeographic histories

We present Ginkgo, a software package for agent-based, forward-time simulations of genealogies of multiple unlinked loci from diploid populations. Ginkgo simulates the evolution of one or more species on a spatially explicit landscape of cells. The user of the software can specify the geographical and environmental characteristics of the landscape, and these properties can change according to a prespecified schedule. The geographical elements modelled include the arrangement of cells and movement rates between particular cells. Each species has a function that can calculate a fitness score for any combination of an individual organism's phenotype and environmental characteristics. The user can control the number of fitness factors (the dimensionality of the cell-specific fitness factors and the individuals phenotypic vectors) and the weighting of each of these dimensions in the fitness calculation. Cell-specific fitness trait optima can be specified across the landscape to mimic differences in habitat. In addition to their differing fitness functions, species can differ in terms of their vagility and fecundity. Genealogies and occurrence data can be produced at any time during the simulation in NEXUS and ESRI Ascii Grid formats, respectively.

Three roads diverged? Routes to phylogeographic inference

Phylogeographic methods facilitate inference of the geographical history of genetic lineages. Recent examples explore human migration and the origins of viral pandemics. There is longstanding disagreement over the use and validity of certain phylogeographic inference methodologies. In this paper, we highlight three distinct frameworks for phylogeographic inference to give a taste of this disagreement. Each of the three approaches presents a different viewpoint on phylogeography, most fundamentally on how we view the relationship between the inferred history of a sample and the history of the population the sample is embedded in. Satisfactory resolution of this relationship between history of the tree and history of the population remains a challenge for all but the most trivial models of phylogeographic processes. Intriguingly, we believe that some recent methods that entirely avoid inference about the history of the population will eventually help to reach a resolution.

Identifying shared genetic structure patterns among Pacific Northwest forest taxa: insights from use of visualization tools and computer simulations

BACKGROUND

Identifying causal relationships in phylogeographic and landscape genetic investigations is notoriously difficult, but can be facilitated by use of multispecies comparisons.

METHODOLOGY/PRINCIPAL FINDINGS

We used data visualizations to identify common spatial patterns within single lineages of four taxa inhabiting Pacific Northwest forests (northern spotted owl: Strix occidentalis caurina; red tree vole: Arborimus longicaudus; southern torrent salamander: Rhyacotriton variegatus; and western white pine: Pinus monticola). Visualizations suggested that, despite occupying the same geographical region and habitats, species responded differently to prevailing historical processes. S. o. caurina and P. monticola demonstrated directional patterns of spatial genetic structure where genetic distances and diversity were greater in southern versus northern locales. A. longicaudus and R. variegatus displayed opposite patterns where genetic distances were greater in northern versus southern regions. Statistical analyses of directional patterns subsequently confirmed observations from visualizations. Based upon regional climatological history, we hypothesized that observed latitudinal patterns may have been produced by range expansions. Subsequent computer simulations confirmed that directional patterns can be produced by expansion events.

CONCLUSIONS/SIGNIFICANCE

We discuss phylogeographic hypotheses regarding historical processes that may have produced observed patterns. Inferential methods used here may become increasingly powerful as detailed simulations of organisms and historical scenarios become plausible. We further suggest that inter-specific comparisons of historical patterns take place prior to drawing conclusions regarding effects of current anthropogenic change within landscapes.

Phylogeography of two moray eels indicates high dispersal throughout the indo-pacific

Reef fishes disperse primarily as oceanic "pelagic" larvae, and debate continues over the extent of this dispersal, with recent evidence for geographically restricted (closed) populations in some species. In contrast, moray eels have the longest pelagic larval stages among reef fishes, possibly providing opportunities to disperse over great distances. We test this prediction by measuring mitochondrial DNA (mtDNA) and nuclear DNA variation in 2 species of moray eels, Gymnothorax undulatus (N = 165) and G. flavimarginatus (N = 124), sampled at 14-15 locations across the Indo-Pacific. The mtDNA data comprise 632 bp of cytochrome b and 596 bp of cytochrome oxidase I. Nuclear markers include 2 recombination-activating loci (421 bp of RAG-1 and 754 bp of RAG-2). Analyses of molecular variance and Mantel tests indicate little or no genetic differentiation, and no isolation by distance, across 22 000 km of the Indo-Pacific. We estimate that mitochondrial genetic variation coalesces within the past about 2.3 million years (My) for G. flavimarginatus and within the past about 5.9 My for G. undulatus. Permutation tests of geographic distance on the mitochondrial haplotype networks indicate recent range expansions for some younger haplotypes (estimated within approximately 600 000 years) and episodic fragmentation of populations at times of low sea level. Our results support the predictions that the extended larval durations of moray eels enable ocean-wide genetic continuity of populations. This is the first phylogeographic survey of the moray eels, and morays are the first reef fishes known to be genetically homogeneous across the entire Indo-Pacific.

Inference of population history by coupling exploratory and model-driven phylogeographic analyses

Understanding the nature, timing and geographic context of historical events and population processes that shaped the spatial distribution of genetic diversity is critical for addressing questions relating to speciation, selection, and applied conservation management. Cladistic analysis of gene trees has been central to phylogeography, but when coupled with approaches that make use of different components of the information carried by DNA sequences and their frequencies, the strength and resolution of these inferences can be improved. However, assessing concordance of inferences drawn using different analytical methods or genetic datasets, and integrating their outcomes, can be challenging. Here we overview the strengths and limitations of different types of genetic data, analysis methods, and approaches to historical inference. We then turn our attention to the potentially synergistic interactions among widely-used and emerging phylogeographic analyses, and discuss some of the ways that spatial and temporal concordance among inferences can be assessed. We close this review with a brief summary and outlook on future research directions.

Matrilinear phylogeography and demographical patterns of Rutilus rutilus: implications for taxonomy and conservation

A phylogeographic analysis of mitochondrial DNA sequence variation was carried out to infer the geographical distribution of the genealogical lineages and the historical demography of roach Rutilus rutilus (L.). A total of 265 individuals from 52 sites covering most of the Eurasian distribution range were sequenced for a 475 bp fragment of the mitochondrial cytochrome b gene. The monophyletic roach contained two deep clades that dated back to the Pliocene. The Ponto-Caspian clade comprised populations from Greece to Siberia with a likely palaeorefugium at the west coast of the Caspian Sea. This clade largely corresponds to individuals with morphological features described as Rutilus heckelii. The west European clade included individuals from central and western Europe with the Danube and Dniester basins as possible palaeorefugia. This clade largely corresponds to individuals with morphological features described as R. rutilus. A suture-zone of the two main lineages was observed along the coastal region of the Black Sea. The neutrality tests and the mismatch distributions indicated a demographic expansion during the Middle-Pleistocene for both clades.

Plastid DNA variation in highly fragmented populations of Microbiota decussata Kom. (Cupressaceae), an endemic to Sikhote Alin Mountains

Microbiota decussata Kom. (Cupressaceae) is a subalpine species endemic to the Sikhote Alin Mountains with populations scattered throughout the range. We used sequence data for four noncoding regions of chloroplast DNA to characterize the genetic diversity in populations sampled from different parts of M. decussata natural range. No variation was observed in the trnT-trnF region, whereas the trnH-psbA, trnS-trnfM, and trnS-trnG regions showed polymorphisms. At the species level, we found a low nucleotide diversity (pi = 0.0009) and high haplotype diversity (h = 0.981) as well as high differentiation (Phi(ST) = 0.420). N(ST) and G(ST) values suggested the existence of a phylogeographic structure in M. decussata. The observed patterns of diversity could be explained in part by ecological features of the species and its long-term persistence throughout the range with population expansion, successive fragmentation and isolation.

The evolution of sexual and parthenogenetic Warramaba: a window onto Plio-Pleistocene diversification processes in an arid biome

Environmental changes over the Plio-Pleistocene have been key drivers of speciation patterns and genetic diversification in high-latitude and mesic environments, yet comparatively little is known about the evolutionary history of species in arid environments. We applied phylogenetic and phylogeographic analyses to understand the evolutionary history of Warramaba grasshoppers from the Australian arid zone, a group including sexual and parthenogenetic lineages. Sequence data (mitochondrial COI) showed that the four major sexual lineages within Warramaba most likely diverged in the Pliocene, around 2-7 million years ago. All sexual lineages exhibited considerable phylogenetic structure. Detailed analyses of the hybrid parthenogenetic species W. virgo and its sexual progenitors showed a pattern of high phylogenetic diversity and phylogeographic structure in northern lineages, and low diversity and evidence for recent expansion in southern lineages. Northern sexual lineages persisted in localized refugia over the Pleistocene, with sustained barriers promoting divergence over this period. Southern parts of the present range became periodically unsuitable during the Pleistocene, and it is into this region that parthenogenetic lineages have expanded. Our results strongly parallel those for sexual and parthenogenetic lineages of the gecko Heteronotia from the same region, indicating a highly general effect of Plio-Pleistocene environmental change on diversification processes in arid Australia.

Faecal DNA analysis enables genetic monitoring of the species recovery program for an arid-dwelling marsupial

The greater bilby, Macrotis lagotis, is a species of conservation significance in the arid and semiarid zones of Australia. A species recovery program has been underway since the mid-1990s but the incorporation of molecular genetic data within the program has been difficult due to the problems of obtaining regular, population-wide samples of this trap-shy and sparsely distributed species. In this study, we demonstrate that faecal pellets collected from around burrows in the dry, arid habitat of western Queensland provide a viable source for DNA extraction and analysis. Faecal DNA was used to generate population-level estimates of microsatellite and mtDNA diversity for comparison with previous estimates for the natural population derived from tissue samples. Data were used to assess both the reliability of faecal-derived genotypes and the extent of any diversity loss since the previous study. Microsatellite diversity recorded from eight polymorphic markers for the natural population (A = 4.31 ± 0.30, HE = 0.76 ± 0.03) was comparable with the previous study, indicating little change in genetic diversity for the natural population in the 10-year interim. Faecal genotypes generated for the recently reintroduced population matched the known number of founders as well as a known genotype, providing support for the reliability of the faecal DNA approach. The captive and reintroduced populations had significantly lower diversity levels than the natural population (A = 3.59 ± 0.28, HE = 0.68 ± 0.03; A = 3.57 ± 0.20, HE = 0.65 ± 0.03 respectively). Mitochondrial control region analysis, incorporating nested clade phylogeographic analysis (NCPA), agrees with earlier findings that populations of bilbies across the arid zone in Australia have only recently become fragmented, but the case for Queensland bilbies being strongly differentiated from other regions is diminished. Implications from this study include the need to further supplement the captive and reintroduced populations with additional out-bred individuals and that faecal DNA can be used effectively for ongoing monitoring and management of this species.

Twenty years of phylogeography: the state of the field and the challenges for the Southern Hemisphere

Phylogeography is a young, vigorous and integrative field of study that uses genetic data to understand the history of populations. This field has recently expanded into many areas of biology and also into several historical disciplines of Earth sciences. In this review, I present a numerical synthesis of the phylogeography literature based on an examination of over 3000 articles published during the first 20 years of the field (i.e. from 1987 to 2006). Information from several topics needed to evaluate the progress, tendencies and deficiencies of the field is summarized for 10 major groups of organisms and at a global scale. The topics include the geography of phylogeographic surveys, comparative nature of studies, temporal scales and major environments investigated, and genetic markers used. I also identify disparities in research productivity between the developing and the developed world, and propose ways to reduce some of the challenges faced by phylogeographers from less affluent countries. Phylogeography has experienced explosive growth in recent years fuelled by developments in DNA technology, theory and statistical analysis. I argue that the intellectual maturation of the field will eventually depend not only on these recent developments, but also on syntheses of comparative information across different regions of the globe. For this to become a reality, many empirical phylogeographic surveys in regions of the Southern Hemisphere (and in developing countries of the Northern Hemisphere) are needed. I expect the information and views presented here will assist in promoting international collaborative work in phylogeography and in guiding research efforts at both regional and global levels.

Complex population genetic and demographic history of the Salangid, Neosalanx taihuensis, based on cytochrome b sequences

BACKGROUND

The Salangid icefish Neosalanx taihuensis (Salangidae) is an economically important fish, which is endemic to China, restricted to large freshwater systems (e.g. lakes, large rivers and estuaries) and typically exhibit low vagility. The continuous distribution ranges from the temperate region of the Huai and Yellow River basins to the subtropical region of the Pearl River basin. This wide ranging distribution makes the species an ideal model for the study of palaeoclimatic effects on population genetic structure and phylogeography. Here, we aim to analyze population genetic differentiation within and between river basins and demographic history in order to understand how this species responded to severe climatic oscillations, decline of the sea levels during the Pleistocene ice ages and tectonic activity.

RESULTS

We obtained the complete mtDNA cytochrome b sequences (1141 bp) of 354 individuals from 13 populations in the Pearl River, the Yangze River and the Huai River basin. Thirty-six haplotypes were detected. Haplotype frequency distributions were strongly skewed, with most haplotypes (n = 24) represented only in single samples each and thus restricted to a single population. The most common haplotype (H36) was found in 49.15% of all individuals. Analysis of molecular variance (AMOVA) revealed a random pattern in the distribution of genetic diversity, which is inconsistent with contemporary hydrological structure. Significant levels of genetic subdivision were detected among populations within basins rather than between the three basins. Demographic analysis revealed that the population size in the Pearl River basin has remained relatively constant whereas the populations in the Yangze River and the Huai River basins expanded about 221 and 190 kyr ago, respectively, with the majority of mutations occurring after the last glacial maximum (LGM).

CONCLUSION

The observed complex genetic pattern of N. taihuensis is coherent with a scenario of multiple unrelated founding events by long-distance colonization and dispersal combined with contiguous population expansion and locally restricted gene flow. We also found that this species was likely severely impacted by past glaciations. More favourable climate and the formation of large suitable habitations together facilitated population expansion after the late Quaternary (especially the LGM). We proposed that all populations should be managed and conserved separately, especially for habitat protection.

Origin and post-colonization evolution of the Chatham Islands skink (Oligosoma nigriplantare nigriplantare)

Island ecosystems provide an opportunity to examine a range of evolutionary and ecological processes. The Chatham Islands are an isolated archipelago situated approximately 800 km east of New Zealand. Geological evidence indicates that the Chatham Islands re-emerged within the last 1-4 million years, following a prolonged period of marine inundation, and therefore the resident flora and fauna is the result of long-distance overwater dispersal. We examine the origin and post-colonization evolution of the Chatham Islands skink, Oligosoma nigriplantare nigriplantare, the sole reptile species occurring on the archipelago. We sampled O. n. nigriplantare from across nine islands within the Chatham Islands group, and representative samples from across the range of its closest relative, the New Zealand mainland common skink (Oligosoma nigriplantare polychroma). Our mitochondrial sequence data indicate that O. n. nigriplantare diverged from O. n. polychroma 5.86-7.29 million years ago. This pre-dates the emergence date for the Chatham Islands, but indicates that O. n. nigriplantare colonized the Chatham Islands via overwater dispersal on a single occasion. Despite the substantial morphological variability evident in O. n. nigriplantare, only relatively shallow genetic divergences (maximum divergence approximately 2%) were found across the Chatham Islands. Our analyses (haplotypic diversity, Phi(ST), analysis of molecular variance, and nested clade phylogeographical analysis) indicated restricted gene flow in O. n. nigriplantare resulting in strong differentiation between islands. However, the restrictions to gene flow might have only arisen recently as there was also a significant pattern of isolation by distance, possibly from when the Chatham Islands were a single landmass during Pleistocene glacial maxima when sea levels were lower. The level of genetic and morphological divergence between O. n. nigriplantare and O. n. polychroma might warrant their recognition as distinct species.

Patterns of diversification in the discus fishes (Symphysodon spp. Cichlidae) of the Amazon basin

We carried out a phylogeograhic and population genetic analysis of fishes of the taxonomically contentious genus Symphysodon from the Amazon basin in order to test hypotheses of relationships among taxonomic units, and potential processes driving diversification within this genus. We sampled 334 individuals of the genus Symphysodon from 24 localities that span the complete geographic distribution of this genus. The sampling scheme included all known phenotypic groups, species and subspecies. Analyses were based on 474 bp of the mitochondrial control region and 1443 bp of the exon 3 of RAG1 gene. We observed 102 mtDNA haplotypes defined by 89 segregating sites, and 5 nuDNA alleles defined by three segregating sites. Maximum-likelihood, Bayesian-inference and statistical parsimony analyses revealed three well defined monophyletic groups. These clades corresponded to the 'green' and 'blue' groups of Symphysodon aequifasciatus, and to a previously morphologically unrecognized clade from the Xingu River drainage. These three clades were nested within a paraphyletic assemblage consisting of the 'brown' group of S. aequifasciatus and of both described subspecies of S. discus, the 'Heckel' and the 'abacaxi' discus. Nuclear allele sharing was observed among groups, but there were significant differences in frequencies. We inferred several processes including past fragmentation among groups, and restricted gene flow with isolation by distance within the paraphyletic 'brown+Heckel+abacaxi' groups, and suggest that differences among the 'blue', 'Heckel' and 'brown' groups are potentially maintained by differences in water chemistry preferences. We further inferred colonization of the western Amazon basin by an ancestor of the 'green' clade. The 'green' group was the only group with a pattern of haplotype distribution consistent of a demographic expansion, and the divergence of this clade from other groups of discus was consistent with recent geologic evidence on the breach of the Purus Arch which separates western Amazon from eastern Amazon. We further hypothesized that the differentiation of the 'Xingu' clade could be due to vicariance events resulting from Pleistocene sea level, and thus Amazon River level fluctuations. We discuss the bearings of our results on the current taxonomy of this group, and on the biological reality of the different forms, subspecies and species of Symphysodon concluding that we are probably observing a process of diversification, and therefore taxonomy will remain contentious.

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.