Genomic landscape of introgression from the ghost lineage in a gobiid fish uncovers the generality of forces shaping hybrid genomes

Extinct lineages can leave legacies in the genomes of extant lineages through ancient introgressive hybridization. The patterns of genomic survival of these extinct lineages provide insight into the role of extinct lineages in current biodiversity. However, our understanding on the genomic landscape of introgression from extinct lineages remains limited due to challenges associated with locating the traces of unsampled 'ghost' extinct lineages without ancient genomes. Herein, we conducted population genomic analyses on the East China Sea (ECS) lineage of Chaenogobius annularis, which was suspected to have originated from ghost introgression, with the aim of elucidating its genomic origins and characterizing its landscape of introgression. By combining phylogeographic analysis and demographic modelling, we demonstrated that the ECS lineage originated from ancient hybridization with an extinct ghost lineage. Forward simulations based on the estimated demography indicated that the statistic γ of the HyDe analysis can be used to distinguish the differences in local introgression rates in our data. Consistent with introgression between extant organisms, we found reduced introgression from extinct lineage in regions with low recombination rates and with functional importance, thereby suggesting a role of linked selection that has eliminated the extinct lineage in shaping the hybrid genome. Moreover, we identified enrichment of repetitive elements in regions associated with ghost introgression, which was hitherto little known but was also observed in the re-analysis of published data on introgression between extant organisms. Overall, our findings underscore the unexpected similarities in the characteristics of introgression landscapes across different taxa, even in cases of ghost introgression.

Phylogeny and spatial distribution of Japanese encephalitis virus vector species in Cambodia

In Southeast Asia, despite the use of Japanese encephalitis vaccines and vaccination coverage, Japanese encephalitis (JE) transmission is still a major public health issue. The main vectors of this virus are mosquitoes from the genus Culex, which diversity and density are important in Southeast Asia. The main vector species of Japanese encephalitis virus (JEV) in Cambodia belong to the Vishnui subgroup. However, their morphological identification solely based on the adult stage remains challenging, making their segregation and detection difficult. In order to identify and describe the distribution of the three main JEV vector species in Cambodia, namely Culex vishnui, Cx. pseudovishnui and Cx. tritaeniorhynchus, mosquito samplings were carried out throughout the country in different environments. Phylogenetic analysis of the cytochrome c oxidase subunit I (coI) gene using maximum-likelihood tree with ultrafast bootstrap and phylogeographic analysis were performed. The three main Culex species are phylogenetically separated, and represent two distinct clades, one with Cx. tritaeniorhynchus and the second with Cx. vishnui and Cx. pseudovishnui, the latter appearing as a subgroup of Cx. vishnui. The phylogeographic analysis shows a distribution of the Vishnui subgroup on the entire Cambodian territory with an overlapped distribution areas leading to a sympatric distribution of these species. The three JEV vector species are geographically well-defined with a strong presence of Cx. pseudovishnui in the forest. Combined with the presence of Cx. tritaeniorhynchus and Cx. vishnui in rural, peri-urban, and urban areas, the presence of JEV-competent vectors is widespread in Cambodia.

Approximate Bayesian computation and ecological niche models elucidate the demographic history and current fragmented population distribution of a Korean endemic shrub

Climatic fluctuations and geological events since the LGM are believed to have significantly impacted the population size, distribution, and mobility of many species that we observe today. In this paper, we determined the processes driving the phylogeographic structure of the Korean endemic white forsythia by combining the use of genome-wide SNPs and predicting paleoclimatic habitats during the LGM (21 kya), Early Holocene (10 kya), Mid-Holocene (6 kya), and Late Holocene (3 kya). Using a maximum of 1897 SNPs retrieved from 124 samples across nine wild populations, five environmental predictors, and the species' natural occurrence records, we aimed to infer the species' demographic history and reconstruct its possible paleodistributions with the use of approximate Bayesian computation and ecological niche models, respectively. Under this integrated framework, we found strong evidence for patterns of range shift and expansion, and population divergence events from the onset of the Holocene, resulting in the formation of its five distinct genetic units. The most highly supported model inferred that after the split of an ancestral population into the southern group and a larger central metapopulation lineage, the latter gave rise to the eastern and northern clusters, before finally dividing into two sub-central groups. While the use of molecular data allowed us to identify and refine the (phylo)genetic relationships of the species' lineages and populations, the use of ecological data helped us infer a past LGM refugium and the directions of post-glacial range dynamics. The time frames of these demographic events were shown to be congruent with climatic and geological events that affected the central Korean Peninsula during these periods. These findings gave us a better understanding of the consequences of past spatiotemporal factors that may have resulted in the current fragmented population distribution of this endangered plant.

Nuclear and chloroplast DNA phylogeography reveals high genetic diversity and postglacial range expansion in Quercus mexicana

PREMISE

Phylogeographical studies are fundamental for understanding factors that influence the spatial distribution of genetic lineages within species. Population expansions and contractions, distribution shifts, and climate changes are among the most important factors shaping the genetic compositions of populations.

METHODS

We investigated the phylogeography of an endemic oak, Quercus mexicana (Fagaceae), which has a restricted distribution in northeastern Mexico along the Sierra Madre Oriental and adjacent areas. Nuclear and chloroplast DNA microsatellite markers were used to describe the genetic diversity and structure of 39 populations of Q. mexicana along its entire distribution area. We tested whether population expansion or contraction events influenced the genetic diversity and structure of the species. We also modeled the historical distributional range of Q. mexicana (for the Mid Holocene, the Last Glacial Maximum, and the Last Interglacial) to estimate the extent to which climate fluctuations have impacted the distribution of this oak species.

RESULTS

Our results revealed high genetic diversity and low genetic structure in Q. mexicana populations. Ecological niche models suggested historical fluctuations in the distributional range of Q. mexicana. Historical range changes, gene flow, and physical barriers seem to have played an important role in shaping the phylogeographic structure of Q. mexicana.

CONCLUSIONS

Our study indicates that the genetic structure of Q. mexicana may have been the result of responses of oak trees not only to heterogeneous environments present in the Sierra Madre Oriental and adjacent areas, but also to elevational and latitudinal shifts in response to climate changes in the past.

Is there a massive glacial-Holocene flora continuity in Central Europe?

The prevailing paradigm about the Quaternary ecological and evolutionary history of Central European ecosystems is that they were repeatedly impoverished by regional extinctions of most species during the glacial periods, followed by massive recolonizations from southern and eastern refugia during interglacial periods. Recent literature partially contradicts this view and provides evidence to re-evaluate this Postglacial Recolonization Hypothesis and develop an alternative one. We examined the long-term history of the flora of the Carpathian (Pannonian) Basin by synthesising recent advances in ecological, phylogeographical, palaeoecological and palaeoclimatological research, and analysing the cold tolerance of the native flora of a test area (Hungary, the central part of the Carpathian Basin). We found that (1) many species have likely occurred there continuously since before the Last Glacial Maximum (LGM); (2) most of the present-day native flora (1404 species, about 80%) can occur in climates as cold as or colder than the LGM (mean annual temperature ≤+3.5°C); and (3) grasslands and forests can be species-rich under an LGM-like cold climate. These arguments support an alternative hypothesis, which we call the Flora Continuity Hypothesis. It states that long-term continuity of much of the flora in the Carpathian Basin is more plausible than regional extinctions during the LGM followed by massive postglacial recolonizations. The long-term continuity of the region's flora may have fundamental implications not only for understanding local biogeography and ecology (e.g. the temporal scale of processes), but also for conservation strategies focusing on protecting ancient species-rich ecosystems and local gene pools.

Genomic epidemiology sheds light on the emergence and spread of Mycobacterium bovis Eu2 Clonal Complex in Portugal

ABSTRACTAnimal tuberculosis (TB) remains a serious concern for animal and human health. Mycobacterium bovis circulates in multi-host systems, dominated by the European 2 clonal complex (Eu2) in Iberia. In this work, we use genomic epidemiology to infer the emergence, spread, and spatiotemporal patterns of Eu2 in the official epidemiological risk area of animal TB in Portugal. Phylogenetic analysis of 144 M. bovis whole-genome sequences from cattle, wild boar, and red deer, representing the 2002-2021 period, distinguished three Eu2 clades that evolved independently. The major Eu2 clade underwent phylodynamic inferences to estimate the time and location of outbreaks, host transitions, and spatial diffusion as well. The origin of this Eu2 clade was attributed to the red deer population in the Castelo Branco district, near the border with Spain. Most host transitions were intraspecific (80%), while interspecific transmissions between wildlife species (wild boar-red deer), and between wild boar and cattle, were highly supported. Phylogeographic reconstruction evidenced that most transitions (82%) occur within municipalities, highlighting local transmission corridors.Our study indicates that M. bovis continues to spread at the cattle-wildlife interface within the animal TB hotspot area, possibly driven by the foraging behaviour of wild boar near agricultural lands. Red deer seems to be an important driver of TB within wildlife hosts, while the wild boar links the multi-host wildlife community and livestock. This work highlights the value of combining genomic epidemiology with phylodynamic inference to resolve host jumps and spatial patterns of M. bovis, providing real-time clues about points of intervention.

