Three roads diverged? Routes to phylogeographic inference

doi:10.1016/j.tree.2010.08.010

Trends in Ecology & Evolution

Volume 25, Issue 11, November 2010, Pages 626-632

https://doi.org/10.1016/j.tree.2010.08.010 Get rights and content

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.

Section snippets

Emerging pathways of phylogeographic inference

The influence of phylogeography is spreading throughout biology. Among other examples, phylogeographic techniques have enabled us to infer the origins of mice [1], modern humans 2, 3, and man's “best friend”, the domesticated dog [4]. Phylogeographic analyses also enable public health officials to understand the origin and spread of emerging infectious diseases 5, 6, 7, 8. In spite of these successes, disagreement and confusion persist regarding the most effective ways to learn about

Comparative approach

For almost a decade, supporters of NCPA [14] and of model-based phylogeographic methods [20] have argued over the merits of each method. Similar to previous debate in phylogenetics, this has generated several positive outcomes; most importantly, investigators now vigorously question assumptions when analyzing geographic, demographic and evolutionary data [21]. Nevertheless, this long-standing discussion has introduced considerable confusion 22, 23. However, the debate now seems to be coming to

Spatial diffusion approach

As an alternative to the comparative method and its Bayesian extensions, we now highlight recent developments towards model-based approaches that take a probabilistic perspective on spatial diffusion 15, 16. Although much of this work was implemented as part of a comprehensive statistical inference package that led to fruitful advances in demographic models 42, 43, we want to emphasize that, unlike spatial coalescent approaches [44], phylogenetic diffusion models do not infer population-based

Population genetics approach

By far the most popular statistical approaches to phylogeography rely on the structured-coalescent framework 9, 44, 61, 62. In general, these methods assume that evolutionary trees are random draws from some underlying population-level process [9]. Essentially, population-level processes fossilize their histories as evolutionary trees that we indirectly view through molecular sequence and other data 10, 61. These processes include selection, migration, population size changes and recombination.

Three routes to the same destination

With the expansion of phylogeography in new Bayesian directions involving NCPA and spatial diffusion, the field might seem to be fragmenting. However, we believe all three approaches address the same basic question posed at the beginning of this article: what are the most effective ways to learn about phylogeographic processes from geospatially identified molecular sequence data? Essentially, we believe that all three frameworks produce effective answers, if we know what questions to ask. Under

Acknowledgments

