Search

Hosting 1460 plant and 126 vertebrate endemic species, the Great Escarpment (hereafter, Escarpment) forms a semi‐circular "amphitheater" of mountains girdling southern Africa from arid west to temperate east. Since arid and temperate biota are usually studied separately, earlier studies overlooked the biogeographical importance of the Escarpment as a whole. Bats disperse more widely than other mammalian taxa, with related species and intraspecific lineages occupying both arid and temperate highlands of the Escarpment, providing an excellent model to address this knowledge gap. We investigated patterns of speciation and micro‐endemism from modeled past, present, and future distributions in six clades of southern African bats from three families (Rhinolophidae, Cistugidae, and Vespertilionidae) having different crown ages (Pleistocene to Miocene) and biome affiliations (temperate to arid). We estimated mtDNA relaxed clock dates of key divergence events across the six clades in relation both to biogeographical features and patterns of phenotypic variation in crania, bacula and echolocation calls. In horseshoe bats (Rhinolophidae), both the western and eastern "arms" of the Escarpment have facilitated dispersals from the Afrotropics into southern Africa. Pleistocene and pre‐Pleistocene "species pumps" and temperate refugia explained observed patterns of speciation, intraspecific divergence and, in two cases, mtDNA introgression. The Maloti‐Drakensberg is a center of micro‐endemism for bats, housing three newly described or undescribed species. Vicariance across biogeographic barriers gave rise to 29 micro‐endemic species and intraspecific lineages whose distributions were congruent with those identified in other phytogeographic and zoogeographic studies. Although Köppen–Geiger climate models predict a widespread replacement of current temperate ecosystems in southern Africa by tropical or arid ecosystems by 2070–2100, future climate Maxent models for 13 bat species (all but one of those analyzed above) showed minimal range changes in temperate species from the eastern Escarpment by 2070, possibly due to the buffering effect of mountains to climate change. [ABSTRACT FROM AUTHOR]

The Black-throated Gray Warbler (Setophaga nigrescens) is a migratory songbird with a widespread and continuous distribution in the North American West. Previously published morphological data support the designation of 2 putative subspecies based on subtle differences in plumage and song; however, whether these patterns are clinal or discrete has yet to be determined. Furthermore, given that patterns of morphological variation occur across multiple distinct habitats and geographical boundaries, underlying population genetic structure may be present. To investigate this question and explore the evolutionary history of S. nigrescens, we conducted the first range-wide phylogeographic analysis of this species using genetic (mitochondrial DNA [mtDNA], ND2 gene sequences), morphological (body measurements), and environmental (niche modeling) data. Overall, our results identified 2 major mtDNA clades. One clade occurs east of the Cascade-Sierra ranges throughout the southwestern, Great Basin, and Rocky Mountain regions in North America (Interior clade), while a second widespread clade is comprised of birds that span the Pacific Northwest, and also across the Great Basin and southwest into the Rocky Mountains (Coastal clade). Divergence time estimates suggest a Pleistocene epoch vicariance event. Environmental niche models predict habitat modification and fragmentation from glacial–interglacial cycles, which were frequent over the Pleistocene. Range expansion northwards following the Last Glacial Maximum was identified from neutrality tests and shifts in suitable habitat from the Mid-Holocene to present-day conditions, where both clades are likely experiencing secondary contact. We found support for nonequivalent environmental niches among both clades based on a spatial analysis of environmental variables, with the Coastal clade showing greater tolerance to both temperate and xeric climatic conditions across the North American West; however, the environmental space occupied by both clades was more similar to each other than expected by chance. High levels of geographic overlap and a lack of discrete patterns in morphology suggest few contemporary reproductive barriers to gene flow. [ABSTRACT FROM AUTHOR]

The bark beetle Ips typographus (L.) is the most important pest in European Norway spruce forests. Especially after abiotic disturbance events, like storms or snowfall, this beetle can utilise vast amounts of trees with impaired defences, leading to mass outbreaks and further mortality of spruce. In 2018, the 'Vaia' storm damaged huge areas of Norway spruce-dominated forests in the Southern Alps. In the following years, the same regions were additionally affected by heavy snowfalls providing high amounts of suitable material for I. typographus brood establishment and development. Since then, these forests have faced unprecedented population outbreaks of I. typographus with severe effects for forest ecosystems, from both an ecological and an economic perspective. A comprehensive understanding of this beetle's genetic structure, including effects of the last ice ages as well as consequences of current population outbreaks, can give insights into the population dynamics of an eruptive species. Here, we studied the genetic structure of five epidemic populations from the Southern Alps and one from a potential glacial refugium in the Apennines. Our results revealed a low genetic diversity and a weak population structure among all studied locations. In the six analysed populations, four already described haplotypes, belonging to three slightly divergent haplogroups, were found. These haplotypes were evenly distributed among populations, including the potential glacial refugium in the Apennines. High gene flow among populations is probably maintained by the beetle's high dispersal capacity and high population densities, resulting in a panmictic population structure. [ABSTRACT FROM AUTHOR]

Life history characteristics of cryptobenthic reef fishes (CRFs), such as their fast growth and reproductive rates, near-shore larval retention and high turnover, predispose these species to rapid diversification and cryptic speciation. This may be particularly true for isolated CRF populations. The sailfin blenny, Emblemaria pandionis, is widespread throughout the Caribbean and is found on the reefs of the Flower Garden Banks, an isolated reef system in the NW Gulf of Mexico. Using DNA barcoding and multi-locus delimitation, we show that E. pandionis is a species complex comprised of at least four distinct taxonomic units in the Caribbean, one at the Flower Garden Banks, a second in eastern Florida, a third in the central Caribbean and a fourth in Curaçao. The phylogenetic relationships within the E. pandionis species complex are defined by well-established phylogeographic barriers to gene flow in the Caribbean. Populations of E. pandionis from eastern Florida are genetically distinct from populations found in the Bahamas, highlighting the role of the Florida Strait as a strong barrier to gene flow and populations in Curaçao are deeply divergent from all other populations sampled in the Caribbean. Results from this study highlight the roles of isolated reefs and ocean currents in the speciation of cryptobenthic reef fishes. [ABSTRACT FROM AUTHOR]

