Molecular characterization of the meq oncogene of Marek's disease virus in vaccinated Brazilian poultry farms reveals selective pressure on prevalent strains

Marek's disease virus (MDV) has become an increasingly virulent pathogen in the poultry industry despite vaccination efforts to control it. Brazil has experienced a significant rise of Marek's disease (MD) outbreaks in recent years. Our study aimed to analyze the complete meq gene sequences to understand the molecular epidemiological basis of MD outbreaks in Brazilian vaccinated layer farms. We detected a high incidence rate of visceral MD (67.74%) and multiple circulating MDV strains. The most prevalent and geographically widespread genotype presented several clinical and molecular characteristics of a highly virulent strain and evolving under positive selective pressure. Phylogenetic and phylogeographic analysis revealed a closer relationship with strains from the USA and Japan. This study sheds light on the circulation of MDV strains capable of infecting vaccinated birds. We emphasize the urgency of adopting preventive measures to manage MDV outbreaks threatening the poultry farming industry.

Spatio-temporal spread of Lassa virus and a new rodent host in the Mano River Union area, West Africa

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

Research Note: Characterization and phylodynamic analysis of new infectious bursal disease virus variants circulating in Argentina

Infectious Bursal Disease is a highly contagious disease that affects young chickens and leads to significant economic losses. Its causal agent is a double-stranded RNA virus that, due to its high error rate during the replication process, gives rise to a constant generation of new virus variants. Until 2014, strains of Infectious Bursal Diseases Virus (IBDV) belonging to genogroup 4 predominated in Argentina, but there have been no reports since then regarding the circulating genogroups in poultry. In this study, 11 recent sequences of Argentine from the hypervariable region of VP2 protein (hvVP2) were analyzed to determine their genogroup, origin, evolution, and amino acid sequence. Samples from chickens showing signs of IBDV infection were collected, and the hvVP2 region was amplified using RT-PCR, followed by sequencing. The results indicated that the analyzed strains belong to genogroup 2, with an estimated evolutionary rate of 1.74 × 10-3 substitutions/site/year. It is speculated that the predominant group of sequences began to spread in Argentina around 2014 and had its origins in China. Another sample is related to strains from South Korea and is not closely linked to the main group. Furthermore, the predicted amino acid sequences show similarity to strains that can evade vaccine-induced immunity. These findings underscore the importance of active surveillance in poultry to mitigate losses caused by IBDV.

Cryptic diversity and speciation in an endemic copepod crustacean Harpacticella inopinata within Lake Baikal

Ancient lakes are hotspots of species diversity, posing challenges and opportunities for exploration of the dynamics of endemic diversification. Lake Baikal in Siberia, the oldest lake in the world, hosts a particularly rich crustacean fauna, including the largest known species flock of harpacticoid copepods with some 70 species. Here, we focused on exploring the diversity and evolution within a single nominal species, Harpacticella inopinata Sars, 1908, using molecular markers (mitochondrial COI, nuclear ITS1 and 28S rRNA) and a set of qualitative and quantitative morphological traits. Five major mitochondrial lineages were recognized, with model-corrected COI distances of 0.20-0.37. A concordant pattern was seen in the nuclear data set, and qualitative morphological traits also distinguish a part of the lineages. All this suggests the presence of several hitherto unrecognized cryptic taxa within the baikalian H. inopinata, with long independent histories. The abundances, distributions and inferred demographic histories were different among taxa. Two taxa, H. inopinata CE and H. inopinata CW, were widespread on the eastern and western coasts, respectively, and were largely allopatric. Patterns in mitochondrial variation, that is, shallow star-like haplotype networks, suggest these taxa have spread through the lake relatively recently. Three other taxa, H. inopinata RE, RW and RW2, instead were rare and had more localized distributions on either coast, but showed deeper intraspecies genealogies, suggesting older regional presence. The rare taxa were often found in sympatry with the others and occasionally introgressed by mtDNA from the common ones. The mitochondrial divergence between and within the H. inopinata lineages is still unexpectedly deep, suggesting an unusually high molecular rate. The recognition of true systematic diversity in the evaluation and management of ecosystems is important in hotspots, as it is everywhere else, while the translation of the diversity into a formal taxonomy remains a challenge.

