Lassa virus is genetically diverse with several lineages circulating in West Africa. This study aimed at describing the sequence variability of Lassa virus across Nigeria and inferring its spatiotemporal evolution. We sequenced and isolated 77 Lassa virus strains from 16 Nigerian states. The final data set, including previous works, comprised metadata and sequences of 219 unique strains sampled between 1969 and 2018 in 22 states. Most of this data originated from Lassa fever patients diagnosed at Irrua Specialist Teaching Hospital, Edo State, Nigeria. The majority of sequences clustered with the main Nigerian lineages II and III, while a few sequences formed a new cluster related to Lassa virus strains from Hylomyscus pamfi Within lineages II and III, seven and five sublineages, respectively, were distinguishable. Phylogeographic analysis suggests an origin of lineage II in the southeastern part of the country around Ebonyi State and a main vector of dispersal toward the west across the Niger River, through Anambra, Kogi, Delta, and Edo into Ondo State. The frontline of virus dispersal appears to be in Ondo. Minor vectors are directed northeast toward Taraba and Adamawa and south toward Imo and Rivers. Lineage III might have spread from northern Plateau State into Kaduna, Nasarawa, Federal Capital Territory, and Bauchi. One sublineage moved south and crossed the Benue River into Benue State. This study provides a geographic mapping of lineages and phylogenetic clusters in Nigeria at a higher resolution. In addition, we estimated the direction and time frame of virus dispersal in the country.IMPORTANCE Lassa virus is the causative agent of Lassa fever, a viral hemorrhagic fever with a case fatality rate of approximately 30% in Africa. Previous studies disclosed a geographical pattern in the distribution of Lassa virus strains and a westward movement of the virus across West Africa during evolution. Our study provides a deeper understanding of the geography of genetic lineages and sublineages of the virus in Nigeria. In addition, we modeled how the virus spread in the country. This knowledge allows us to predict into which geographical areas the virus might spread in the future and prioritize areas for Lassa fever surveillance. Our study not only aimed to generate Lassa virus sequences from across Nigeria but also to isolate and conserve the respective viruses for future research. Both isolates and sequences are important for the development and evaluation of medical countermeasures to treat and prevent Lassa fever, such as diagnostics, therapeutics, and vaccines.

Characterizing the evolutionary relationships and population structure of parasites can provide important insights into the epidemiology of human disease.

Climatic fluctuations in the Pleistocene caused glacial expansion-contraction cycles in Eurasia and other parts of the world. Consequences of these cycles, such as population expansion and subsequent subdivision, have been studied in many taxa at intraspecific population level across much of the Northern Hemisphere. However, the consequences for the potential of hybridization and polyploidization are poorly understood. Here, we investigated the phylogeographic structure of two widespread, closely related species, Veronica spicata and Veronica longifolia, across their European distribution ranges. We assessed the extent and the geographic pattern of polyploidization in both species and hybridization between them. We used genome-scale SNP data to clarify phylogenetic relationships and detect possible hybridization/introgression events. In addition, crossing experiments were performed in different combination between V. spicata and V. longifolia individuals of two ploidy levels and of different geographic origins. Finally, we employed ecological niche modeling to infer macroclimatic differences between both species and both ploidy levels. We found a clear genetic structure reflecting the geographical distribution patterns in both species, with V. spicata showing higher genetic differentiation than V. longifolia. We retrieved significant signals of hybridization and introgression in natural populations from the genetic data and corroborated this with crossing experiments. However, there were no clear phylogeographic patterns and unequivocal macroclimatic niche differences between diploid and tetraploid lineages. This favors the hypothesis, that autopolyploidization has happened frequently and in different regions. The crossing experiments produced viable hybrids when the crosses were made between plants of the same ploidy levels but not in the interploidy crosses. The results suggest that hybridization occurs across the overlapping areas of natural distribution ranges of both species, with apparently directional introgression from V. spicata to V. longifolia. Nevertheless, the two species maintain their species-level separation due to their adaptation to different habitats and spatial isolation rather than reproductive isolation.