Genetic diversity and microevolution in clinical Cryptococcus isolates from Cameroon

Cryptococcal meningitis is the second most common cause of death in people living with HIV/AIDS, yet we have a limited understanding of how cryptococcal isolates change over the course of infection. Cryptococcal infections are environmentally acquired, and the genetic diversity of these infecting isolates can also be geographically linked. Here, we employ whole genome sequences for 372 clinical Cryptococcus isolates from 341 patients with HIV-associated cryptococcal meningitis obtained via a large clinical trial, across both Malawi and Cameroon, to enable population genetic comparisons of isolates between countries. We see that isolates from Cameroon are highly clonal, when compared to those from Malawi, with differential rates of disruptive variants in genes with roles in DNA binding and energy use. For a subset of patients (22) from Cameroon, we leverage longitudinal sampling, with samples taken at days 7 and 14 post-enrollment, to interrogate the genetic changes that arise over the course of infection, and the genetic diversity of isolates within patients. We see disruptive variants arising over the course of infection in several genes, including the phagocytosis-regulating transcription factor GAT204. In addition, in 13% of patients sampled longitudinally, we see evidence for mixed infections. This approach identifies geographically linked genetic variation, signatures of microevolution, and evidence for mixed infections across a clinical cohort of patients affected by cryptococcal meningitis in Central Africa.

Complete genome reconstruction of the global and European regional dispersal history of the lumpy skin disease virus

IMPORTANCE

Lumpy skin disease virus (LSDV) has a complex epidemiology involving multiple strains, recombination, and vaccination. Its DNA genome provides limited genetic variation to trace outbreaks in space and time. Sequencing of LSDV whole genomes has also been patchy at global and regional scales. Here, we provide the first fine-grained whole genome sequence sampling of a constrained LSDV outbreak (southeastern Europe, 2015-2017), which we analyze along with global publicly available genomes. We formally evaluate the past occurrence of recombination events as well as the temporal signal that is required for calibrating molecular clock models and subsequently conduct a time-calibrated spatially explicit phylogeographic reconstruction. Our study further illustrates the importance of accounting for recombination events before reconstructing global and regional dynamics of DNA viruses. More LSDV whole genomes from endemic areas are needed to obtain a comprehensive understanding of global LSDV dispersal dynamics.

Speciation in Coastal Basins Driven by Staggered Headwater Captures: Dispersal of a Species Complex, Leporinus bahiensis, as Revealed by Genome-wide SNP Data

Past sea level changes and geological instability along watershed boundaries have largely influenced fish distribution across coastal basins, either by dispersal via palaeodrainages now submerged or by headwater captures, respectively. Accordingly, the South American Atlantic coast encompasses several small and isolated drainages that share a similar species composition, representing a suitable model to infer historical processes. Leporinus bahiensis is a freshwater fish species widespread along adjacent coastal basins over narrow continental shelf with no evidence of palaeodrainage connections at low sea level periods. Therefore, this study aimed to reconstruct its evolutionary history to infer the role of headwater captures in the dispersal process. To accomplish this, we employed molecular-level phylogenetic and population structure analyses based on Sanger sequences (5 genes) and genome-wide SNP data. Phylogenetic trees based on Sanger data were inconclusive, but SNPs data did support the monophyletic status of L. bahiensis. Both COI and SNP data revealed structured populations according to each hydrographic basin. Species delimitation analyses revealed from 3 (COI) to 5 (multilocus approach) MOTUs, corresponding to the sampled basins. An intricate biogeographic scenario was inferred and supported by Approximate Bayesian Computation (ABC) analysis. Specifically, a staggered pattern was revealed and characterized by sequential headwater captures from basins adjacent to upland drainages into small coastal basins at different periods. These headwater captures resulted in dispersal throughout contiguous coastal basins, followed by deep genetic divergence among lineages. To decipher such recent divergences, as herein represented by L. bahiensis populations, we used genome-wide SNPs data. Indeed, the combined use of genome-wide SNPs data and ABC method allowed us to reconstruct the evolutionary history and speciation of L. bahiensis. This framework might be useful in disentangling the diversification process in other neotropical fishes subject to a reticulate geological history.

Genetic signatures of lineage fusion closely resemble population decline

Accurate interpretation of the genetic signatures of past demographic events is crucial for reconstructing evolutionary history. Lineage fusion (complete merging, resulting in a single panmictic population) is a special case of secondary contact that is seldom considered. Here, the circumstances under which lineage fusion can be distinguished from population size constancy, growth, bottleneck, and decline were investigated. Multi-locus haplotype data were simulated under models of lineage fusion with different divergence versus sampling lag times (D:L ratios). These pseudo-observed datasets also differed in their allocation of a fixed amount of sequencing resources (number of sampled alleles, haplotype length, number of loci). Distinguishability of lineage fusion versus each of 10 untrue non-fusion scenarios was quantified based on six summary statistics (neutrality tests). Some datasets were also analyzed using extended Bayesian skyline plots. Results showed that signatures of lineage fusion very closely resemble those of decline-high distinguishability was generally limited to the most favorable scenario (D:L = 9), using the most sensitive summary statistics (F S and Z nS), coupled with the optimal sequencing resource allocation (maximizing number of loci). Also, extended Bayesian skyline plots often erroneously inferred population decline. Awareness of the potential for lineage fusion to carry the hallmarks of population decline is critical.

Evolution of horn length and lifting strength in the Japanese rhinoceros beetle Trypoxylus dichotomus

What limits the size of nature's most extreme structures? For weapons like beetle horns, one possibility is a tradeoff associated with mechanical levers: as the output arm of the lever system-the beetle horn-gets longer, it also gets weaker. This "paradox of the weakening combatant" could offset reproductive advantages of additional increases in weapon size. However, in contemporary populations of most heavily weaponed species, males with the longest weapons also tend to be the strongest, presumably because selection drove the evolution of compensatory changes to these lever systems that ameliorated the force reductions of increased weapon size. Therefore, we test for biomechanical limits by reconstructing the stages of weapon evolution, exploring whether initial increases in weapon length first led to reductions in weapon force generation that were later ameliorated through the evolution of mechanisms of mechanical compensation. We describe phylogeographic relationships among populations of a rhinoceros beetle and show that the "pitchfork" shaped head horn likely increased in length independently in the northern and southern radiations of beetles. Both increases in horn length were associated with dramatic reductions to horn lifting strength-compelling evidence for the paradox of the weakening combatant-and these initial reductions to horn strength were later ameliorated in some populations through reductions to horn length or through increases in head height (the input arm for the horn lever system). Our results reveal an exciting geographic mosaic of weapon size, weapon force, and mechanical compensation, shedding light on larger questions pertaining to the evolution of extreme structures.

Molecular Epidemiology of Rotavirus A in Calves: Evolutionary Analysis of a Bovine G8P[11] Strain and Spatio-Temporal Dynamics of G6 Lineages in the Americas

Rotavirus A (RVA) causes diarrhea in calves and frequently possesses the G6 and P[5]/P[11] genotypes, whereas G8 is less common. We aimed to compare RVA infections and G/P genotypes in beef and dairy calves from major livestock regions of Argentina, elucidate the evolutionary origin of a G8 strain and analyze the G8 lineages, infer the phylogenetic relationship of RVA field strains, and investigate the evolution and spatio-temporal dynamics of the main G6 lineages in American countries. Fecal samples (n = 422) from diarrheic (beef, 104; dairy, 137) and non-diarrheic (beef, 78; dairy, 103) calves were analyzed by ELISA and semi-nested multiplex RT-PCR. Sequencing, phylogenetic, phylodynamic, and phylogeographic analyses were performed. RVA infections were more frequent in beef (22.0%) than in dairy (14.2%) calves. Prevalent genotypes and G6 lineages were G6(IV)P[5] in beef (90.9%) and G6(III)P[11] (41.2%) or mixed genotypes (23.5%) in dairy calves. The only G8 strain was phylogenetically related to bovine and artiodactyl bovine-like strains. Re-analyses inside the G8 genotype identified G8(I) to G8(VIII) lineages. Of all G6 strains characterized, the G6(IV)P[5](I) strains from "Cuenca del Salado" (Argentina) and Uruguay clustered together. According to farm location, a clustering pattern for G6(IV)P[5] strains of beef farms was observed. Both G6 lineage strains together revealed an evolutionary rate of 1.24 × 10-3 substitutions/site/year, and the time to the most recent common ancestor was dated in 1853. The most probable ancestral locations were Argentina in 1981 for G6(III) strains and the USA in 1940 for G6(IV) strains. The highest migration rates for both G6 lineages together were from Argentina to Brazil and Uruguay. Altogether, the epidemiology, genetic diversity, and phylogeny of RVA in calves can differ according to the production system and farm location. We provide novel knowledge about the evolutionary origin of a bovine G8P[11] strain. Finally, bovine G6 strains from American countries would have originated in the USA nearly a century before its first description.