Context: Striostrea prismatica is a valuable native oyster found along the Eastern Pacific coast from Mexico to Peru. There is evidence of overexploitation of this natural resource. Thus, population genetic data can help in designing management strategies. Aims: This study aimed to assess the genetic diversity, population structure and demographic history of S. prismatica using nuclear and mitochondrial DNA markers. Methods: Oysters were collected from seven natural fishing banks in two biogeographic realms of the eastern Pacific Ocean. Inter-simple sequence repeat markers and partial sequences of the mitochondrial cytochrome c oxidase subunit I (COI), 12S rRNA (12S) and 16S rRNA (16S) genes were amplified by polymerase chain reaction and analysed. Key results: Nuclear data revealed high intraspecific genetic diversity, moderate to strong genetic differentiation among fishing banks and a subtle genetic structure among marine ecoregions. The analyses supported strong isolation by distance, but chaotic genetic patchiness was also detected. Mitochondrial data showed high haplotype diversity and evidence of a population expansion c. 25,000 years ago. Conclusions: Dispersal capacity, oceanographic features and historical factors can explain the levels of genetic diversity and the pattern of population differentiation and connectivity observed in S. prismatica. Implications: The information generated in this study will support future management plans for this fishing resource. This study aimed to assess the genetic diversity, population structure and demographic history of Striostrea prismatica , an ecologically and commercially important native oyster in Mexico, Ecuador and Peru, using nuclear and mitochondrial DNA markers. Nuclear data revealed high intraspecific genetic diversity and strong genetic differentiation. Mitochondrial data showed high haplotypic diversity, no genetic structure and a population expansion c. 25,000 years ago. [ABSTRACT FROM AUTHOR]

Despite the current disjoint distribution of Amazonian and Atlantic forests, evidence suggests past historical connections. Here we investigated the historical connections between three didelphid mammal species from Amazonian and Atlantic forests (Caluromys philander, Marmosa demerarae, and M. murina) using comparative phylogeography and paleodistribution models. We generated ecological niche models from the present until the Pliocene to evaluate changes in distributions over time, focusing on the previously suggested connection routes. We inferred divergence times between populations of each species using the mitochondrial Cytochrome b, estimated the phylogeographic relationships with haplotype networks, and calculated the genetic distances. All species exhibited populations that were separated between the Amazonian and Atlantic forests around 1–2 million years ago, while some populations showed more recent divergences. Paleodistributions were more extensive for periods predating the estimated times of divergence for all species, becoming narrower after this period. We suggest that the climatic oscillations during the onset of the Pleistocene largely influenced the phylogenetic structuring of these forest-dwelling species. This led to the separation of populations currently restricted to Amazonian or Atlantic forests. Our findings also point to repeated connections over time, with both North- and South-Eastern routes compatible with ancient connections. This represents the first simultaneous evaluation of past connections between Amazonian and Atlantic forests combining phylogeographic inferences with paleodistribution models for didelphid marsupials. We suggested the relation of the responses to past climate change, such as increases in the paleodistribution, and species resilience. [ABSTRACT FROM AUTHOR]

Iriomote Island is a small continental island hosting approximately 40 rivers characterized by short, steep-gradient streams and small waterfalls. We examined the population genetic structure and diversity of the freshwater prawn Macrobrachium shokitai Fujino & Baba, 1973 (Palaemonidae) using the mitochondrial cytochrome c oxidase subunit I (COI) gene sequence. Specimens were collected from 19 sites along ten rivers on the island, covering the known geographic distribution of the species. Haplotype and nucleotide diversities were very low within each river on the island. Nearly all haplotypes were exclusive to a given river, confirming distinct population structures among rivers at a small geographical scale. Slight genetic differentiation was discerned between the upper and lower sites of waterfalls in only the Yuchin River. Phylogenetic tree reconstruction, haplotype network, and principle coordinate analysis classified five clades roughly correspond to geographical groups in the western, southern, central, and northern regions of the island. Extreme caution should therefore be taken when translocating individuals to other rivers to avoid introducing genetic disturbance, even during conservation and mitigation efforts. [ABSTRACT FROM AUTHOR]

Determining the genetic diversity and source rookeries of sea turtles collected from feeding grounds can facilitate effective conservation initiatives. To ascertain the genetic composition and source rookery, we examined a partial sequence of the mitochondrial control region (CR, 796 bp) of 40 green turtles (Chelonia mydas) collected from feeding grounds around the Korean Peninsula between 2014 and 2022. We conducted genetic and mixed-stock analyses (MSA) and identified 10 CR haplotypes previously reported in Japanese populations. In the haplotype network, six, three, and one haplotype(s) grouped with the Japan, Indo-Pacific, and Central South Pacific clades, respectively. The primary rookeries of the green turtles were two distantly remote sites, Ogasawara (OGA) and Central Ryukyu Island (CRI), approximately 1,300 km apart from each other. Comparing three parameters (season, maturity, and specific feeding ground), we noted that OGA was mainly associated with summer and the Jeju Sea, whereas CRI was with fall and the East (Japan) Sea ground. The maturity did not show a distinct pattern. Our results indicate that green turtles in the feeding grounds around the Korean Peninsula originate mainly from the Japan MU and have genetic origins in the Japan, Indo-Pacific, and Central South Pacific clades. Our results provide crucial insights into rookeries and MUs, which are the focus of conservation efforts of the Republic of Korea and potential parties to collaborate for green turtle conservation. [ABSTRACT FROM AUTHOR]