The cephalopod genus Nautilus is considered a "living fossil" with a contested number of extant and extinct species, and a benthic lifestyle that limits movement of animals between isolated seamounts and landmasses in the Indo-Pacific. Nautiluses are fished for their shells, most heavily in the Philippines, and these fisheries have little monitoring or regulation. Here, we evaluate the hypothesis that multiple species of Nautilus (e.g., N. belauensis, N. repertus and N. stenomphalus) are in fact one species with a diverse phenotypic and geologic range. Using mitochondrial markers, we show that nautiluses from the Philippines, eastern Australia (Great Barrier Reef), Vanuatu, American Samoa, and Fiji fall into distinct geographical clades. For phylogenetic analysis of species complexes across the range of nautilus, we included sequences of Nautilus pompilius and other Nautilus species from GenBank from localities sampled in this study and others. We found that specimens from Western Australia cluster with samples from the Philippines, suggesting that interbreeding may be occurring between those locations, or that there is limited genetic drift due to large effective population sizes. Intriguingly, our data also show that nautilus identified in other studies as N. belauensis, N. stenomphalus, or N. repertus are likely N. pompilius displaying a diversity of morphological characters, suggesting that there is significant phenotypic plasticity within N. pompilius.

The Mississippi Embayment is a prominent physiographic feature of eastern North America consisting of primarily lowland aquatic habitats and a fish fauna that is largely distinct from nearby highland regions. Numerous studies have demonstrated that both pre-Pleistocene and Pleistocene events have had a strong influence on the distributions and relationships of highland fishes in eastern North America. However, the extent to which these same events affected Embayment distributed taxa remains largely unexplored. The purpose of this study was to investigate the relative roles of pre-Pleistocene and Pleistocene events in shaping phylogeographic relationships of four stream dwelling fishes in the Mississippi Embayment. Molecular genetic analyses of the mitochondrial gene cytochrome b were performed for three ictalurid catfish species (Noturus miurus, n = 67; Noturus hildebrandi, n = 93, and Noturus phaeus, n = 44) and one minnow species (Cyprinella camura, n = 78), all distributed in tributary streams of the Mississippi Embayment. Phylogenetic relationships and divergence times among haplotypes for each species were estimated using maximum likelihood and Bayesian methods. Phylogenetic analyses recovered 6 major haplotype clades within N. miurus, 5 within N. hildbrandi, 8 within N. phaeus, and 8 within C. camura. All three Noturus species show a high degree of isolation by drainage, which is less evident in C. camura. A clade of haplotypes from tributaries in the southern portion of the Mississippi Embayment was consistently recovered in all four species. Divergence times among clades spanned the Pleistocene, Pliocene, and Miocene. Novel relationships presented here for C. camura and N. phaeus suggest the potential for cryptic species. Pre-Pleistocene and Pleistocene era sea level fluctuations coincide with some divergence events, but no single event explains any common divergence across all taxa. Like their highland relatives, a combination of both pre-Pleistocene and Pleistocene era events have driven divergences among Embayment lineages.

For a remote oceanic archipelago of up to 8 Myr age, the Azores have a comparatively low level of endemism. We present an analysis of phylogeographic patterns of endemic Azorean island arthropods aimed at testing patterns of diversification in relation to the ontogeny of the archipelago, in order to distinguish between alternative models of evolutionary dynamics on islands. We collected individuals of six species (representing Araneae, Hemiptera and Coleoptera) from 16 forest fragments from 7 islands. Using three mtDNA markers, we analysed the distribution of genetic diversity within and between islands, inferred the differentiation time-frames and investigated the inter-island migration routes and colonization patterns.

Organisms that live at the freshwater surface layer (the neuston) occupy a high energy habitat that is threatened by human activities. Daphniids of the genera Scapholeberis and Megafenestra are adapted to the neuston but are poorly studied for biogeography and diversity. Here we assess the global phylogeography of neustonic daphniids. We obtained 402 new multigene sequences from the 12S rRNA, 16S rRNA, and tRNA (val) regions of the mitochondrial genomes of daphniids from 186 global sites. We assessed the intercontinental origins and boundaries of mitochondrial lineages and the relative rates of evolution in neustonic and planktonic daphniids. We identified 17 divergent lineages in the neustonic daphniids that were associated with biogeographic regions. Six of these lineages had intercontinental ranges - four of these were Transberingian. Patagonian populations of Scapholeberis rammneri were monophyletic and nested within a closely related clade of western North American haplotypes, suggesting an introduction from the Western Nearctic to South America. The Eastern Palearctic was more diverse than other regions, containing eight of the major lineages detected in the Scapholeberinae. The Genus Scapholeberis had high levels of divergence compared to non-neustonic daphniids. Neustonic daphniids have more divergent biogeographic lineages than previously appreciated.

The water buffalo- Bubalus bubalis holds tremendous potential in livestock sector in many Asian countries, particularly India. The origin, domestication and genetic structure of the Indian river buffalo are poorly understood. Therefore, to understand the relationship among the maternal lineages of Indian river buffalo breeds and their domestication process, we analysed mitochondrial D-loop region of 217 animals representing eight breeds from eight different locations in India along with published sequences of Mediterranean buffalo.

