Genetic differentiation and phylogeography of Erythroneurini (Hemiptera, Cicadellidae, Typhlocybinae) in the southwestern karst area of China

Erythroneurini is the largest tribe of the microleafhopper subfamily Typhlocybinae. Most prior research on this tribe has focused on traditional classification, phylogeny, and control of agricultural pests, and the phylogeography of the group remains poorly understood. In this study, the mitochondrial genomes of 10 erythroneurine species were sequenced, and sequences of four genes were obtained for 12 geographical populations of Seriana bacilla. The new sequence data were combined with previously available mitochondrial DNA sequence data and analyzed using Bayesian and Maximum-Likelihood-based phylogenetic methods to elucidate relationships among genera and species and estimate divergence times. Seriana was shown to be derived from within Empoascanara. Phylogeographic and population genetic analysis of the endemic Chinese species Seriana bacilla suggest that the species diverged about 54.85 Mya (95% HPD: 20.76-66.23 million years) in the Paleogene period and that population divergence occurred within the last 14 million years. Ancestral area reconstruction indicates that Seriana bacilla may have originated in the central region of Guizhou, and geographical barriers are the main factors affecting gene flow among populations. Ecological niche modeling using the MaxEnt model suggests that the distribution of the species was more restricted in the past but is likely to expand in the future years 2050 and 2070.

Phylogenetics and Population Genetics of the Asian House Shrew, Suncus murinus-S. montanus Species Complex, Inferred From Whole-Genome and Mitochondrial DNA Sequences, with Special Reference to the Ryukyu Archipelago, Japan

The house shrew (Suncus murinus-S. montanus species complex) colonized regions across southern Asia and the Indian Ocean following human activity. The house shrew is distributed on islands of the Ryukyu Archipelago, the southernmost part of Japan, but the evolutionary history of the shrew on those islands and possible associations between these populations and humans remain unknown. In this study, we conducted phylogenetic and population genetic analyses based on both nuclear and mitochondrial genome sequences of house shrews. Phylogenetic analyses based on mitochondrial cytochrome b (cytb) sequences revealed that shrews from the Ryukyu Archipelago showed strong genetic affinity to Vietnamese and southern Chinese shrews. Demographic analyses of cytb sequences indicated a rapid population expansion event affecting the haplotype group in Vietnam, southern China, and the Ryukyu Archipelago 3300-7900 years ago. Furthermore, gene flow between Ryukyu (Yonaguni Island) and Taiwan and between Ryukyu and Vietnam inferred from f4 statistics of the nuclear genomes suggested repeated immigration to Ryukyu in recent years. The present study demonstrates that the Nagasaki population has a different origin from the Ryukyu population. These findings elucidate the complex pattern of genetic admixture in house shrews and provide insights into their evolutionary history.

Multiple transitions between realms shape relict lineages of Proteus cave salamanders

In comparison to biodiversity on Earth's surface, subterranean biodiversity has largely remained concealed. The olm (Proteus anguinus) is one of the most enigmatic extant cave inhabitants, and until now little was known regarding its genetic structure and evolutionary history. Olms inhabit subterranean waters throughout the Dinaric Karst of the western Balkans, with a seemingly uniform phenotypic appearance of cave-specialized traits: an elongate body, snout and limbs, degenerated eyes and loss of pigmentation ("white olm"). Only a single small region in southeastern Slovenia harbours olms with a phenotype typical of surface animals: pigmented skin, eyes, a blunt snout and short limbs ("black olm"). We used a combination of mitochondrial DNA and genome-wide single nucleotide polymorphism data to investigate the molecular diversity, evolutionary history and biogeography of olms along the Dinaric Karst. We found nine deeply divergent species-level lineages that separated between 17 and 4 million years ago, while molecular diversity within lineages was low. We detected no signal of recent admixture between lineages and only limited historical gene flow. Biogeographically, the contemporaneous distribution of lineages mostly mirrors hydrologically separated subterranean environments, while the historical separation of olm lineages follows microtectonic and climatic changes in the area. The reconstructed phylogeny suggests at least four independent transitions to the cave phenotype. Two of the species-level lineages have miniscule ranges and may represent Europe's rarest amphibians. Their rarity and the decline in other lineages call for protection of their subterranean habitats.

Population Variation and Phylogeography of Cherry Blossom (Prunus conradinae) in China

