The narrow endemic Norwegian peat moss Sphagnum troendelagicum originated before the last glacial maximum

It is commonly found that individual hybrid, polyploid species originate recurrently and that many polyploid species originated relatively recently. It has been previously hypothesized that the extremely rare allopolyploid peat moss Sphagnum troendelagicum has originated multiple times, possibly after the last glacial maximum in Scandinavia. This conclusion was based on low linkage disequilibrium in anonymous genetic markers within natural populations, in which sexual reproduction has never been observed. Here we employ microsatellite markers and chloroplast DNA (cpDNA)-encoded trnG sequence data to test hypotheses concerning the origin and evolution of this species. We find that S. tenellum is the maternal progenitor and S. balticum is the paternal progenitor of S. troendelagicum. Using various Bayesian approaches, we estimate that S. troendelagicum originated before the Holocene but not before c. 80 000 years ago (median expected time since speciation 40 000 years before present). The observed lack of complete linkage disequilibrium in the genome of this species suggests cryptic sexual reproduction and recombination. Several lines of evidence suggest multiple origins for S. troendelagicum, but a single origin is supported by approximate Bayesian computation analyses. We hypothesize that S. troendelagicum originated in a peat-dominated refugium before last glacial maximum, and subsequently immigrated to central Norway by means of spore flow during the last thousands of years.     

Phylogeography and population genetic structure of the European roe deer in Switzerland following recent recolonization

In the early 1800s, the European roe deer (Capreolus capreolus) was probably extirpated from Switzerland, due to overhunting and deforestation. After a federal law was enacted in 1875 to protect lactating females and young, and limiting the hunting season, the roe deer successfully recovered and recolonized Switzerland. In this study, we use mitochondrial DNA and nuclear DNA markers to investigate the recolonization and assess contemporary genetic structure in relation to broad topographic features, in order to understand underlying ecological processes, inform future roe deer management strategies, and explore the opportunity for development of forensic traceability tools. The results concerning the recolonization origin support natural, multidirectional immigration from neighboring countries. We further demonstrate that there is evidence of weak genetic differentiation within Switzerland among topographic regions. Finally, we conclude that the genetic data support the recognition of a single roe deer management unit within Switzerland, within which there is a potential for broad‐scale geographic origin assignment using nuclear markers to support law enforcement.     

Phylogeography: retrospect and prospect

Phylogeography has grown explosively in the two decades since the word was coined and the discipline was outlined in 1987. Here I summarize the many achievements and novel perspectives that phylogeography has brought to population genetics, phylogenetic biology and biogeography. I also address future directions for the field. From the introduction of mitochondrial DNA assays in the late 1970s, to the key distinction between gene trees and species phylogenies, to the ongoing era of multi-locus coalescent theory, phylogeographic perspectives have consistently challenged conventional genetic and evolutionary paradigms, and they have forged empirical and conceptual bridges between the formerly separate disciplines of population genetics (microevolutionary analysis) and phylogenetic biology (in macroevolution).     

Mitochondrial DNA diversity in North American and European Atlantic salmon with emphasis on the Downeast rivers of Maine

The displacement loop and NADH-1 dehydrogenase regions of mitochondrial DNA (mtDNA) were amplified by the polymerase chain reaction in 954 Atlantic salmon and digested with 40 restriction endonucleases. Variation was detected with 10 enzymes, resulting in 21 composite haplotypes which were strongly patterned geographically with a major discontinuity observed between most North American (NA) and European salmon. Significant heterogeneity of haplotype frequencies was found within and among all classification levels (continent, country, and river). Haplotype frequencies were significantly different across continents, within European samples, within NA samples, within Canadian samples, within wild Maine samples, within captive Maine strains, and between captive and wild Maine strains. Nine haplotypes occurred only in NA, seven in Maine, three only in Maine, and 11 occurred only in Europe. Some Maine rivers had only a single haplotype, suggesting that effective population sizes may be low. The second most frequent European haplotype occurred in tributaries to one Newfoundland river. Gene trees based on parsimony and genetic distance suggest that the haplotypes are monophyletic within each continent, and that the haplotype found on both continents is intermediate between those of Europe and NA, suggesting common ancestry of all haplotypes.     

Life‐history predicts past and present population connectivity in two sympatric sea stars