The Andes as a semi-permeable geographical barrier: Genetic connectivity between structured populations in a widespread spider

Geographical barriers like mountain ranges impede genetic exchange among populations, promoting diversification. The effectiveness of these barriers in limiting gene flow varies between lineages due to each species' dispersal modes and capacities. Our understanding of how the Andes orogeny contributes to species diversification comes from well-studied vertebrates and a few arthropods and plants, neglecting organisms unable to fly or walk long distances. Some arachnids, such as Gasteracantha cancriformis, have been hypothesized to disperse long distances via ballooning (i.e. using their silk to interact with the wind). Yet, we do not know how the environment and geography shape its genetic diversity. Therefore, we tested whether the Andes contributed to the diversification of G. cancriformis acting as an absolute or semi-permeable barrier to genetic connectivity between populations of this spider at opposite sides of the mountain range. We sampled thousands of loci across the distribution of the species and implemented population genetics, phylogenetic, and landscape genetic analyses. We identified two genetically distinct groups structured by the Central Andes, and a third less structured group in the Northern Andes that shares ancestry with the previous two. This structure is largely explained by the altitude along the Andes, which decreases in some regions, possibly facilitating cross-Andean dispersal and gene flow. Our findings support that altitude in the Andes plays a major role in structuring populations in South America, but the strength of this barrier can be overcome by organisms with long-distance dispersal modes together with altitudinal depressions.

Analysis of Genetic Structure and Genetic Diversity in Japanese Grey-Headed Lapwing Population Using mtDNA

The grey-headed lapwing (Vanellus cinereus) is a wading species in East Asia. However, examples of regional population dynamics and genetic research are limited. To reconsider the natural history and current status of the grey-headed lapwing in Japan, we analyzed the genetic diversity of the Japanese grey-headed lapwing population. We collected 77 grey-headed lapwing samples from 12 locations across Japan during the breeding season and three individuals during the wintering season and extracted DNA; 496-bp sequences of the ND2, which form part of the mitochondrial DNA, were determined for genetic analysis of the population. Consequently, 10 haplotypes were detected in 80 individuals, and 67 individuals, 84% of the total, shared two haplotypes, namely Vc1 and Vc2. Furthermore, the results showed that the prevalence of Vc1 was higher mainly in northern Japan, while that of Vc2 was higher mainly in southern Japan. Genetic diversity analysis showed that the overall haplotype diversity in Japan was 0.617, which is not particularly low. The sequence of Vc1 was exactly the same as that of grey-headed lapwing in China. Our study revealed the genetic structure of the grey-headed lapwing, suggesting that as the grey-headed lapwing expanded its distribution area into southern Japan, many Vc2-positive individuals migrated southward, resulting in a higher detection rate of Vc2 in southern Japan.

Deep Learning and Likelihood Approaches for Viral Phylogeography Converge on the Same Answers Whether the Inference Model Is Right or Wrong

Analysis of phylogenetic trees has become an essential tool in epidemiology. Likelihood-based methods fit models to phylogenies to draw inferences about the phylodynamics and history of viral transmission. However, these methods are often computationally expensive, which limits the complexity and realism of phylodynamic models and makes them ill-suited for informing policy decisions in real-time during rapidly developing outbreaks. Likelihood-free methods using deep learning are pushing the boundaries of inference beyond these constraints. In this paper, we extend, compare, and contrast a recently developed deep learning method for likelihood-free inference from trees. We trained multiple deep neural networks using phylogenies from simulated outbreaks that spread among 5 locations and found they achieve close to the same levels of accuracy as Bayesian inference under the true simulation model. We compared robustness to model misspecification of a trained neural network to that of a Bayesian method. We found that both models had comparable performance, converging on similar biases. We also implemented a method of uncertainty quantification called conformalized quantile regression that we demonstrate has similar patterns of sensitivity to model misspecification as Bayesian highest posterior density (HPD) and greatly overlap with HPDs, but have lower precision (more conservative). Finally, we trained and tested a neural network against phylogeographic data from a recent study of the SARS-Cov-2 pandemic in Europe and obtained similar estimates of region-specific epidemiological parameters and the location of the common ancestor in Europe. Along with being as accurate and robust as likelihood-based methods, our trained neural networks are on average over 3 orders of magnitude faster after training. Our results support the notion that neural networks can be trained with simulated data to accurately mimic the good and bad statistical properties of the likelihood functions of generative phylogenetic models.