We thank Chris Simon, Allen Rodrigo, Brian C. Carstens, H. Lisle Gibbs, Laura S. Kubatko, Peter Beerli and Ioanna Manolopoulou for their comments and suggestions. We also thank the National Evolutionary Synthesis Center (NSF #EF-0423641) for fostering our collaboration. Alexei J. Drummond contributed greatly to discussions on an earlier version of this paper. EWB is supported in part by the National Science Foundation under Agreement No. 0635561. PL is supported by a postdoctoral fellowship

Glossary

Ancestral history: any information about the direct ancestors of a sample of molecular sequences. This term can refer, for example, to inferred sequence composition or phenotype, such as geography, and is often associated with a time scale.
Approximate Bayesian computation (ABC): simulation technique used to draw statistical inference based on data summaries, often used when computation of the full data likelihood is impractical.
Bayes factor: ratio of the marginal likelihoods of a given data set

References (86)

J. Searle
Of mice and (Viking?) men: phylogeography of British and Irish house mice
Proc. R. Soc. B Biol. Sci.
(2009)
N. Fagundes
Statistical evaluation of alternative models of human evolution
Proc. Natl. Acad. Sci. U. S. A
(2007)
J. Li
Worldwide human relationships inferred from genome-wide patterns of variation
Science
(2008)
B. von Holdt
Genome-wide SNP and haplotype analyses reveal a rich history underlying dog domestication
Nature
(2010)
R. Biek
A virus reveals population structure and recent demographic history of its carnivore host
Science
(2006)
R. Biek
A high-resolution genetic signature of demographic and spatial expansion in epizootic rabies virus
Proc. Natl. Acad. Sci. U. S. A.
(2007)
G. Smith
Origins and evolutionary genomics of the 2009 swine-origin H1N1 influenza A epidemic
Nature
(2009)
P. Lemey
Reconstructing the initial global spread of a human influenza pandemic: a Bayesian spatial--temporal model for the global spread of H1N1
PLoS Curr. Influenza
(2009)
R. Nielsen et al.
Statistical inferences in phylogeography
Mol. Ecol.
(2009)
L. Knowles
Statistical phylogeography
Annu. Rev. Ecol. Evol. Syst.
(2009)

J. Avise

Phylogeography: retrospect and prospect

J. Biogeog.

(2009)

M. Hickerson

Phylogeography's past, present, and future: 10 years after Avise, 2000

Mol. Phylogenet. Evol.

(2010)

S. Brooks

Assessing the effect of genetic mutation: a Bayesian framework for determining population history from DNA sequence data

A. Templeton

Coalescent-based, maximum likelihood inference in phylogeography

Mol. Ecol.

(2010)

P. Lemey

Bayesian phylogeography finds its roots

PLoS Comput. Biol.

(2009)

P. Lemey

Phylogeography takes a relaxed random walk in continuous space and time

Mol. Biol. Evol.

(2010)

J. Hey

Isolation with migration models for more than two populations

Mol. Biol. Evol.

(2010)

J. Heled et al.

Bayesian inference of species trees from multilocus data

Mol. Biol. Evol.

(2010)

P. Beerli et al.

Unified framework to evaluate panmixia and migration direction among multiple sampling locations

Genetics

(2010)

M. Beaumont

In defence of model-based inference in phylogeography

Mol. Ecol.

(2010)

A. Camargo

Phylogeography of the frog Leptodactylus validus (Amphibia: Anura): patterns and timing of colonization events in the Lesser Antilles

Mol. Phylogenet. Evol.

(2009)

I. Phillipsen et al.

Phylogeography of a stream-dwelling frog ( Pseudacris cadaverina) in southern California

Mol. Phylogenet. Evol.

(2009)

H. Gante

Diversification within glacial refugia: tempo and mode of evolution of the polytypic fish Barbus sclateri

Mol. Ecol.

(2009)

J. Felsenstein

Phylogenies and the comparative method

Am. Nat.

(1985)

P. Harvey et al.

The Comparative Method in Evolutionary Biology

(1991)

A. Templeton

Nested clade analyses of phylogeographic data: testing hypotheses about gene flow and population history

Mol. Ecol.

(1998)

M. Panchal et al.

The automation and evaluation of nested clade phylogeographic analysis

Evolution

(2007)

A. Templeton

A cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping. III. Cladogram estimation

Genetics

(1992)

M. Clement

TCS: a computer program to estimate gene genealogies

Mol. Ecol.

(2000)

D. Posada et al.

Intraspecific gene genealogies: trees grafting into networks

Trends Ecol. Evol.

(2001)

A. Templeton

A cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping. I. Basic theory and an analysis of alcohol dehydrogenase activity in Drosophilia

Genetics

(1987)

A. Templeton

A cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping. IV. Nested analyses with cladogram uncertainty and recombination

Genetics

(1993)

I. Cassens

Evaluating intraspecific “network” construction methods using simulated sequence data: do existing algorithms outperform the global maximum parsimony approach?

Syst. Biol.

(2005)

L. Knowles

Why does a method that fails continue to be used?

Evolution

(2008)

A. Templeton

A maximum likelihood framework for cross validation of phylogeographic hypotheses

M. Panchal et al.

Evaluating nested clade phylogenetic analysis under models of restricted gene flow

Syst. Biol.

(2010)

K. Wong

Alignment uncertainty and genomic analysis

Science

(2008)

B. Redelings et al.

Robust inferences from ambiguous alignments

M. Suchard

Bayesian selection of continuous-time Markov chain evolutionary models

Mol. Biol. Evol.

(2001)

B. Redelings et al.

Joint Bayesian estimation of alignment and phylogeny

Syst. Biol.

(2005)

A. Novák

StatAlign: an extendable software package for joint Bayesian estimation of alignments and evolutionary trees

Bioinformatics

(2008)

A. Drummond

Bayesian coalescent inference of past population dynamics from molecular sequences

Mol. Biol. Evol.

(2005)

V. Minin

Smooth skyride through a rough skyline: Bayesian coalescent-based inference of population dynamics

Mol. Biol. Evol.

(2008)

Cited by (88)

Population structure, genetic diversity and evolutionary history of the narrow endemic squirrel Sciurus meridionalis
2024, Global Ecology and Conservation
Endemic species and their genetic diversity are critical components of global biodiversity, yet they usually face a higher risk of decline and extinction than species with broader distributions. Loss of genetic diversity within geographically restricted and fragmented ranges plays a significant role in driving endemic species to extinction, threatening their ability to adapt to environmental challenges. Hence, considering genetic diversity and its spatial distribution within populations is mandatory in planning conservation actions for narrow endemic species. Here, we investigated the population genetic diversity, structure and evolutionary history of the Calabrian black squirrel Sciurus meridionalis, a species endemic to the Calabrian region in the south of the Italian peninsula. A fine-scale sampling procedure was carried out throughout the whole species range, and genetic variation was assessed from multiple mitochondrial and nuclear markers. Phylogeographic and population structure analyses showed three evolutionary significant units within the current species range, each characterised by extremely low levels of genetic diversity. Bayesian phylogeographic analyses suggested that Calabrian black squirrel populations experienced a range fragmentation during the Last Glacial Maximum (∼20 ka), which geographically mirrors the three main mountain massifs within the species range. Historical demographic reconstructions revealed a recent and still ongoing demographic expansion for Calabrian black squirrel populations, likely promoted by a concurrent expansion of the Calabrian black pine in this area. Overall, our findings clarify the designation of the evolutionary significant units for a narrow endemic squirrel and underscore the importance of analysing range-wide population genetic data and considering species interactions when defining management strategies for narrow endemic species’ conservation.
Emergence of an early SARS-CoV-2 epidemic in the United States
2021, Cell
Citation Excerpt :
Prior to Mardi Gras, our analyses demonstrated that Texas is more than twice as likely as the next most probable state to be the source of SARS-CoV-2 lineages in New Orleans, while SARS-CoV-2 in Shreveport likely originated from New Orleans itself (Figures 4A and 4B). Although these analyses point to Texas as a likely source of the Louisiana clade, our phylogeographic inference is limited by geographic and temporal sampling (Bloomquist et al., 2010). Therefore, we also investigated movement between New Orleans, Shreveport, and other U.S. states by analyzing human mobility patterns.
The emergence of the COVID-19 epidemic in the United States (U.S.) went largely undetected due to inadequate testing. New Orleans experienced one of the earliest and fastest accelerating outbreaks, coinciding with Mardi Gras. To gain insight into the emergence of SARS-CoV-2 in the U.S. and how large-scale events accelerate transmission, we sequenced SARS-CoV-2 genomes during the first wave of the COVID-19 epidemic in Louisiana. We show that SARS-CoV-2 in Louisiana had limited diversity compared to other U.S. states and that one introduction of SARS-CoV-2 led to almost all of the early transmission in Louisiana. By analyzing mobility and genomic data, we show that SARS-CoV-2 was already present in New Orleans before Mardi Gras, and the festival dramatically accelerated transmission. Our study provides an understanding of how superspreading during large-scale events played a key role during the early outbreak in the U.S. and can greatly accelerate epidemics.
Unequal sisters – Past and potential future range development of Anatolian and Hyrcanian brown frogs
2021, Zoology
Citation Excerpt :
We also changed the a priori distribution of the parameter ucld.mean into a uniform prior ranging from 0 to 1. In Beast, the phylogeographic history is reconstructed by the inference of the geographic history of genetic lineages (Bloomquist et al., 2010), so we linked each sequence to a geographic coordinate (Table S1). With the RRW option in Beast, dispersal rates are allowed to vary across branches, so we used this model to perform a RRW phylogeographic analysis to reconstruct the ancestral ranges of the Anatolian and Hyrcanian brown frog lineages.
Phylogeography can reconstruct historical evolutionary processes by comparing historical patterns of gene flow, divergence among species and by using species distribution models (SDM) upon geographic distribution. We investigate the phylogeographic patterns of Anatolian brown frogs including R. macrocnemis and R. tavasensis as well as the Hyrcanian brown frog, R. pseudodalmatina, using a fragment of the mitochondrial 16S rRNA gene for 145 specimens across the entire range of these frogs. We calculate parameters of molecular diversity, such as the number of variable sites (S), the number of haplotypes (h), haplotype diversity (Hd) and nucleotide diversity (π). We generated a haplotype network and used three methods (Neutrality tests, mismatch distributions and Bayesian skyline plots) to reconstruct the demographic histories of R. macrocnemis and R. pseudodalmatina. Finally, we used SDMs to predict the habitat suitability for three periods: The Present Day, the Last Glacial Maximum (LGM) and the future until 2070 for R. macrocnemis and R. pseudodalmatina. Our phylogenetic analyses support a late Miocene origin of Anatolian and Hyrcanian lineages. Hyrcanian brown frogs were enclosed in lowlands of the southern coast of the Caspian Sea after the uplift of the Elburz range and the Armenian plateau. The formation of a salinity belt from the north Aegean corridor (the south western Turkey) to northward during the Late Tortonian led to the subdivision of ancestor of the Anatolian lineage into today isolated western and eastern populations. The salinity belt had a considerable impact on the divergence of R. tavasensis from R. macrocnemis. Combined historical demographic analyses and SDMs revealed a rapid expansion occurring during the Pleistocene in R. macrocnemis and R. pseudodalmatina. Currently, suitable habitat for R. macrocnemis has declined compared to the LGM, and the species is predicted to do even worse under future climatic conditions. In contrast, R. pseudodalmatina found suitable habitat from the LGM to present within its restricted distribution area; it is predicted to do fine even under future climate.
HIV Rebound Is Predominantly Fueled by Genetically Identical Viral Expansions from Diverse Reservoirs
2019, Cell Host and Microbe
Citation Excerpt :
Visualizing genetic similarity using haplotype networks further confirmed the diverse origins of the rebound viruses (Figure S2). To further unravel to what extent specific compartments can act as the source of rebound viruses, we employed a phylogeographic approach that is commonly used in molecular epidemiological research (Bloomquist et al., 2010; Faria et al., 2011) and quantified the rebound virus emergence events from each cell subset. The estimated number of rebound virus emergence events for 3 participants are summarized in a radar plot in Figure 3D and across participants in Figure S4.
Viral rebound upon stopping combined antiretroviral therapy poses a major barrier toward an HIV cure. Cellular and anatomical sources responsible for reinitiating viral replication remain a subject of ardent debate, despite extensive research efforts. To unravel the source of rebounding viruses, we conducted a large-scale HIV-STAR (HIV-1 sequencing before analytical treatment interruption to identify the anatomically relevant HIV reservoir) clinical trial. We collected samples from 11 participants and compared the genetic composition of (pro)viruses collected under treatment from different cellular and anatomical compartments with that of plasma viruses sampled during analytical treatment interruption. We found a remarkably heterogeneous source of viral rebound. In addition, irrespective of the compartment or cell subset, genetically identical viral expansions played a significant role in viral rebound. Our study suggests that although there does not seem to be a primary source for rebound HIV, cellular proliferation is an important driver of HIV persistence and should therefore be considered in future curative strategies.
Time-scaled phylogeography of complete Zika virus genomes using discrete and continuous space diffusion models
2019, Infection, Genetics and Evolution
Citation Excerpt :
These models allow a more realistic reconstruction of spatial movements because, unlike discrete models, they do not necessarily need the ancestral location to be represented in the sampling location set (Lemey et al., 2010). The differences between discrete and continuous phylogeographic models have been efficiently described in previous reviews (Bloomquist et al., 2010; Faria et al., 2017). The aim of this study was to infer the origin and dispersion routes of ZIKV in the world using a classical discrete method and to reconstruct the recent epidemic in the Americas using a continuous phylogeographical method to better describe the local spread of ZIKV and to make hypothesis about the eco/epidemiology of the virus.
Zika virus (ZIKV), a vector-borne infectious agent that has recently been associated with neurological diseases and congenital microcephaly, was first reported in the Western hemisphere in early 2015.
A number of authors have reconstructed its epidemiological history using advanced phylogenetic approaches, and the majority of Zika phylogeography studies have used discrete diffusion models. Continuous space diffusion models make it possible to infer the possible origin of the virus in real space by reconstructing its ancestral location on the basis of geographical coordinates deduced from the latitude and longitude of the sampling locations. We analysed all the ZIKV complete genome isolates whose sampling times and localities were available in public databases at the time the study began, using a Bayesian approach for discrete and continuous phylogeographic reconstruction.
The discrete phylogeographic analysis suggested that ZIKV emerged to become endemic/epidemic in the first decade of the 1900s in the Ugandan rainforests, and then reached Western Africa and Asia between the 1930s and 1950s. After a long period of about 40 years, it spread to the Pacific islands and reached Brazil from French Polynesia. Continuous phylogeography of the American epidemic showed that the virus entered in north-eastern Brazil in late 2012 and started to spread in early 2013 from two high probability regions: one corresponding to the entire north-east Brazil and the second surrounding the city of Rio de Janeiro, in a mainly northwesterly direction to Central America, the north-western countries of south America and the Caribbean islands. Our data suggest its cryptic circulation in both French Polynesia and Brazil, thus raising questions about the mechanisms underlying its undetected persistence in the absence of a known animal reservoir, and underline the importance of continuous diffusion models in making more reliable phylogeographic reconstructions of emerging viruses.
The origin and spread of the southern black ant, a widely distributed leaf-cutting ant
2023, Journal of Biogeography

View all citing articles on Scopus

View full text

Trends in Ecology & Evolution

OpinionThree roads diverged? Routes to phylogeographic inference

Section snippets

Emerging pathways of phylogeographic inference

Comparative approach

Spatial diffusion approach

Population genetics approach

Three routes to the same destination

Acknowledgments

Glossary

Of mice and (Viking?) men: phylogeography of British and Irish house mice

Proc. R. Soc. B Biol. Sci.

Statistical evaluation of alternative models of human evolution

Proc. Natl. Acad. Sci. U. S. A

Worldwide human relationships inferred from genome-wide patterns of variation

Science

Genome-wide SNP and haplotype analyses reveal a rich history underlying dog domestication

Nature

A virus reveals population structure and recent demographic history of its carnivore host

Science

A high-resolution genetic signature of demographic and spatial expansion in epizootic rabies virus

Proc. Natl. Acad. Sci. U. S. A.

Origins and evolutionary genomics of the 2009 swine-origin H1N1 influenza A epidemic

Nature

Reconstructing the initial global spread of a human influenza pandemic: a Bayesian spatial--temporal model for the global spread of H1N1

PLoS Curr. Influenza

Statistical inferences in phylogeography

Mol. Ecol.

Statistical phylogeography

Annu. Rev. Ecol. Evol. Syst.

Phylogeography: retrospect and prospect

J. Biogeog.

Phylogeography's past, present, and future: 10 years after Avise, 2000

Mol. Phylogenet. Evol.

Assessing the effect of genetic mutation: a Bayesian framework for determining population history from DNA sequence data

Coalescent-based, maximum likelihood inference in phylogeography

Mol. Ecol.

Bayesian phylogeography finds its roots

PLoS Comput. Biol.

Phylogeography takes a relaxed random walk in continuous space and time

Mol. Biol. Evol.

Isolation with migration models for more than two populations

Mol. Biol. Evol.

Bayesian inference of species trees from multilocus data

Mol. Biol. Evol.

Unified framework to evaluate panmixia and migration direction among multiple sampling locations

Genetics

In defence of model-based inference in phylogeography

Mol. Ecol.

Phylogeography of the frog Leptodactylus validus (Amphibia: Anura): patterns and timing of colonization events in the Lesser Antilles

Mol. Phylogenet. Evol.

Phylogeography of a stream-dwelling frog ( Pseudacris cadaverina) in southern California

Mol. Phylogenet. Evol.

Diversification within glacial refugia: tempo and mode of evolution of the polytypic fish Barbus sclateri

Mol. Ecol.

Phylogenies and the comparative method

Am. Nat.

The Comparative Method in Evolutionary Biology

Nested clade analyses of phylogeographic data: testing hypotheses about gene flow and population history

Mol. Ecol.

The automation and evaluation of nested clade phylogeographic analysis

Evolution

A cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping. III. Cladogram estimation

Genetics

TCS: a computer program to estimate gene genealogies

Mol. Ecol.

Intraspecific gene genealogies: trees grafting into networks

Trends Ecol. Evol.

A cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping. I. Basic theory and an analysis of alcohol dehydrogenase activity in Drosophilia

Genetics

A cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping. IV. Nested analyses with cladogram uncertainty and recombination

Genetics

Evaluating intraspecific “network” construction methods using simulated sequence data: do existing algorithms outperform the global maximum parsimony approach?

Syst. Biol.

Why does a method that fails continue to be used?

Evolution

A maximum likelihood framework for cross validation of phylogeographic hypotheses

Evaluating nested clade phylogenetic analysis under models of restricted gene flow

Syst. Biol.

Alignment uncertainty and genomic analysis

Opinion
Three roads diverged? Routes to phylogeographic inference