Phylogeographic Inference Research Articles

Erratum to: DNA barcoding of aquatic insects reveals unforeseen diversity and recurrent population divergence patterns through broad-scale sampling in northern Canada

Arctic and subarctic environments are among the most inaccessible regions in the world, making biological surveys difficult to conduct. Thus, the insect fauna of these regions has remained inadequately surveyed. The aquatic insect orders Ephemeroptera, Plecoptera, and Trichoptera (EPTs) are particularly abundant and diverse at high latitudes, playing key roles in trophic chains where nutrients are scarce. However, particular aspects of their life cycle make them difficult to study. Specifically, species-level identification requires last-instar larvae or adults which, because they are short lived, are typically not available for all taxa during a particular collecting event. With the initial goal of surveying the biodiversity of these insect groups, we sampled ca.10,000 EPT specimens from 12 locations across northern Canada over two years. Approximately 800 of these were subsequently selected for COI DNA barcoding. Overall, we identified 155 EPT species (58 Ephemeroptera, 41 Plecoptera, 56 Trichoptera) based on a 2% divergence criterion. Compared to other similar studies on EPTs we found higher (particularly among the Plecoptera and Ephemeroptera) and more even diversity, potentially reflecting environmental differences in sampling localities. We further assessed phylogeographic divergence patterns among seven species, finding that eastern and western populations diverged during the Pleistocene Epoch (<2.5 Ma), with overlapping time frames. This finding highlights the role of potential glacial refugia and subsequent recolonization, as well as the dispersal potential of some EPT species. This study exemplifies how large-scale DNA-based surveys can be combined with phylogeographic inference to better understand the biodiversity and natural history of northern aquatic insect communities.

Tracking the Origin and Deciphering the Phylogenetic Relationship of Porcine Epidemic Diarrhea Virus in Ecuador.

In 2010, new Chinese strains of porcine epidemic diarrhea virus (PEDV), clinically more severe than the classical strains, emerged. These strains were spread to United States in 2013 through an intercontinental transmission from China with further spreading across the world, evidencing the emergent nature of these strains. In the present study, an analysis of PEDV field sequences from Ecuador was conducted by comparing all the PEDV S gene sequences available in the GenBank database. Phylogenetic comparisons and Bayesian phylogeographic inference based on complete S gene sequences were also conducted to track the origin and putative route of PEDV. The sequence from the PED-outbreak in Ecuador was grouped into the clade II of PEDV genogroup 2a together with other sequences of isolates from Mexico, Canada, and United States. The phylogeographic study revealed the emergence of the Chinese PEDV strains, followed by spreading to US in 2013, from US to Korea, and later the introduction of PEDV to Canada, Mexico, and Ecuador directly from the US. The sources of imports of live swine in Ecuador in 2014 were mainly from Chile and US. Thus, this movement of pigs is suggested as the main way for introducing PEDV to Ecuador.

Species Delimitation with Gene Flow

Species are commonly thought to be evolutionarily independent in a way that populations within a species are not. In recent years, studies that seek to identify evolutionarily independent lineages (i.e., to delimit species) using genetic data have typically adopted multispecies coalescent approaches that assume that evolutionary independence is formed by the differential sorting of ancestral alleles due to genetic drift. However, gene flow appears to be common among populations and nascent species, and while this process may inhibit lineage divergence (and thus independence), it is usually not explicitly considered when delimiting species. In this article, we apply Phylogeographic Inference using Approximate Likelihoods (PHRAPL), a recently described method for phylogeographic model selection, to species delimitation. We describe an approach to delimiting species using PHRAPL that attempts to account for both genetic drift and gene flow, and we compare the method's performance to that of a popular delimitation approach (BPP) using both simulated and empirical datasets. PHRAPL generally infers the correct demographic-delimitation model when the generating model includes gene flow between taxa, given a sufficient amount of data. When the generating model includes only isolation in the recent past, PHRAPL will in some cases fail to differentiate between gene flow and divergence, leading to model misspecification. Nevertheless, the explicit consideration of gene flow by PHRAPL is an important complement to existing delimitation approaches, particularly in systems where gene flow is likely important. [approximate likelihoods; coalescent simulations; genealogical divergence index; Homo sapiens; isolation-with-migration; multispecies coalescent; Sarracenia; Scincella.].

Genetics, Gene Flow, and Glaciation: The Case of the South American Limpet Nacella mytilina.