Geogenomic predictors of genetree heterogeneity explain phylogeographic and introgression history: a case study in an Amazonian bird (Thamnophilus aethiops)

Can knowledge about genome architecture inform biogeographic and phylogenetic inference? Selection, drift, recombination, and gene flow interact to produce a genomic landscape of divergence wherein patterns of differentiation and genealogy vary nonrandomly across the genomes of diverging populations. For instance, genealogical patterns that arise due to gene flow should be more likely to occur on smaller chromosomes, which experience high recombination, whereas those tracking histories of geographic isolation (reduced gene flow caused by a barrier) and divergence should be more likely to occur on larger and sex chromosomes. In Amazonia, populations of many bird species diverge and introgress across rivers, resulting in reticulated genomic signals. Herein, we used reduced representation genomic data to disentangle the evolutionary history of four populations of an Amazonian antbird, Thamnophilus aethiops, whose biogeographic history was associated with the dynamic evolution of the Madeira River Basin. Specifically, we evaluate whether a large river capture event ca. 200 Ka, gave rise to reticulated genealogies in the genome by making spatially explicit predictions about isolation and gene flow based on knowledge about genomic processes. We first estimated chromosome-level phylogenies and recovered two primary topologies across the genome. The first topology (T1) was most consistent with predictions about population divergence and was recovered for the Z chromosome. The second (T2), was consistent with predictions about gene flow upon secondary contact. To evaluate support for these topologies, we trained a convolutional neural network to classify our data into alternative diversification models and estimate demographic parameters. The best-fit model was concordant with T1 and included gene flow between non-sister taxa. Finally, we modeled levels of divergence and introgression as functions of chromosome length and found that smaller chromosomes experienced higher gene flow. Given that (1) gene-trees supporting T2 were more likely to occur on smaller chromosomes and (2) we found lower levels of introgression on larger chromosomes (and especially the Z-chromosome), we argue that T1 represents the history of population divergence across rivers and T2 the history of secondary contact due to barrier loss. Our results suggest that a significant portion of genomic heterogeneity arises due to extrinsic biogeographic processes such as river capture interacting with intrinsic processes associated with genome architecture. Future phylogeographic studies would benefit from accounting for genomic processes, as different parts of the genome reveal contrasting, albeit complementary histories, all of which are relevant for disentangling the intricate geogenomic mechanisms of biotic diversification.

The Temporal and Geographical Dynamics of Potato Virus Y Diversity in Russia

Potato virus Y, an important viral pathogen of potato, has several genetic variants and geographic distributions which could be affected by environmental factors, aphid vectors, and reservoir plants. PVY is transmitted to virus-free potato plants by aphids and passed on to the next vegetative generations through tubers, but the effects of tuber transmission in PVY is largely unknown. By using high-throughput sequencing, we investigated PVY populations transmitted to potato plants by aphids in different climate zones of Russia, namely the Moscow and Astrakhan regions. We analyzed sprouts from the tubers produced by field-infected plants to investigate the impact of tuber transmission on PVY genetics. We found a significantly higher diversity of PVY isolates in the Astrakhan region, where winters are shorter and milder and summers are warmer compared to the Moscow region. While five PVY types, NTNa, NTNb, N:O, N-Wi, and SYR-I, were present in both regions, SYRI-II, SYRI-III, and 261-4 were only found in the Astrakhan region. All these recombinants were composed of the genome sections derived from PVY types O and N, but no full-length sequences of such types were present. The composition of the PVY variants in the tuber sprouts was not always the same as in their parental plants, suggesting that tuber transmission impacts PVY genetics.

Phylodynamic approaches to studying avian influenza virus

Avian influenza viruses can cause severe disease in domestic and wild birds and are a pandemic threat. Phylodynamics is the study of how epidemiological, evolutionary, and immunological processes can interact to shape viral phylogenies. This review summarizes how phylodynamic methods have and could contribute to the study of avian influenza viruses. Specifically, we assess how phylodynamics can be used to examine viral spread within and between wild or domestic bird populations at various geographical scales, identify factors associated with virus dispersal, and determine the order and timing of virus lineage movement between geographic regions or poultry production systems. We discuss factors that can complicate the interpretation of phylodynamic results and identify how future methodological developments could contribute to improved control of the virus.

Population genetic structure of a Chihuahuan Desert endemic mammal, the desert pocket gopher, Geomys arenarius

The biogeographic history of the Chihuahuan Desert is complex, driven by numerous physiographic events and climatic changes. This dynamic history would have influenced the flora and fauna of the region including the desert pocket gopher, Geomys arenarius, a subterranean rodent endemic to the northern Chihuahuan Desert. G. arenarius is restricted to sandy soils and are considered to have a disjunct distribution. Two subspecies are recognized: G. a. arenarius and G. a. brevirostris. We used multilocus nuclear (amplified fragment length polymorphisms) and mitochondrial DNA (ND2) sequence data to uncover patterns of genetic diversity within and among populations of G. arenarius. We evaluated correspondence of genetic patterns to traditionally accepted subspecies boundaries, mapped the distribution of potentially suitable soils to identify barriers or corridors to dispersal and to guide future survey efforts, provided evidence that could be used to recognize distinct population segments, and quantified genetic diversity within populations. Both datasets were largely concordant and demonstrated hierarchical patterns of genetic divergence. The greatest divergence was consistent with the two recognized subspecies. Mapping of potentially habitable soils revealed likely barriers to dispersal contributing to the allopatric pattern of geographic distribution and areas, which may be occupied by G. arenarius but not yet documented. Because G. arenarius is restricted to soils with high sand content, and these habitable soils are disjunct within the region occupied by this species, historical factors that impacted soil deposition and deflation likely contributed to the observed patterns of genetic divergence. Genetic diversity was higher within populations of the southern subspecies (G. a. arenarius) compared to G. a. brevirostris. This may be due to a greater availability of continuous suitable soils within the range of G. a. arenarius or higher density due to greater food availability (currently or historically)-both of which could allow for a higher effective population size.

Loss of species and genetic diversity during colonization: Insights from acanthocephalan parasites in northern European seals

Studies on host-parasite systems that have experienced distributional shifts, range fragmentation, and population declines in the past can provide information regarding how parasite community richness and genetic diversity will change as a result of anthropogenic environmental changes in the future. Here, we studied how sequential postglacial colonization, shifts in habitat, and reduced host population sizes have influenced species richness and genetic diversity of Corynosoma (Acanthocephala: Polymorphidae) parasites in northern European marine, brackish, and freshwater seal populations. We collected Corynosoma population samples from Arctic, Baltic, Ladoga, and Saimaa ringed seal subspecies and Baltic gray seals, and then applied COI barcoding and triple-enzyme restriction-site associated DNA (3RAD) sequencing to delimit species, clarify their distributions and community structures, and elucidate patterns of intraspecific gene flow and genetic diversity. Our results showed that Corynosoma species diversity reflected host colonization histories and population sizes, with four species being present in the Arctic, three in the Baltic Sea, two in Lake Ladoga, and only one in Lake Saimaa. We found statistically significant population-genetic differentiation within all three Corynosoma species that occur in more than one seal (sub)species. Genetic diversity tended to be high in Corynosoma populations originating from Arctic ringed seals and low in the landlocked populations. Our results indicate that acanthocephalan communities in landlocked seal populations are impoverished with respect to both species and intraspecific genetic diversity. Interestingly, the loss of genetic diversity within Corynosoma species seems to have been less drastic than in their seal hosts, possibly due to their large local effective population sizes resulting from high infection intensities and effective intra-host population mixing. Our study highlights the utility of genomic methods in investigations of community composition and genetic diversity of understudied parasites.

Two Lineages of Oribatid Mites Morphologically Correspond to the Circumpolar Species Ameronothrus nigrofemoratus (Acari, Oribatida) but Differ Genetically as Distinct Species Are Revealed on the Kolguev Island

Data reported from Northern Canada were until recently the only available data on the genetic characteristics of the oribatid mite Ameronothrus nigrofemoratus, which has a circumpolar distribution on the coasts of Arctic seas. A partial cytochrome oxidase I gene (COI) mtDNA sequence was examined in mites morphologically assigned to this species from the Kolguev Island. Two highly divergent phylogenetic lineages of A. nigrofemoratus (7% divergence) were revealed, neither of which was found on the Canadian coast. Four COI amino acid substitutions distinguished one of the lineages from North American A. nigrofemoratus, corresponding to the degree of difference between A. nigrofemoratus and its sister species A. lineatus.