The biota of Sulawesi is noted for its high degree of endemism and for its substantial levels of in situ biological diversification. While the island's long period of isolation and dynamic tectonic history have been implicated as drivers of the regional diversification, this has rarely been tested in the context of an explicit geological framework. Here, we provide a tectonically informed biogeographical framework that we use to explore the diversification history of Sulawesi flying lizards (the Draco lineatus Group), a radiation that is endemic to Sulawesi and its surrounding islands. We employ a framework for inferring cryptic speciation that involves phylogeographic and genetic clustering analyses as a means of identifying potential species followed by population demographic assessment of divergence-timing and rates of bi-directional migration as means of confirming lineage independence (and thus species status). Using this approach, phylogenetic and population genetic analyses of mitochondrial sequence data obtained for 613 samples, a 50-SNP data set for 370 samples, and a 1249-locus exon-capture data set for 106 samples indicate that the current taxonomy substantially understates the true number of Sulawesi Draco species, that both cryptic and arrested speciations have taken place, and that ancient hybridization confounds phylogenetic analyses that do not explicitly account for reticulation. The Draco lineatus Group appears to comprise 15 species-9 on Sulawesi proper and 6 on peripheral islands. The common ancestor of this group colonized Sulawesi ~11 Ma when proto-Sulawesi was likely composed of two ancestral islands, and began to radiate ~6 Ma as new islands formed and were colonized via overwater dispersal. The enlargement and amalgamation of many of these proto-islands into modern Sulawesi, especially during the past 3 Ma, set in motion dynamic species interactions as once-isolated lineages came into secondary contact, some of which resulted in lineage merger, and others surviving to the present. [Genomics; Indonesia; introgression; mitochondria; phylogenetics; phylogeography; population genetics; reptiles.].

Currents are unique drivers of oceanic phylogeography and thus determine the distribution of marine coastal species, along with past glaciations and sea-level changes. Here we reconstruct the worldwide colonization history of eelgrass (Zostera marina L.), the most widely distributed marine flowering plant or seagrass from its origin in the Northwest Pacific, based on nuclear and chloroplast genomes. We identified two divergent Pacific clades with evidence for admixture along the East Pacific coast. Two west-to-east (trans-Pacific) colonization events support the key role of the North Pacific Current. Time-calibrated nuclear and chloroplast phylogenies yielded concordant estimates of the arrival of Z. marina in the Atlantic through the Canadian Arctic, suggesting that eelgrass-based ecosystems, hotspots of biodiversity and carbon sequestration, have only been present there for ~243 ky (thousand years). Mediterranean populations were founded ~44 kya, while extant distributions along western and eastern Atlantic shores were founded at the end of the Last Glacial Maximum (~19 kya), with at least one major refuge being the North Carolina region. The recent colonization and five- to sevenfold lower genomic diversity of the Atlantic compared to the Pacific populations raises concern and opportunity about how Atlantic eelgrass might respond to rapidly warming coastal oceans.

ABSTRACT Giant clams (Tridacna and Hippopus) are large marine bivalves occupying tropical and subtropical reefs in the Indo‐Pacific. Giant clam populations have declined in many areas of the Indo‐Pacific and continue to be threatened by harvesting and environmental change. The small giant clam (Tridacna maxima) occurs throughout the Indo‐Pacific and has been subject to several phylogeographic studies across its range. However, given its broad range, there are several areas where the genetic diversity and connectivity of T. maxima populations has not been characterised. Here, we analyse the mitochondrial marker cytochrome oxidase 1 (CO1) to examine the genetic diversity and connectivity of T. maxima in two regions: Australia's Coral Sea Marine Park and the Cook Islands. Samples were collected from 13 reefs within the Coral Sea Marine Park and ten islands within the Cook Islands archipelago. Tridacna maxima across the sampled region of the Coral Sea did not display any population structure, whereas significant population structure was detected for T. maxima within the Cook Islands. For the Cook Islands, most pairwise comparisons involving an island in the northern group (Manihiki) were significant, as were comparisons for Palmerston (a more centrally located island) and the southern islands, Rarotonga and Mangaia. Both regions displayed high haplotype diversities (> 0.90), indicating that they are important repositories of genetic diversity. Additional CO1 data from throughout T. maxima's distribution showed that the Coral Sea clams belonged to the clade occurring in the South‐Western Pacific Ocean, whilst those from the Cook Islands belonged to a unique clade found in the Central Pacific Ocean. This clade extended from Fiji in the west to French Polynesia in the east and the atolls of Palmyra and Tarawa (Kiribati) in the north. Our assessment of genetic diversity and population structure in these regions will assist with management decisions for the species.

IntroductionPorcine circovirus 2 (PCV2) emerged more than three decades ago as one of the most impactful pathogens in the swine industry. Despite being a DNA virus, one of the hallmarks of PCV2 is its high evolutionary rate, which has led to the emergence of different genotypes, each exhibiting varying degrees of evolutionary success. Current knowledge suggests the occurrence of three main waves of genotype dominance, alternating over time (i.e., PCV2a, PCV2b, and PCV2d), alongside less prevalent genotypes. However, although PCV2d is currently the most common genotype nowadays, the others continue being circulating in the pig population.MethodsThe present study reconsidered the epidemiological and evolutionary patterns of PCV2 genotypes using phylodynamic analyses, benefiting from an almost 10-fold increase in ORF2 sequence availability compared to previous studies. Additionally, a phylogeographic analysis was performed to investigate viral dispersal patterns and frequency, and the selective pressures acting on the capsid protein were estimated and compared among genotypes.ResultsWhile successive emergence of major genotypes was confirmed, this study extends previous findings by revealing subsequent prevalence fluctuations of PCV2a and PCV2b after their initial decline. This evolutionary process may represent an example of balancing selection, specifically negative frequency-dependent selection, where a genotype fitness diminishes as it becomes more common, and vice versa. Variations in genotype- or clade-specific immunity—affected by the local prevalence of viral groups—combined with the periodic introduction of strains that have independently evolved in different regions, may have led to fluctuations in the population dynamics of major genotypes over time. These fluctuations were associated with ongoing evolution and variations in the capsid amino acid profile.DiscussionThese findings have profound implications for future control strategies. Although PCV2d remains the most prevalent and widespread genotype, other genotypes should not be neglected. Control strategies should thus target the entire PCV2 population, with a focus on fostering broader and more cross-protective immunity.