Glacial episodes of the Quaternary, and particularly the Last Glacial Maximum (LGM) drastically altered the distribution of the Southern-Hemisphere biota, principally at higher latitudes. The irregular coastline of Patagonia expanding for more than 84.000 km constitutes a remarkable area to evaluate the effect of Quaternary landscape and seascape shifts over the demography of near-shore marine benthic organisms. Few studies describing the biogeographic responses of marine species to the LGM have been conducted in Patagonia, but existing data from coastal marine species have demonstrated marked genetic signatures of post-LGM recolonization and expansion. The kelp-dweller limpet Nacella mytilina is broadly distributed along the southern tip of South America and at the Falkland/Malvinas Islands. Considering its distribution, abundance, and narrow bathymetry, N. mytilina represents an appropriate model to infer how historical and contemporary processes affected the distribution of intraspecific genetic diversity and structure along the southern tip of South America. At the same time, it will be possible to determine how life history traits and the ecology of the species are responsible for the current pattern of gene flow and connectivity across the study area. We conducted phylogeographic and demographic inference analyses in N. mytilina from 12 localities along Pacific Patagonia (PP) and one population from the Falkland/Malvinas Islands (FI). Analyses of the mitochondrial gene COI in 300 individuals of N. mytilina revealed low levels of genetic polymorphism and the absence of genetic differentiation along PP. In contrast, FI showed a strong and significant differentiation from Pacific Patagonian populations. Higher levels of genetic diversity were also recorded in the FI population, together with a more expanded genealogy supporting the hypothesis of glacial persistence of the species in these islands. Haplotype genealogy, and mismatch analyses in the FI population recognized an older and more complex demographic history than in PP. Demographic reconstructions along PP suggest a post-LGM expansion process (7.5 ka), also supported by neutrality tests, mismatch distribution and maximum parsimony haplotype genealogies. Migration rate estimations showed evidence of asymmetrical gene flow from PP to FI. The absence of genetic differentiation, the presence of a single dominant haplotype, high estimated migration rates, and marked signal of recent demographic growth, support the hypothesis of rapid post-glacial expansion in N. mytilina along PP. This expansion could have been sustained by larval and rafting-mediated dispersal of adults from northernmost populations following the Cape Horn Current System. Marked genetic differentiation between PP and FI could be explained through differences in their respective glacial histories. During the LGM, Pacific Patagonia (PP) was almost fully covered by the Patagonian Ice Sheet, while sheet coverage in the FI ice was restricted to small cirques and valleys. As previously recorded in the sister-species N. magellanica, the FI rather than represent a classical glacial refugium for N. mytilina, seems to represent a sink area and/or a secondary contact zone. Accordingly, historical and contemporary processes, contrasting glacial histories between the analyzed sectors, as well as life history traits constitute the main factors explaining the current biogeographical patterns of most shallow Patagonian marine benthic organisms.

Phylogeography by diffusion on a sphere: whole world phylogeography.

BackgroundTechniques for reconstructing geographical history along a phylogeny can answer many questions of interest about the geographical origins of species. Bayesian models based on the assumption that taxa move through a diffusion process have found many applications. However, these methods rely on diffusion processes on a plane, and do not take the spherical nature of our planet in account. Performing an analysis that covers the whole world thus does not take in account the distortions caused by projections like the Mercator projection.ResultsIn this paper, we introduce a Bayesian phylogeographical method based on diffusion on a sphere. When the area where taxa are sampled from is small, a sphere can be approximated by a plane and the model results in the same inferences as with models using diffusion on a plane. For taxa sampled from the whole world, we obtain substantial differences. We present an efficient algorithm for performing inference in a Markov Chain Monte Carlo (MCMC) algorithm, and show applications to small and large samples areas. We compare results between planar and spherical diffusion in a simulation study and apply the method by inferring the origin of Hepatitis B based on sequences sampled from Eurasia and Africa.ConclusionsWe describe a framework for performing phylogeographical inference, which is suitable when the distortion introduced by map projections is large, but works well on a smaller scale as well. The framework allows sampling tips from regions, which is useful when the exact sample location is unknown, and placing prior information on locations of clades in the tree. The method is implemented in the GEO_SPHERE package in BEAST 2, which is open source licensed under LGPL and allows joint tree and geography inference under a wide range of models.

Phylogeographic model selection leads to insight into the evolutionary history of four-eyed frogs