Evolutionary dynamics of Euphorbia carniolica suggest a complex Plio-Pleistocene history of understorey species of deciduous forest in southeastern Europe

Deciduous forests form the dominant natural vegetation of Europe today, but were restricted to small refugia during Pleistocene cold stages, implying an evolutionary past shaped by recurrent range contractions and expansions. Cold-stage forest refugia were probably widespread in southern and central Europe, with the northwestern Balkan Peninsula being of particular importance. However, the actual number and location of deciduous forest refugia, as well as the connections between them, remain disputed. Here, we address the evolutionary dynamics of the deciduous forest understorey species Euphorbia carniolica as a proxy for past forest dynamics. To do so, we obtained genomic and morphometric data from populations representing the species' entire range, investigated phylogenetic position and intraspecific genetic variation, tested explicit demographic scenarios and applied species distribution models. Our data support two disjoint groups linked to separate refugia on the northwestern and central Balkan Peninsula. We find that genetic differentiation between groups started in the early Pleistocene via vicariance, suggesting a larger distribution in the past. Both refugia acted as sources for founder events to the southeastern Alps and the Carpathians; the latter were likely colonised before the last cold stage. In line with traditional views on the pre-Pleistocene origin of many southeastern European deciduous forest species, the origin of E. carniolica was dated to the late Pliocene. The fact that E. carniolica evolved at a time when a period of continuous forestation was ending in much of Eurasia provides an interesting biogeographical perspective on the past links between Eurasian deciduous forests and their biota.

Multiple introductions and country-wide spread of DENV-2 genotype II (Cosmopolitan) in Brazil

Dengue virus serotype 2, genotype Cosmopolitan (DENV-2-GII), is one of the most widespread DENV strains globally. In the USA, DENV-2 epidemics have been dominated by DENV-2 genotype Asian-American (DENV-2-GIII), and the first cases of DENV-2-GII were only described in 2019, in Peru, and in 2021 in Brazil. To gain new information about the circulation of DENV-2-GII in Brazil, we sequenced 237 DENV-2 confirmed cases sampled between March 2021 and March 2023 and revealed that DENV-2-GII is already present in all geographic regions of Brazil. The phylogeographic analysis inferred that DENV-2-GII was introduced at least four times in Brazil, between May 2020 and August 2022, generating multiple clades that spread throughout the country with different success. Despite multiple introductions of DENV-2-GII, analysis of the country-wide laboratory surveillance data showed that the Brazilian dengue epidemic in 2022 was dominated by DENV-1 in most states. We hypothesize that massive circulation of DENV-2-GIII in previous years in Brazil might have created a population immune barrier against symptomatic homotypic reinfections by DENV-2-GII, leading to sustained cryptic circulation in asymptomatic cases and localized outbreaks of this new genotype. In summary, our study stresses the importance of arboviral genomic surveillance to close monitoring and better understanding the potential impact of DENV-2-GII in the coming years.

A multifaceted interplay between virulence, drug resistance, and the phylogeographic landscape of Mycobacterium tuberculosis

Latin-American Mediterranean (LAM) family is one of the most significant and global genotypes of Mycobacterium tuberculosis. Here, we used the murine model to study the virulence and lethality of the genetically and epidemiologically distinct LAM strains. The pathobiological characteristics of the four LAM strains (three drug resistant and one drug susceptible) and the susceptible reference strain H37Rv were studied in the C57BL/6 mouse model. The whole-genome sequencing was performed using the HiSeq Illumina platform, followed by bioinformatics and phylogenetic analysis. The susceptible strain H37Rv showed the highest virulence. Drug-susceptible LAM strain (spoligotype SIT264) was more virulent than three multidrug-resistant (MDR) strains (SIT252, SIT254, and SIT266). All three MDR isolates were low lethal, while the susceptible isolate and H37Rv were moderately/highly lethal. Putting the genomic, phenotypic, and virulence features of the LAM strains/spoligotypes in the context of their dynamic phylogeography over 20 years reveals three types of relationships between virulence, resistance, and transmission. First, the most virulent and more lethal drug-susceptible SIT264 increased its circulation in parts of Russia. Second, moderately virulent and pre-XDR SIT266 was prevalent in Belarus and continues to be visible in North-West Russia. Third, the low virulent and MDR strain SIT252 previously considered as emerging has disappeared from the population. These findings suggest that strain virulence impacts the transmission, irrespective of drug resistance properties. The increasing circulation of susceptible but more virulent and lethal strains implies that personalized TB treatment should consider not only resistance but also the virulence of the infecting M. tuberculosis strains. IMPORTANCE The study is multidisciplinary and investigates the epidemically/clinically important and global lineage of Mycobacterium tuberculosis, named Latin-American-Mediterranean (LAM), yet insufficiently studied with regard to its pathobiology. We studied different LAM strains (epidemic vs endemic and resistant vs susceptible) in the murine model and using whole-genome analysis. We also collected long-term, 20-year data on their prevalence in Eurasia. The findings are both expected and unexpected. (i) We observe that a drug-susceptible but highly virulent strain increased its prevalence. (ii) By contrast, the multidrug-resistant (MDR) but low-virulent, low-lethal strain (that we considered as emerging 15 years ago) has almost disappeared. (iii) Finally, an intermediate case is the MDR strain with moderate virulence that continues to circulate. We conclude that (i) the former and latter strains are the most hazardous and require close epidemiological monitoring, and (ii) personalized TB treatment should consider not only drug resistance but also the virulence of the infecting strains and development of anti-virulence drugs is warranted.

Genetic features of the Puumala virus (Hantaviridae: Orthohantavirus) identified in the Moscow region

INTRODUCTION

Puumala virus (family Hantaviridae, genus Orthohantavirus) is distributed in most regions of the European part of Russia. However, information about its genetic variants circulating on the territory of the Central Federal District is extremely scarce.

MATERIALS AND METHODS

Rodents' tissue samples were tested after reverse transcription by PCR for the presence of hantaviral RNA. The amplified fragments of the L segment were sequenced by the Sanger method. For two samples, sequences of all three segments were obtained using the NGS method. Phylogenetic trees were built in the MEGA-X software.

RESULTS

Puumala virus was found in six samples. Based on the phylogenetic analysis of sequences of three segments, the obtained genetic variants belong to the sublineage previously designated as W-RUS.

CONCLUSION

A genetic variant of the Puumala virus, belonging to the subline W-RUS, circulates on the territory of the Volokolamsk district of Moscow region.

Effects of Climatic Change on Phylogeography and Ecological Niche of the Endemic Herb Elymus breviaristatus on the Qinghai-Tibet Plateau

Past climatic and topographic variations have created strong biogeographic barriers for alpine species and are key drivers of the distribution of genetic variation and population dynamics of species on the Qinghai-Tibet Plateau (QTP). Therefore, to better conserve and use germplasm resources, it is crucial to understand the distribution and differentiation of genetic variation within species. Elymus breviaristatus, an ecologically important rare grass species with strong resistance, is restricted to a limited area of the QTP. In this study, we investigated the phylogeography of E. breviaristatus using five chloroplast genes and spacer regions in natural populations distributed along the eastern QTP. We identified a total of 25 haplotypes among 216 individuals from 18 E. breviaristatus populations, which were further classified into four haplogroups based on geographical distribution and haplotype network analysis. Notably, we did not observe any signs of population expansion. High genetic diversity was exhibited at both species and population levels, with precipitation being the main limiting factor for population genetic diversity levels. Higher genetic diversity was exhibited by populations located near the Mekong-Salween Divide genetic barrier, suggesting that they may have served as a glacial refuge. The significant pattern of genetic differentiation by environmental isolation highlights the influence of heterogeneous environments on the genetic structure of E. breviaristatus populations. Additionally, the results of ecological niche models indicated that the geographic distribution of E. breviaristatus populations has decreased rapidly since the Last Glacial Maximum but is not threatened by future global warming.

The Role of the Hercynian Mountains of Central Europe in Shaping Plant Migration Patterns in the Pleistocene-A Review

The climatic changes that took place in Europe during the Quaternary period influenced plant habitats as well as their species and vegetation composition. In this article, biogeographical studies on Hercynian mountain plants that include data for the Alps, Carpathians, and European lowlands are reviewed in order to discuss the phylogeographical structure and divergence of the Hercynian populations from those in other European mountain ranges, Scandinavia, and lowlands. The analyzed studies show specific phylogeographical relations between the Hercynian mountains, Alps, Scandinavia, Carpathians, and European lowlands. The results also indicate that the genetic patterns of plant populations in the Hercynian Mountains may differ significantly in terms of origin. The main migration routes of species to the Hercynian ranges began in the Alps or Carpathians. Some species, such as Rubus chamaemorus L., Salix lapponum L., and Salix herbacea L., are glacial relics that may have arrived and settled in the Hercynian Mountains during the Ice Age and that survived in isolated habitats. The Hercynian Mountains are composed of various smaller mountain ranges and are a crossroads of migration routes from different parts of Europe; thus, intensive hybridization has occurred between the plant populations therein, which is indicated by the presence of several divergent genetic lines.