The bacterium Burkholderia pseudomallei causes melioidosis, which is mainly associated with tropical areas. We analyzed single-nucleotide polymorphisms (SNPs) among genome sequences from isolates of B. pseudomallei that originated in the Western Hemisphere by comparing them with genome sequences of isolates that originated in the Eastern Hemisphere. Analysis indicated that isolates from the Western Hemisphere form a distinct clade, which supports the hypothesis that these isolates were derived from a constricted seeding event from Africa. Subclades have been resolved that are associated with specific regions within the Western Hemisphere and suggest that isolates might be correlated geographically with cases of melioidosis. One isolate associated with a former World War II prisoner of war was believed to represent illness 62 years after exposure in Southeast Asia. However, analysis suggested the isolate originated in Central or South America.

Ectromelia virus (ECTV), the causative agent of mousepox, has threatened laboratory mouse colonies worldwide for almost a century. Mousepox has been valuable for the understanding of poxvirus pathogenesis and immune evasion. Here, we have monitored in parallel the pathogenesis of nine ECTVs in BALB/cJ mice and report the full-length genome sequence of eight novel ECTV isolates or strains, including the first ECTV isolated from a field mouse, ECTV-MouKre. This approach allowed us to identify several genes, absent in strains attenuated through serial passages in culture, that may play a role in virulence and a set of putative genes that may be involved in enhancing viral growth in vitro. We identified a putative strong inhibitor of the host inflammatory response in ECTV-MouKre, an isolate that did not cause local foot swelling and developed a moderate virulence. Most of the ECTVs, except ECTV-Hampstead, encode a truncated version of the P4c protein that impairs the recruitment of virions into the A-type inclusion bodies, and our data suggest that P4c may play a role in viral dissemination and transmission. This is the first comprehensive report that sheds light into the phylogenetic and geographic relationship of the worldwide outbreak dynamics for the ECTV species.

Puumala virus is an RNA virus hosted by the bank vole (Myodes glareolus) and is today present in most European countries. Whilst it is generally accepted that hantaviruses have been tightly co-evolving with their hosts, Puumala virus (PUUV) evolutionary history is still controversial and so far has not been studied at the whole European level. This study attempts to reconstruct the phylogeographical spread of modern PUUV throughout Europe during the last postglacial period in the light of an upgraded dataset of complete PUUV small (S) segment sequences and by using most recent computational approaches. Taking advantage of the knowledge on the past migrations of its host, we identified at least three potential independent dispersal routes of PUUV during postglacial recolonization of Europe by the bank vole. From the Alpe-Adrian region (Balkan, Austria, and Hungary) to Western European countries (Germany, France, Belgium, and Netherland), and South Scandinavia. From the vicinity of Carpathian Mountains to the Baltic countries and to Poland, Russia, and Finland. The dissemination towards Denmark and North Scandinavia is more hypothetical and probably involved several independent streams from south and north Fennoscandia.

To further probe into whether swamp buffaloes were domesticated once or multiple times in China, this survey examined the mitochondrial DNA (mtDNA) Control Region (D-loop) diversity of 471 individuals representing 22 populations of 455 Chinese swamp buffaloes and 16 river buffaloes. Phylogenetic analysis revealed that Chinese swamp buffaloes could be divided into two distinct lineages, A and B, which were defined previously. Of the two lineages, lineage A was predominant across all populations. For predominant lineage A, Southwestern buffalo populations possess the highest genetic diversity among the three hypothesized domestication centers (Southeastern, Central, and Southwestern China), suggesting Southwestern China as the most likely location for the domestication of lineage A. However, a complex pattern of diversity is detected for the lineage B, preventing the unambiguous pinpointing of the exact place of domestication center and suggesting the presence of a long-term, strong gene flow among swamp buffalo populations caused by extensive migrations of buffaloes and frequent human movements along the Yangtze River throughout history. Our current study suggests that Southwestern China is the most likely domestication center for lineage A, and may have been a primary center of swamp buffalo domestication. More archaeological and genetic evidence is needed to show the process of domestication.

Statistical phylogeography provides useful tools to characterize and quantify the spread of organisms during the course of evolution. Analyzing georeferenced genetic data often relies on the assumption that samples are preferentially collected in densely populated areas of the habitat. Deviation from this assumption negatively impacts the inference of the spatial and demographic dynamics. This issue is pervasive in phylogeography. It affects analyses that approximate the habitat as a set of discrete demes as well as those that treat it as a continuum. The present study introduces a Bayesian modeling approach that explicitly accommodates for spatial sampling strategies. An original inference technique, based on recent advances in statistical computing, is then described that is most suited to modeling data where sequences are preferentially collected at certain locations, independently of the outcome of the evolutionary process. The analysis of georeferenced genetic sequences from the West Nile virus in North America along with simulated data shows how assumptions about spatial sampling may impact our understanding of the forces shaping biodiversity across time and space.