Phylogeographic research investigates biodiversity at the interface between populations and species, in a temporal and geographic context. Phylogeography has benefited from analytical approaches that allow empiricists to estimate parameters of interest from the genetic data (e.g., θ = 4Neμ, population divergence, gene flow), and the widespread availability of genomic data allow such parameters to be estimated with greater precision. However, the actual inferences made by phylogeographers remain dependent on qualitative interpretations derived from these parameters’ values and as such may be subject to overinterpretation and confirmation bias. Here we argue in favor of using an objective approach to phylogeographic inference that proceeds by calculating the probability of multiple demographic models given the data and the subsequent ranking of these models using information theory. We illustrate this approach by investigating the diversification of two sister species of four-eyed frogs of northeastern Brazil using single nucleotide polymorphisms obtained via restriction-associated digest sequencing. We estimate the composite likelihood of the observed data given nine demographic models and then rank these models using Akaike information criterion. We demonstrate that estimating parameters under a model that is a poor fit to the data is likely to produce values that lead to spurious phylogeographic inferences. Our results strongly imply that identifying which parameters to estimate from a given system is a key step in the process of phylogeographic inference and is at least as important as being able to generate precise estimates of these parameters. They also illustrate that the incorporation of model uncertainty should be a component of phylogeographic hypothesis tests.

A new amplicon based approach of whole mitogenome sequencing for phylogenetic and phylogeographic analysis: An example of East African white-eyes (Aves, Zosteropidae)

Classical Sanger sequencing is still frequently used to generate sequence data for phylogenetic and phylogeographic inference. In this contribution we present a novel approach to genotype whole mitogenomic haplotypes using Illumina MiSeq reads from indexed amplicons. Our new approach reduces preparation time by multiplexing loci within a single or few PCR reactions and by plate format library construction. The use of paired-end reads allows covering amplicons of about 0.5kb and thus no nebulisation and assembly are necessary. We tested the power and effectiveness of this technique by analysing the mitogenomic diversity of East African white-eye bird species (Zosteropidae), a taxonomically highly diverse and complex species flock found in various ecosystems spread across major parts of Africa. We compare the newly generated mitogenomic data set with published data of three mitochondrial genes for a similar set of populations and taxa. The comparison demonstrates that our new procedure represents a cost effective use of NGS for medium throughput phylogenetic analyses. Using this method, we were able to increase the amount of phylogenetic information significantly, while reducing the costs and effort in the laboratory. The mitogenomic data show a higher resolution than previous studies providing higher support and new insights in the relationships of Zosterops species. Our data suggest to split Z. poliogaster into four distinct species, three of which had previously been proposed: Z. silvanus, Z. mbulensis, Z. kikyuensis and Z. kulalensis. Our approach allows the genotyping of whole mitogenomes for a large number of individuals and thus allows more reliable reconstruction of phylogenetic and phylogeographic relationships – also for non-model organisms.

Phylogenetic Analysis of Human Immunodeficiency Virus Type 2 Group B.

Context:Human immunodeficiency virus type 2 (HIV-2) infections are mainly restricted to West Africa; however, in the recent years, the prevalence of HIV-2 is a growing concern in some European countries and the Southwestern region of India. Despite the presence of different HIV-2 groups, only A and B Groups have established human-to-human transmission chains.Aims:This work aimed to evaluate the phylogeographic inference of HIV-2 Group B worldwide to estimate their data of origin and the population dynamics.Materials and Methods:The evolutionary rates, the demographic history for HIV-2 Group B dataset, and the phylogeographic analysis were estimated using a Bayesian approach. The viral gene flow analysis was used to count viral gene out/in flow among different locations.Results:The root of the Bayesian maximum clade credibility tree of HIV-2 Group B dated back to 1957. The demographic history of HIV-2 Group B showed that the epidemic remained constant up to 1970 when started an exponential growth. From 1985 to early 2000s, the epidemic reached a plateau, and then it was characterized by two bottlenecks and a new plateau at the end of 2000s. Phylogeographic reconstruction showed that the most probable location for the root of the tree was Ghana. Regarding the viral gene flow of HIV-2 Group B, the only observed viral gene flow was from Africa to France, Belgium, and Luxembourg.Conclusions:The study gives insights into the origin, history, and phylogeography of HIV-2 Group B epidemic. The growing number of infections of HIV-2 worldwide indicates the need for strengthening surveillance.

Alpine endemic spiders shed light on the origin and evolution of subterranean species.