Uncovering the endemic circulation of rabies in Cambodia

In epidemiology, endemicity characterizes sustained pathogen circulation in a geographical area, which involves a circulation that is not being maintained by external introductions. Because it could potentially shape the design of public health interventions, there is an interest in fully uncovering the endemic pattern of a disease. Here, we use a phylogeographic approach to investigate the endemic signature of rabies virus (RABV) circulation in Cambodia. Cambodia is located in one of the most affected regions by rabies in the world, but RABV circulation between and within Southeast Asian countries remains understudied. Our analyses are based on a new comprehensive data set of 199 RABV genomes collected between 2014 and 2017 as well as previously published Southeast Asian RABV sequences. We show that most Cambodian sequences belong to a distinct clade that has been circulating almost exclusively in Cambodia. Our results thus point towards rabies circulation in Cambodia that does not rely on external introductions. We further characterize within-Cambodia RABV circulation by estimating lineage dispersal metrics that appear to be similar to other settings, and by performing landscape phylogeographic analyses to investigate environmental factors impacting the dispersal dynamic of viral lineages. The latter analyses do not lead to the identification of environmental variables that would be associated with the heterogeneity of viral lineage dispersal velocities, which calls for a better understanding of local dog ecology and further investigations of the potential drivers of RABV spread in the region. Overall, our study illustrates how phylogeographic investigations can be performed to assess and characterize viral endemicity in a context of relatively limited data.

Genomic analysis reveals complex population structure within the smooth newt, Lissotriton vulgaris, in Central Europe

Species with wide-range distributions usually display high genetic variation. This variation can be partly explained by historical lineages that were temporally isolated from each other and are back into secondary reproductive contact, and partly by local adaptations. The smooth newt (Lissotriton vulgaris) is one of the most widely distributed amphibians species across Eurasia and forms a species complex with a partially overlapping distribution and morphology. In the present study, we explored the population genomic structure of smooth newt lineages in the Carpathian Basin (CB) relying on single-nucleotide polymorphisms. Our dataset included new and previously published data to study the secondary contact zone between lineages in the CB and also tested for the barrier effect of rivers to gene flow between these lineages. We confirmed the presence of the South L. v. vulgaris Lineage distributed in Transdanubia and we provided new distribution records of L. v. ampelensis inhabiting the eastern territories of the CB. High genetic diversity of smooth newts was observed, especially in the North Hungarian Mountains and at the interfluves of the main rivers in the South with four distinct lineages of L. v. vulgaris and one lineage of L. v. ampelensis showing a low level of admixture with the spatially closest L. v. vulgaris lineage. Moreover, admixture detected at the interfluve of the main rivers (i.e. Danube and Tisza) suggested a secondary contact zone in the area. Finally, we found that the river Danube has a very weak effect on population divergence, while the river Tisza is a geographical barrier limiting gene flow between smooth newt lineages. As the range boundaries of L. v. vulgaris and L. v. ampelensis in the CB coincide with the river Tisza, our study underpins the influence of rivers in lineage diversification.

Parsing variance by marker type: Testing biogeographic hypotheses and differential contribution of historical processes to population structure in a desert lizard

A fundamental goal of population genetic studies is to identify historical biogeographic patterns and understand the processes that generate them. However, localized demographic events can skew population genetic inference. Assessing populations with multiple types of genetic markers, each with unique mutation rates and responses to changes in population size, can help to identify potentially confounding population-specific demographic processes. Here, we compared population structure and connectivity inferred from microsatellites and restriction site-associated DNA loci among 17 populations of an arid-specialist lizard, the desert night lizard, Xantusia vigilis, in central California to test among historical processes structuring population genetic diversity. We found that both marker types yielded generally concordant insights into population genetic structure including a major phylogenetic break maintained between two populations separated by less than 10 km, suggesting that either marker type could be used to understand generalized demographic patterns across the region for management purposes. However, we also found that the effects of demography on marker discordance could be used to elucidate population histories and distinguish among competing biogeographic hypotheses. Our results suggest that comparisons of within-population diversity across marker types provide powerful opportunities for leveraging marker discordance, particularly for understanding the creation and maintenance of contact zones among clades.

A Customized Human Mitochondrial DNA Database (hMITO DB v1.0) for Rapid Sequence Analysis, Haplotyping and Geo-Mapping

The field of mitochondrial genomics has advanced rapidly and has revolutionized disciplines such as molecular anthropology, population genetics, and medical genetics/oncogenetics. However, mtDNA next-generation sequencing (NGS) analysis for matrilineal haplotyping and phylogeographic inference remains hindered by the lack of a consolidated mitogenome database and an efficient bioinformatics pipeline. To address this, we developed a customized human mitogenome database (hMITO DB) embedded in a CLC Genomics workflow for read mapping, variant analysis, haplotyping, and geo-mapping. The database was constructed from 4286 mitogenomes. The macro-haplogroup (A to Z) distribution and representative phylogenetic tree were found to be consistent with published literature. The hMITO DB automated workflow was tested using mtDNA-NGS sequences derived from Pap smears and cervical cancer cell lines. The auto-generated read mapping, variants track, and table of haplotypes and geo-origins were completed in 15 min for 47 samples. The mtDNA workflow proved to be a rapid, efficient, and accurate means of sequence analysis for translational mitogenomics.

Insights on the Evolutionary History and Genetic Patterns of Octopus vulgaris Cuvier, 1797 in the Northeastern Atlantic Using Mitochondrial DNA

Octopus vulgaris is one of the most harvested octopus species in the world. In the Iberian Peninsula, there are several small-scale fisheries that have a long-term tradition of harvesting octopus. The Asturias fleet (in Northern Spain) has an internationally recognized MSC label for its exploitation. Of concern, genetic assessments of exploited stocks are currently scarce, which could prevent the implementation of adequate managing strategies. We use two mitochondrial regions (cytochrome oxidase subunit 1 and control region) to analyze the genetic status and evolutionary events that conditioned octopus populations' characteristics in the Northeastern Atlantic. A total of 90 individuals were sampled from three different localities in the Iberian Peninsula as well as a location in Macaronesia. Temporal genetic analyses on Asturias and Algarve populations were also performed. Results indicated the absence of fine spatial genetic structuring but showed the Canary Islands (in Macaronesia) as the most distinct population. Our analyses detected two distinct clades, already described in the literature, but, for the first time, we confirmed the presence of the α-southern haplogroup in the Northern Iberian Peninsula. This result indicates a more continuous cline for the distribution of these two haplogroups than previously reported. Temporal changes in the distribution of both haplogroups in contact zones were also detected.

Deep Sequencing to Reveal Phylo-Geographic Relationships of Juquitiba Virus in Paraguay

Several hantaviruses result in zoonotic infections of significant public health concern, causing hemorrhagic fever with renal syndrome (HFRS) or hantavirus cardiopulmonary syndrome (HCPS) in the Old and New World, respectively. Given a 35% case fatality rate, disease-causing New World hantaviruses require a greater understanding of their biology, genetic diversity, and geographical distribution. Juquitiba hantaviruses have been identified in Oligoryzomys nigripes in Brazil, Paraguay, and Uruguay. Brazil has reported the most HCPS cases associated with this virus. We used a multiplexed, amplicon-based PCR strategy to screen and deep-sequence the virus harbored within lung tissues collected from Oligoryzomys species during rodent field collections in southern (Itapúa) and western (Boquerón) Paraguay. No Juquitiba-like hantaviruses were identified in Boquerón. Herein, we report the full-length S and M segments of the Juquitiba hantaviruses identified in Paraguay from O. nigripes. We also report the phylogenetic relationships of the Juquitiba hantaviruses in rodents collected from Itapúa with those previously collected in Canindeyú. We showed, using the TN93 nucleotide substitution model, the coalescent (constant-size) population tree model, and Bayesian inference implemented in the Bayesian evolutionary analysis by sampling trees (BEAST) framework, that the Juquitiba virus lineage in Itapúa is distinct from that in Canindeyú. Our spatiotemporal analysis showed significantly different time to the most recent ancestor (TMRA) estimates between the M and S segments, but a common geographic origin. Our estimates suggest the additional geographic diversity of the Juquitiba virus within the Interior Atlantic Forest and highlight the need for more extensive sampling across this biome.

Helmet Shape and Phylogeography of the Treehopper Membracis mexicana