Abstract Identification of the geographic origin of invasive species can be critical to effective management and amelioration of negative impacts in the introduced range. Liriomyza huidobrensis is a polyphagous leafmining fly that is a devastating pest of many vegetable and floriculture crops around the world. Considered native to South and possibly Central America, L. huidobrensis became invasive in the 1980s and has since spread to at least 30 countries on five continents. We used phylogeographic analysis of over 2 kb of mitochondrial cytochrome oxidase I and II sequence data from 403 field‐collected specimens from both native and introduced populations to investigate the geographic origins of invasive L. huidobrensis worldwide. Within South America, there was substantial genetic variation, as well as the strong phylogeographic structure typical of a native range. In contrast, leafminers from the introduced range and Central America all contained little genetic variation and shared the same small set of haplotypes. These haplotypes trace to Peru as the ultimate geographic origin of invasive populations. Central America is rejected as part of the original geographic range of L. huidobrensis. Within Peru, the primary export region of Lima shared an extremely similar pattern of reduced haplotype variation to the invasive populations. An additional 18 specimens collected at US ports of entry did not share the same haplotype profile as contemporary invasive populations, raising perplexing questions on global pathways and establishment success in this species.

ABSTRACT Bartonella grahamii is one of the most prevalent Bartonella species in wild rodents and has been associated with human cases of neuroretinitis. The structure and distribution of genomic diversity in natural B. grahamii is largely unexplored. Here, we have applied a comprehensive population genomic and phylogenomic analysis to 172 strains of B. grahamii to unravel the genetic differences and influencing factors that shape its populations. The findings reveal a remarkable genomic diversity within the species, primarily in the form of single- nucleotide polymorphisms. The open pangenome of B. grahamii indicates a dynamic genomic evolution influenced by its ecological niche. Whole-genome data allowed us to decompose B. grahamii diversity into six phylogroups, each characterized by a unique “mosaic pattern” of hosts and biogeographic regions. This suggests a complex interplay between host specificity and biogeography. In addition, our study suggests a possible origin of European strains from Asian lineages, and host factors have a more significant impact on the genetic differentiation of B. grahamii than geographical factors. These insights contribute to understanding the evolutionary history of this pathogen and provide a foundation for future epidemiological research and public health strategies.IMPORTANCEBartonella grahamii has been reported worldwide and shown to infect humans. Up to now, an effective transmission route of B. grahamii to humans has not been confirmed. The genetic evolution of B. grahamii and the relationship between B. grahamii and its host need to be further studied. The factors driving the genetic diversity of B. grahamii are still controversial. The results showed that the European isolates shared a common ancestor with the Chinese isolates. Host factors were shown to play an important role in driving the genetic diversity of B. grahamii. When host factors were fixed, geographic barriers drove B. grahamii microevolution. Our study emphasizes the importance of characterizing isolate genomes derived from hosts and geographical locations and provides a new reference for the origin of B. grahamii.

The prevalence of avian-derived Escherichia coli (E. coli) carrying mcr-1 poses a significant threat to the development of the poultry industry and public health safety. Despite ongoing in-depth epidemiological research worldwide, a comprehensive macroscopic study based on genomics is still lacking. In response, this study collected 1104 genomic sequences of avian-derived mcr-1-positive E. coli (MCRPEC) from the NCBI public database, covering 31 countries. The majority of sequences originated from China (48.82 %), followed by the Netherlands (10.41 %). In terms of avian hosts, chicken accounted for the largest proportion (44.11 %), followed by gallus (24.09 %). Avian-derived MCRPEC also serves as a reservoir for other antibiotic resistance genes (ARGs), with 179 ARGs coexisting with mcr-1 identified. A total of 206 virulence-associated genes were also identified, revealing the pathogenic risks of MCRPEC. Pan-genome analysis revealed that avian-derived MCRPEC from different hosts, countries of origin, and serotypes exhibit minor SNP differences, indicating a high risk of cross-regional and cross-host transmission. The ST types of MCRPRC are diverse, with ST10 being the most prevalent (n=70). Spearman analysis showed a significant correlation between the number of ARGs and the insertion sequences (ISs) as well as plasmid replicon in ST10 strains. Furthermore, ST10 strains share a similar genetic basis with human-derived MCRPEC, suggesting the possibility of clonal dissemination. Pan-genome-wide association studies (pan-GWAS) indicated that the differential genes of MCRPEC from different countries and host sources are significantly different, mainly related to genes encoding type IV secretion systems and mobile genetic elements (MGEs). Plasmid mapping of showed that the prevalent plasmid types vary by country and host, with IncI2 and IncX4 being the main mcr-1-positive plasmids. Among the 12 identified mcr-1 genetic contexts with ISs, the Tn6330 transposon was the predominant carrier of mcr-1. In summary, the potential threat of avian-derived MCRPEC cannot be ignored, and long-term and comprehensive monitoring are essential.

The puma (Puma concolor Linnaeus, 1771) is the top predator with the widest distribution in America. Since the establishment of European settlers on the American continent, puma populations have experienced significant contractions and reductions in their original distribution. In Argentina, the management of the conflict between humans and pumas (direct persecution and habitat modification) focused on reduction or elimination methods, leading to a drastic contraction, even total eradication, of puma populations as seen in Patagonia and the eastern part of the country. Despite the lack of knowledge about puma population demographic trends, there are taxonomic issues that remain controversial and need to be resolved to implement appropriate management and conservation measures. Therefore, the aim of this study was to genetically characterize puma populations in the central-southern region of Argentina using two mitochondrial markers, evaluate their demographic history, compare our results at a macro-regional level, and discuss our findings in a conservation and management context. A total of 203 individuals were used, and a fragment of the control region and another of ND5 were sequenced. The genetic variability obtained was moderate. Substitution rates for each locus and the tMRCA were calculated from calibrated trees. In a concatenated tree, two main phylogenetic clades were identified (posterior probability = 1), although a reciprocal monophyly was not observed, with a divergence time of 228 thousand years and a 95% CI [117–363 thousand years]. When evaluating population structuring, three genetic clusters were found, one characteristic of the Patagonian region and the others in the central part of the country. Calculating the ФST values for pairs resulted in significant structuring between Patagonia and the rest of the populations, suggesting the arid diagonal as a possible barrier to gene flow. When evaluating the demographic history, neutrality tests would support a recent expansion in Patagonia. These findings are crucial in defining two distinct Management Units (MUs) in the southern part of puma distribution and providing valuable information for management and conservation measures for the species.