Phylogenetic and phylogeographic insights into Sri Lankan killifishes (Teleostei: Aplocheilidae)

Three nominal species of the killifish genus Aplocheilus are reported from the lowlands of Sri Lanka. Two of these, Aplocheilus dayi and Aplocheilus werneri, are considered endemic to the island, whereas Aplocheilus parvus is reported from both Sri Lanka and Peninsular India. Here, based on a collection from 28 locations in Sri Lanka, also including a dataset of Asian Aplocheilus downloaded from GenBank, we present a phylogeny constructed from the mitochondrial cytochrome b (cytb), mitochondrial cytochrome c oxidase subunit 1 (cox1), and nuclear recombination activating protein 1 (rag1), and investigate the interrelationships of the species of Aplocheilus in Sri Lanka. The endemic Sri Lankan aplocheilid clade comprising A. dayi and A. werneri is recovered as the sister group to the clade comprising A. parvus from Sri Lanka and Aplocheilus blockii from Peninsular India. The reciprocal monophyly of A. dayi and A. werneri is not supported in our molecular phylogeny. A. dayi and A. werneri display strong sexual dimorphism, but species-level differences are subtle, explained mostly by pigmentation patterns. Their phenotypes exhibit a parapatric distribution and may represent locally adapted forms of a single species. Alternatively, the present study does not rule out the possibility that A. dayi and A. werneri may represent an incipient species pair or that they have undergone introgression or hybridization in their contact zones. We provide evidence that the Nilwala-Gin region of southwestern Sri Lanka may have acted as a drought refugium for these fishes.

Genomic Diversity and Geographic Distribution of Newcastle Disease Virus Genotypes in Africa: Implications for Diagnosis, Vaccination, and Regional Collaboration

The emergence of new virulent genotypes and the continued genetic drift of Newcastle disease virus (NDV) implies that distinct genotypes of NDV are simultaneously evolving in different geographic locations across the globe, including throughout Africa, where NDV is an important veterinary pathogen. Expanding the genomic diversity of NDV increases the possibility of diagnostic and vaccine failures. In this review, we systematically analyzed the genetic diversity of NDV genotypes in Africa using the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines. Information published between 1999 and 2022 were used to obtain the genetic background of different genotypes of NDV and their geographic distributions in Africa. The following genotypes were reported in Africa: I, II, III, IV, V, VI, VII, VIII, XI, XIII, XIV, XVII, XVIII, XX, and XXI. A new putative genotype has been detected in the Democratic Republic of the Congo. However, of 54 African countries, only 26 countries regularly report information on NDV outbreaks, suggesting that this number may be vastly underestimated. With eight different genotypes, Nigeria is the country with the greatest genotypic diversity of NDV among African countries. Genotype VII is the most prevalent group of NDV in Africa, which was reported in 15 countries. A phylogeographic analysis of NDV sequences revealed transboundary transmission of the virus in Eastern Africa, Western and Central Africa, and in Southern Africa. A regional and continental collaboration is recommended for improved NDV risk management in Africa.

Phylogeographic Substructuring in the Southernmost Refugium of the European Common Frog Rana temporaria

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

Environmental DNA unveils deep phylogeographic structure of a freshwater fish