Membracis mexicana (Hemiptera: Membracidae) is distributed in four biogeographic provinces of Mexico. Field observations indicate that there are different forms of this species, but the distribution of the phenotype and the genetic variation of this species have not been clarified. The aim of this study was to quantify the phenotypic and genetic variation of M. mexicana and determine whether the configuration of biogeographic provinces impacts the distribution of this variation. To achieve this, we analyzed 307 photographs using 19 landmarks and geometric morphometrics to quantify the phenotypic variation in helmets. We sequenced five molecular markers for 205 individuals to describe the phylogeographic pattern. As a result, we identified three morphological configurations of the helmet of M. mexicana and two genetic lineages. The morphotypes are (1) a large and wide helmet with small dorsal spots, (2) a small and narrow helmet with large dorsal spots, and (3) a small and narrow helmet with small spots. Genetic lineages are distributed in southeast and western Mexico. The western lineage corresponds to two helmet morphotypes (1 and 2) and the southeast lineage to morphotype 3. We found that the larger helmets correspond to the western lineage and are distributed in Trans-Mexican Volcanic Belt and Pacific lowlands provinces, whereas the smallest helmets correspond to the southeast lineage and are present in the Veracruzan and Yucatan Peninsula provinces.

Systematics, biogeography and phylogenomics of northern bog lemmings (Cricetidae), cold-temperate rodents of conservation concern under global change

Northern bog lemmings, Mictomys (Synaptomys) borealis, are currently being assessed for protections under the U.S. Endangered Species Act. A major impediment to comprehensive evaluation is a deficiency of data towards understanding the biology of these rodents. Inherent rarity and scarce specimen sampling, despite a continent-wide distribution, has precluded our ability to implement modern methods for resolving taxonomy, evolutionary history, and investigating multiple other species traits. Here we use a maternally inherited locus (mitochondrial cytochrome b) and between 5,939 and 11,513 nuclear loci from reduced representation sequencing (ddRADseq) to investigate the evolutionary history of northern bog lemmings. We 1) qualify evidence based on morphological and early molecular studies for the genus assignment of Mictomys, 2) test the validity of multiple sub-species designations, 3) provide spatial and temporal historical biogeographic perspectives, and 4) discuss how incomplete sampling might influence conservation efforts. Both mitochondrial and nuclear datasets exhibit deep divergence and paraphyly between two recognized species, the northern (Mictomys borealis) and southern (Synaptomys cooperi) bog lemmings. Based on mtDNA, the geographically isolated subspecies (M. b. sphagnicola) was found to be divergent from all other specimens. The remainder of the species exhibited shallow intra-specific differentiation in mtDNA, however nuclear data supports genetic distinction consistent with four geographic subspecies. Recent coalescence of all northern bog lemmings (except for M. b. sphagnicola) reflects divergence in multiple refugia through the last glacial cycle, including a well-known coastal center of endemism and multiple regions south of continental ice-sheets. Regional lineages across North America suggest strong latitudinal displacement with global climate change, coupled with isolation-reconnection dynamics. This taxon suffers from a lack of modern samples through most of its distribution, severely limiting interpretation of ongoing evolutionary processes, particularly in southern portions of the species' range. Limited voucher specimen sampling of vulnerable populations could aid in rigorous conservation decision-making.

Population dynamics in newts of the Carpathian Mountains

Rarity, range restriction, and narrow endemism tend to carry dire and urgent conservation implications for imperilled species. What is also clear is that human-associated extinction risk factors such as urbanization and deforestation pose overwhelming threats to range-restricted species. In this issue of Molecular Ecology, Antunes et al. (2022) demonstrate that these threats can also impact widespread species. By comparing newts in the genus Lissotriton that co-occur in the same geographical region, they expose the distinctness of risks facing species with different habitat preferences. Their study emphasizes the importance of local-scale landscape genetics to reveal the nuances of population connectivity that might otherwise be missed by studying a broader spatial scale.

A temporal and spatial study of genetic structure in four species of bladed Bangiales (Rhodophyta) from the southeastern Pacific coast of Chile

The coastline is a heterogeneous and highly dynamic environment influenced by abiotic and biotic variables affecting the temporal stability of genetic diversity and structure of marine organisms. The aim of this study was to determine how much the genetic structure of four species of marine Bangiales vary in time and space. Partial sequences of the cytochrome oxidase I (COI) gene obtained from two Pyropia (Py. sp. CHJ and Py. orbicularis) and two Porphyra (P. mumfordii and P. sp. FIH) species were used to compare the effect of the 40° S/41° S biogeographic break (spatial-regional scale) and the one of the Valdivia River discharges (spatial-local scale) and determine their temporal stability. Four seasonal samplings were taken during 1 year at five sites, one site located in Melinka (Magallanes province) and four sites along the coast of Valdivia (Intermediate area), on both sides of the river mouth. Results showed a strong genetic spatial structure at regional scale (ΦST > 0.4) in Py. sp. CHJ, Py. orbicularis, and P. mumfordii, congruent with the 41° S/42° S biogeographic break. A potential barrier to gene flow, related to the Valdivia River discharge, was detected only in P. mumfordii. In P. sp. FIH, spatial genetic structure was not detected at any scale. The genetic structure of all four species is stable throughout the year. The potential effect of main currents and river discharge in limiting the transport of Bangiales spores are discussed. We propose that both a restricted propagule dispersal and the formation potential for persistent banks of microscopic stages could lead to a temporally stable spatial partitioning of genetic variation in bladed Bangiales.

Phylogeography of the Sunda pangolin, Manis javanica: Implications for taxonomy, conservation management and wildlife forensics

The Sunda pangolin (Manis javanica) is the most widely distributed Asian pangolin species, occurring across much of Southeast Asia and in southern China. It is classified as Critically Endangered and is one of the most trafficked mammals in the world, which not only negatively impacts wild Sunda pangolin populations but also poses a potential disease risk to other species, including humans and livestock. Here, we aimed to investigate the species' phylogeography across its distribution to improve our understanding of the species' evolutionary history, elucidate any taxonomic uncertainties and enhance the species' conservation genetic management and potential wildlife forensics applications. We sequenced mtDNA genomes from 23 wild Sunda pangolins of known provenance originating from Malaysia to fill sampling gaps in previous studies, particularly in Borneo. To conduct phylogenetic and population genetic analyses of Sunda pangolins across their range, we integrated these newly generated mitochondrial genomes with previously generated mtDNA and nuclear DNA data sets (RAD-seq SNP data). We identified an evolutionarily distinct mtDNA lineage in north Borneo, estimated to be ~1.6 million years divergent from lineages in west/south Borneo and the mainland, comparable to the divergence time from the Palawan pangolin. There appeared to be mitonuclear discordance, with no apparent genetic structure across Borneo based on analysis of nuclear SNPs. These findings are consistent with the 'out of Borneo hypothesis', whereby Sunda pangolins diversified in Borneo before subsequently migrating throughout Sundaland, and/or a secondary contact scenario between mainland and Borneo. We have elucidated possible taxonomic issues in the Sunda/Palawan pangolin complex and highlight the critical need for additional georeferenced samples to accurately apportion its range-wide genetic variation into appropriate taxonomic and conservation units. Additionally, these data have improved forensic identification testing involving these species and permit the implementation of geographic provenance testing in some scenarios.

Variable rates of hybridization among contact zones between a pair of topminnow species, Fundulus notatus and F. olivaceus

Pairs of species that exhibit broadly overlapping distributions, and multiple geographically isolated contact zones, provide opportunities to investigate the mechanisms of reproductive isolation. Such naturally replicated systems have demonstrated that hybridization rates can vary substantially among populations, raising important questions about the genetic basis of reproductive isolation. The topminnows, Fundulus notatus and F. olivaceus, are reciprocally monophyletic, and co-occur in drainages throughout much of the central and southern United States. Hybridization rates vary substantially among populations in isolated drainage systems. We employed genome-wide sampling to investigate geographic variation in hybridization, and to assess the possible importance of chromosome fusions to reproductive isolation among nine separate contact zones. The species differ by chromosomal rearrangements resulting from Robertsonian (Rb) fusions, so we hypothesized that Rb fusion chromosomes would serve as reproductive barriers, exhibiting steeper genomic clines than the rest of the genome. We observed variation in hybridization dynamics among drainages that ranged from nearly random mating to complete absence of hybridization. Contrary to predictions, our use of genomic cline analyses on mapped species-diagnostic SNP markers did not indicate consistent patterns of variable introgression across linkage groups, or an association between Rb fusions and genomic clines that would be indicative of reproductive isolation. We did observe a relationship between hybridization rates and population phylogeography, with the lowest rates of hybridization tending to be found in populations inferred to have had the longest histories of drainage sympatry. Our results, combined with previous studies of contact zones between the species, support population history as an important factor in explaining variation in hybridization rates.

An elevational shift facilitated the Mesoamerican diversification of Azure-hooded Jays (Cyanolyca cucullata) during the Great American Biotic Interchange