Researchers studying the evolution of viral pathogens and other organisms increasingly encounter and use large and complex data sets from multiple different sources. Statistical research in Bayesian phylogenetics has risen to this challenge. Researchers use phylogenetics not only to reconstruct the evolutionary history of a group of organisms, but also to understand the processes that guide its evolution and spread through space and time. To this end, it is now the norm to integrate numerous sources of data. For example, epidemiologists studying the spread of a virus through a region incorporate data including genetic sequences (e.g. DNA), time, location (both continuous and discrete) and environmental covariates (e.g. social connectivity between regions) into a coherent statistical model. Evolutionary biologists routinely do the same with genetic sequences, location, time, fossil and modern phenotypes, and ecological covariates. These complex, hierarchical models readily accommodate both discrete and continuous data and have enormous combined discrete/continuous parameter spaces including, at a minimum, phylogenetic tree topologies and branch lengths. The increased size and complexity of these statistical models have spurred advances in computational methods to make them tractable. We discuss both the modeling and computational advances below, as well as unsolved problems and areas of active research.

The objective is to determine the complete nucleotide sequence and conduct a phylogenetic analysis of genome variants of the Puumala virus isolated in the Saratov region. Materials and methods. The samples for the study were field material collected in the Gagarinsky (formerly Saratovsky), Engelssky, Novoburassky and Khvalynsky districts of the Saratov region in the period from 2019 to 2022. To specifically enrich the Puumala virus genome in the samples, were used PCR and developed a specific primer panel. Next, the resulting PCR products were sequenced and the fragments were assembled into one sequence for each segment of the virus genome. To construct phylogenetic trees, the maximum parsimony algorithm was used. Results. Genetic variants of the Puumala virus isolated in the Saratov region have a high degree of genome similarity to each other, which indicates their unity of origin. According to phylogenetic analysis, they all form a separate branch in the cluster formed by hantaviruses from other subjects of the Volga Federal District. The virus variants from the Republics of Udmurtia and Tatarstan, as well as from the Samara and Ulyanovsk regions, are closest to the samples from the Saratov region. Conclusion. The data obtained show the presence of a pronounced territorial confinement of strains to certain regions or areas that are the natural biotopes of their carriers. This makes it possible to fairly accurately determine the territory of possible infection of patients and/or the circulation of carriers of these virus variants based on the sequence of individual segments of their genome.

Abstract Background There is no consensus as to the origin of the domestic yak (Bos grunniens). Previous studies on yak mitochondria mainly focused on mitochondrial displacement loop (D-loop), a region with low phylogenetic resolution. Here, we analyzed the entire mitochondrial genomes of 509 yaks to obtain greater phylogenetic resolution and a comprehensive picture of geographical diversity. Results A total of 278 haplotypes were defined in 509 yaks from 21 yak breeds. Among them, 28 haplotypes were shared by different varieties, and 250 haplotypes were unique to specific varieties. The overall haplotype diversity and nucleotide diversity of yak were 0.979 ± 0.0039 and 0.00237 ± 0.00076, respectively. Phylogenetic tree and network analysis showed that yak had three highly differentiated genetic branches with high support rate. The differentiation time of clades I and II were about 0.4328 Ma, and the differentiation time of clades (I and II) and III were 0.5654 Ma. Yushu yak is shared by all haplogroups. Most (94.70%) of the genetic variation occurred within populations, and only 5.30% of the genetic variation occurred between populations. The classification showed that yaks and wild yaks were first clustered together, and yaks were clustered with American bison as a whole. Altitude had the highest impact on the distribution of yaks. Conclusions Yaks have high genetic diversity and yak populations have experienced population expansion and lack obvious phylogeographic structure. During the glacial period, yaks had at least three or more glacial refugia.

Background. The study of the phylogeography of the tularemia microbe is not only practical, but also of scientific interest, as it allows us to evaluate the evolutionary processes of Francisella tularensis in modern conditions.Aim. Phylogenetic analysis of strains of the causative agent of tularemia in the Rostov region and identification of erythromycin-sensitive genovariants.Results. A representative collection of 122 cultures, isolated from various sources from 1945 to 2022, was studied. As it was established, all the studied strains of the causative agent of tularemia had typical cultural, morphological, and biological characteristics. The results of the study of bacteria antibiograms showed that all strains were resistant to β-lactams, polymyxin, clindamycin, and sensitive to aminoglycosides, rifampicin, and fluoroquinolones. In relation to erythromycin, only 2 erythromycin-sensitive strains were found. Phylogenetic analysis of 34 European and Russian strains of F. tularensis made it possible to determine whether these strains belong to subgroup B.7 of the main genetic group B.6.Conclusion. For the first time, in the Tselinsky district of the Rostov region, the presence of «atypical» genovariants of the tularemia microbe belonging to subgroup B.7 of the main group B.6 of the subspecies holarctica was registered, according to the modern scheme of genetic typing of F. tularensis subsp. holarctica.