We designed a comparative study to unravel the phylogeography of two Alpine endemic spiders characterized by a different degree of adaptation to subterranean life: Troglohyphantes vignai (Araneae, Linyphiidae) and Pimoa rupicola (Araneae, Pimoidae), the latter showing minor adaptation to hypogean life. We sampled populations of the model species in caves and other subterranean habitats across their known geographical range in the Western Alps. By combining phylogeographic inferences and Ecological Niche Modeling techniques, we inferred the biogeographic scenario that led to the present day population structure of the two species. According to our divergent time estimates and relative uncertainties, the isolation of T. vignai and P. rupicola from their northern sister groups was tracked back to Middle–Late Miocene. Furthermore, the fingerprint left by Pleistocene glaciations on the population structure revealed by the genetic data, led to the hypothesis that a progressive adaptation to subterranean habitats occurred in T. vignai, followed by strong population isolation. On the other hand, P. rupicola underwent a remarkable genetic bottleneck during the Pleistocene glaciations, that shaped its present population structure. It seems likely that such shallow population structure is both the result of the minor degree of specialization to hypogean life and the higher dispersal ability characterizing this species. The simultaneous study of overlapping spider species showing different levels of adaptation to hypogean life, disclosed a new way to clarify patterns of biological diversification and to understand the effects of past climatic shift on the subterranean biodiversity.

Phylogeographic inference using Bayesian model comparison across a fragmented chorus frog species complex.

Fragmented species complexes provide an interesting system for investigating biogeographic history and the present distribution of genetic variation. Recent advances in sequencing technology and statistical phylogeography enable the collection and rigorous analysis of large multilocus data sets, but designing studies that produce meaningful phylogeographic inferences remains challenging. We implemented a Bayesian model comparison approach to investigate previous biogeographic hypotheses while simultaneously inferring the presence of genetic structure in a chorus frog species complex. The Illinois chorus frog (Pseudacris illinoensis), originally described as a subspecies of the broadly distributed Strecker's chorus frog (Pseudacris streckeri), occurs in small, disjunct regions associated with scarce sand prairie habitats that have been impacted by human development. We used high-throughput sequencing to develop and collect a multitiered genetic data set comprised of three different marker types (23 anonymous nuclear sequence loci, four mitochondrial genes and 14 microsatellite loci) designed to address questions across different evolutionary timescales. Phylogenetic analyses uncovered a deep divergence between populations in the Edwards Plateau of central Texas and all other P.streckeri/P. illinoensis populations, but suggest the disjunct distribution of P.illinoensis occurred more recently. Our best-supported migration model is consistent with the hypothesis that central Texas represented a refugium from which populations expanded via multiple routes. This model also indicates that disjunct northern and southern regions of P.illinoensis should be considered genetically distinct management units. Our study provides an evolutionary context for future studies and conservation efforts in P.illinoensis and demonstrates the utility of model-based approaches for phylogeographic inference.

Phylogeography of the Chinese endemic freshwater crab Sinopotamon acutum (Brachyura, Potamidae)

Phylogeographical inferences based on multiple molecular markers from different modes of inheritance (i.e. mtDNA and nucDNA) can help to explore drivers of current species divergence over different timescales and allow a deeper understanding of evolutionary processes. River systems and mountains, owing to their drainage networks and complex topography, are often a high priority when inferring external causes of phylogeographical patterns. Here, we selected the Chinese endemic freshwater crab Sinopotamon acutum inhabiting drainage networks across the Qinling–Daba Mountains, a watershed that divides northern from southern fauna in China, as a model species (1) to investigate whether river networks and mountain systems act as barriers in shaping the phylogeography of freshwater crabs and (2) to test the impact of historic environments and ecological habitats on crab distribution. Phylogenetic tree and network analyses based on mtDNA revealed divergence between eastern and western populations, but microsatellite DNA did not recover this. Frequent and male‐biased gene flow is evident. The adjacent tributaries, such as Weihe River 1 and Hanjiang River 1, although isolated by Qinling Mts, have low genetic variation in mtDNA. In addition, microsatellite DNA showed low variation between all rivers. Frequent gene flow, genetic variation and mitochondrial genetic structure in S. acutum indicated that the Qinling–Daba Mountains have had no effect on divergence in this species; instead, hydrological networks across these mountains serve as a continuous migration corridor for these crabs. Glacial refugia are probably responsible for the initial isolation of the eastern and western clades, and local adaptations after dispersal have further driven divergence. Cytonuclear discordance in phylogeographical patterns may be attributable to ongoing male‐biased gene flow.