Life‐history traits, especially the mode and duration of larval development, are expected to strongly influence the population connectivity and phylogeography of marine species. Comparative analysis of sympatric, closely related species with differing life histories provides the opportunity to specifically investigate these mechanisms of evolution but have been equivocal in this regard. Here, we sample two sympatric sea stars across the same geographic range in temperate waters of Australia. Using a combination of mitochondrial DNA sequences, nuclear DNA sequences, and microsatellite genotypes, we show that the benthic‐developing sea star, Parvulastra exigua, has lower levels of within‐ and among‐population genetic diversity, more inferred genetic clusters, and higher levels of hierarchical and pairwise population structure than Meridiastra calcar, a species with planktonic development. While both species have populations that have diverged since the middle of the second glacial period of the Pleistocene, most P. exigua populations have origins after the last glacial maxima (LGM), whereas most M. calcar populations diverged long before the LGM. Our results indicate that phylogenetic patterns of these two species are consistent with predicted dispersal abilities; the benthic‐developing P. exigua shows a pattern of extirpation during the LGM with subsequent recolonization, whereas the planktonic‐developing M. calcar shows a pattern of persistence and isolation during the LGM with subsequent post‐Pleistocene introgression.     

Reproductive mode and genetic variation suggest a North American origin of European Letharia vulpina

Our data on the intercontinental population biology of Letharia vulpina show an unexpected shift from a recombining North American population with unique haplotypes to genetically depauperate Swedish and Italian populations, each with many representatives of a single repeated haplotype. Analysis of eight loci in 47 individuals supported recombination in North American populations and showed almost no variation among European populations. We infer that a genetic bottleneck caused by limited long-distance dispersal accounts for the lack of genetic variation found in marginal populations. This lack of variation in the European populations makes it impossible to use population genetic means to distinguish clonal reproduction from self-fertilization or even outcrossing, but phenotype indicates that reproduction in the marginal populations is by clonal spread, via soredia and isidioid soredia.     

In defence of model‐based inference in phylogeography

Recent papers have promoted the view that model‐based methods in general, and those based on Approximate Bayesian Computation (ABC) in particular, are flawed in a number of ways, and are therefore inappropriate for the analysis of phylogeographic data. These papers further argue that Nested Clade Phylogeographic Analysis (NCPA) offers the best approach in statistical phylogeography. In order to remove the confusion and misconceptions introduced by these papers, we justify and explain the reasoning behind model‐based inference. We argue that ABC is a statistically valid approach, alongside other computational statistical techniques that have been successfully used to infer parameters and compare models in population genetics. We also examine the NCPA method and highlight numerous deficiencies, either when used with single or multiple loci. We further show that the ages of clades are carelessly used to infer ages of demographic events, that these ages are estimated under a simple model of panmixia and population stationarity but are then used under different and unspecified models to test hypotheses, a usage the invalidates these testing procedures. We conclude by encouraging researchers to study and use model‐based inference in population genetics.     

Divergence,gene flow,and the origin of leapfrog geographic distributions: The history of colour pattern variation in Phyllobates poison‐dart frogs

The geographic distribution of phenotypic variation among closely related populations is a valuable source of information about the evolutionary processes that generate and maintain biodiversity. Leapfrog distributions, in which phenotypically similar populations are disjunctly distributed and separated by one or more phenotypically distinct populations, represent geographic replicates for the existence of a phenotype, and are therefore especially informative. These geographic patterns have mostly been studied from phylogenetic perspectives to understand how common ancestry and divergent evolution drive their formation. Other processes, such as gene flow between populations, have not received as much attention. Here, we investigate the roles of divergence and gene flow between populations in the origin and maintenance of a leapfrog distribution in Phyllobates poison frogs. We found evidence for high levels of gene flow between neighbouring populations but not over long distances, indicating that gene flow between populations exhibiting the central phenotype may have a homogenizing effect that maintains their similarity, and that introgression between ‘leapfroging’ taxa has not played a prominent role as a driver of phenotypic diversity in Phyllobates. Although phylogenetic analyses suggest that the leapfrog distribution was formed through independent evolution of the peripheral (i.e. leapfrogging) populations, the elevated levels of gene flow between geographically close populations poise alternative scenarios, such as the history of phenotypic change becoming decoupled from genome‐averaged patterns of divergence, which we cannot rule out. These results highlight the importance of incorporating gene flow between populations into the study of geographic variation in phenotypes, both as a driver of phenotypic diversity and as a confounding factor of phylogeographic inferences.     

Selection promotes organ compartmentalization in HIV-1: evidence from gag and pol genes