The Great American Biotic Interchange (GABI) was a key biogeographic event in the history of the Americas. The rising of the Panamanian land bridge ended the isolation of South America and ushered in a period of dispersal, mass extinction, and new community assemblages, which sparked competition, adaptation, and speciation. Diversification across many bird groups, and the elevational zonation of others, ties back to events triggered by the GABI. But the exact timing of these events is still being revealed, with recent studies suggesting a much earlier time window for faunal exchange, perhaps as early as 20 million years ago (Mya). Using a time-calibrated phylogenetic tree, we show that the jay genus Cyanolyca is emblematic of bird dispersal trends, with an early, pre-land bridge dispersal from Mesoamerica to South America 6.3-7.3 Mya, followed by a back-colonization of C. cucullata to Mesoamerica 2.3-4.8 Mya, likely after the land bridge was complete. As Cyanolyca species came into contact in Mesoamerica, they avoided competition due to a prior shift to lower elevation in the ancestor of C. cucullata. This shift allowed C. cucullata to integrate itself into the Mesoamerican highland avifauna, which our time-calibrated phylogeny suggests was already populated by higher-elevation, congeneric dwarf-jays (C. argentigula, C. pumilo, C. mirabilis, and C. nanus). The outcome of these events and fortuitous elevational zonation was that C. cucullata could continue colonizing new highland areas farther north during the Pleistocene. Resultingly, four C. cucullata lineages became isolated in allopatric, highland regions from Panama to Mexico, diverging in genetics, morphology, plumage, and vocalizations. At least two of these lineages are best described as species (C. mitrata and C. cucullata). Continued study will further document the influence of the GABI and help clarify how dispersal and vicariance shaped modern-day species assemblages in the Americas.

Impacts of Quaternary glaciation, geological history and geography on animal species history in continental East Asia: A phylogeographic review

Continental East Asia has a mild Pleistocene climate and a complex recent geological history. Phylogeographic studies of animals over the last 30 years have produced several distinctive patterns. Glaciation refugia are numerous and are not restricted to any particular regions. Most of them are localized and species-specific, although several large refugia, for example the mountains of SW China, are shared by multiple species and have refugia-within-refugia. Furthermore, postglaciation range expansion events vary greatly in time, scale and direction. Large-scale south-to-north post-LGM expansions are few and mostly occurred in the northern regions. Additionally, several unique geographic features, including the three-step terrain of China and the northern arid belt, have significant impacts on many species histories. Overall, the impacts of Pleistocene glaciations, particularly the LGM, on species history vary drastically from nondetectable to significant. The impacts are the least for species from the southwestern region and are most dominant for species from the north. Geological events play a more significant role in shaping species history than Pleistocene climatic changes. Phylogeographic patterns among animals species are highly consistent with those of plants. Future phylogeographic endeavour in East Asia should be hypothesis-driven and seek processes that underlie common patterns. The wide use of genomic data allow accurate estimates of historical population processes and exploration of older history beyond the Pleistocene.

Mitochondrial DNA in Human Diversity and Health: From the Golden Age to the Omics Era

Mitochondrial DNA (mtDNA) is a small fraction of our hereditary material. However, this molecule has had an overwhelming presence in scientific research for decades until the arrival of high-throughput studies. Several appealing properties justify the application of mtDNA to understand how human populations are-from a genetic perspective-and how individuals exhibit phenotypes of biomedical importance. Here, I review the basics of mitochondrial studies with a focus on the dawn of the field, analysis methods and the connection between two sides of mitochondrial genetics: anthropological and biomedical. The particularities of mtDNA, with respect to inheritance pattern, evolutionary rate and dependence on the nuclear genome, explain the challenges of associating mtDNA composition and diseases. Finally, I consider the relevance of this single locus in the context of omics research. The present work may serve as a tribute to a tool that has provided important insights into the past and present of humankind.

Dispersal patterns and influence of air travel during the global expansion of SARS-CoV-2 variants of concern

The Alpha, Beta, and Gamma SARS-CoV-2 variants of concern (VOCs) co-circulated globally during 2020 and 2021, fueling waves of infections. They were displaced by Delta during a third wave worldwide in 2021, which, in turn, was displaced by Omicron in late 2021. In this study, we use phylogenetic and phylogeographic methods to reconstruct the dispersal patterns of VOCs worldwide. We find that source-sink dynamics varied substantially by VOC and identify countries that acted as global and regional hubs of dissemination. We demonstrate the declining role of presumed origin countries of VOCs in their global dispersal, estimating that India contributed <15% of Delta exports and South Africa <1%-2% of Omicron dispersal. We estimate that >80 countries had received introductions of Omicron within 100 days of its emergence, associated with accelerated passenger air travel and higher transmissibility. Our study highlights the rapid dispersal of highly transmissible variants, with implications for genomic surveillance along the hierarchical airline network.

Complete chloroplast genome of two nutmeg species, Myristica argentea Warb. 1891 and Myristica fatua Houtt. 1774 (Myristicaceae)

Myristica argentea Warb. 1891 and M. fatua Houtt. 1774 are two South-East Asian food tree species. They are harvested from the wild or cultivated for local uses as a condiment (nutmeg and mace), medicine, and source of wood. In this study, we reconstructed the complete chloroplast (cp) genomes of these two species from whole genome sequencing data using the Illumina NovaSeq platform. The genome sizes of M. argentea and M. fatua were respectively 155,871 base pairs (bp) and 155,898 bp, including 126 genes and an overall GC content of 39.20% in both species. Our study provides useful resources for future evolutionary research and diversity analysis of Myristica species.

Diversity, phylogeny and intraspecific variability of Paradiplozoon species (Monogenea: Diplozoidae) parasitizing endemic cyprinoids in the Middle East

Diplozoidae are common monogenean ectoparasites of cyprinoid fish, with the genus Paradiplozoon being the most diversified. Despite recent studies on Diplozoidae from Europe, Africa and Asia, the diversity, distribution and phylogeny of this parasite group appears to be still underestimated in the Middle East. The objective of this study was to investigate the diversity, endemism and host specificity of diplozoids parasitizing cyprinoid fish from the Middle East, considering this region as an important historical interchange of fish fauna, and to elucidate the phylogenetic position of Middle Eastern Paradiplozoon species within Diplozoidae. Four Paradiplozoon species were collected from 48 out of 94 investigated cyprinoid species. Three known species, Paradiplozoon homoion, Paradiplozoon bliccae and Paradiplozoon bingolensis, were recorded on new cyprinoid host species, and a new species, Paradiplozoon koubkovae n. sp., was recorded on Luciobarbus capito and Capoeta capoeta from the Caspian Sea basin in Iran and Turkey. Paradiplozoon bliccae, exhibiting a wide host range in the Middle East, expressed both morphological and genetic intraspecific variabilities. The four Paradiplozoon species collected in the Middle East were placed in divergent clades, showing the rich evolutionary history of diplozoid parasites in the Middle East. Our study also revealed that two lineages of African diplozoids have a Middle Eastern origin. We stress the importance of applying an integrative approach combining morphological, ecological and molecular methods to reveal the real diversity of diplozoids.

Genetic characteristics of classical astroviruses in Shenzhen, China, 2016-2019

Human astrovirus (HAstV) is a single-stranded, positive-sense RNA virus and is the leading cause of viral gastroenteritis. However, despite its prevalence, astroviruses still remain one of the least studied enteroviruses. In this study, we sequenced 11 classical astrovirus strains from clinical samples collected in Shenzhen, China from 2016 to 2019, analyzed their genetic characteristics, and deposited them into GenBank. We conducted phylogenetic analysis using IQ-TREE software, with references to astrovirus sequences worldwide. The phylogeographic analysis was performed using the Bayesian Evolutionary Analysis Sampling Trees program, through Bayesian Markov Chain Monte Carlo sampling. We also conducted recombination analysis with the Recombination Detection Program. The newly sequenced strains were categorized as HAstV genotype 1, which is the predominant genotype in Shenzhen. Phylogeographic reconstruction indicated that HAstV-1 may have migrated from the United States to China, followed by frequent transmission between China and Japan. The recombination analysis revealed recombination events within and across genotypes, and identified a recombination-prone region that produced relatively uniform recombination breakpoints and fragment lengths. The genetic analysis of HAstV strains in Shenzhen addresses the current lack of astrovirus data in the region of Shenzhen and provides key insights to the evolution and transmission of astroviruses worldwide. These findings highlight the importance of improving surveillance of astroviruses.

Parallel anagenetic patterns in endemic Artemisia species from three Macaronesian archipelagos