To address the phylogenetic and phylogeographic relationship between different lineages of Anopheles (An.) subpictus species complex in most parts of the Asian continent by maximum utilization of Internal Transcriber Spacer 2 (ITS2) and cytochrome C oxidase I (COI) sequences deposited at the GenBank. Seventy-five ITS2, 210 COI and 26 concatenated sequences available in the NCBI database were used. Phylogenetic analysis was performed using Bayesian likelihood trees, whereas median-joining haplotype networks and time-scale divergence trees were generated for phylogeographic analysis. Genetic diversity indices and genetic differentiation were also calculated. Two genetically divergent molecular forms of An. subpictus species complex corresponding to sibling species A and B are established. Species A evolved around 37-82 million years ago in Sri Lanka, India, and the Netherlands, and species B evolved around 22-79 million years ago in Sri Lanka, India, and Myanmar. Vietnam, Thailand, and Cambodia have two molecular forms: one is phylogenetically similar to species B. Other forms differ from species A and B and evolved recently in the above mentioned countries, Indonesia and the Philippines. Genetic subdivision among Sri Lanka, India, and the Netherlands is almost absent. A substantial genetic differentiation was obtained for some populations due to isolation by large geographical distances. Genetic diversity indices reveal the presence of a long-established stable mosquito population, at mutation-drift equilibrium, regardless of population fluctuations. An. subpictus species complex consists of more than two genetically divergent molecular forms. Species A is highly divergent from the rest. Sri Lanka and India contain only species A and B.

Abstract Background The eastern edge of the Qinghai‒Tibet Plateau (QTP) and subtropical China have various regions where plant species originate and thrive, but these regions have been the focus of very few integrative studies. Here, we elucidated the phylogeographic structure of a continuous and widespread Akebia trifoliata population across these two regions. Results Sixty-one populations consisting of 391 genotypes were examined to assess population diversity and structure via network distribution analysis, maximum likelihood phylogenetic tree reconstruction, divergence time estimation, demographic history inference, and ancestral area reconstruction of both conserved internal transcribed spacer (ITS) and chloroplast (rps16) DNA sequences. The results showed that the ITS region was more variable than the rps16 region and could be suitable for studying intraspecific phylogeography. The A. trifoliata population displayed high genetic diversity, genetic differentiation and obvious phylogeographical structure, possibly originating on the eastern QTP, expanding during the last glacial-interglacial cycle, diverging in the early Pleistocene and middle Pleistocene, and extensively migrating thereafter. The migration route from west to east along rivers could be largely responsible for the long-distance dispersal of this species, while three main refuges (Qinba Mountains, Nanling Mountains and Yunnan-Guizhou Plateau) with multiple ice shelters facilitated its wide distribution. Conclusions Our results suggested that the from west to east long migration accompanying with the minor short reciprocal migration in the south-north direction, and the three main refuges (the Qinba Mountains, Nanling Mountains and Yunnan-Guizhou Plateau) contributed to the extant geographical distribution of A. trifoliata. In addition, this finding also strongly reduced the discrepancy between glacial contraction and postglacial expansion and the in situ survival hypothesis by simultaneously considering the existence of many similar climate-related ecological niches and migration influences.

Abstract Background The sustainable supply of medicinal plants is important, and cultivating and domesticating them has been suggested as an optimal strategy. However, this can lead to a loss of genetic diversity. Tripterygium wilfordii Hook. f. is a medicinal plant commonly used in traditional Chinese medicine, but its wild populations are dwindling due to excessive harvesting. To protect the species and meet the increasing demand, it is urgent to cultivate it on a large scale. However, distinguishing between T. wilfordii and T. hypoglaucum, two similar species with different medicinal properties, is challenging. Therefore, it is crucial to understand the genetic diversity and population structure of these species for their sustainable utilization. Results In this study, we investigated the genetic diversity and population structure of the two traditional medicinal semiwoody vines plant species, Tripterygium wilfordii and T. hypoglaucum, including wild and cultivated populations using chloroplast DNA (cpDNA) sequences and microsatellite loci. Our results indicated that the two species maintain a high level of genetic divergence, indicating possible genetic bases for the different contents of bioactive compounds of the two species. T. wilfordii showed lower genetic diversity and less subdivided population structures of both markers than T. hypoglaucum. The potential factors in shaping these interesting differences might be differentiated pollen-to-seed migration rates, interbreeding, and history of population divergence. Analyses of cpDNA and microsatellite loci supported that the two species are genetically distinct entities. In addition, a significant reduction of genetic diversity was observed for cultivated populations of the two species, which mainly resulted from the small initial population size and propagated vegetative practice during their cultivation. Conclusion Our findings indicate significant genetic divergence between T. wilfordii and T. hypoglaucum. The genetic diversity and population structure analyses provide important insights into the sustainable cultivation and utilization of these medicinal plants. Accurate identification and conservation efforts are necessary for both species to ensure the safety and effectiveness of crude drug use. Our study also highlighted the importance of combined analyses of different DNA markers in addressing population genetics of medicinal plants because of the contrasts of inheritance and rates of gene flow. Large-scale cultivation programs should consider preserving genetic diversity to enhance the long-term sustainability of T. wilfordii and T. hypoglaucum. Our study proposed that some populations showed higher genetic diversity and distinctness, which can be considered with priority for conservation and as the sources for future breeding and genetic improvement.

Abstract Background Lassa fever is a hemorrhagic disease caused by Lassa virus (LASV), which has been classified by the World Health Organization as one of the top infectious diseases requiring prioritized research. Previous studies have provided insights into the classification and geographic characteristics of LASV lineages. However, the factor of the distribution and evolution characteristics and phylodynamics of the virus was still limited. Methods To enhance comprehensive understanding of LASV, we employed phylogenetic analysis, reassortment and recombination detection, and variation evaluation utilizing publicly available viral genome sequences. Results The results showed the estimated the root of time of the most recent common ancestor (TMRCA) for large (L) segment was approximately 634 (95% HPD: [385879]), whereas the TMRCA for small (S) segment was around 1224 (95% HPD: [10301401]). LASV primarily spread from east to west in West Africa through two routes, and in route 2, the virus independently spread to surrounding countries through Liberia, resulting in a wider spread of LASV. From 1969 to 2018, the effective population size experienced two significant increased, indicating the enhanced genetic diversity of LASV. We also found the evolution rate of L segment was faster than S segment, further results showed zinc-binding protein had the fastest evolution rate. Reassortment events were detected in multiple lineages including sub-lineage IIg, while recombination events were observed within lineage V. Significant amino acid changes in the glycoprotein precursor of LASV were identified, demonstrating sequence diversity among lineages in LASV. Conclusion This study comprehensively elucidated the transmission and evolution of LASV in West Africa, providing detailed insights into reassortment events, recombination events, and amino acid variations.