The existence of organ-specific human immunodeficiency virus type 1 (HIV-1) populations within infected hosts has been long lasting studied. Previous work established that population subdivision by organs occurs at the envelope env gene, but less is known about other genomic regions. Here, we used a population genetics approach to detect organ compartmentalization in proviral sequences of HIV-1 gag and pol genes. Significant population structure was found in pol (100% of cases) and gag (33%) pair-wise organ comparisons. The degree of compartmentalization positively correlated with the ratio of nonsynonymous to synonymous substitutions, and codons showing organ compartmentalization were more likely to be under significantly positive selection. This suggests that HIV-1 populations dynamically adapt to locally variable intra-host environments. In the case of pol gene, differential penetration of antiretroviral drugs might account for the observed pattern, whereas for gag gene, local selective pressures remain unexplored.     

Mapping the evolutionary twilight zone: molecular markers, populations and geography

Since evolutionary processes, such as dispersal, adaptation and drift, occur in a geographical context, at multiple hierarchical levels, biogeography provides a central and important unifying framework for understanding the patterns of distribution of life on Earth. However, the advent of molecular markers has allowed a clearer evaluation of the relationships between microevolutionary processes and patterns of genetic divergence among populations in geographical space, triggering the rapid development of many research programmes. Here we provide an overview of the interpretation of patterns of genetic diversity in geographical and ecological space, using both implicit and explicit spatial approaches. We discuss the actual or potential interaction of phylogeography, molecular ecology, ecological genetics, geographical genetics, landscape genetics and conservation genetics with biogeography, identifying their respective roles and their ability to deal with ecological and evolutionary processes at different levels of the biological hierarchy. We also discuss how each of these research programmes can improve strategies for biodiversity conservation. A unification of these research programmes is needed to better achieve their goals, and to do this it is important to develop cross‐disciplinary communication and collaborations among geneticists, ecologists, biogeographers and spatial statisticians.     

Phylogeographic breaks in low-dispersal species: the emergence of concordance across gene trees

Computer simulations were used to investigate population conditions under which phylogeographic breaks in gene genealogies can be interpreted with confidence to infer the existence and location of historical barriers to gene flow in continuously distributed, low-dispersal species. We generated collections of haplotypic gene trees under a variety of demographic scenarios and analyzed them with regard to salient genealogical breaks in their spatial patterns. In the first part of the analysis, we estimated the frequency in which the spatial location of the deepest phylogeographic break between successive pairs of populations along a linear habitat coincided with a spatial physical barrier to dispersal. Results confirm previous reports that individual gene trees can show ‘haphazard’ phylogeographic discontinuities even in the absence of historical barriers to gene flow. In the second part of the analysis, we assessed the probability that pairs of gene genealogies from a set of population samples agree upon the location of a geographical barrier. Our findings extend earlier reports by demonstrating that spatially concordant phylogeographic breaks across independent neutral loci normally emerge only in the presence of longstanding historical barriers to gene flow. Genealogical concordance across multiple loci thus becomes a deciding criterion by which to distinguish between stochastic and deterministic causation in accounting for phylogeographic discontinuities in continuously distributed species.     

Refugial isolation and secondary contact in the dyeing poison frog Dendrobates tinctorius

Recent palaeoclimactic research suggests that fluctuating environmental conditions throughout the Pleistocene of Amazonia occurred with previously unrecognized frequency. This has resulted in a theoretical shift from glacially influenced fluctuations to those driven by precessional rhythms. This theoretical revolution has a profound impact on expectations of biotic diversity within biogeographical regions that have long been based on the idea of large-scale landscape fragmentation associated with increased aridity and glacial cycles. Generally speaking, this shifts phylogeographical expectations from that of (i) large areas of sympatry of closely related (but not sister) species whose origins lie in separate refugia, and current distributions are the results of cyclic connectivity of those two refugia (refuge hypothesis), to that of (ii) fine scale genetic structure, often associated with elevation, and divergence well below expected speciation levels [disturbance-vicariance (DV) hypothesis]. To date there have been few tests of the expectations of the DV hypothesis based on empirical studies of Neotropical floral and faunal communities. Herein we examine phylogeographical structure of Dendrobates tinctorius, an amphibian species endemic to the uplands of the eastern Guiana Shield, based on sampling of 114 individuals from 24 localities. Phylogenetic, nested clade, and dispersal-vicariance (DIVA) analyses of cytochrome b sequence data reveal the presence of two mitochondrial lineages that are associated with previously identified western and eastern uplands of this area. The geographical distribution of mitochondrial haplotypes and the results of DIVA and coalescent analyses suggest that there has been extensive secondary contact between these lineages indicating a complex history of connectivity between these western and eastern highlands, supporting the predictions of the DV hypothesis.     