New Routes to Phylogeography: A Bayesian Structured Coalescent Approximation.

Phylogeographic methods aim to infer migration trends and the history of sampled lineages from genetic data. Applications of phylogeography are broad, and in the context of pathogens include the reconstruction of transmission histories and the origin and emergence of outbreaks. Phylogeographic inference based on bottom-up population genetics models is computationally expensive, and as a result faster alternatives based on the evolution of discrete traits have become popular. In this paper, we show that inference of migration rates and root locations based on discrete trait models is extremely unreliable and sensitive to biased sampling. To address this problem, we introduce BASTA (BAyesian STructured coalescent Approximation), a new approach implemented in BEAST2 that combines the accuracy of methods based on the structured coalescent with the computational efficiency required to handle more than just few populations. We illustrate the potentially severe implications of poor model choice for phylogeographic analyses by investigating the zoonotic transmission of Ebola virus. Whereas the structured coalescent analysis correctly infers that successive human Ebola outbreaks have been seeded by a large unsampled non-human reservoir population, the discrete trait analysis implausibly concludes that undetected human-to-human transmission has allowed the virus to persist over the past four decades. As genomics takes on an increasingly prominent role informing the control and prevention of infectious diseases, it will be vital that phylogeographic inference provides robust insights into transmission history.

Disentangling the complex evolutionary history of the Western Palearctic blue tits (Cyanistes spp.) - phylogenomic analyses suggest radiation by multiple colonization events and subsequent isolation.

Isolated islands and their often unique biota continue to play key roles for understanding the importance of drift, genetic variation and adaptation in the process of population differentiation and speciation. One island system that has inspired and intrigued evolutionary biologists is the blue tit complex (Cyanistes spp.) in Europe and Africa, in particular the complex evolutionary history of the multiple genetically distinct taxa of the Canary Islands. Understanding Afrocanarian colonization events is of particular importance because of recent unconventional suggestions that these island populations acted as source of the widespread population in mainland Africa. We investigated the relationship between mainland and island blue tits using a combination of Sanger sequencing at a population level (20 loci; 12 500 nucleotides) and next-generation sequencing of single population representatives (>3 200 000 nucleotides), analysed in coalescence and phylogenetic frameworks. We found (i) that Afrocanarian blue tits are monophyletic and represent four major clades, (ii) that the blue tit complex has a continental origin and that the Canary Islands were colonized three times, (iii) that all island populations have low genetic variation, indicating low long-term effective population sizes and (iv) that populations on La Palma and in Libya represent relicts of an ancestral North African population. Further, demographic reconstructions revealed (v) that the Canary Islands, conforming to traditional views, hold sink populations, which have not served as source for back colonization of the African mainland. Our study demonstrates the importance of complete taxon sampling and an extensive multimarker study design to obtain robust phylogeographical inferences.

Comparative phylogeography highlights the double-edged sword of climate change faced by arctic- and alpine-adapted mammals.

Recent studies suggest that alpine and arctic organisms may have distinctly different phylogeographic histories from temperate or tropical taxa, with recent range contraction into interglacial refugia as opposed to post-glacial expansion out of refugia. We use a combination of phylogeographic inference, demographic reconstructions, and hierarchical Approximate Bayesian Computation to test for phylodemographic concordance among five species of alpine-adapted small mammals in eastern Beringia. These species (Collared Pikas, Hoary Marmots, Brown Lemmings, Arctic Ground Squirrels, and Singing Voles) vary in specificity to alpine and boreal-tundra habitat but share commonalities (e.g., cold tolerance and nunatak survival) that might result in concordant responses to Pleistocene glaciations. All five species contain a similar phylogeographic disjunction separating eastern and Beringian lineages, which we show to be the result of simultaneous divergence. Genetic diversity is similar within each haplogroup for each species, and there is no support for a post-Pleistocene population expansion in eastern lineages relative to those from Beringia. Bayesian skyline plots for four of the five species do not support Pleistocene population contraction. Brown Lemmings show evidence of late Quaternary demographic expansion without subsequent population decline. The Wrangell-St. Elias region of eastern Alaska appears to be an important zone of recent secondary contact for nearctic alpine mammals. Despite differences in natural history and ecology, similar phylogeographic histories are supported for all species, suggesting that these, and likely other, alpine- and arctic-adapted taxa are already experiencing population and/or range declines that are likely to synergistically accelerate in the face of rapid climate change. Climate change may therefore be acting as a double-edged sword that erodes genetic diversity within populations but promotes divergence and the generation of biodiversity.