Anagenetic speciation is an important mode of evolution in oceanic islands, yet relatively understudied compared to adaptive radiation. In the Macaronesian region, three closely related species of Artemisia (i.e. A. argentea, A. thuscula and A. gorgonum) are each endemic from a single archipelago (i.e. Madeira, Canary Islands and Cape Verde, respectively), representing a perfect opportunity to study three similar but independent anagenetic speciation processes. By analysing plastid and nuclear DNA sequences, as well as nuclear DNA amount data, generated from a comprehensive sampling in all the islands and archipelagos where these species are currently distributed, we intend to find common evolutionary patterns that help us explain the limited taxonomic diversification experienced by endemic Macaronesian Artemisia. Our time-calibrated phylogenetic reconstruction suggested that divergence among the three lineages occurred in a coincidental short period of time during the Pleistocene. Haplotype and genetic differentiation analyses showed similar diversity values among A. argentea, A. thuscula and A. gorgonum. Clear phylogeographic patterns-showing comparable genetic structuring among groups of islands-were also found within the three archipelagos. Even from the cytogenetic point of view, the three species presented similarly lower genome size values compared to the mainland closely related species A. arborescens. We hypothesize that the limited speciation experienced by the endemic Artemisia in Madeira, Canary Islands and Cape Verde archipelagos could be related to their recent parallel evolutionary histories as independent lineages, combined with certain shared characteristics of seed dispersal, pollen transport and type of habitat.

Diversification processes in Gerp's mouse lemur demonstrate the importance of rivers and altitude as biogeographic barriers in Madagascar's humid rainforests

Madagascar exhibits exceptionally high levels of biodiversity and endemism. Models to explain the diversification and distribution of species in Madagascar stress the importance of historical variability in climate conditions which may have led to the formation of geographic barriers by changing water and habitat availability. The relative importance of these models for the diversification of the various forest-adapted taxa of Madagascar has yet to be understood. Here, we reconstructed the phylogeographic history of Gerp's mouse lemur (Microcebus gerpi) to identify relevant mechanisms and drivers of diversification in Madagascar's humid rainforests. We used restriction site associated DNA (RAD) markers and applied population genomic and coalescent-based techniques to estimate genetic diversity, population structure, gene flow and divergence times among M. gerpi populations and its two sister species M. jollyae and M. marohita. Genomic results were complemented with ecological niche models to better understand the relative barrier function of rivers and altitude. We show that M. gerpi diversified during the late Pleistocene. The inferred ecological niche, patterns of gene flow and genetic differentiation in M. gerpi suggest that the potential for rivers to act as biogeographic barriers depended on both size and elevation of headwaters. Populations on opposite sides of the largest river in the area with headwaters that extend far into the highlands show particularly high genetic differentiation, whereas rivers with lower elevation headwaters have weaker barrier functions, indicated by higher migration rates and admixture. We conclude that M. gerpi likely diversified through repeated cycles of dispersal punctuated by isolation to refugia as a result of paleoclimatic fluctuations during the Pleistocene. We argue that this diversification scenario serves as a model of diversification for other rainforest taxa that are similarly limited by geographic factors. In addition, we highlight conservation implications for this critically endangered species, which faces extreme habitat loss and fragmentation.

The evolutionary history of hepaciviruses

In the search for natural reservoirs of hepatitis C virus (HCV), a broad diversity of non-human viruses within the Hepacivirus genus has been uncovered. However, the evolutionary dynamics that shaped the diversity and timescale of hepaciviruses evolution remain elusive. To gain further insights into the origins and evolution of this genus, we screened a large dataset of wild mammal samples (n = 1,672) from Africa and Asia, and generated 34 full-length hepacivirus genomes. Phylogenetic analysis of these data together with publicly available genomes emphasizes the importance of rodents as hepacivirus hosts and we identify 13 rodent species and 3 rodent genera (in Cricetidae and Muridae families) as novel hosts of hepaciviruses. Through co-phylogenetic analyses, we demonstrate that hepacivirus diversity has been affected by cross-species transmission events against the backdrop of detectable signal of virus-host co-divergence in the deep evolutionary history. Using a Bayesian phylogenetic multidimensional scaling approach, we explore the extent to which host relatedness and geographic distances have structured present-day hepacivirus diversity. Our results provide evidence for a substantial structuring of mammalian hepacivirus diversity by host as well as geography, with a somewhat more irregular diffusion process in geographic space. Finally, using a mechanistic model that accounts for substitution saturation, we provide the first formal estimates of the timescale of hepacivirus evolution and estimate the origin of the genus to be about 22 million years ago. Our results offer a comprehensive overview of the micro- and macroevolutionary processes that have shaped hepacivirus diversity and enhance our understanding of the long-term evolution of the Hepacivirus genus.

Range extensions of Pacific bone-eating worms (Annelida, Siboglinidae, Osedax)

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.

Evolutionarily diverse origins of deformed wing viruses in western honey bees

Novel transmission routes can allow infectious diseases to spread, often with devastating consequences. Ectoparasitic varroa mites vector a diversity of RNA viruses, having switched hosts from the eastern to western honey bees (Apis cerana to Apis mellifera). They provide an opportunity to explore how novel transmission routes shape disease epidemiology. As the principal driver of the spread of deformed wing viruses (mainly DWV-A and DWV-B), varroa infestation has also driven global honey bee health declines. The more virulent DWV-B strain has been replacing the original DWV-A strain in many regions over the past two decades. Yet, how these viruses originated and spread remains poorly understood. Here, we use a phylogeographic analysis based on whole-genome data to reconstruct the origins and demography of DWV spread. We found that, rather than reemerging in western honey bees after varroa switched hosts, as suggested by previous work, DWV-A most likely originated in East Asia and spread in the mid-20th century. It also showed a massive population size expansion following the varroa host switch. By contrast, DWV-B was most likely acquired more recently from a source outside East Asia and appears absent from the original varroa host. These results highlight the dynamic nature of viral adaptation, whereby a vector's host switch can give rise to competing and increasingly virulent disease pandemics. The evolutionary novelty and rapid global spread of these host-virus interactions, together with observed spillover into other species, illustrate how increasing globalization poses urgent threats to biodiversity and food security.

Genomic Surveillance of Recent Dengue Outbreaks in Colombo, Sri Lanka

All four serotypes of the dengue virus (DENV1-4) cause a phenotypically similar illness, but serial infections from different serotypes increase the risk of severe disease. Thus, genomic surveillance of circulating viruses is important to detect serotype switches that precede community outbreaks of disproportionate magnitude. A phylogenetic analysis was conducted on near full length DENV genomes sequenced from serum collected from a prospective cohort study from the Colombo district, Sri Lanka during a 28-month period using Oxford nanopore technology, and the consensus sequences were analyzed using maximum likelihood and Bayesian evolutionary analysis. From 523 patients, 328 DENV sequences were successfully generated (DENV1: 43, DENV2: 219, DENV3:66). Most circulating sequences originated from a common ancestor that was estimated to have existed from around 2010 for DENV2 and around 2015/2016 for DENV1 and DENV3. Four distinct outbreaks coinciding with monsoon rain seasons were identified during the observation period mostly driven by DENV2 cosmopolitan genotype, except for a large outbreak in 2019 contributed by DENV3 genotype I. This serotype switch did not result in a more clinically severe illness. Phylogeographic analyses showed that all outbreaks started within Colombo city and then spread to the rest of the district. In 2019, DENV3 genotype I, previously, rarely reported in Sri Lanka, is likely to have contributed to a disease outbreak. However, this did not result in more severe disease in those infected, probably due to pre-existing DENV3 immunity in the community. Targeted vector control within Colombo city before anticipated seasonal outbreaks may help to limit the geographic spread of outbreaks.

The Principles of SARS-CoV-2 Intervariant Competition Are Exemplified in the Pre-Omicron Era of the Colombian Epidemic

The first 18 months of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections in Colombia were characterized by three epidemic waves. During the third wave, from March through August 2021, intervariant competition resulted in Mu replacing Alpha and Gamma. We employed Bayesian phylodynamic inference and epidemiological modeling to characterize the variants in the country during this period of competition. Phylogeographic analysis indicated that Mu did not emerge in Colombia but acquired increased fitness there through local transmission and diversification, contributing to its export to North America and Europe. Despite not having the highest transmissibility, Mu's genetic composition and ability to evade preexisting immunity facilitated its domination of the Colombian epidemic landscape. Our results support previous modeling studies demonstrating that both intrinsic factors (transmissibility and genetic diversity) and extrinsic factors (time of introduction and acquired immunity) influence the outcome of intervariant competition. This analysis will help set practical expectations about the inevitable emergences of new variants and their trajectories. IMPORTANCE Before the appearance of the Omicron variant in late 2021, numerous SARS-CoV-2 variants emerged, were established, and declined, often with different outcomes in different geographic areas. In this study, we considered the trajectory of the Mu variant, which only successfully dominated the epidemic landscape of a single country: Colombia. We demonstrate that Mu competed successfully there due to its early and opportune introduction time in late 2020, combined with its ability to evade immunity granted by prior infection or the first generation of vaccines. Mu likely did not effectively spread outside of Colombia because other immune-evading variants, such as Delta, had arrived in those locales and established themselves first. On the other hand, Mu's early spread within Colombia may have prevented the successful establishment of Delta there. Our analysis highlights the geographic heterogeneity of early SARS-CoV-2 variant spread and helps to reframe the expectations for the competition behaviors of future variants.