Lineage specific histories of Mycobacterium tuberculosis dispersal in Africa and Eurasia

Mycobacterium tuberculosis (M.tb) is a globally distributed, obligate pathogen of humans that can be divided into seven clearly defined lineages. An emerging consensus places the origin and global dispersal of M.tb within the past 6,000 years: identifying how the ancestral clone of M.tb spread and differentiated within this timeframe is important for identifying the ecological drivers of the current pandemic. We used Bayesian phylogeographic inference to reconstruct the migratory history of M.tb in Africa and Eurasia and to investigate lineage specific patterns of spread from a geographically diverse sample of 552 M.tb genomes. Applying evolutionary rates inferred with ancient M.tb genome calibration, we estimated the timing of major events in the migratory history of the pathogen. Inferred timings contextualize M.tb dispersal within historical phenomena that altered patterns of connectivity throughout Africa and Eurasia: trans‐Indian Ocean trade in spices and other goods, the Silk Road and its predecessors, the expansion of the Roman Empire, and the European Age of Exploration. We found that Eastern Africa and Southeast Asia have been critical in the dispersal of M.tb. Our results further reveal that M.tb populations have grown through range expansion, as well as in situ, and delineate the independent evolutionary trajectories of bacterial subpopulations underlying the current pandemic.     

Pan‐African phylogeography of a model organism,the African clawed frog ‘Xenopus laevis’

The African clawed frog Xenopus laevis has a large native distribution over much of sub‐Saharan Africa and is a model organism for research, a proposed disease vector, and an invasive species. Despite its prominent role in research and abundance in nature, surprisingly little is known about the phylogeography and evolutionary history of this group. Here, we report an analysis of molecular variation of this clade based on 17 loci (one mitochondrial, 16 nuclear) in up to 159 individuals sampled throughout its native distribution. Phylogenetic relationships among mitochondrial DNA haplotypes were incongruent with those among alleles of the putatively female‐specific sex‐determining gene DM‐W, in contrast to the expectation of strict matrilineal inheritance of both loci. Population structure and evolutionarily diverged lineages were evidenced by analyses of molecular variation in these data. These results further contextualize the chronology, and evolutionary relationships within this group, support the recognition of X. laevis sensu stricto, X. petersii, X. victorianus and herein revalidated X. poweri as separate species. We also propose that portions of the currently recognized distributions of X. laevis (north of the Congo Basin) and X. petersii (south of the Congo Basin) be reassigned to X. poweri.     

Migratory flyway and geographical distance are barriers to the gene flow of influenza virus among North American birds

Despite the importance of migratory birds in the ecology and evolution of avian influenza virus (AIV), there is a lack of information on the patterns of AIV spread at the intra-continental scale. We applied a variety of statistical phylogeographic techniques to a plethora of viral genome sequence data to determine the strength, pattern and determinants of gene flow in AIV sampled from wild birds in North America. These analyses revealed a clear isolation-by-distance of AIV among sampling localities. In addition, we show that phylogeographic models incorporating information on the avian flyway of sampling proved a better fit to the observed sequence data than those specifying homogeneous or random rates of gene flow among localities. In sum, these data strongly suggest that the intra-continental spread of AIV by migratory birds is subject to major ecological barriers, including spatial distance and avian flyway.     

Plant self-incompatibility in natural populations: a critical assessment of recent theoretical and empirical advances

Self-incompatibility systems in plants are genetic systems that prevent self-fertilization in hermaphrodites through recognition and rejection of pollen expressing the same allelic specificity as that expressed in the pistils. The evolutionary properties of these self-recognition systems have been revealed through a fascinating interplay between empirical advances and theoretical developments. In 1939, Wright suggested that the main evolutionary force driving the genetic and molecular properties of these systems was strong negative frequency-dependent selection acting on pollination success. The empirical observation of high allelic diversity at the self-incompatibility locus in several species, followed by the discovery of very high molecular divergence among alleles in all plant families where the locus has been identified, supported Wright's initial theoretical predictions as well as many of its later developments. In the last decade, however, advances in the molecular characterization of the incompatibility reaction and in the analysis of allelic frequencies and allelic divergence from natural populations have stimulated new theoretical investigations that challenged some important assumptions of Wright's model of gametophytic self-incompatibility. We here review some of these recent empirical and theoretical advances that investigated: (i) the hypothesis that S-alleles are selectively equivalent, and the evolutionary consequences of genetic interactions between alleles; (ii) the occurrence of frequency-dependent selection in female fertility; (iii) the evolutionary genetics of self-incompatibility systems in subdivided populations; (iv) the evolutionary implications of the self-incompatibility locus's genetic architecture; and (v) of its interactions with the genomic environment.     