Extracting transmission networks from phylogeographic data for epidemic and endemic diseases: Ebola virus in Sierra Leone, 2009 H1N1 pandemic influenza and polio in Nigeria.

BackgroundPhylogeography improves our understanding of spatial epidemiology. However, application to practical problems requires choices among computational tools to balance statistical rigor, computational complexity, sensitivity to sampling strategy and interpretability.MethodsWe introduce a fast, heuristic algorithm to reconstruct partially-observed transmission networks (POTN) that combines features of phylogenetic and transmission tree approaches. We compare the transmission network generated by POTN with existing algorithms (BEAST and SeqTrack), and discuss the benefits and challenges of phylogeographic analysis on examples of epidemic and endemic diseases: Ebola virus, H1N1 pandemic influenza and polio.ResultsFor the 2014 Sierra Leone Ebola virus outbreak and the 2009 H1N1 outbreak, all three methods provide similarly plausible transmission histories but differ in detail. For polio in northern Nigeria, we discuss performance trade-offs between the POTN and discrete phylogeography in BEAST and conclude that spatial history reconstruction is limited by under-sampling.ConclusionsPOTN is complementary to available tools on densely-sampled data, fails gracefully on under-sampled data and is scalable to accommodate larger datasets. We provide further evidence for the utility of phylogeography for understanding transmission networks of rapidly evolving epidemics. We propose simple heuristic criteria to identify how sampling rates and disease dynamics interact to determine fundamental limitations of phylogeographic inference.

The spatiotemporal expansion of human rabies and its probable explanation in mainland China, 2004-2013.

BackgroundHuman rabies is a significant public health concern in mainland China. However, the neglect of rabies expansion and scarce analyses of the dynamics have made the spatiotemporal spread pattern of human rabies and its determinants being poorly understood.MethodsWe collected geographic locations and timeline of reported human rabies cases, rabies sequences and socioeconomic variables for the years 2004-2013, and integrated multidisciplinary approaches, including epidemiological characterization, hotspots identification, risk factors analysis and phylogeographic inference, to explore the spread pattern of human rabies in mainland China during the last decade.ResultsThe results show that human rabies distribution and hotspots were expanding from southeastern regions to north or west regions, which could be associated with the evolution of the virus, especially the clade I-G. A Panel Poisson Regression analysis reveals that human rabies incidences had significant correlation with the education level, GDP per capita, temperature at one-month lag and canine rabies outbreak at two-month lag.ConclusionsThe reduction in the overall human rabies incidence was accompanied by a westward and northward expansion of the circulating region in mainland China. Higher risk of human rabies was associated with lower level of education and economic status. New clades of rabies, especial Clade I-G, played an important role in recent spread. Our findings provide valuable information for rabies control and prevention in the future.

Host ecology determines the dispersal patterns of a plant virus.

Since its isolation in 1966 in Kenya, rice yellow mottle virus (RYMV) has been reported throughout Africa resulting in one of the economically most important tropical plant emerging diseases. A thorough understanding of RYMV evolution and dispersal is critical to manage viral spread in tropical areas that heavily rely on agriculture for subsistence. Phylogenetic analyses have suggested a relatively recent expansion, perhaps driven by the intensification of agricultural practices, but this has not yet been examined in a coherent statistical framework. To gain insight into the historical spread of RYMV within Africa rice cultivations, we analyse a dataset of 300 coat protein gene sequences, sampled from East to West Africa over a 46-year period, using Bayesian evolutionary inference. Spatiotemporal reconstructions date the origin of RMYV back to 1852 (1791–1903) and confirm Tanzania as the most likely geographic origin. Following a single long-distance transmission event from East to West Africa, separate viral populations have been maintained for about a century. To identify the factors that shaped the RYMV distribution, we apply a generalised linear model (GLM) extension of discrete phylogenetic diffusion and provide strong support for distances measured on a rice connectivity landscape as the major determinant of RYMV spread. Phylogeographic estimates in continuous space further complement this by demonstrating more pronounced expansion dynamics in West Africa that are consistent with agricultural intensification and extensification. Taken together, our principled phylogeographic inference approach shows for the first time that host ecology dynamics have shaped the historical spread of a plant virus.