ABSTRACTThe spread of Lassa virus (LASV) in Guinea, Liberia and Sierra Leone, which together are named the Mano River Union (MRU) area, was examined phylogeographically. To provide a reliable evolutionary scenario, new rodent-derived, whole LASV sequences were included. These were generated by metatranscriptomic next-generation sequencing from rodents sampled between 2003 and 2020 in 21 localities of Guinea and Sierra Leone. An analysis was performed using BEAST to perform continuous phylogeographic inference and EvoLaps v36 to visualize spatio-temporal spread. LASV was identified as expected in its primary host reservoir, the Natal multimammate mouse (Mastomys natalensis), and also in two Guinean multimammate mice (Mastomys erythroleucus) in northern Sierra Leone and two rusty-bellied brush-furred mice (Lophuromys sikapusi) in southern Sierra Leone. This finding is consistent with the latter two species being secondary host reservoirs. The strains in these three species were very closely related in LASV lineage IV. Phylogenetic analysis indicated that the most recent common ancestor of lineage IV existed 316–374 years ago and revealed distinct, well-supported clades from Sierra Leone (Bo, Kabala and Kenema), Guinea (Faranah, Kissidougou-Guekedou and Macenta) and Liberia (Phebe-Ganta). The phylogeographic scenario suggests southern Guinea as the point of origin of LASV in the MRU area, with subsequent spread to towards Mali, Liberia and Sierra Leone at a mean speed of 1.6 to 1.1 km/year.

Getah virus (GETV) mainly causes disease in livestock and may pose an epidemic risk due to its expanding host range and the potential of long-distance dispersal through animal trade. Here, we used metagenomic next-generation sequencing (mNGS) to identify GETV as the pathogen responsible for reemerging swine disease in China and subsequently estimated key epidemiological parameters using phylodynamic and spatially-explicit phylogeographic approaches. The GETV isolates were able to replicate in a variety of cell lines, including human cells, and showed high pathogenicity in a mouse model, suggesting the potential for more mammal hosts. We obtained 16 complete genomes and 79 E2 gene sequences from viral strains collected in China from 2016 to 2021 through large-scale surveillance among livestock, pets, and mosquitoes. Our phylogenetic analysis revealed that three major GETV lineages are responsible for the current epidemic in livestock in China. We identified three potential positively selected sites and mutations of interest in E2, which may impact the transmissibility and pathogenicity of the virus. Phylodynamic inference of the GETV demographic dynamics identified an association between livestock meat consumption and the evolution of viral genetic diversity. Finally, phylogeographic reconstruction of GETV dispersal indicated that the sampled lineages have preferentially circulated within areas associated with relatively higher mean annual temperature and pig population density. Our results highlight the importance of continuous surveillance of GETV among livestock in southern Chinese regions associated with relatively high temperatures. IMPORTANCE Although livestock is known to be the primary reservoir of Getah virus (GETV) in Asian countries, where identification is largely based on serology, the evolutionary history and spatial epidemiology of GETV in these regions remain largely unknown. Through our sequencing efforts, we provided robust support for lineage delineation of GETV and identified three major lineages that are responsible for the current epidemic in livestock in China. We further analyzed genomic and epidemiological data to reconstruct the recent demographic and dispersal history of GETV in domestic animals in China and to explore the impact of environmental factors on its genetic diversity and its diffusion. Notably, except for livestock meat consumption, other pig-related factors such as the evolution of live pig transport and pork production do not show a significant association with the evolution of viral genetic diversity, pointing out that further studies should investigate the potential contribution of other host species to the GETV outbreak. Our analysis of GETV demonstrates the need for wider animal species surveillance and provides a baseline for future studies of the molecular epidemiology and early warning of emerging arboviruses in China.

First described in 2004 off California, Osedax worms are now known from many of the world's oceans, ranging from 10 to over 4000 m in depth. Currently, little is known about species ranges, since most descriptions are from single localities. In this study, we used new sampling in the north-eastern Pacific and available GenBank data from off Japan and Brazil to report expanded ranges for five species: Osedaxfrankpressi, O.knutei, O.packardorum, O.roseus and O.talkovici. We also provided additional DNA sequences from previously reported localities for two species: Osedaxpriapus and O.randyi. To assess the distribution of each species, we used cytochrome c oxidase subunit I (COI) sequences to generate haplotype networks and assess connectivity amongst localities where sampling permitted. Osedaxfrankpressi, O.packardorum, O.priapus, O.roseus and O.talkovici all had one or more dominant COI haplotypes shared by individuals at multiple localities, suggesting high connectivity throughout some or all of their ranges. Low ΦST values amongst populations for O.packardorum, O.roseus and O.talkovici confirmed high levels of gene flow throughout their known ranges. High ΦST values for O.frankpressi between the eastern Pacific and the Brazilian Atlantic showed little gene flow, reflected by the haplotype network, which had distinct Pacific and Atlantic haplotype clusters. This study greatly expands the ranges and provides insights into the phylogeography for these nine species.

The two major rainforests of the neotropics, the Amazon and Atlantic forests, show maximum expansion during the warm and wet conditions of interglacial periods, including the current Holocene. They are connected by a network of gallery forests through the Cerrado biome. However, the extent of their expansion during glacial periods, when they were more disjunct, is unknown. During glacial periods, a pollen assemblage comprising Podocarpus–Ilex–Hedyosmum–Myrsine displays higher frequencies in marine, continental and coastal Brazilian pollen records. This assemblage is observed today in the high-elevation grasslands of the Cerrado and Atlantic Forest biomes and in the coastal vegetation, the restinga, of southern Brazil. We therefore reviewed the possible migration routes for these species by tracking glacial period Podocarpus–Ilex–Hedyosmum–Myrsine assemblages in published pollen records. The marine pollen records provide evidence of a glacial expansion of restinga, its floristic composition being continuous with the dominant regional vegetation, a cold type of shrubby grassland. There appear to be two migration routes, one involving the expansion of high-elevation grassland taxa within the lowlands, and the other low-sea level stands of coastal restinga. We conclude that the Cerrado was a node of migration between the Amazon and Atlantic rainforests, linking the Andes to the central and coastal mountains of Brazil. The Brazilian mountain ranges represent refugia of ancient taxa that colonized the continent up to the Andes and modified the floristic composition of the two rainforests during the Pleistocene glacial periods.