Prunus conradinae (subgenus Cerasus, Rosaceae) is a significant germplasm resource of wild cherry blossom in China. To ensure the comprehensiveness of this study, we used a large sample size (12 populations comprising 244 individuals) which involved the fresh leaves of P. conradinae in Eastern, Central, and Southwestern China. We combined morphological and molecular evidence (three chloroplast DNA (cpDNA) sequences and one nuclear DNA (nr DNA) sequence) to examine the population of P. conradinae variation and differentiation. Our results revealed that Central, East, and Southwest China are important regions for the conservation of P. conradinae to ensure adequate germplasm resources in the future. We also found support for a new variant, P. conradinae var. rubrum. We observed high genetic diversity within P. conradinae (haplotype diversity [Hd] = 0.830; ribotype diversity [Rd] = 0.798), with novel genetic variation and a distinct genealogical structure among populations. There was genetic variation among populations and phylogeographic structure among populations and three geographical groups (Central, East, and Southwest China). The genetic differentiation coefficient was the lowest in the Southwest region and the gene exchange was obvious, while the differentiation was obvious in Central China. In the three geographic groups, we identified two distinct lineages: an East China lineage (Central China and East China) and a Southwest China lineage ((Central China and Southwest China) and East China). These two lineages originated approximately 4.38 million years ago (Mya) in the early Pliocene due to geographic isolation. P. conradinae expanded from Central China to East China at 3.32 Mya (95% HPD: 1.12-5.17 Mya) in the Pliocene. The population of P. conradinae spread from East China to Southwest China, and the differentiation time was 2.17 Mya (95% (HPD: 0.47-4.54 Mya), suggesting that the population of P. conradinae differentiated first in Central and East China. The population of P. conradinae experienced differentiation from Central China to Southwest China around 1.10 Mya (95% HPD: 0.11-2.85 Mya) during the early Pleistocene of the Quaternary period. The southeastern region of East China, near Mount Wuyi, likely serves as a refuge for P. conradinae. This study establishes a theoretical foundation for the classification, identification, conservation, and exploitation of germplasm resources of P. conradinae.

Tracing evolutionary trajectories in the presence of gene flow in South American temperate lizards (Squamata: Liolaemus kingii group)

Evolutionary processes behind lineage divergence often involve multidimensional differentiation. However, in the context of recent divergences, the signals exhibited by each dimension may not converge. In such scenarios, incomplete lineage sorting, gene flow, and scarce phenotypic differentiation are pervasive. Here, we integrated genomic (RAD loci of 90 individuals), phenotypic (linear and geometric traits of 823 and 411 individuals, respectively), spatial, and climatic data to reconstruct the evolutionary history of a speciation continuum of liolaemid lizards (Liolaemus kingii group). Specifically, we (a) inferred the population structure of the group and contrasted it with the phenotypic variability; (b) assessed the role of postdivergence gene flow in shaping phylogeographic and phenotypic patterns; and (c) explored ecogeographic drivers of diversification across time and space. We inferred eight genomic clusters exhibiting leaky genetic borders coincident with geographic transitions. We also found evidence of postdivergence gene flow resulting in transgressive phenotypic evolution in one species. Predicted ancestral niches unveiled suitable areas in southern and eastern Patagonia during glacial and interglacial periods. Our study underscores integrating different data and model-based approaches to determine the underlying causes of diversification, a challenge faced in the study of recently diverged groups. We also highlight Liolaemus as a model system for phylogeographic and broader evolutionary studies.

Genomic analysis supports Cape Lion population connectivity prior to colonial eradication and extinction

Cape lions (Panthera leo melanochaitus) formerly ranged throughout the grassland plains of the "Cape Flats" in what is today known as the Western Cape Province, South Africa. Cape lions were likely eradicated because of overhunting and habitat loss after European colonization. European naturalists originally described Cape lions as "black-maned lions" and claimed that they were phenotypically distinct. However, other depictions and historical descriptions of lions from the Cape report mixed or light coloration and without black or extensively developed manes. These findings suggest that, rather than forming a distinct population, Cape lions may have had phenotypic and genotypic variation similar to other African lions. Here we investigate Cape lion genome characteristics, population dynamics, and genetic distinctiveness prior to their extinction. We generated genomic data from 2 historic Cape lions to compare to 118 existing high-coverage mitogenomes, and low-coverage nuclear genomes of 53 lions from 13 African countries. We show that, before their eradication, lions from the Cape Flats had diverse mitogenomes and nuclear genomes that clustered with lions from both southern and eastern Africa. Cape lions had high genome-wide heterozygosity and low inbreeding coefficients, indicating that populations in the Cape Flats went extinct so rapidly that genomic effects associated with long-term small population size and isolation were not detectable. Our findings do not support the characterization of Cape lions as phylogeographically distinct, as originally put forth by some European naturalists, and illustrates how alternative knowledge systems, for example, Indigenous perspectives, could potentially further inform interpretations of species histories.

Phylogeographic analysis reveals extensive genetic variation of native grass Elymus nutans (Poaceae) on the Qinghai-Tibetan plateau

Introduction

Elymus nutans holds ecological and pastoral significance due to its adaptability and nutritional value, the Qinghai-Tibet Plateau (QTP) is a key hub for its genetic diversity. To conserve and harness its genetic resources in highland ecosystems, a thorough assessment is vital. However, a comprehensive phylogeographic exploration of E. nutans is lacking. The objective of this study was to unravel the genetic diversity, adaptation, and phylogenetics of E. nutans populations.

Methods

Encompassing 361 individuals across 35 populations, the species' genetic landscape and dynamic responses to diverse environments were decoded by using four chloroplast DNA (cpDNA) sequences and nine microsatellite markers derived from the transcriptome.

Results and discussion

This study unveiled a notable degree of genetic diversity in E. nutans populations at nuclear (I = 0.46, He = 0.32) and plastid DNA levels (Hd = 0.805, π = 0.67). Analysis via AMOVA highlighted genetic variation predominantly within populations. Despite limited isolation by distance (IBD), the Mekong-Salween Divide (MSD) emerged as a significant factor influencing genetic differentiation and conserving diversity. Furthermore, correlations were established between external environmental factors and effective alleles of three EST-SSRs (EN5, EN57 and EN80), potentially linked to glutathione S-transferases T1 or hypothetical proteins, affecting adaptation. This study deepens the understanding of the intricate relationship between genetic diversity, adaptation, and environmental factors within E. nutans populations on the QTP. The findings shed light on the species' evolutionary responses to diverse ecological conditions and contribute to a broader comprehension of plant adaptation mechanisms.

Spatiotemporal dynamics and epidemiological impact of SARS-CoV-2 XBB lineage dissemination in Brazil in 2023

The SARS-CoV-2 XBB is a group of highly immune-evasive lineages of the Omicron variant of concern that emerged by recombining BA.2-descendent lineages and spread worldwide during 2023. In this study, we combine SARS-CoV-2 genomic data (n = 11,065 sequences) with epidemiological data of severe acute respiratory infection (SARI) cases collected in Brazil between October 2022 and July 2023 to reconstruct the space-time dynamics and epidemiologic impact of XBB dissemination in the country. Our analyses revealed that the introduction and local emergence of lineages carrying convergent mutations within the Spike protein, especially F486P, F456L, and L455F, propelled the spread of XBB* lineages in Brazil. The average relative instantaneous reproduction numbers of XBB* + F486P, XBB* + F486P + F456L, and XBB* + F486P + F456L + L455F lineages in Brazil were estimated to be 1.24, 1.33, and 1.48 higher than that of other co-circulating lineages (mainly BQ.1*/BE*), respectively. Despite such a growth advantage, the dissemination of these XBB* lineages had a reduced impact on Brazil's epidemiological scenario concerning previous Omicron subvariants. The peak number of SARI cases from SARS-CoV-2 during the XBB wave was approximately 90%, 80%, and 70% lower than that observed during the previous BA.1*, BA.5*, and BQ.1* waves, respectively. These findings revealed the emergence of multiple XBB lineages with progressively increasing growth advantage, yet with relatively limited epidemiological impact in Brazil throughout 2023. The XBB* + F486P + F456L + L455F lineages stand out for their heightened transmissibility, warranting close monitoring in the months ahead.

IMPORTANCE

Brazil was one the most affected countries by the SARS-CoV-2 pandemic, with more than 700,000 deaths by mid-2023. This study reconstructs the dissemination of the virus in the country in the first half of 2023, a period characterized by the dissemination of descendants of XBB.1, a recombinant of Omicron BA.2 lineages evolved in late 2022. The analysis supports that XBB dissemination was marked by the continuous emergence of indigenous lineages bearing similar mutations in key sites of their Spike protein, a process followed by continuous increments in transmissibility, and without repercussions in the incidence of severe cases. Thus, the results suggest that the epidemiological impact of the spread of a SARS-CoV-2 variant is influenced by an intricate interplay of factors that extend beyond the virus's transmissibility alone. The study also underlines the need for SARS-CoV-2 genomic surveillance that allows the monitoring of its ever-shifting composition.

A renewed glance at the Palearctic golden eagle: Genetic variation in space and time

Anthropogenic pressures on nature have been causing population declines for centuries. Intensified persecution of apex predators, like the golden eagle, resulted in population bottlenecks during the 19th and 20th centuries. To study population genetics and demographic history of the golden eagle throughout its distribution, we collected museum samples from previously underrepresented regions, such as Russia and Central Asia. We used 12 microsatellite loci and a fragment of the mitochondrial DNA control region to re-evaluate phylogeography of Eurasian golden eagles and study the impacts of the population bottleneck. Our results revealed a north-south genetic gradient, expressed by the difference between Mediterranean and Holarctic lineages, as well as genetically distinct Northern Europe and Central Asia and Caucasus regions. Furthermore, Northern Europe exhibited the lowest, whereas Central Asia and Caucasus had the highest genetic diversity. Although golden eagles maintained relatively high genetic diversity, we detected genetic signatures of the recent bottleneck, including reduced genetic diversity and a decline in the effective female population size around the year 1975. Our study improves the knowledge of the genetic composition of Eurasian golden eagles and highlights the importance of understanding their historical population dynamics in the face of ongoing and future conservation efforts.

Population genomic evidence of structured and connected Plasmodium vivax populations under host selection in Latin America

Pathogen genomic epidemiology has the potential to provide a deep understanding of population dynamics, facilitating strategic planning of interventions, monitoring their impact, and enabling timely responses, and thereby supporting control and elimination efforts of parasitic tropical diseases. Plasmodium vivax, responsible for most malaria cases outside Africa, shows high genetic diversity at the population level, driven by factors like sub-patent infections, a hidden reservoir of hypnozoites, and early transmission to mosquitoes. While Latin America has made significant progress in controlling Plasmodium falciparum, it faces challenges with residual P. vivax. To characterize genetic diversity and population structure and dynamics, we have analyzed the largest collection of P. vivax genomes to date, including 1474 high-quality genomes from 31 countries across Asia, Africa, Oceania, and America. While P. vivax shows high genetic diversity globally, Latin American isolates form a distinctive population, which is further divided into sub-populations and occasional clonal pockets. Genetic diversity within the continent was associated with the intensity of transmission. Population differentiation exists between Central America and the North Coast of South America, vs. the Amazon Basin, with significant gene flow within the Amazon Basin, but limited connectivity between the Northwest Coast and the Amazon Basin. Shared genomic regions in these parasite populations indicate adaptive evolution, particularly in genes related to DNA replication, RNA processing, invasion, and motility - crucial for the parasite's survival in diverse environments. Understanding these population-level adaptations is crucial for effective control efforts, offering insights into potential mechanisms behind drug resistance, immune evasion, and transmission dynamics.

Piecing the barcoding puzzle of Palearctic water frogs (Pelophylax) sheds light on amphibian biogeography and global invasions

Palearctic water frogs (genus Pelophylax) are an outstanding model in ecology and evolution, being widespread, speciose, either threatened or threatening to other species through biological invasions, and capable of siring hybrid offspring that escape the rules of sexual reproduction. Despite half a century of genetic research and hundreds of publications, the diversity, systematics and biogeography of Pelophylax still remain highly confusing, in no small part due to a lack of correspondence between studies. To provide a comprehensive overview, we gathered >13,000 sequences of barcoding genes from >1700 native and introduced localities and built multigene mitochondrial (~17 kb) and nuclear (~10 kb) phylogenies. We mapped all currently recognized taxa and their phylogeographic lineages (>40) to get a grasp on taxonomic issues, cyto-nuclear discordances, the genetic makeup of hybridogenetic hybrids, and the origins of introduced populations. Competing hypotheses for the molecular calibration were evaluated through plausibility tests, implementing a new approach relying on predictions from the anuran speciation continuum. Based on our timetree, we propose a new biogeographic paradigm for the Palearctic since the Paleogene, notably by attributing a prominent role to the dynamics of the Paratethys, a vast paleo-sea that extended over most of Europe. Furthermore, our results show that distinct marsh frog lineages from Eastern Europe, the Balkans, the Near East, and Central Asia (P. ridibundus ssp.) are naturally capable of inducing hybridogenesis with pool frogs (P. lessonae). We identified 14 alien lineages (mostly of P. ridibundus) over ~20 areas of invasions, especially in Western Europe, with genetic signatures disproportionally pointing to the Balkans and Anatolia as the regions of origins, in line with exporting records of the frog leg industry and the stocks of pet sellers. Pelophylax thus emerges as one of the most invasive amphibians worldwide, and deserves much higher conservation concern than currently given by the authorities fighting biological invasions.

Molecular evolution and phylogeographic analysis of wheat dwarf virus

Wheat dwarf virus (WDV) has caused considerable economic loss in the global production of grain crops. Knowledge of the evolutionary biology and population history of the pathogen remain poorly understood. We performed molecular evolution and worldwide phylodynamic analyses of the virus based on the genes in the protein-coding region of the entire viral genome. Our results showed that host-driven and geography-driven adaptation are major factors that affects the evolution of WDV. Bayesian phylogenetic analysis estimates that the average WDV substitution rate was 4.240 × 10-4 substitutions/site/year (95% credibility interval, 2.828 × 10-4-5.723 × 10-4), and the evolutionary rates of genes encoding proteins with virion-sense transcripts and genes encoding proteins with complementary-sense transcripts were different. The positively selected sites were detected in only two genes encoding proteins with complementary-sense, and WDV-barley are subject to stronger purifying selection than WDV-wheat. The time since the most recent common WDV ancestor was 1746 (95% credibility interval, 1517-1893) CE. Further analyses identified that the WDV-barley population and WDV-wheat population experienced dramatic expansion-decline episodes, and the expansion time of the WDV-barley population was earlier than that of the WDV-wheat population. Our phylogeographic analysis showed that the WDV population originating in Iran was subsequently introduced to Europe, and then spread from Eastern Europe to China.

Regional environmental differences significantly affect the genetic structure and genetic differentiation of Carpinus tientaiensis Cheng, an endemic and extremely endangered species from China

Differences in topography and environment greatly affect the genetic structure and genetic differentiation of species, and endemic or endangered species with limited geographic ranges seem to be more sensitive to changes in climate and other environmental factors. The complex topography of eastern China is likely to affect genetic differentiation of plants there. Carpinus tientaiensis Cheng is a native and endangered plants from China, and exploring its genetic diversity has profound significance for protection and the collection of germplasm resources. Based on AFLP markers, this study found that C. tientaiensis has low genetic diversity, which mainly came from within populations, while Shangshantou and Tiantai Mountain populations have relatively high genetic diversity. The Nei genetic distance was closely related to geographical distance, and temperature and precipitation notablely affected the genetic variation and genetic differentiation of C. tientaiensis. Based on cpDNA, this study indicated that C. tientaiensis exhibits a moderate level of genetic diversity, and which mainly came from among populations, while Tiantai Mountain population have the highest genetic diversity. It demonstrated that there was genetic differentiation between populations, which can be divided into two independent geographical groups, but there was no significant phylogeographic structure between them. The MaxEnt model showed that climate change significantly affects its distribution, and the suitable distribution areas in Zhejiang were primarily divided into two regions, eastern Zhejiang and southern Zhejiang, and there was niche differentiation in its suitable distribution areas. Therefore, this study speculated that the climate and the terrain of mountains and hills in East China jointly shape the genetic structure of C. tientaiensis, which gived rise to an obvious north-south differentiation trend of these species, and the populations located in the hilly areas of eastern Zhejiang and the mountainous areas of southern Zhejiang formed two genetic branches respectively.

Evolutionary dynamics and geographical dispersal of Borrelia lusitaniae

Background

Borrelia lusitaniae is a species within the complex Borrelia burgdorferi sensu lato, associated with lizards as reservoirs and Ixodes ricinus as its main vector. Borrelia lusitaniae is predominantly distributed in Central and Southeast Europe, and in countries of the Mediterranean basin, such as Portugal, Morocco, Tunisia, and Italy where this spirochete appears to infect vector ticks more frequently than other genospecies. Evolution of this zoonotic tick-borne microparasite is shaped by different environmental factors. Comprehensive phylogenetic analysis may give insight into how B. lusitaniae spreads to new geographic locations.

Aim

We applied Bayesian statistical methods to B. lusitaniae multilocus sequence typing (MLST) data to study the migration routes of B. lusitaniae and its potential for further spread.

Results

The discrete phylogeographic analysis placed origins of B. lusitaniae in Southeast Europe and identified at least two introductions of B. lusitaniae from Europe to North Africa. Estimated effective reproductive potential (Re), as a key indicator for a pathogen spread, suggested potential for further spread.

Conclusion

The results of this study can provide beneficial information about the potential for further spread of B. lusitaniae in Europe and North Africa and estimation of necessity for the development of strategies to monitor and control Lyme borreliosis.

Repeated patterns of reptile diversification in Western North America supported by the Northern Alligator Lizard (Elgaria coerulea)

Understanding the processes that shape genetic diversity by either promoting or preventing population divergence can help identify geographic areas that either facilitate or limit gene flow. Furthermore, broadly distributed species allow us to understand how biogeographic and ecogeographic transitions affect gene flow. We investigated these processes using genomic data in the Northern Alligator Lizard (Elgaria coerulea), which is widely distributed in Western North America across diverse ecoregions (California Floristic Province and Pacific Northwest) and mountain ranges (Sierra Nevada, Coastal Ranges, and Cascades). We collected single-nucleotide polymorphism data from 120 samples of E. coerulea. Biogeographic analyses of squamate reptiles with similar distributions have identified several shared diversification patterns that provide testable predictions for E. coerulea, including deep genetic divisions in the Sierra Nevada, demographic stability of southern populations, and recent post-Pleistocene expansion into the Pacific Northwest. We use genomic data to test these predictions by estimating the structure, connectivity, and phylogenetic history of populations. At least 10 distinct populations are supported, with mixed-ancestry individuals situated at most population boundaries. A species tree analysis provides strong support for the early divergence of populations in the Sierra Nevada Mountains and recent diversification into the Pacific Northwest. Admixture and migration analyses detect gene flow among populations in the Lower Cascades and Northern California, and a spatial analysis of gene flow identified significant barriers to gene flow across both the Sierra Nevada and Coast Ranges. The distribution of genetic diversity in E. coerulea is uneven, patchy, and interconnected at population boundaries. The biogeographic patterns seen in E. coerulea are consistent with predictions from co-distributed species.

Ancient mitochondrial genomes unveil the origins and evolutionary history of New Zealand's enigmatic takahē and moho

Many avian species endemic to Aotearoa New Zealand were driven to extinction or reduced to relict populations following successive waves of human arrival, due to hunting, habitat destruction and the introduction of mammalian predators. Among the affected species were the large flightless South Island takahē (Porphyrio hochstetteri) and the moho (North Island takahē; P. mantelli), with the latter rendered extinct and the former reduced to a single relictual population. Little is known about the evolutionary history of these species prior to their decline and/or extinction. Here we sequenced mitochondrial genomes from takahē and moho subfossils (12 takahē and 4 moho) and retrieved comparable sequence data from takahē museum skins (n = 5) and contemporary individuals (n = 17) to examine the phylogeny and recent evolutionary history of these species. Our analyses suggest that prehistoric takahē populations lacked deep phylogeographic structure, in contrast to moho, which exhibited significant spatial genetic structure, albeit based on limited sample sizes (n = 4). Temporal genetic comparisons show that takahē have lost much of their mitochondrial genetic diversity, likely due to a sudden demographic decline soon after human arrival (~750 years ago). Time-calibrated phylogenetic analyses strongly support a sister species relationship between takahē and moho, suggesting these flightless taxa diverged around 1.5 million years ago, following a single colonisation of New Zealand by a flighted Porphyrio ancestor approximately 4 million years ago. This study highlights the utility of palaeogenetic approaches for informing the conservation and systematic understanding of endangered species whose ranges have been severely restricted by anthropogenic impacts.

The influence of climate oscillations and geological events on population differentiation of Camponotus japonicus in the Chinese mainland

Camponotus japonicus (Hymenoptera: Formicidae) is an omnivorous social insect which builds sizable colonies in sparse woodlands or cropland and spreads across multiple climatic zones in the Chinese mainland. This study aims to reveal the role of climate changes and geological events in driving the genetic structure of social insect populations by investigating the phylogenetics and historical demography of C. japonicus in the Chinese mainland. Phylogenetic analyses were conducted based on the mitochondria DNA dataset using MrBayes and IQ-TREE. We constructed a haplotype network, calculated analyses of molecular variance, estimated the divergence time, and reconstructed the maximum clade credibility tree. Mismatch distribution and Bayesian skyline plots were used to infer historical population fluctuations. Additionally, ecological niche modeling was employed to predict the potential distribution of the species during the present, mid-holocene, and last glacial maximum periods in the Chinese mainland. The phylogenetic tree and median-joining network analyses support the presence of four distinct lineages in C. japonicus. These lineages exhibit significant genetic differentiation and limited gene flow. The divergence among the four lineages began in the early Pleistocene, approximately 1.41 million years ago (Ma). Subsequently, the central lineage diverged from both the northern and southern lineages around 1.16 Ma, while the northern and southern lineages diverged from each other at approximately 1.07 Ma. Population expansion was observed in the southern, central, and northern lineages prior to the last glacial maximum, while the Yunnan-Sichuan lineage experienced a slight increase in population size in more recent times. The predicted distribution of the species corresponds well with the actual distribution. Furthermore, the current suitable habitat areas in northern Xinjiang, southern Tibet, and the southeast coastal regions have significantly decreased compared to the last glacial maximum and the mid-holocene periods. Our results suggest that climate oscillations and geological events play an important role in driving genetic patterns and differentiation of C. japonicus. Mountain barriers isolate populations from each other, hinder the flow of genes, and effectively prevent the spread of this species. But at the same time, it also formed refugia at low altitudes areas such as Qinling-Bashan Mountains and Yanshan-Taihang Mountains and provide suitable habitats during glaciation. This study provides a good model for understanding how complex climate changes and geological events affect population genetic differentiation of social insects in the Chinese mainland.

Phylogenetic Relationships of Avian Cestodes from Brine Shrimp and Congruence with Larval Morphology

Determining molecular markers for parasites provides a useful tool for their identification, particularly for larval stages with few distinguishable diagnostic characters. Avian cestodes play a key role in the food webs and biodiversity of hypersaline wetlands, yet they remain understudied. Using naturally infected Artemia, we identified cestode larvae (cysticercoids), assessed their genetic diversity, and explored phylogenetic relationships in relation to larval morphology and waterbird final hosts. We obtained partial 18S rDNA sequences for 60 cysticercoids of the family Hymenolepidae infecting Artemia spp. from seven localities and three countries (Spain, the USA, and Chile). We present the first DNA sequences for six taxa: Confluaria podicipina, Fimbriarioides sp., Flamingolepis liguloides, Flamingolepis sp. 1, Flamingolepis sp. 2, and Hymenolepis californicus. Intraspecific sequence variation (0.00-0.19% diversity) was lower than intergroup genetic distance (0.7-14.75%). Phylogenetic analysis revealed three main clades: 1-Flamingolepis, 2-Fimbriarioides, 3-Confluaria and Hymenolepis, all of which separated from hymenolepidids from mammals and terrestrial birds. This clear separation among taxa is congruent with previous morphological identification, validating the 18S gene as a useful marker to discriminate at generic/species level. Working with intermediate hosts allows the expansion of knowledge of taxonomic and genetic diversity of cestodes in wildlife, as well as elucidation of their life cycles.

Low Genetic Variability of the Tundra Vole in Lithuania

The distribution and spread of the tundra vole (Alexandromys oeconomus) in Lithuania have been documented over the last 70 years, but the genetic diversity of the species has not been studied. In this study, we examined A. oeconomus trapped in three sites in northern and western Lithuania using mtDNA sequence analysis of the cytb and control region. The western and northern sites are separated by anthropogenic landscape barriers. The western site is subject to regular spring flooding. Phylogenetic analyses of the studied individuals placed them in the Central European phylogroup, suggesting that Lithuanian A. oeconomus originated from northeastern Poland. In Lithuania, the genetic diversity of A. oeconomus at both mtDNA loci was relatively low (Hd < 0.6, π < 0.002) compared to that found in other European samples (Hd = 0.833-0.958; π = 0.00402-0.01552). Individuals analyzed in Lithuania were genetically different from samples collected in Poland and Northern Europe (ΦST > 0.15, p < 0.05). The genetic divergence between the western and northern samples of A. oeconomus in Lithuania, together with the low genetic variability among the voles studied, provides new insights into the phylogeography of the species and the influence of barriers on the colonization of the country.

Assessing genome-wide adaptations associated with range expansion in the pink rice borer, Sesamia inferens

Understanding the genetic basis of adaptive evolution following habitat expansion can have important implications for pest management. The pink rice borer (PRB), Sesamia inferens (Walker), is a destructive pest of rice that was historically restricted to regions south of 34° N latitude in China. However, with changes in global climate and farming practices, the distribution of this moth has progressively expanded, encompassing most regions in North China. Here, 3 highly differentiated subpopulations were discovered using high-quality single-nucleotide polymorphism and structural variant datasets across China, corresponding to northern, southern China regions, and the Yunnan-Guizhou Plateau, with significant patterns of isolation by geographic and environmental distances. Our estimates of evolutionary history indicate asymmetric migration with varying population sizes across the 3 subpopulations. Selective sweep analyses estimated strong selection at insect cuticle glycine-rich cuticular protein genes which are associated with enhanced desiccation adaptability in the northern group, and at the histone-lysine-N-methyltransferase gene associated with range expansion and local adaptation in the Shandong population. Our findings have significant implications for the development of effective strategies to control this pest.

Genetic Structure and Population History of the Zaisan Toad-Headed Agama (Phrynocephalus melanurus) Inferred from Mitochondrial DNA

The agamid lizard Phrynocephalus melanurus is restricted to Northwest China (Dzungar Basin) and the adjacent Eastern Kazakhstan (Zaisan and Alakol basins). To elucidate the phylogeography of P. melanurus, we obtained the mitochondrial DNA COI segments of 175 sampled lizards from 44 localities across the whole distribution. Phylogenetic analyses revealed two main Clades comprising five geographically structured lineages (I, IIa, IIb1, IIb2, and IIb3) that fit an isolation-by-distance (IBD) model. The divergence from the most recent common ancestor was dated to ~1.87 million years ago (Ma). Demographic analyses demonstrated lineage-specific response to past climate change: stable population for Clade I, Subclade IIb1; past population expansion for IIb3 since 0.18 Ma, respectively. Bayesian phylogeographic diffusion analyses detected initial spreading at the Saur Mount vicinity, approximately 1.8 Ma. Historical species distribution model (SDM) projected expansion of the suitable habitat in the last interglacial and shift and contraction in the last glacial maximum and Holocene epochs. The SDM predicted a drastic reduction in suitable area throughout the range as a response to future climate change. Our findings suggest that the evolution of P. melanurus followed a parapatric divergence with subsequent dispersal and adaptation to cold and dry environments during the Quaternary. Overall, this work improves our understanding of the lineage diversification and population dynamics of P. melanurus, providing further insights into the evolutionary processes that occurred in Northwest China and adjacent Eastern Kazakhstan.

Gene Flow and Isolation in the Arid Nearctic Revealed by Genomic Analyses of Desert Spiny Lizards

The opposing forces of gene flow and isolation are two major processes shaping genetic diversity. Understanding how these vary across space and time is necessary to identify the environmental features that promote diversification. The detection of considerable geographic structure in taxa from the arid Nearctic has prompted research into the drivers of isolation in the region. Several geographic features have been proposed as barriers to gene flow, including the Colorado River, Western Continental Divide, and a hypothetical Mid-Peninsular Seaway in Baja California. However, recent studies suggest that the role of barriers in genetic differentiation may have been overestimated when compared to other mechanisms of divergence. In this study, we infer historical and spatial patterns of connectivity and isolation in Desert Spiny Lizards (Sceloporus magister) and Baja Spiny Lizards (S. zosteromus), which together form a species complex composed of parapatric lineages with wide distributions in arid western North America. Our analyses incorporate mitochondrial sequences, genomic-scale data, and past and present climatic data to evaluate the nature and strength of barriers to gene flow in the region. Our approach relies on estimates of migration under the multispecies coalescent to understand the history of lineage divergence in the face of gene flow. Results show that the S. magister complex is geographically structured, but we also detect instances of gene flow. The Continental Divide is a strong barrier to gene flow, while the Colorado River is more permeable. Analyses yield conflicting results for the catalyst of differentiation of peninsular lineages in S. zosteromus. Our study shows how large-scale genomic data for thoroughly sampled species can shed new light on biogeography. Furthermore, our approach highlights the need for the combined analysis of multiple sources of evidence to adequately characterize the drivers of divergence.

Contrasting genetic and morphological differentiation among geographical lineages of a stenotopic miniature rasborine, Boraras maculatus, in Peninsular Malaysia

The variability in the stenotopic miniature rasborine Boraras maculatus (Cypriniformes: Danionidae: Rasborinae) across acidic-water habitats of Peninsular Malaysia (PM) was investigated using two molecular markers (the mitochondrial cytochrome c oxidase subunit I [COI] gene and the nuclear rhodopsin gene), as well as morphological evidence. Molecular phylogenetic analyses revealed differentiation among populations of B. maculatus in PM with the distinction of four allopatric lineages. Each of them was recognized as a putative species by automatic species delimitation methods. These lineages diverged from each other between 7.4 and 1.9 million years ago. A principal component analysis (PCA) was conducted to examine the multivariate variation in 11 morphometric measurements among three of these lineages. PCA results showed a significant overlap in morphological characteristics among these lineages. Additionally, a photograph-based machine learning approach failed to fully differentiate these lineages, suggesting limited morphological differentiation. B. maculatus represents a case of morphological stasis in a stenotopic miniature species. Strong habitat preference, coupled with long-term habitat fragmentation, may explain why each lineage of B. maculatus has a restricted distribution and did not disperse to other regions within and outside of PM, despite ample possibilities when the Sunda shelf was emerged and drained by large paleodrainages for most of the past 7 million years. The conservation status of B. maculatus and its peat swamp habitats are discussed, and it is concluded that peat swamps comprise several evolutionary units. Each of these units is considered a conservation unit and deserves appropriate protection.

EvoLaps 2: Advanced phylogeographic visualization

EvoLaps is a user-friendly web application designed to visualize the spatial and temporal spread of pathogens. It takes an annotated tree as entry, such as a maximum clade credibility tree obtained through continuous phylogeographic inference. By following a 'Top-Down' reading of a tree recursively, transitions (latitude/longitude changes from a node to its children) are represented on a cartographic background using graphical paths. The complete set of paths forms the phylogeographic scenario. EvoLaps offers several features to analyze complex scenarios: (1) enhanced path display using multiple graphical variables with time-dependent gradients, (2) cross-highlighting and selection capabilities between the phylogeographic scenario and the phylogenetic tree, (3) production of specific spatio-temporal scales and synthetic views through dynamic and iterative clustering of localities into spatial clusters, (4) animation of the phylogeographic scenario using tree brushing, which can be done manually or automatically, gradually over time or at specific time intervals, and for the entire tree or a specific clade, and (5) an evolving library of additional tools. EvoLaps is freely available for use at evolaps.org.