Genealogy of the nuclear β-fibrinogen intron 7 in Lissotriton boscai (Caudata, Salamandridae): concordance with mtDNA and implications for phylogeography and speciation

The power of phylogeographic inference resides in its ability to integrate information from multiple sources in an iterative hypothesis- testing framework. In this paper, we build upon previous mtDNA-based hypotheses about the evolutionary history of the Iberian newt Lissotriton boscai using sequences of the highly variable nuclear β-fibrinogen intron 7. In addition to the nuclear sequences, we produced new mtDNA data across the species range to delineate contact zones and test the congruence between nuclear and mitochondrial datasets at the same level of spatial organization. Through a combination of phylogenetic, phylogeographic continuous diffusion, and genetic landscape modelling analyses, we infer the evolutionary history of the species. We found notable congruence between nuclear and mtDNA datasets, which confirms deep and consistent differentiation between two major lineages that originated in the Miocene. Additionally, we found a new nuclear haplogroup with no mtDNA counterpart, roughly circumscribed to the Iberian Sistema Central mountains, and extensive areas of nuclear admixture across mtDNA lineages. We describe potential historical dispersal routes from an ancestral hypothetical refugium in the western end of the Sistema Central in central Portugal and highlight how deep phylogeographic breaks do not necessarily indicate cryptic speciation events.

Chloroplast DNA sequence utility for the lowest phylogenetic and phylogeographic inferences in angiosperms: the tortoise and the hare IV.

Noncoding chloroplast DNA (NC-cpDNA) sequences are the staple data source of low-level phylogeographic and phylogenetic studies of angiosperms. We followed up on previous papers (tortoise and hare II and III) that sought to identify the most consistently variable regions of NC-cpDNA. We used an exhaustive literature review and newly available whole plastome data to assess applicability of previous conclusions at low taxonomic levels. We aligned complete plastomes of 25 species pairs from across angiosperms, comparing the number of genetic differences found in 107 NC-cpDNA regions and matK. We surveyed Web of Science for the plant phylogeographic literature between 2007 and 2013 to assess how NC-cpDNA has been used at the intraspecific level. Several regions are consistently the most variable across angiosperm lineages: ndhF-rpl32, rpl32-trnL((UAG)), ndhC-trnV((UAC)), 5'rps16-trnQ((UUG)), psbE-petL, trnT((GGU))-psbD, petA-psbJ, and rpl16 intron. However, there is no universally best region. The average number of regions applied to low-level studies is ∼2.5, which may be too little to access the full discriminating power of this genome. Plastome sequences have been used successfully at lower and lower taxonomic levels. Our findings corroborate earlier works, suggesting that there are regions that are most likely to be the most variable. However, while NC-cpDNA sequences are commonly used in plant phylogeographic studies, few of the most variable regions are applied in that context. Furthermore, it appears that in most studies too few NC-cpDNAs are used to access the discriminating power of the cpDNA genome.

Understanding the formation of Mediterranean-African-Asian disjunctions: evidence for Miocene climate-driven vicariance and recent long-distance dispersal in the Tertiary relict Smilax aspera (Smilacaceae).

Tethyan plant disjunctions, including Mediterranean-African-Asian disjunctions, are thought to be vicariant, but their temporal origin and underlying causes remain largely unknown. To address this issue, we reconstructed the evolutionary history of Smilax aspera, a hypothesized component of the European Tertiary laurel forest flora. Thirty-eight populations and herbarium specimens representing 57 locations across the species range were sequenced at seven plastid regions and the nuclear ribosomal internal transcribed spacer region. Time-calibrated phylogenetic and phylogeographic inferences were used to trace ancestral areas and biogeographical events. The deep intraspecific split between Mediterranean and African-Asian lineages is attributable to range fragmentation of a southern Tethyan ancestor, as colder and more arid climates developed shortly after the mid-Miocene. In the Mediterranean, climate-induced vicariance has shaped regional population structure since the Late Miocene/Early Pliocene. At around the same time, East African and South Asian lineages split by vicariance, with one shared haplotype reflecting long-distance dispersal. Our results support the idea that geographic range formation and divergence of Tertiary relict species are more or less gradual (mostly vicariant) processes over long time spans, rather than point events in history. They also highlight the importance of the Mediterranean Basin as a centre of intraspecific divergence for Tertiary relict plants.