In this study of evolutionary relationships in the subfamily Rubioideae (Rubiaceae), we take advantage of the off-target proportion of reads generated via previous target capture sequencing projects based on nuclear genomic data to build a plastome phylogeny and investigate cytonuclear discordance. The assembly of off-target reads resulted in a comprehensive plastome dataset and robust inference of phylogenetic relationships, where most intratribal and intertribal relationships are resolved with strong support. While the phylogenetic results were mostly in agreement with previous studies based on plastome data, novel relationships in the plastid perspective were also detected. For example, our analyses of plastome data provide strong support for the SCOUT clade and its sister relationship to the remaining members of the subfamily, which differs from previous results based on plastid data but agrees with recent results based on nuclear genomic data. However, several instances of highly supported cytonuclear discordance were identified across the Rubioideae phylogeny. Coalescent simulation analysis indicates that while ILS could, by itself, explain the majority of the discordant relationships, plastome introgression may be the better explanation in some cases. Our study further indicates that plastomes across the Rubioideae are, with few exceptions, highly conserved and mainly conform to the structure, gene content, and gene order present in the majority of the flowering plants. [ABSTRACT FROM AUTHOR]

Simple Summary: Rana temporaria is one of the most widespread Palearctic brown frogs, with two main genetic clades in Europe, geographic spatial pattern of which is insufficiently known. We analyzed samples from the understudied western and central Balkans to evaluate the haplotype diversity of widely used 16SrRNA and MT-CYTB mitochondrial gene sequences and to recognize the contour of a possible contact zone between the main clades. The results revealed a suture zone between the Western and Eastern Clades in the central part of the Balkan Peninsula. Overall, haplotype diversity in the western and central Balkans sample is high. Harboring both main genetic clades of R. temporaria qualifies the Balkan Peninsula as another important center of species' genetic diversity, as well as rich in unique haplotypes. Rana temporaria is one of the most widespread Palearctic brown frogs. We aimed to clarify distribution pattern of two main genetic clades in the understudied Balkan peninsula by using 16SrRNA and MT-CYTB sequences, already widely applied in analyses of populations from other parts of Europe, while focusing on the broad area along the Morava river (central Balkans) as a known gap in the species distribution. Additionally, we were interested in revealing the extent of haplotype diversity within the main genetic clades in the Balkans, particularly around the supposed suture zone. The results revealed a suture zone between the Western and Eastern Clades in the central part of the Balkan Peninsula. This indicated the existence of a historical barrier between the Balkan Mountain Belt and geographically close mountains surrounding the Vlasina Plateau (Rhodope/Serbian–Macedonian Massif). The overall observed haplotype diversity in populations of R. temporaria from the Balkan Peninsula seems high. Harboring both main genetic clades of R. temporaria qualifies the Balkan Peninsula as another important center of species' genetic diversity, as well as rich in unique haplotypes. This points out the necessity of applying conservation measures focused on the common European frog populations and habitats in this part of the species' distribution area. [ABSTRACT FROM AUTHOR]

Traditional methods of phylogenetic reconstruction and species delimitation may be impeded by frequent hybridization among lineages. In this study, we conducted phylogenetic and clustering analyses of ddRAD genomic data on the entire genus Salamandra, which includes six species and over 25 subspecies of terrestrial salamanders. We expanded previous datasets to include missing subspecies and incorporated new samples, with an emphasis on secondary contact zones. Results obtained from a full dataset of 392 individuals (356,874 bp; 24,192 SNPs) were compared with those obtained after excluding substantially admixed individuals (n = 95; 835,467 bp; 51,557 SNPs) to explore the consequences of introgression on phylogenetic inference and taxonomic arrangement of subspecies. We found conflicting phylogenetic placements for taxa represented by many admixed individuals (identified by clustering ancestries). In contrast, a time-calibrated tree constructed without hybrids largely agrees with previous phylogenetic hypotheses. Within S. atra, we found paraphyly of S. atra atra, suggesting an additional candidate subspecies. Within S. infraimmaculata, two lineages are assignable to known subspecies and we additionally identified a third, deeply diverged lineage sampled near the Turkish/Syrian border. In S. algira, we found limited admixture between the subspecies S. a. tingitana and S. a. splendens despite their geographic proximity. Finally, within S. salamandra, we detected significant levels of hybridization between subspecies, which blurred their phylogenetic relationships, although the removal of admixed samples in subset analyses clarified the situation in most cases. Monophyly was recovered for subspecies that were previously found paraphyletic, including S. s. salamandra, S. s. gallaica, and S. s. fastuosa. Salamandra s. "alfredschmidti" was confirmed to be a junior synonym of S. s. bernardezi. Previously disputed subspecies, like S. s. "molleri" and S. s. "hispanica", correspond to separated lineages but are affected by admixture with other lineages. Further newly identified candidate subspecies in S. salamandra included a southern lineage within S. s. wer-neri and a western lineage within S. s. bernardezi. Finally, we re-evaluate the status of recognized subspecies in Salamandra, based on evidence from multiple delimitation criteria. Given that the evolutionary history could not be resolved for all subspecies, we highlight taxa within Salamandra that warrant further molecular examination and taxonomic revision, notably within the S. s. gallaica/"molleri"/bejarae complex. This study illustrates the impact of hybridization in phylogenetic analyses and its downstream effects in the identification of conservation units and their naming in the Linnean classification. [ABSTRACT FROM AUTHOR]