Little is known about the phylogeography of norovirus (NoV) in China. In norovirus, a clear understanding for the characteristics of tree topology, migration patterns and its demographic dynamics in viral circulation are needed to identify its prevalence trends, which can help us better prepare for its epidemics as well as develop useful control strategies. The aim of this study was to explore the genetic diversity, temporal distribution, demographic dynamics and migration patterns of NoV that circulated in China.

Daphnia pulex is quickly becoming an attractive model species in the field of ecological genomics due to the recent release of its complete genome sequence, a wide variety of new genomic resources, and a rich history of ecological data. Sequences of the mitochondrial NADH dehydrogenase subunit 5 and cytochrome c oxidase subunit 1 genes were used to assess the global phylogeography of this species, and to further elucidate its phylogenetic relationship to other members of the Daphnia pulex species complex. Using both newly acquired and previously published data, we analyzed 398 individuals from collections spanning five continents. Eleven strongly supported lineages were found within the D. pulex complex, and one lineage in particular, panarctic D. pulex, has very little phylogeographical structure and a near worldwide distribution. Mismatch distribution, haplotype network, and population genetic analyses are compatible with a North American origin for this lineage and subsequent spatial expansion in the Late Pleistocene. In addition, our analyses suggest that dispersal between North and South America of this and other species in the D. pulex complex has occurred multiple times, and is predominantly from north to south. Our results provide additional support for the evolutionary relationships of the eleven main mitochondrial lineages of the D. pulex complex. We found that the well-studied panarctic D. pulex is present on every continent except Australia and Antarctica. Despite being geographically very widespread, there is a lack of strong regionalism in the mitochondrial genomes of panarctic D. pulex--a pattern that differs from that of most studied cladocerans. Moreover, our analyses suggest recent expansion of the panarctic D. pulex lineage, with some continents sharing haplotypes. The hypothesis that hybrid asexuality has contributed to the recent and unusual geographic success of the panarctic D. pulex lineage warrants further study.

Eight traditional subspecies of tiger (Panthera tigris),of which three recently became extinct, are commonly recognized on the basis of geographic isolation and morphological characteristics. To investigate the species' evolutionary history and to establish objective methods for subspecies recognition, voucher specimens of blood, skin, hair, and/or skin biopsies from 134 tigers with verified geographic origins or heritage across the whole distribution range were examined for three molecular markers: (1) 4.0 kb of mitochondrial DNA (mtDNA) sequence; (2) allele variation in the nuclear major histocompatibility complex class II DRB gene; and (3) composite nuclear microsatellite genotypes based on 30 loci. Relatively low genetic variation with mtDNA,DRB,and microsatellite loci was found, but significant population subdivision was nonetheless apparent among five living subspecies. In addition, a distinct partition of the Indochinese subspecies P. t. corbetti in to northern Indochinese and Malayan Peninsula populations was discovered. Population genetic structure would suggest recognition of six taxonomic units or subspecies: (1) Amur tiger P. t. altaica; (2) northern Indochinese tiger P. t. corbetti; (3) South China tiger P. t. amoyensis; (4) Malayan tiger P. t. jacksoni, named for the tiger conservationist Peter Jackson; (5) Sumatran tiger P. t. sumatrae; and (6) Bengal tiger P. t. tigris. The proposed South China tiger lineage is tentative due to limited sampling. The age of the most recent common ancestor for tiger mtDNA was estimated to be 72,000-108,000 y, relatively younger than some other Panthera species. A combination of population expansions, reduced gene flow, and genetic drift following the last genetic diminution, and the recent anthropogenic range contraction, have led to the distinct genetic partitions. These results provide an explicit basis for subspecies recognition and will lead to the improved management and conservation of these recently isolated but distinct geographic populations of tigers.