Phylogeography bears an important part in ecology and evolution. However, current phylogeographic studies are largely constrained by limited numbers of individual samples. Using an environmental DNA (eDNA) assay for phylogeographic analyses, this study provides detailed information regarding the history of Siberian stone loach Barbatula toni, a primary freshwater fish across the whole range of Hokkaido, Japan. Based on an eDNA metabarcoding on 293 river water samples, we detected eDNA from B. toni in 189 rivers. A total of 51 samples, representing the entire island, were then selected from the B. toni eDNA-positive sample set for the subsequent analyses. To elucidate the phylogeographic structure of B. toni, newly developed eDNA metabarcoding primers (Barba-cytb-F/R) were applied to these samples, specifically targeting their haplotypic variation in cytochrome b. After a bioinformatic processing to mitigate haplotypic false positives, a total of 50 eDNA haplotypes were identified. Two regionally restricted, genetically distinct lineages of the species were revealed as a result of phylogeographic analyses on the haplotypes and tissue-derived DNA from B. toni. According to a molecular clock analysis, they have been genetically isolated for at least 1.5 million years, suggesting their ancient origin and colonisation of Hokkaido, presumably in the glacial periods. These results demonstrate how freshwater fishes can alter their distributions over evolutionary timescales and how eDNA assay can deepen our understanding of phylogeography.

Survival in nunatak and peripheral glacial refugia of three alpine plant species is partly predicted by altitudinal segregation

Mountain biota survived the Quaternary cold stages most probably in peripheral refugia and/or ice-free peaks within ice-sheets (nunataks). While survival in peripheral refugia has been broadly demonstrated, evidence for nunatak refugia is still scarce. We generated RADseq data from three mountain plant species occurring at different elevations in the southeastern European Alps to investigate the role of different glacial refugia during the Last Glacial Maximum (LGM). We tested the following hypotheses. (i) The deep Piave Valley forms the deepest genetic split in the species distributed across it, delimiting two peripheral refugia. (ii) The montane to alpine species Campanula morettiana and Primula tyrolensis survived the LGM in peripheral refugia, while high-alpine to subnival Saxifraga facchinii likely survived in several nunatak refugia. (iii) The lower elevation species suffered a strong population decline during the LGM. By contrast, the higher elevation species shows long-term stability of population sizes due to survival on permanently ice-free peaks and small population sizes at present. We found peripheral refugia on both sides of the Piave Valley, which acted as a major genetic barrier. Demographic modelling confirmed nunatak survival not only for S. facchinii but also for montane to alpine C. morettiana. Altitudinal segregation influenced the species' demographic fluctuations, with the lower elevation species showing a significant population increase at the end of the LGM, and the higher elevation species either showing decrease towards the present or stable population sizes with a short bottleneck. Our results highlight the role of nunatak survival and species ecology in the demographic history of mountain species.

mapmixture: An R package and web app for spatial visualisation of admixture and population structure