Divergent and linked selection shape patterns of genomic differentiation between European and North American Atlantic salmon (Salmo salar)

As populations diverge many processes can shape genomic patterns of differentiation. Regions of high differentiation can arise due to divergent selection acting on selected loci, genetic hitchhiking of nearby loci, or through repeated selection against deleterious alleles (linked background selection); this divergence may then be further elevated in regions of reduced recombination. Atlantic salmon (Salmo salar) from Europe and North America diverged >600,000 years ago and despite some evidence of secondary contact, the majority of genetic data indicate substantial divergence between lineages. This deep divergence with potential gene flow provides an opportunity to investigate the role of different mechanisms that shape the genomic landscape during early speciation. Here, using 184,295 single nucleotide polymorphisms (SNPs) and 80 populations, we investigate the genomic landscape of differentiation across the Atlantic Ocean with a focus on highly differentiated regions and the processes shaping them. We found evidence of high (mean F_ST = 0.26) and heterogeneous genomic differentiation between continents. Genomic regions associated with high trans‐Atlantic differentiation ranged in size from single loci (SNPs) within important genes to large regions (1–3 Mbp ) on four chromosomes (Ssa06, Ssa13, Ssa16 and Ssa19). These regions showed signatures consistent with selection, including high linkage disequilibrium, despite no significant reduction in recombination. Genes and functional enrichment of processes associated with differentiated regions may highlight continental differences in ocean navigation and parasite resistance. Our results provide insight into potential mechanisms underlying differences between continents, and evidence of near‐fixed and potentially adaptive trans‐Atlantic differences concurrent with a background of high genome‐wide differentiation supports subspecies designation in Atlantic salmon.     

Differential patterns of male and female mtDNA exchange across the Atlantic Ocean in the blue mussel, Mytilus edulis

Comparisons among loci with differing modes of inheritance can reveal unexpected aspects of population history. We employ a multilocus approach to ask whether two types of independently assorting mitochondrial DNAs (maternally and paternally inherited: F- and M-mtDNA) and a nuclear locus (ITS) yield concordant estimates of gene flow and population divergence. The blue mussel, Mytilus edulis, is distributed on both North American and European coastlines and these populations are separated by the waters of the Atlantic Ocean. Gene flow across the Atlantic Ocean differs among loci, with F-mtDNA and ITS showing an imprint of some genetic interchange and M-mtDNA showing no evidence for gene flow. Gene flow of F-mtDNA and ITS causes trans-Atlantic population divergence times to be greatly underestimated for these loci, although a single trans-Atlantic population divergence time (1.2 MYA) can be accommodated by considering all three loci in combination in a coalescent framework. The apparent lack of gene flow for M-mtDNA is not readily explained by different dispersal capacities of male and female mussels. A genetic barrier to M-mtDNA exchange between North American and European mussel populations is likely to explain the observed pattern, perhaps associated with the double uniparental system of mitochondrial DNA inheritance.     

Phylogeography of a pan-Atlantic abyssal protobranch bivalve: implications for evolution in the Deep Atlantic

The deep sea is a vast and essentially continuous environment with few obvious barriers to gene flow. How populations diverge and new species form in this remote ecosystem is poorly understood. Phylogeographical analyses have begun to provide some insight into evolutionary processes at bathyal depths (<3000 m), but much less is known about evolution in the more extensive abyssal regions (>3000 m). Here, we quantify geographical and bathymetric patterns of genetic variation (16S rRNA mitochondrial gene) in the protobranch bivalve Ledella ultima, which is one of the most abundant abyssal protobranchs in the Atlantic with a broad bathymetric and geographical distribution. We found virtually no genetic divergence within basins and only modest divergence among eight Atlantic basins. Levels of population divergence among basins were related to geographical distance and were greater in the South Atlantic than in the North Atlantic. Ocean‐wide patterns of genetic variation indicate basin‐wide divergence that exceeds what others have found for abyssal organisms, but considerably less than bathyal protobranchs across similar geographical scales. Populations on either side of the Mid‐Atlantic Ridge in the North Atlantic differed, suggesting the Ridge might impede gene flow at abyssal depths. Our results indicate that abyssal populations might be quite large (cosmopolitan), exhibit only modest genetic structure and probably provide little potential for the formation of new species.