Tardigrades constitute a micrometazoan phylum usually considered as taxonomically challenging and therefore difficult for biogeographic analyses. The genus Pseudechiniscus, the second most speciose member of the family Echiniscidae, is commonly regarded as a particularly difficult taxon for studying due to its rarity and homogenous sculpturing of the dorsal plates. Recently, wide geographic ranges for some representatives of this genus and a new hypothesis on the subgeneric classification have been suggested. In order to test these hypotheses, we sequenced 65 Pseudechiniscus populations extracted from samples collected in 19 countries distributed on 5 continents, representing the Neotropical, Afrotropical, Holarctic, and Oriental realms. The deep subdivision of the genus into the cosmopolitan suillus-facettalis clade and the mostly tropical-Gondwanan novaezeelandiae clade is demonstrated. Meridioniscus subgen. nov. is erected to accommodate the species belonging to the novaezeelandiae lineage characterised by dactyloid cephalic papillae that are typical for the great majority of echiniscids (in contrast to pseudohemispherical papillae in the suillus-facettalis clade, corresponding to the subgenus Pseudechiniscus). Moreover, the evolution of morphological traits (striae between dorsal pillars, projections on the pseudosegmental plate IV', ventral sculpturing pattern) crucial in the Pseudechiniscus taxonomy is reconstructed. Furthermore, broad distributions are emphasised as characteristic of some taxa. Finally, the Malay Archipelago and Indochina are argued to be the place of origin and extensive radiation of Pseudechiniscus.

Rhizophora is a key genus for revealing the formation process of the pantropical distribution of mangroves. In this study, in order to fully understand the historical scenario of Rhizophora that achieved pantropical distribution, we conducted phylogeographic analyses based on nucleotide sequences of chloroplast and nuclear DNA as well as microsatellites for samples collected worldwide. Phylogenetic trees suggested the monophyly of each AEP and IWP lineages respectively except for R. samoensis and R. × selala. The divergence time between the two lineages was 10.6 million years ago on a dated phylogeny, and biogeographic stochastic mapping analyses supported these lineages separated following a vicariant event. These data suggested that the closure of the Tethys Seaway and the reduction in mangrove distribution followed by Mid-Miocene cooling were key factors that caused the linage diversification. Phylogeographic analyses also suggested the formation of the distinctive genetic structure at the AEP region across the American continents around Pliocene. Furthermore, long-distance trans-pacific dispersal occurred from the Pacific coast of American continents to the South Pacific and formed F1 hybrid, resulting in gene exchange between the IWP and AEP lineages after 11 million years of isolation. Considering the phylogeny and phylogeography with divergence time, a comprehensive picture of the historical scenario behind the pantropical distribution of Rhizophora is updated.

Rabies is a major fatal zoonotic disease in Indonesia. This study was conducted to determine the recent dynamics of rabies virus (RABV) in various areas and animal species throughout Indonesia. A total of 27 brain samples collected from rabid animals of various species in Bali, Sumatra, Kalimantan, Sulawesi, Java, and Flores in 2008 to 2010 were investigated. The cDNA of the nucleoprotein gene from each sample was generated and amplified by one-step reverse transcription-PCR, after which the products were sequenced and analyzed. The symmetric substitution model of a Bayesian stochastic search variable selection extension of the discrete phylogeographic model of the social network was applied in BEAST ver. 1.7.5 software. The spatial dispersal was visualized in Cartographica using Spatial Phylogenetic Reconstruction of Evolutionary Dynamics. We demonstrated inter-island introduction and reintroduction, and dog was found to be the only source of infection of other animals. Ancestors of Indonesian RABVs originated in Java and its descendants were transmitted to Kalimantan, then further to Sumatra, Flores, and Bali. The Flores descendent was subsequently transmitted to Sulawesi and back to Kalimantan. The viruses found in various animal species were transmitted by the dog.

Papaya ringspot virus (PRSV) is the most important virus affecting papaya and cucurbit plants in tropical and subtropical areas. PRSV isolates are divided into biotypes P and W: both the P and W types naturally infect plants in the family Cucurbitaceae, whereas the P type naturally infects papaya (Carica papaya). Understanding the origin and nature of the PRSV genetic diversity and evolution is critical for the implementation of control strategies based on cross-protection and the deployment of transgenic plants that show resistance to virus isolates highly similar to the transgene. The molecular epidemiology of PRSV was evaluated by analyzing the nucleotide sequence of the capsid protein (CP) and helper component-proteinase (HC-Pro) genes of isolates from around the world, including newly characterized ones from Colombia and Venezuela, using a relaxed molecular clock-based approach and a phylogeographic study. Our results confirm previous estimates on the origin of PRSV around 400 years ago and suggest distinct dispersion events from the Indian Peninsula to the rest of Asia, via Thailand, and subsequently to the Americas. A historical reconstruction of the P- and W-type characters in the phylogenetic study supports the need to revise the hypothesis that PRSV-P derives from PRSV-W since our results suggest that the ancestral state could be either of the two biotypes. Moreover, estimates of epidemic growth predict an increasing genetic diversity of the virus over time that has direct implications for control strategies of PRSV based on cross-protection and the use of transgenic plants.