The mapmixture R package and interactive web app are tools to aid visualisation of admixture and population structure in geographic space. The purpose of mapmixture is to enable and encourage molecular ecologists, and in particular population geneticists and phylogeneticists, to plot their admixture, ancestry or assignment results on a map when location information is available. mapmixture accepts data in the format typically generated by admixture analyses and visualises proportions to each genetic cluster per site as pie charts on a projected (optional) map. Combining this site-based map presentation approach with the routine individual-based presentation of admixture (structure) barplots will enhance interpretation of genetic-geographic patterns. Additionally, in the context of science communication, this enables clearer transfer of spatial genetic information to readers or listeners, and especially to audiences that do not have a background in genetics but who are able to use the genetic information as evidence in conservation management. The latest version of mapmixture is available on GitHub (https://github.com/tom-jenkins/mapmixture), which details installation instructions and examples of how to use the package and interactive web app.

Investigation of CaMV-host co-evolution through synonymous codon pattern

Cauliflower mosaic virus (CaMV) has a double-stranded DNA genome and is globally distributed. The phylogeny tree of 121 CaMV isolates was categorized into two primary groups, with Iranian isolates showing the greatest genetic variations. Nucleotide A demonstrated the highest percentage (36.95%) in the CaMV genome and the dinucleotide odds ratio analysis revealed that TC dinucleotide (1.34 ≥ 1.23) and CG dinucleotide (0.63 ≤ 0.78) are overrepresented and underrepresented, respectively. Relative synonymous codon usage (RSCU) analysis confirmed codon usage bias in CaMV and its hosts. Brassica oleracea and Brassica rapa, among the susceptible hosts of CaMV, showed a codon adaptation index (CAI) value above 0.8. Additionally, relative codon deoptimization index (RCDI) results exhibited the highest degree of deoptimization in Raphanus sativus. These findings suggest that the genes of CaMV underwent codon adaptation with its hosts. Among the CaMV open reading frames (ORFs), genes that produce reverse transcriptase and virus coat proteins showed the highest CAI value of 0.83. These genes are crucial for the creation of new virion particles. The results confirm that CaMV co-evolved with its host to ensure the optimal expression of its genes in the hosts, allowing for easy infection and effective spread. To detect the force behind codon usage bias, an effective number of codons (ENC)-plot and neutrality plot were conducted. The results indicated that natural selection is the primary factor influencing CaMV codon usage bias.

Intraspecific divergence in a coastal plant, Euphorbia jolkinii, at a major biogeographic boundary in East Asia

PREMISE

Quaternary climatic fluctuations and long-distance seed dispersal across the sea are critical factors affecting the distribution of coastal plants, but the spatiotemporal nature of population expansion and distribution change of East Asian coastal plants during this period are rarely examined. To explore this process, we investigated the genome-wide phylogenetic patterns of Euphorbia jolkinii Boiss. (Euphorbiaceae), which grows widely on littoral areas of Japan, Korea, and Taiwan.

METHODS

We used plastome sequences and genome-wide single nucleotide polymorphisms in samples across the species range to reveal phylogeographic patterns and spatiotemporal distributional changes. We conducted ecological niche modeling for the present and the last glacial maximum (LGM).

RESULTS

Genetic differentiation was observed between the northern and southern populations of E. jolkinii, separated by the major biogeographic boundary, the Tokara Gap. These two groups of populations differentiated during the glacial period and subsequently intermingled in the intermorainic areas of the central Ryukyu Islands after the LGM. Ecological niche models suggested that the potential range of E. jolkinii was restricted to southern Kyushu; however, it was widespread in the southern Ryukyu Islands and Taiwan during the LGM.

CONCLUSIONS

This study provides evidence of genetic differentiation among coastal plant populations separated by the prominent biogeographical boundary. Although coastal plants are typically expected to maintain population connectivity through sea-drifted seed dispersal, our findings suggest that genetic differences may arise because of a combination of limited gene flow and changes in climate during the glacial period.

Phylogeographic Pattern of the Assassin Bug Sycanus bifidus Inferred from Mitochondrial Genomes and Nuclear Genes

The assassin bug Sycanus bifidus has a wide distribution across southern China. This study explored its distribution and evolution by analyzing mitochondrial and nuclear ribosomal RNA genes, revealing how Pleistocene climate and geological changes shaped its phylogeography. We identified two main clades, A and B, that diverged in the Middle Pleistocene. Hainan Island's populations form a unique group within Clade A, suggesting that the Qiongzhou Strait served as a dispersal corridor during glaciation. Rising sea levels likely separated the Hainan population afterward. Ecological niche modeling showed that both populations have been viable since the last interglacial period, with demographic analyses indicating possible expansions during the Middle and Late Pleistocene, driven by favorable climates. This study highlights the significant effects of Pleistocene sea-level and climatic changes on the distribution and evolution of S. bifidus in China.

Phylodynamic analysis of foot-and-mouth disease virus evolution in Mar Chiquita, Argentina

Foot-and-mouth disease is a highly contagious disease affecting cloven-hoofed animals, resulting in considerable economic losses. Its causal agent is foot-and-mouth disease virus (FMDV), a picornavirus. Due to its error-prone replication and rapid evolution, the transmission and evolutionary dynamics of FMDV can be studied using genomic epidemiological approaches. To analyze FMDV evolution and identify possible transmission routes in an Argentinean region, field samples that tested positive for FMDV by PCR were obtained from 21 farms located in the Mar Chiquita district. Whole FMDV genome sequences were obtained by PCR amplification in seven fragments and sequencing using the Sanger technique. The genome sequences obtained from these samples were then analyzed using phylogenetic, phylogeographic, and evolutionary approaches. Three local transmission clusters were detected among the sampled viruses. The dataset was analyzed using Bayesian phylodynamic methods with appropriate coalescent and relaxed molecular clock models. The estimated mean viral evolutionary rate was 1.17 × 10- 2 substitutions/site/year. No significant differences in the rate of viral evolution were observed between farms with vaccinated animals and those with unvaccinated animals. The most recent common ancestor of the sampled sequences was dated to approximately one month before the first reported case in the outbreak. Virus transmission started in the south of the district and later dispersed to the west, and finally arrived in the east. Different transmission routes among the studied herds, such as non-replicating vectors and close contact contagion (i.e., aerosols), may be responsible for viral spread.

Divergence of the freshwater sleeper, Neodontobutishainanensis (Chen, 1985) (Teleostei, Odontobutidae), in the Pearl River basin and on Hainan Island of southern China

Study of divergence of freshwater fish populations between island and adjacent mainland areas can shed light on the phylogeographical relationships of these regions. Neodontobutishainanensis is a freshwater fish species restricted to Hainan Island and in Guangdong and Guangxi provinces in the southern mainland China. We examine the phylogenetic relationship and population structure of N.hainanensis based on 3,176 nuclear loci using a gene-capture method. STRUCTURE analysis and principal coordinate analyses (PCA) indicate that populations from Guangdong, Guangxi and Hainan are each distinct, except that some individuals of the Guangdong population share minor genetic components with individuals of the Guangxi population. In the concatenated gene tree, the Hainan population is grouped with the Guangdong population, but the coalescent tree groups the Hainan population as the sister to the Guangxi population. Finally, coalescent simulations confirmed the divergence pattern supported by the coalescent tree and revealed a one-way introgression from the Guangxi population to the Guangdong population, which can explain the discordant results supported by the concatenated and coalescent phylogenetic analyses. Due to recent decline of N.hainanensis populations and the genetic patterns in this species, as revealed in this study, the populations in the three areas should be treated as separate conservation units.

The impact of climate change and the conservation of the keystone Asian honeybee using niche models and systematic prioritization

Global warming has seriously disturbed the Earth's ecosystems, and in this context, Asian honeybee (Apis cerana) has experienced a dramatic decline in recent decades. Here, we examined both direct and indirect effects of climate change on A. cerana through ecological niche modeling of A. cerana, and its disease pathogens (i.e., Chinese sacbrood virus and Melissococcus plutonius) and enemies (i.e., Galleria mellonella and Vespa mandarinia). Ecological niche modeling predicts that climate change will increase the potential suitability of A. cerana, but it will also cause some of the original habitat areas to become unsuitable. Outbreak risks of Chinese sacbrood disease and European Foulbrood will increase dramatically, while those of G. mellonella and V. mandarinia will decrease only slightly. Thus, climate change will produce an unfavorable situation for even maintaining some A. cerana populations in China in the future. Genetic structure analyses showed that the A. cerana population from Hainan Island had significant genetic differentiation from that of the mainland, and there was almost no gene flow between the 2, suggesting that urgent measures are needed to protect the unique genetic resources there. Through taking an integrated planning technique with the Marxan approach, we optimized conservation planning, and identified potential nature reserves (mainly in western Sichuan and southern Tibet) for conservation of A. cerana populations. Our results can provide insights into the potential impact of climate change on A. cerana, and will help to promote the conservation of the keystone honeybee in China and the long-term sustainability of its ecosystem services.

Emergence and spread of Mycobacterium ulcerans at different geographic scales

The classical lineage of Mycobacterium ulcerans is the most prevalent clonal group associated with Buruli ulcer in humans. Its reservoir is strongly associated with the environment. We analyzed together 1,045 isolates collected from 13 countries on two continents to define the evolutionary history and population dynamics of this lineage. We confirm that this lineage spread over 7,000 years from Australia to Africa with the emergence of outbreaks in distinct waves in the 18th and 19th centuries. In sharp contrast with its global spread over the last century, transmission chains are now mostly local, with little or no dissemination between endemic areas. This study provides new insights into the phylogeography and population dynamics of M. ulcerans, highlighting the importance of comparative genomic analyses to improve our understanding of pathogen transmission.

IMPORTANCE

Mycobacterium ulcerans is an environmental mycobacterial pathogen that can cause Buruli ulcer, a severe cutaneous infection, mostly spread in Africa and Australia. We conducted a large genomic study of M. ulcerans, combining genomic and evolutionary approaches to decipher its evolutionary history and pattern of spread at different geographic scales. At the scale of villages in an endemic area of Benin, the circulating genotypes have been introduced in recent decades and are not randomly distributed along the river. On a global scale, M. ulcerans has been spreading for much longer, resulting in distinct and compartmentalized endemic foci across Africa and Australia.

Emergence and dissemination of equine-like G3P[8] rotavirus A in Brazil between 2015 and 2021

Rotavirus A (RVA) is a major cause of acute gastroenteritis globally that is classically genotyped by its two immunodominant outer capsid proteins, VP7 (G-) and VP4 (P-). Recent evidence suggests that the reassortant equine-like G3P[8] strain played a substantial role in RVA transmission in Brazil since 2015. To understand its global emergence and dissemination in Brazilian territory, stool samples collected from 11 Brazilian states (n = 919) were genotyped by RT-qPCR and proceeded to sequence the VP7 gene (n = 102, 79 being newly generated) of the G3P[8] samples with pronounced viral loads. Our phylogenetic genotyping showed that G3P[8] became the dominant strain in Brazil between 2017 and 2020, with equine-like variants representing 75%-100% of VP7 samples in this period. A Bayesian discrete phylogeographic analysis strongly suggests that the equine-like G3P[8] strain originated in Asia during the early 2010s and subsequently spread to Europe, the Caribbean, and South America. Multiple introductions were detected in Brazil between 2014 and 2017, resulting in five national clusters. The reconstruction of the effective population size of the largest Brazilian cluster showed an expansion until 2017, followed by a plateau phase until 2019 and subsequent contraction. Our study also supports that most mutations fixed during equine-like G3P[8] evolution were synonymous, suggesting that adaptive evolution was not an important driving force during viral dissemination in humans, potentially increasing its susceptibility to acquired immunity. This research emphasizes the need for comprehensive rotavirus genomic surveillance that allows close monitoring of its ever-shifting composition and informs more effective public health policies.IMPORTANCEOur original article demonstrated the origin and spread in a short time of equine-like G3P[8] in Brazil and the world. Due to its segmented genome, it allows numerous mechanisms including genetic drift and reassortment contribute substantially to the genetic diversity of rotavirus. Although the effectiveness and increasing implementation of vaccination have not been questioned, a matter of concern is its impact on the emergence of escape mutants or even the spread of unusual strains of zoonotic transmission that could drive epidemic patterns worldwide. This research emphasizes the need for comprehensive rotavirus genomic surveillance, which could facilitate the formulation of public policies aimed at preventing and mitigating its transmission.

Speciation across biomes: Rapid diversification with reproductive isolation in the Australian delicate mice

Phylogeographic studies of continental clades, especially when combined with palaeoclimate modelling, provide powerful insight into how environment drives speciation across climatic contexts. Australia, a continent characterized by disparate modern biomes and dynamic climate change, provides diverse opportunity to reconstruct the impact of past and present environments on diversification. Here, we use genomic-scale data (1310 exons and whole mitogenomes from 111 samples) to investigate Pleistocene diversification, cryptic diversity, and secondary contact in the Australian delicate mice (Hydromyini: Pseudomys), a recent radiation spanning almost all Australian environments. Across northern Australia, we find no evidence for introgression between cryptic lineages within Pseudomys delicatulus sensu lato, with palaeoclimate models supporting contraction and expansion of suitable habitat since the last glacial maximum. Despite multiple contact zones, we also find little evidence of introgression at a continental scale, with the exception of a potential hybrid zone in the mesic biome. In the arid zone, combined insights from genetic data and palaeomodels support a recent expansion in the arid specialist P. hermannsburgensis and contraction in the semi-arid P. bolami. In the face of repeated secondary contact, differences in sperm morphology and chromosomal rearrangements are potential mechanisms that maintain species boundaries in these recently diverged species. Additionally, we describe the western delicate mouse as a new species and recommend taxonomic reinstatement of the eastern delicate mouse. Overall, we show that speciation in an evolutionarily young and widespread clade has been driven by environmental change, and potentially maintained by divergence in reproductive morphology and chromosome rearrangements.

Phylogenetics and phylogeography of Euphorbia canariensis reveal an extreme Canarian-Asian disjunction but limited inter-island colonization

Euphorbia canariensis is an iconic endemic species representative of the lowland xerophytic communities of the Canary Islands. It is widely distributed in the archipelago despite having diasporas unspecialized for long-distance dispersal. Here, we reconstructed the evolutionary history of E. canariensis at two levels: a time-calibrated phylogenetic analysis aimed at clarifying interspecific relationships and large-scale biogeographic patterns; and a phylogeographic study focused on the history of colonization across the Canary Islands. For the phylogenetic study, we sequenced the ITS region for E. canariensis and related species of Euphorbia sect. Euphorbia. For the phylogeographic study, we sequenced two cpDNA regions for 28 populations representing the distribution range of E. canariensis. The number of inter-island colonization events was explored using PAICE, a recently developed method that includes a sample size correction. Additionally, we used species distribution modelling (SDM) to evaluate environmental suitability for E. canariensis through time. Phylogenetic results supported a close relationship between E. canariensis and certain Southeast Asian species (E. epiphylloides, E. lacei, E. sessiliflora). In the Canaries, E. canariensis displayed a west-to-east colonization pattern, not conforming to the "progression rule", i.e. the concordance between phylogeographic patterns and island emergence times. We estimated between 20 and 50 inter-island colonization events, all of them in the Quaternary, and SDM suggested a late Quaternary increase in environmental suitability for E. canariensis. The extreme biogeographic disjunction between Macaronesia and Southeast Asia (ca. 11,000 km) parallels that found in a few other genera (Pinus, Dracaena). We hypothesize that these disjunctions are better explained by extinction across north Africa and southwest Asia rather than long-distance dispersal. The relatively low number of inter-island colonization events across the Canaries is congruent with the low dispersal capabilities of E. canariensis.

Population genomics of flat-tailed horned lizards (Phrynosoma mcallii) informs conservation and management across a fragmented Colorado Desert landscape

Phrynosoma mcallii (flat-tailed horned lizards) is a species of conservation concern in the Colorado Desert of the United States and Mexico. We analysed ddRADseq data from 45 lizards to estimate population structure, infer phylogeny, identify migration barriers, map genetic diversity hotspots, and model demography. We identified the Colorado River as the main geographic feature contributing to population structure, with the populations west of this barrier further subdivided by the Salton Sea. Phylogenetic analysis confirms that northwestern populations are nested within southeastern populations. The best-fit demographic model indicates Pleistocene divergence across the Colorado River, with significant bidirectional gene flow, and a severe Holocene population bottleneck. These patterns suggest that management strategies should focus on maintaining genetic diversity on both sides of the Colorado River and the Salton Sea. We recommend additional lands in the United States and Mexico that should be considered for similar conservation goals as those in the Rangewide Management Strategy. We also recommend periodic rangewide genomic sampling to monitor ongoing attrition of diversity, hybridization, and changing structure due to habitat fragmentation, climate change, and other long-term impacts.

Phylogeography and population structure of Lagocephalus spadiceus (Richardson, 1845) (Tetraodontiformes, Tetraodontidae) in the South China Sea

The climate fluctuations during the Late Pleistocene significantly influenced the phylogeographic structure and historical dynamics of marine fishes in the marginal seas of the western Pacific Ocean. The puffer fish, Lagocephalus spadiceus, holds substantial nutritional and economic value in the South China Sea. To investigate the demographic history and population structure of the L. spadiceus, the mitochondrial DNA COI and Cyt b gene datasets from 300 individuals across eight populations in the South China Sea were sequenced. Our findings revealed high haplotype diversity (0.874 ± 0.013) and low nucleotide diversity (0.00075 ± 0.00058). The phylogenetic tree and haplotype networks revealed no significant genetic differentiation along the northern coast of South China Sea. Neutrality tests, mismatch distribution analyses, and Bayesian skyline plots suggested that L. spadiceus underwent population expansion during the Late Pleistocene. Both ocean currents and climate change significantly influenced the geographical distribution and genetic population structure of L. spadiceus.