Bayesian phylogenetics with BEAUti and the BEAST 1.7

Evolution of H6N6 viruses in China between 2014 and 2019 involves multiple reassortment events

H6N6 avian influenza viruses (AIVs) have been widely detected in wild birds, poultry, and even mammals. Recently, H6N6 viruses were reported to be involved in the generation of H5 and H7 subtype viruses. To investigate the emergence, evolutionary pattern, and potential for an epidemic of H6N6 viruses, the complete genomes of 198 H6N6 viruses were analyzed, including 168 H6N6 viruses deposited in the NCBI and GISAID databases from inception to January 2019 and 30 isolates collected from China between November 2014 and January 2019. Using phylogenetic analysis, the 198 strains of H6N6 viruses were identified as 98 genotypes. Molecular clock analysis indicated that the evolution of H6N6 viruses in China was constant and not interrupted by selective pressure. Notably, the laboratory isolates reassorted with six subtype viruses: H6N2, H5N6, H7N9, H5N2, H4N2, and H6N8, resulting in nine novel H6N6 reassortment events. These results suggested that H6N6 viruses can act as an intermediary in the evolution of H5N6, H6N6, and H7N9 viruses. Animal experiments demonstrated that the 10 representative H6N6 viruses showed low pathogenicity in chickens and were capable of infecting mice without prior adaptation. Our findings suggest that H6N6 viruses play an important role in the evolution of AIVs, and it is necessary to continuously monitor and evaluate the potential epidemic of the H6N6 subtype viruses.

Emergence and circulation of enterovirus B species in infants in southern China: A multicenter retrospective analysis

BACKGROUND

Enteroviruses (EV) are common and can cause severe diseases, particularly in young children. However, the information of EV infection in infants in China is limited due to the vast population size and extensive geographical area of the country. Here, we conducted a retrospective multicenter analysis of available EV data to assess the current epidemiological situation in the infant population in southern China.

METHODS

The study enrolled infants with suspected EV infection from 34 hospitals across 12 cities in southern China between 2019 to 2022, and the confirmation of EV was done using RT-PCR and VP1 gene sequencing.

RESULTS

Out of 1221 infants enrolled, 330 (27.03%) were confirmed as EV-infected. Of these, 260 (78.79%) were newborns aged 0-28 days. The EV belonged to three species: EV-B (80.61%), EV-A (11.82%), and human rhinovirus (7.58%). Newborns were more susceptible to EV-B than older infants (p < 0.001). Within EV-B, we identified 15 types, with coxsackievirus (CV) B3 (20.91%), echovirus (E) 11 (19.70%), and E18 (16.97%) being the most common. The predominant EV types changed across different years. EV infection in infants followed a seasonal pattern, with a higher incidence from May to August. Furthermore, perinatal mother-to-child EV transmission in 12 mother-newborn pairs were observed.

CONCLUSION

Our study is the first to demonstrate the emergence and widespread circulation of EV-B species, mainly CVB3, E11, and E18, in southern China, primarily affecting young infants. This research provides valuable insights for future epidemic assessment, prediction, as well as the elimination of mother-to-child transmission.

Climate, demography, immunology, and virology combine to drive two decades of dengue virus dynamics in Cambodia

The incidence of dengue virus disease has increased globally across the past half-century, with highest number of cases ever reported in 2019 and again in 2023. We analyzed climatological, epidemiological, and phylogenomic data to investigate drivers of two decades of dengue in Cambodia, an understudied endemic setting. Using epidemiological models fit to a 19-y dataset, we first demonstrate that climate-driven transmission alone is insufficient to explain three epidemics across the time series. We then use wavelet decomposition to highlight enhanced annual and multiannual synchronicity in dengue cycles between provinces in epidemic years, suggesting a role for climate in homogenizing dynamics across space and time. Assuming reported cases correspond to symptomatic secondary infections, we next use an age-structured catalytic model to estimate a declining force of infection for dengue through time, which elevates the mean age of reported cases in Cambodia. Reported cases in >70-y-old individuals in the 2019 epidemic are best explained when also allowing for waning multitypic immunity and repeat symptomatic infections in older patients. We support this work with phylogenetic analysis of 192 dengue virus (DENV) genomes that we sequenced between 2019 and 2022, which document emergence of DENV-2 Cosmopolitan Genotype-II into Cambodia. This lineage demonstrates phylogenetic homogeneity across wide geographic areas, consistent with invasion behavior and in contrast to high phylogenetic diversity exhibited by endemic DENV-1. Finally, we simulate an age-structured, mechanistic model of dengue dynamics to demonstrate how expansion of an antigenically distinct lineage that evades preexisting multitypic immunity effectively reproduces the older-age infections witnessed in our data.

Speciation in savanna birds in South America: The case of the Least Nighthawk Chordeiles pusillus (Aves: Caprimulgidae) in and out of the Amazon

The Least Nighthawk Chordeiles pusillus is widespread wherever there are savannas in the South American tropics, often in isolated patches, such as white-sands savannas in the Amazon rainforest realm. Here, we investigate genetic relationships between populations of the Least Nighthawk to understand historical processes leading to its diversification and to determine dispersal routes between northern and southern savannas by way of three hypothesized dispersal corridors by comparing samples from white-sand savannas to samples from other savannas outside of the Amazon rainforest region. We use 32 mtDNA samples from the range of C. pusillus to infer a dated phylogeny. In a subset of 17 samples, we use shotgun sequences to infer a distance-based phylogeny and to estimate individual admixture proportions. We calculate gene flow and shared alleles between white-sand and non-Amazonian populations using the ABBA-BABA test (D statistics), and Principal Component Analysis (PCA) to examine genetic structure within and between lineages. Finally, we use species distribution modelling (SDM) of conditions during the Last Glacial Maximum (LGM), currently, and in the future (2050-2080) to predict potential species occurrence under a climate change scenario. Two main clades (estimated to have diverged around 1.07 million years ago) were recovered with mtDNA sequences and Single Nucleotide Polymorphism (SNPs) and were supported by NGSadmix and PCA: one in the Amazon basin white-sand savannas, the other in the non-Amazonian savannas. Possible allele sharing between these clades was indicated by the D-statistics between northern non-Amazonian populations and the white-sand savanna population, but this was not corroborated by the admixture analyses. Dispersal among northern non-Amazonian populations may have occurred in a dry corridor between the Guianan and the Brazilian Shield, which has since moved eastward. Our data suggest that the lineages separated well before the Last Glacial Maximum, consequently dispersal could have happened at any earlier time during similar climatic conditions. Subsequently, non-Amazonian lineages became more divergent among themselves, possibly connecting and dispersing across the mouth of the Amazon River across Marajó island during favourable climatic conditions in the Pleistocene.

Comparative phylogeography reveals dissimilar genetic differentiation patterns in two sympatric amphibian species

Global climate change is expected to have a profound effect on species distribution. Due to the temperature constraints, some narrow niche species could shift their narrow range to higher altitudes or latitudes. In this study, we explored the correlation between species traits, genetic structure, and geographical range size. More specifically, we analyzed how these variables are affected by differences in fundamental niche breadth or dispersal ability in the members of two sympatrically distributed stream-dwelling amphibian species (frog, Quasipaa yei; salamander, Pachyhynobius shangchengensis), in Dabie Mountains, East China. Both species showed relatively high genetic diversity in most geographical populations and similar genetic diversity patterns (JTX, low; BYM, high) correlation with habitat changes and population demography. Multiple clustering analyses were used to disclose differentiation among the geographical populations of these two amphibian species. Q. yei disclosed the relatively shallow genetic differentiation, while P. shangchengensis showed an opposite pattern. Under different historical climatic conditions, all ecological niche modeling disclosed a larger suitable habitat area for Q. yei than for P. shangchengensis; these results indicated a wider environment tolerance or wider niche width of Q. yei than P. shangchengensis. Our findings suggest that the synergistic effects of environmental niche variation and dispersal ability may help shape genetic structure across geographical topology, particularly for species with extremely narrow distribution.

An expanded phylogeny of the genus Pseudamnicola (Gastropoda; Truncatelloidea; Hydrobiidae) across the Mediterranean Basin

The genus Pseudamnicola Paulucci, 1878, is commonly found throughout the Mediterranean region. The genus displays considerable levels of endemism, accompanied by notable systematic and taxonomic ambiguity. However, the application of molecular data has proven highly effective in clarifying taxonomy and unveiling the diversity of cryptic species within the genus. Therefore, we employed all cytochrome c oxidase subunit I sequence data available and generated new ones from Greece to infer the phylogeny of the genus throughout its Mediterranean range and estimate the divergence times as well as the ancestral area of diversification. Our phylogenetic and time-estimate analyses demonstrate that with 36 to 38 extant Pseudamnicola species and genetic divergences across species ranging from 0.5% to 11.9% on average, the genus underwent relatively recent diversification during late Miocene (6.53 Ma), and the primary speciation events occurred during Plio-Pleistocene. The Italian Peninsula and Islands and the Ionian Drainages as defined by the Freshwater Ecoregions of the World are the ancestral regions of the genus following two different dispersal routes. Our study contributes to deepening our understanding of Pseudamnicola phylogeny by using data from throughout its range for the first time. This phylogeny provides evidence and confirms previous studies that relatively recent habitat isolation, followed by founder and dispersal events, has been one of the primary reasons for the evolution of the genus Pseudamnicola in the Mediterranean basin.

Phylogenetic evidence clarifies the history of the extrusion of Indochina

The southeastward extrusion of Indochina along the Ailao Shan-Red River shear zone (ARSZ) is one of two of the most prominent consequences of the India-Asia collision. This plate-scale extrusion has greatly changed Southeast Asian topography and drainage patterns and effected regional climate and biotic evolution. However, little is known about how Indochina was extruded toward the southeast over time. Here, we sampled 42 plant and animal clades (together encompassing 1,721 species) that are distributed across the ARSZ and are not expected to disperse across long distances. We first assess the possible role of climate on driving the phylogenetic separations observed across the ARSZ. We then investigate the temporal dynamics of the extrusion of Indochina through a multitaxon analysis. We show that the lineage divergences across the ARSZ were most likely associated with the Indochinese extrusion rather than climatic events. The lineage divergences began at ~53 Ma and increased sharply ~35 Ma, with two peaks at ~19 Ma and ~7 Ma, and one valley at ~13 Ma. Our results suggest a two-phase model for the extrusion of Indochina, and in each phase, the extrusion was subject to periods of acceleration and decrease, in agreement with the changes of the India-Asia convergence rate and angle from the early Eocene to the late Miocene. This study highlights that a multitaxon analysis can illuminate the timing of subtle historical events that may be difficult for geological data to pinpoint and can be used to explore other tectonic events.

Comparative chloroplast genome analyses of cultivated and wild Capsicum species shed light on evolution and phylogeny

The chloroplast (cp.) genome, also known as plastome, plays crucial roles in plant survival, adaptation, and evolution. The stable genetic structure of cp. genomes provides an ideal system for investigating species evolution. We sequenced three complete cp. genome sequences of Capsicum species and analyzed them using sequences of various Capsicum species retrieved from the NCBI database. The cp. genome of Capsicum species maintains a well-preserved quadripartite structure consisting of two inverted repeats (IRs) flanked by a large single copy (LSC) region and a small single copy (SSC) region. The sizes of cp. genome sequences ranged from 156,583 bp (C. lycianthoides) to 157,390 bp (C.pubescens). A total of 127-132 unique genes, including 83-87 protein-coding, 36-37 tRNA, and eight rRNA genes, were predicted. Comparison of cp. genomes of 10 Capsicum species revealed high sequence similarity in genome-wide organization and gene arrangements. Fragments of trnT-UGU/trnL-UAA, ccsA, ndhD, rps12, and ycf1 were identified as variable regions, and nucleotide variability of LSC and SSC was higher than that of IR. Phylogenetic speciation analysis showed that the major domesticated C. annuum species were the most extensively divergent species and closely related to C. tovarii and C. frutescens. Analysis of divergent times suggested that a substantial range of speciation events started occurring ~ 25.79 million years ago (Mya). Overall, comparative analysis of cp. genomes of Capsicum species not only offers new insights into their genetic variation and phylogenetic relationships, but also lays a foundation for evolutionary history, genetic diversity, conservation, and biological breeding of Capsicum species.

Shifts in reproductive strategies in the evolutionary trajectory of plant lineages

Understanding the maintenance and shift in reproductive strategies is a fundamental question in evolutionary research. Although many efforts have been made to compare different reproductive strategies, the association between reproductive strategies and lineage divergence is largely unknown. To explore the impact of different reproductive strategies on lineage divergence, we investigated the evolution of clonality in Saxifraga sect. Irregulares+Heterisia. By integrating several lines of evidence, we found that the loss of clonality in Irregulares+Heterisia was associated with a progressive increase in diversification rate and intraspecific morphological diversity but with a reduction in species distribution range. Our findings provide insights into the ecological and evolutionary effects of different reproductive strategies, suggesting the necessity of integrating clonality into ecological and evolutional research.

Phylogenetic assessment within a species complex of a subterranean rodent (Geomys bursarius) with conservation implications for isolated subspecies

Range contraction and expansion from glaciation have led to genetic divergence that may be particularly pronounced in fossorial species with low dispersal. The plains pocket gopher (Geomys bursarius) is a fossorial species that ranges widely across North America but has a poorly understood phylogeny. We used mitogenomes (14,996 base pairs) from 56 individuals across seven subspecies, plus two outgroup species, to assess genetic divergence from minimum spanning trees, measure genetic distances, and infer phylogenetic trees using BEAST. We found G. b. wisconsinensis was monophyletic with recent divergence. Further assessment is needed for G. b. major because it was paraphyletic and exhibited inconsistent groupings with other clades. Importantly, we identified G. b. illinoensis as being genetically distinct and monophyletic likely due to a unique colonization event eastward across the Mississippi River. Because G. b. illinoensis faces continued pressures from niche reduction and habitat loss, we recommend that G. b. illinoensis be considered an evolutionary significant unit warranting conservation actions to promote connectivity and restore suitable habitat. Such conservation efforts should benefit other grassland species including those originating from clades west of the Mississippi River that may also be evolutionary significant units.

Recent speciation and adaptation to aridity in the ecologically diverse Pilbara region of Australia enabled the native tobaccos (Nicotiana; Solanaceae) to colonize all Australian deserts

Over the last 6 million years, the arid Australian Eremaean Zone (EZ) has remained as dry as it is today. A widely accepted hypothesis suggests that the flora and fauna of arid regions were more broadly distributed before aridification began. In Australia, this process started around 20 million years ago (Ma), leading to gradual speciation as the climate became increasingly arid. Here, we use genomic data to investigate the biogeography and timing of divergence of native allotetraploid tobaccos, Nicotiana section Suaveolentes (Solanaceae). The original allotetraploid migrants from South America were adapted to mesic areas of Australia and recently radiated in the EZ, including in sandy dune fields (only 1.2 Ma old), after developing drought adaptations. Coalescent and maximum likelihood analyses suggest that Nicotiana section Suaveolentes arrived on the continent around 6 Ma, with the ancestors of the Pilbara (Western Australian) lineages radiating there at the onset of extreme aridity 5 Ma by locally adapting to these various ancient, highly stable habitats. The Pilbara thus served as both a mesic refugium and cradle for adaptations to harsher conditions, due to its high topographical diversity, providing microhabitats with varying moisture levels and its proximity to the ocean, which buffers against extreme aridity. This enabled species like Nicotiana to survive in mesic refugia and subsequently adapt to more arid conditions. These results demonstrate that initially poorly adapted plant groups can develop novel adaptations in situ, permitting extensive and rapid dispersal despite the highly variable and unpredictable extreme conditions of the EZ.

Genetic diversity of a recovering European roller (Coracias garrulus) population from Serbia

The European Roller (Coracias garrulus), a long-distance migratory bird, faced a considerable decline in breeding pairs throughout Europe at the end of the 20th century. Due to conservation efforts and the installation of nesting boxes, the population of the European Roller in Serbia has made a remarkable recovery. Here, we used the variability of nucleotide sequences of the mitochondrial DNA (mtDNA) control region and 10 microsatellite loci to assess the genetic diversity and structuring, phylogeographic patterns and demographic history of this species using 224 individuals from Serbia. Our results showed moderate level of genetic diversity (HO = 0.392) and a slightly elevated level of inbreeding and homozygosity (FIS = 0.393). Genetic structuring based on microsatellite data indicated three genetic clusters, but without a clear spatial pattern. High haplotype diversity (Hd = 0.987) of the mtDNA control region sequences was detected, and neutrality tests indicated a recent demographic expansion. The phylogeographic analysis, which also included previously published sequences of the mtDNA control region, supported the subdivision into two distinct European and Asian haplogroups (ΦST = 0.712). However, the results of our study showed that a larger number of haplotypes sampled in Serbia are clustered in the Asian haplogroup as compared to previous studies, indicating a historically continuous distribution of this species and possibly a wider distribution of the subspecies Coracias garrulus semenovwi. Our results suggest that the European Roller population in Serbia is genetically stable, with no evidence of recent bottlenecks, and emphasize the importance of artificial nest boxes for promoting and maintaining population dynamics of European Rollers.

Genome sequencing and CAZymes repertoire analysis of Diaporthe eres P3-1W causing postharvest fruit rot of 'Hongyang' kiwifruit in China

Postharvest rot caused by various fungal pathogens is a damaging disease affecting kiwifruit production and quality, resulting in significant annual economic losses. This study focused on isolating the strain P3-1W, identified as Diaporthe eres, as the causal agent of 'Hongyang' postharvest rot disease in China. The investigation highlighted cell wall degrading enzymes (CWDEs) as crucial pathogenic factors. Specially, the enzymatic activities of cellulase, β-galactosidase, polygalacturonase, and pectin methylesterases peaked significantly on the second day after infection of D. eres P3-1W. To gain a comprehensive understanding of these CWDEs, the genome of this strain was sequenced using PacBio and Illumina sequencing technologies. The analysis revealed that the genome of D. eres P3-1W spans 58,489,835 bp, with an N50 of 5,939,879 bp and a GC content of 50.7%. A total of 15,407 total protein-coding genes (PCGs) were predicted and functionally annotated. Notably, 857 carbohydrate-active enzymes (CAZymes) were identified in D. eres P3-1W, with 521 CWDEs consisting of 374 glycoside hydrolases (GHs), 108 carbohydrate esterase (CEs) and 91 polysaccharide lyases (PLs). Additionally, 221 auxiliary activities (AAs), 91 glycosyltransferases (GTs), and 108 carbohydrate binding modules (CBMs) were detected. These findings offer valuable insights into the CAZymes of D. eres P3-1W.

Combining distribution modelling and phylogeography to understand present, past and future of an endangered spider

BACKGROUND

Understanding how endangered species respond to climatic changes is fundamental for their conservation. Due to its restricted geographic range, its sensitivity to the ongoing global warming and its continuing decline, the Southwestern-Alpine endemic wolf spider Vesubia jugorum is currently classified as Endangered in the IUCN Red List. Here, we combined species distribution modelling (SDM) and phylogeographic inference to describe the present, the past and the future of this species in light of the mtDNA genetic structure of extant populations.

RESULTS

Phylogenetic and network analyses show a high level of genetic differentiation and a strong genetic structure of the populations, likely explicable by a long history of isolation and survival in separate refugia. The SDM projection into past climatic conditions supports these results by showing a smaller distribution range compared to present, mostly restricted to the Maritime and Ligurian Alps, which possibly served as main refugium. Future forecast shows a significant shift in the bioclimatic range towards higher altitudes and latitudes, with a drastic decrease of habitat suitability in the central and south-eastern parts of the range, with consequent general loss of haplotype diversity.

CONCLUSION

SDM and phylogeographic inference support the hypothesis that the current distribution and the genetic structure of the extant populations mirror the survival in situ of Vesubia jugorum across repeated glacial and interglacial phases, in line with the 'long-term stability hypothesis'. Future predictions show a significant shift in the bioclimatic range that V. jugorum will be likely unable to track, with profound impact on its long-term survival and its genetic diversity. Our considerations have implication for conservation genetics, highlighting the pivotal role of the transboundary protected areas of the SW-Alps in promoting conservation efforts for this species.

Bayesian tip-dated phylogeny and biogeography of Cissampelideae (Menispermaceae): Mitigating the effects of homoplastic morphological characters

The integration of morphological and molecular data is essential to understand the affinities of fossil taxa and spatio-temporal evolutionary processes of organisms. However, homoplastic morphological characters can mislead the placement of fossil taxa and impact downstream analyses. Here, we provide an example of how to mitigate effectively the effect of morphological homoplasy on the placement of fossil taxa and biogeographic inferences of Cissampelideae. We assembled three data types, morphological data only, morphological data with a molecular scaffold and combined morphological and molecular data. By removing high-level homoplastic morphological data or reweighting the morphological characters, we conducted 15 parsimony, 12 undated Bayesian and four dated Bayesian analyses. Our results show that the 14 selected Cissampelideae fossil taxa are placed poorly when based only on morphological data, but the addition of molecular scaffold and combination of morphological and molecular data greatly improve the resolution of fossil nodes. We raise the monotypic Stephania subg. Botryodiscia to generic status and discover that three fossils previously assigned to Stephania should be members of Diploclisia. The Bayesian tip-dated tree recovered by removing homoplastic morphological characters with a Rescaled Consistency Index <0.25 has the highest stratigraphic fit and consequently generates more reasonable biogeographic reconstruction for Cissampelideae. Cissampelideae began to diversify in Asia in the latest Cretaceous and subsequently dispersed to South America around the Cretaceous-Palaeogene boundary. Two dispersal events from Asia to Africa occurred in the Early Eocene and the Late Eocene-Late Oligocene, respectively. These findings provide guidelines and practical methods for mitigating the effects of homoplastic morphological characters on fossil placements and Bayesian tip-dating, as well as insights into the past tropical floristic exchanges among different continents.

Annual temperature, body size, and sexual size dimorphism in the evolution of Pyrgomorphidae

In many animal species, larger body size is positively correlated with male mating success and female fecundity. However, in the case of insects, in high seasonality environments, natural selection favors a faster maturation that decreases the risk of pre-reproductive death. However, this advantageous adaptation comes at a tradeoff, resulting in a reduction in body size. Maturation time is influenced by environmental factors, such as temperature, season length, and food availability during the rains. The geographic variation in these parameters provides an opportunity to study their impact on the adaptive evolution of body size in Pyrgomorphidae grasshoppers. These grasshoppers exhibit remarkable variation in body size and wing development and can be found in diverse plant communities across Africa, Asia, Australia, and tropical America. In this study, we utilized a phylogenetic approach to examine the evolution of body size, considering climatic factors, and the influence of sexual selection on size differences between males and females. We found a positive correlation between mean annual temperature and sexual size dimorphism (SSD). Remarkably, species exhibiting a strong bias toward larger females were found to be adapted to regions with higher temperatures. In the Pyrgomorphidae family, an intermediate body size was identified as the ancestral trait. Additionally, winged male and female grasshoppers were observed to be larger than their wingless counterparts. Despite the potential conflicting pressures on body size in males and females, these grasshoppers adhere to Rench's Rule, suggesting that sexual selection on males' body size may explain the evolution of SSD.

Evolution, recombination and geographic spreading of global Coxsackievirus A6

BACKGROUND

The increasing incidence of hand, foot, and mouth disease (HFMD) associated with Coxsackievirus A6 (CVA6) has become a very significant public health problem. The aim of this study is to investigate the recombination, geographic transmission, and evolutionary characteristics of the global CVA6.

METHODS

From 2019 to 2022, 73 full-length CVA6 sequences were obtained from HFMD patients in China and analyzed in combination with 1032 published whole genome sequences. Based on this dataset, the phylogenetic features, recombinant diversity, Bayesian phylodynamic characteristics, and key amino acid variations in CVA6 were analyzed.

RESULTS

The four genotypes of CVA6, A, D, E, and F, are divided into 24 recombinant forms (RFs, RF-A - RF-X) based on differences in the P3 coding region. The eastern China region plays a key role in the dissemination of CVA6 in China. VP1-137 and VP1-138 are located in the DE loop on the surface of the CVA6 VP1 protein, with the former being a highly variable site and the latter having more non-synonymous substitutions.

CONCLUSIONS

Based on whole genome sequences, this study contributes to the CVA6 monitoring, early warning, and the pathogenic mechanism by studying recombination diversity, geographical transmission characteristics, and the variation of important amino acid sites.

Endless forms most frustrating: disentangling species boundaries in the Ramalina decipiens group (Lecanoromycetes, Ascomycota), with the description of six new species and a key to the group

Oceanic islands have been recognized as natural laboratories in which to study a great variety of evolutionary processes. One such process is evolutionary radiations, the diversification of a single ancestor into a number of species that inhabit different environments and differ in the traits that allow them to exploit those environments. The factors that drive evolutionary radiations have been studied for decades in charismatic organisms such as birds or lizards, but are lacking in lichen-forming fungi, despite recent reports of some lineages showing diversification patterns congruent with radiation. Here we propose the Ramalina decipiens group as a model system in which to carry out such studies. This group is currently thought to be comprised of five saxicolous species, all of them endemic to the Macaronesian region (the Azores, Madeira, Selvagens, Canary and Cape Verde islands). Three species are single-island endemics (a rare geographic distribution pattern in lichens), whereas two are widespread and show extreme morphological variation. The latter are suspected to harbor unrecognized species-level lineages. In order to use the Ramalina decipiens group as a model system it is necessary to resolve the group's phylogeny and to clarify its species boundaries. In this study we attempt to do so following an integrative taxonomy approach. We constructed a phylogenetic tree based on six molecular markers, four of which are newly developed and generated competing species hypotheses based on molecular (species discovery strategies based on both single locus and multilocus datasets) and phenotypic data (unsupervised clustering algorithms based on morphology, secondary chemistry and geographic origin). We found that taxonomic diversity in the Ramalina decipiens group has been highly underestimated in previous studies. In consequence, we describe six new species, most of them single-island endemics and provide a key to the group. Phylogenetic relationships among species have been reconstructed with almost full support which, coupled with the endemic character of the group, makes it an excellent system for the study of island radiations in lichen-forming fungi. Citation: Blázquez M, Pérez-Vargas I, Garrido-Benavent I, et al. 2024. Endless forms most frustrating: disentangling species boundaries in the Ramalina decipiens group (Lecanoromycetes, Ascomycota), with the description of six new species and a key to the group. Persoonia 52: 44-93. https://doi.org/10.3767/persoonia.2024.52.03 .

Immune imprinting in early life shapes cross-reactivity to influenza B virus haemagglutinin

Influenza exposures early in life are believed to shape future susceptibility to influenza infections by imprinting immunological biases that affect cross-reactivity to future influenza viruses. However, direct serological evidence linked to susceptibility is limited. Here we analysed haemagglutination-inhibition titres in 1,451 cross-sectional samples collected between 1992 and 2020, from individuals born between 1917 and 2008, against influenza B virus (IBV) isolates from 1940 to 2021. We included testing of 'future' isolates that circulated after sample collection. We show that immunological biases are conferred by early life IBV infection and result in lineage-specific cross-reactivity of a birth cohort towards future IBV isolates. This translates into differential estimates of susceptibility between birth cohorts towards the B/Yamagata and B/Victoria lineages, predicting lineage-specific birth-cohort distributions of observed medically attended IBV infections. Our data suggest that immunological measurements of imprinting could be important in modelling and predicting virus epidemiology.

Mitochondrial Genomes of Streptopelia decaocto: Insights into Columbidae Phylogeny

In this research, the mitochondrial genome of the Streptopelia decaocto was sequenced and examined for the first time to enhance the comprehension of the phylogenetic relationships within the Columbidae. The complete mitochondrial genome of Streptopelia decaocto (17,160 bp) was structurally similar to the recognized members of the Columbidae family, but with minor differences in gene size and arrangement. The structural AT content was 54.12%. Additionally, 150 mitochondrial datasets, representing valid species, were amassed in this investigation. Maximum likelihood (ML) and Bayesian inference (BI) phylogenetic trees and evolutionary time relationships of species were reconstructed based on cytb gene sequences. The findings from the phylogenetic evaluations suggest that the S. decaocto was classified under the Columbinae subfamily, diverging from the Miocene approximately 8.1 million years ago, indicating intricate evolutionary connections with its close relatives, implying a history of species divergence and geographic isolation. The diversification of the Columbidae commenced during the Late Oligocene and extended into the Miocene. This exploration offers crucial molecular data for the S. decaocto, facilitating the systematic taxonomic examination of the Columbidae and Columbiformes, and establishing a scientific foundation for species preservation and genetic resource management.

Genotype Analysis of Respiratory Syncytial Virus Before and After the COVID-19 Pandemic Using Whole-Genome Sequencing: A Prospective, Single-Center Study in Korea From 2019 to 2022

BACKGROUND

Respiratory syncytial virus (RSV), a highly transmissible virus, is the leading cause of lower respiratory tract infections. We examined molecular changes in the RSV genome before and after the coronavirus disease 2019 (COVID-19) pandemic in Korea, and investigated whether drug-resistant mutations were present.

METHODS

In this prospective, single-center study, RSV-positive respiratory samples were collected between September 2019 and December 2022. Long-read whole-genome sequencing (WGS) was performed, and the presence of known drug-resistant substitutions for palivizumab, nirsevimab, and suptavumab was investigated.

RESULTS

Overall, 288 respiratory samples were collected from 276 children. WGS data were available for 133 samples (71 and 62 samples from the pre- and post-pandemic periods, respectively). All RSV-A strains (n = 56) belonged to the GA2.3.5 (ON1) genotype, whereas all RSV-B strains (n = 77) belonged to the GB5.0.5a (BA) genotype. No significant differences in genotypes were observed between the pre- and post-pandemic periods. In addition, no notable mutations related to nirsevimab or palivizumab resistance were detected in the F gene. However, the L172Q and S173L substitutions, which are known to confer resistance to suptavumab, were present in all RSV-B samples.

CONCLUSION

Despite the unprecedented interruption of RSV seasonality, there were no significant molecular changes in circulating RSV strains in Korea related to nirsevimab or palivizumab resistance before and after the COVID-19 pandemic. However, RSV-specific drug-resistance substitutions for suptavumab were identified.

Identification and Characterization of an Alphacoronavirus in Rhinolophus sinicus and a Betacoronavirus in Apodemus ilex in Yunnan, China

Coronaviruses (CoVs), the largest positive-sense RNA viruses, have caused infections in both humans and animals. The cross-species transmission of CoVs poses a serious threat to public health. Rodents and bats, the two largest orders of mammals, serve as significant natural reservoirs for CoVs. It is important to monitor the CoVs carried by bats and rodents. In this study, we collected 410 fecal samples from bats and 74 intestinal samples from rats in Yunnan Province, China. Using RT-PCR, we identified one positive sample for alphacoronavirus (TC-14) from Rhinolophus sinicus (Chinese rufous horseshoe bat) and two positive samples for betacoronavirus (GS-53, GS-56) from Apodemus ilex (Rodentia: Muridae). We successfully characterized the complete genomes of TC-14 and GS-56. Phylogenetic analysis revealed that TC-14 clustered with bat CoV HKU2 and SADS-CoV, while GS-56 was closely related to rat CoV HKU24. The identification of positive selection sites and estimation of divergence dates further helped characterize the genetic evolution of TC-14 and GS-56. In summary, this research reveals the genetic evolution characteristics of TC-14 and GS-56, providing valuable references for the study of CoVs carried by bats and rodents in Yunnan Province.

Convergent evolution toward a slow pace of life predisposes insular endotherms to anthropogenic extinctions

Island vertebrates have evolved a number of morphological, physiological, and life history characteristics that set them apart from their mainland relatives. However, to date, the evolution of metabolism and its impact on the vulnerability to extinction of insular vertebrates remains poorly understood. This study used metabolic data from 2813 species of tetrapod vertebrates, including 695 ectothermic and 2118 endothermic species, to reveal that island mammals and birds evolved convergent metabolic strategies toward a slow pace of life. Insularity was associated with shifts toward slower metabolic rates and greater generation lengths in endotherms, while insularity just drove the evolution of longer generation lengths in ectotherms. Notably, a slow pace of life has exacerbated the extinction of insular endemic species in the face of anthropogenic threats. These findings have important implications for understanding physiological adaptations associated with the island syndrome and formulating conservation strategies across taxonomic groups with different metabolic modes.

Molecular evolutionary analysis of novel NSP4 mono-reassortant G1P[8]-E2 rotavirus strains that caused a discontinuous epidemic in Japan in 2015 and 2018

In the 2010s, several unusual rotavirus strains emerged, causing epidemics worldwide. This study reports a comprehensive molecular epidemiological study of rotaviruses in Japan based on full-genome analysis. From 2014 to 2019, a total of 489 rotavirus-positive stool specimens were identified, and the associated viral genomes were analyzed by next-generation sequencing. The genotype constellations of those strains were classified into nine patterns (G1P[8] (Wa), G1P[8]-E2, G1P[8] (DS-1), G2P[4] (DS-1), G3P[8] (Wa), G3P[8] (DS-1), G8P[8] (DS-1), G9P[8] (Wa), and G9P[8]-E2). The major prevalent genotype differed by year, comprising G8P[8] (DS-1) (37% of that year's isolates) in 2014, G1P[8] (DS-1) (65%) in 2015, G9P[8] (Wa) (72%) in 2016, G3P[8] (DS-1) (66%) in 2017, G1P[8]-E2 (53%) in 2018, and G9P[8] (Wa) (26%) in 2019. The G1P[8]-E2 strains (G1-P[8]-I1-R1-C1-M1-A1-N1-T1-E2-H1) isolated from a total of 42 specimens in discontinuous years (2015 and 2018), which were the newly-emerged NSP4 mono-reassortant strains. Based on the results of the Bayesian evolutionary analyses, G1P[8]-E2 and G9P[8]-E2 were hypothesized to have been generated from distinct independent inter-genogroup reassortment events. The G1 strains detected in this study were classified into multiple clusters, depending on the year of detection. A comparison of the predicted amino acid sequences of the VP7 epitopes revealed that the G1 strains detected in different years encoded VP7 epitopes harboring distinct mutations. These mutations may be responsible for immune escape and annual changes in the prevalent strains.

Palms of the past: can morphometric phytolith analysis inform deep time evolution and palaeoecology of Arecaceae?

BACKGROUND AND AIMS

Palm fossils are often used as evidence for warm and wet palaeoenvironments, reflecting the affinities of most modern palms. However, several extant palm lineages tolerate cool and/or arid climates, making a clear understanding of the taxonomic composition of ancient palm communities important for reliable palaeoenvironmental inference. However, taxonomically identifiable palm fossils are rare and often confined to specific facies. Although the resolution of taxonomic information they provide remains unclear, phytoliths (microscopic silica bodies) provide a possible solution because of their high preservation potential under conditions where other plant fossils are scarce. We thus evaluate the taxonomic and palaeoenvironmental utility of palm phytoliths.

METHODS

We quantified phytolith morphology of 97 modern palm and other monocot species. Using this dataset, we tested the ability of five common discriminant methods to identify nine major palm clades. We then compiled a dataset of species' climate preferences and tested if they were correlated with phytolith morphology using a phylogenetic comparative approach. Finally, we reconstructed palm communities and palaeoenvironmental conditions at six fossil sites.

KEY RESULTS

Best-performing models correctly identified phytoliths to their clade of origin only 59 % of the time. Although palms were generally distinguished from non-palms, few palm clades were highly distinct, and phytolith morphology was weakly correlated with species' environmental preferences. Reconstructions at all fossil sites suggested that palm communities were dominated by Trachycarpeae and Areceae, with warm, equable climates and high, potentially seasonal rainfall. However, fossil site reconstructions had high uncertainty and often conflicted with other climate proxies.

CONCLUSIONS

While phytolith morphology provides some distinction among palm clades, caution is warranted. Unlike prior spatially restricted studies, our geographically and phylogenetically broad study indicates phytolith morphology may not reliably differentiate most palm taxa in deep time. Nevertheless, it reveals distinct clades, including some likely to be palaeoenvironmentally informative.

Mitochondrial Genome Characteristics Reveal Evolution of Acanthopsetta nadeshnyi (Jordan and Starks, 1904) and Phylogenetic Relationships

In the present study, the mitochondrial genomic characteristics of Acanthopsetta nadeshnyi have been reported and have depicted the phylogenetic relationship among Pleuronectidae. Combined with a comparative analysis of 13 PCGs, the TN93 model was used to review the neutral evolution and habitat evolution catalysis of the mitogenome to verify the distancing and purification selectivity of the mitogenome in Pleuronectidae. At the same time, a species differentiation and classification model based on mitogenome analysis data was established. This study is expected to provide a new perspective on the phylogenetic relationship and taxonomic status of A. nadeshnyi and lay a foundation for further exploration of environmental and biological evolutionary mechanisms.

Complete Chloroplast Genome of Megacarpaea megalocarpa and Comparative Analysis with Related Species from Brassicaceae

Megacarpaea megalocarpa, a perennial herbaceous species belonging to the Brassicaceae family, has potential medicinal value. We isolated and characterized the chloroplast (cp) genome of M. megalocarpa and compared it with closely related species. The chloroplast genome displayed a typical quadripartite structure, spanning 154,877 bp, with an overall guanine-cytosine (GC) content of 36.20%. Additionally, this genome contained 129 genes, 105 simple sequence repeats (SSRs), and 48 long repeat sequences. Significantly, the ycf1 gene exhibited a high degree of polymorphism at the small single copy (SSC) region and the inverted repeat a (IRa) boundary. Despite this polymorphism, relative synonymous codon usage (RSCU) values were found to be similar across species, and no large segment rearrangements or inversions were detected. The large single copy (LSC) and SSC regions showed higher sequence variations and nucleotide polymorphisms compared to the IR region. Thirteen distinct hotspot regions were identified as potential molecular markers. Our selection pressure analysis revealed that the protein-coding gene rpl20 is subjected to different selection pressures in various species. Phylogenetic analysis positioned M. megalocarpa within the expanded lineage II of the Brassicaceae family. The estimated divergence time suggests that M. megalocarpa diverged approximately 4.97 million years ago. In summary, this study provides crucial baseline information for the molecular identification, phylogenetic relationships, conservation efforts, and utilization of wild resources in Megacarpaea.

N-glycosylation on hemagglutinin head reveals inter-branch antigenic variability of avian influenza virus H5-subtypes

H5-subtype avian influenza virus (AIV) is globally prevalent and undergoes frequent antigenic drift, necessitating regular updates to vaccines. One of the many influencing elements that cause incompatibility between vaccinations and epidemic strains is the dynamic alteration of glycosylation sites. However, the biological significance of N-glycosylation in the viral evolution and antigenic changes is unclear. Here, we performed a systematic analysis of glycosylation sites on the HA1 subunit of H5N1, providing insights into the changes of primary glycosylation sites, including 140 N, 156 N, and 170 N within the antigenic epitopes of HA1 protein. Multiple recombinant viruses were then generated based on HA genes of historical vaccine strains and deactivated for immunizing SPF chickens. Inactivated recombinant strains showed relatively closer antigenicity compared to which has identical N-glycosylation patterns. The N-glycosylation modification discrepancy highlights the inter-branch antigenic diversity of H5-subtype viruses in avian influenza and serves as a vital foundation for improving vaccination tactics.

The first complete genome of the extinct European wild ass (Equus hemionus hydruntinus)

We present palaeogenomes of three morphologically unidentified Anatolian equids dating to the first millennium BCE, sequenced to a coverage of 0.6-6.4×. Mitochondrial DNA haplotypes of the Anatolian individuals clustered with those of Equus hydruntinus (or Equus hemionus hydruntinus), the extinct European wild ass, secular name 'hydruntine'. Further, the Anatolian wild ass whole genome profiles fell outside the genomic diversity of other extant and past Asiatic wild ass (E. hemionus) lineages. These observations suggest that the three Anatolian wild asses represent hydruntines, making them the latest recorded survivors of this lineage, about a millennium later than the latest observations in the zooarchaeological record. Our mitogenomic and genomic analyses indicate that E. h. hydruntinus was a clade belonging to ancient and present-day E. hemionus lineages that radiated possibly between 0.6 and 0.8 Mya. We also find evidence consistent with recent gene flow between hydruntines and Middle Eastern wild asses. Analyses of genome-wide heterozygosity and runs of homozygosity suggest that the Anatolian wild ass population may have lost genetic diversity by the mid-first millennium BCE, a possible sign of its eventual demise.

Exploring selection signatures in the divergence and evolution of lipid droplet (LD) associated genes in major oilseed crops

BACKGROUND

Oil bodies or lipid droplets (LDs) in the cytosol are the subcellular storage compartments of seeds and the sites of lipid metabolism providing energy to the germinating seeds. Major LD-associated proteins are lipoxygenases, phospholipaseD, oleosins, TAG-lipases, steroleosins, caleosins and SEIPINs; involved in facilitating germination and enhancing peroxidation resulting in off-flavours. However, how natural selection is balancing contradictory processes in lipid-rich seeds remains evasive. The present study was aimed at the prediction of selection signatures among orthologous clades in major oilseeds and the correlation of selection effect with gene expression.

RESULTS

The LD-associated genes from the major oil-bearing crops were analyzed to predict natural selection signatures in phylogenetically close-knit ortholog clusters to understand adaptive evolution. Positive selection was the major force driving the evolution and diversification of orthologs in a lineage-specific manner. Significant positive selection effects were found in 94 genes particularly in oleosin and TAG-lipases, purifying with excess of non-synonymous substitution in 44 genes while 35 genes were neutral to selection effects. No significant selection impact was noticed in Brassicaceae as against LOX genes of oil palm. A heavy load of deleterious mutations affecting selection signatures was detected in T-lineage oleosins and LOX genes of Arachis hypogaea. The T-lineage oleosin genes were involved in mainly anther, tapetum and anther wall morphogenesis. In Ricinus communis and Sesamum indicum > 85% of PLD genes were under selection whereas selection pressures were low in Brassica juncea and Helianthus annuus. Steroleosin, caleosin and SEIPINs with large roles in lipid droplet organization expressed mostly in seeds and were under considerable positive selection pressures. Expression divergence was evident among paralogs and homeologs with one gene attaining functional superiority compared to the other. The LOX gene Glyma.13g347500 associated with off-flavor was not expressed during germination, rather its paralog Glyma.13g347600 showed expression in Glycine max. PLD-α genes were expressed on all the tissues except the seed,δ genes in seed and meristem while β and γ genes expressed in the leaf.

CONCLUSIONS

The genes involved in seed germination and lipid metabolism were under strong positive selection, although species differences were discernable. The present study identifies suitable candidate genes enhancing seed oil content and germination wherein directional selection can become more fruitful.

Phylogeographic patterns of Cyphocharax from trans-Andean rivers and northward expansion to lower Central America (Teleostei, Curimatidae)

Phylogenetic analyses of mitochondrial and nuclear data of 31 specimens of Cyphocharax from trans-Andean rivers support the presence of one lineage of Cyphocharax aspilos in Lago Maracaibo and three cryptic lineages of Cyphocharax magdalenae: (1) Cauca-Magdalena and Ranchería, (2) León and Atrato, and (3) Chucunaque-Tuira, Santa María, and Chiriquí basins of Central America. Results suggest that the Serranía del Perijá facilitated Late Miocene cladogenetic events, whereas post-Isthmian C. magdalenae expansion was enabled by gene flow across the lower Magdalena valley and Central American lowlands. Time-calibrated phylogenetics indicate that the C. magdalenae colonized lower Central America in the Pliocene (3.7 MYA; Ma), the divergence Atrato-Magdalena occurred in Late Pliocene (3.0 Ma) and the split Ranchería-Magdalena during the Middle Pleistocene (1.3 Ma). Updated geographic distribution data support the hypothesis that the Cordillera de Talamanca functions as a barrier to northward expansion of C. magdalenae in Central America.

Bacteriophage-driven emergence and expansion of Staphylococcus aureus in rodent populations

Human activities such as agriculturalization and domestication have led to the emergence of many new pathogens via host-switching events between humans, domesticated and wild animals. Staphylococcus aureus is a multi-host opportunistic pathogen with a global healthcare and economic burden. Recently, it was discovered that laboratory and wild rodents can be colonised and infected with S. aureus, but the origins and zoonotic potential of rodent S. aureus is unknown. In order to trace their evolutionary history, we employed a dataset of 1249 S. aureus genome sequences including 393 of isolates from rodents and other small mammals (including newly determined sequences for 305 isolates from 7 countries). Among laboratory mouse populations, we identified multiple widespread rodent-specific S. aureus clones that likely originated in humans. Phylogeographic analysis of the most common murine lineage CC88 suggests that it emerged in the 1980s in laboratory mouse facilities most likely in North America, from where it spread to institutions around the world, via the distribution of mice for research. In contrast, wild rodents (mice, voles, squirrels) were colonized with a unique complement of S. aureus lineages that are widely disseminated across Europe. In order to investigate the molecular basis for S. aureus adaptation to rodent hosts, genome-wide association analysis was carried out revealing a unique complement of bacteriophages associated with a rodent host ecology. Of note, we identified novel prophages and pathogenicity islands in rodent-derived S. aureus that conferred the potential for coagulation of rodent plasma, a key phenotype of abscess formation and persistence. Our findings highlight the remarkable capacity of S. aureus to expand into new host populations, driven by the acquisition of genes promoting survival in new host-species.

The evolutionary and molecular history of a chikungunya virus outbreak lineage

In 2018-2019, Thailand experienced a nationwide spread of chikungunya virus (CHIKV), with approximately 15,000 confirmed cases of disease reported. Here, we investigated the evolutionary and molecular history of the East/Central/South African (ECSA) genotype to determine the origins of the 2018-2019 CHIKV outbreak in Thailand. This was done using newly sequenced clinical samples from travellers returning to Sweden from Thailand in late 2018 and early 2019 and previously published genome sequences. Our phylogeographic analysis showed that before the outbreak in Thailand, the Indian Ocean lineage (IOL) found within the ESCA, had evolved and circulated in East Africa, South Asia, and Southeast Asia for about 15 years. In the first half of 2017, an introduction occurred into Thailand from another South Asian country, most likely Bangladesh, which subsequently developed into a large outbreak in Thailand with export to neighbouring countries. Based on comparative phylogenetic analyses of the complete CHIKV genome and protein modelling, we identified several mutations in the E1/E2 spike complex, such as E1 K211E and E2 V264A, which are highly relevant as they may lead to changes in vector competence, transmission efficiency and pathogenicity of the virus. A number of mutations (E2 G205S, Nsp3 D372E, Nsp2 V793A), that emerged shortly before the outbreak of the virus in Thailand in 2018 may have altered antibody binding and recognition due to their position. This study not only improves our understanding of the factors contributing to the epidemic in Southeast Asia, but also has implications for the development of effective response strategies and the potential development of new vaccines.

Molecular epidemiology of HIV among people who inject drugs after the HIV-outbreak in Athens, Greece: Evidence for a 'slow burn' outbreak

BACKGROUND

New diagnoses of HIV-1 infection among people who inject drugs (PWID) in Athens, Greece, saw a significant increase in 2011 and a subsequent decline after 2013. Despite this, ongoing HIV-1 transmission persisted from 2014 to 2020 within this population. Our objective was to estimate the time of infection for PWID in Athens following the HIV-1 outbreak, explore the patterns of HIV-1 dispersal over time, and determine the duration from infection to diagnosis.

METHODS

Time from HIV-1 infection to diagnosis was estimated for 844 individuals infected within 4 PWID-specific clusters and for 8 PWID infected with sub-subtype A6 diagnosed during 2010-2019. Phylogeny reconstruction was performed using the maximum-likelihood method. HIV-1 infection dates were based on molecular clock calculations.

RESULTS

In total 86 of 92 (93.5%) sequences from PWID diagnosed during 2016-2019 were either related to the previously identified PWID-specific clusters (n = 81) or belonged to a new A6 cluster (n = 5). The median time between infection and diagnosis was 0.42 years during the outbreak period and 0.70 years during 2016-2019 (p < 0.001). The proportion of clustered sequences from PWID was very low at 5.3% during the pre-outbreak period (1998-2009), saw an increase to 41.7% one year before the outbreak in 2010, and consistently remained high during the whole period after 2011, spanning the post-outbreak period (2016-2019) with a range from 92.9% to 100%.

CONCLUSIONS

The substantial proportion of clustered infections (93.5%) during 2016-2019 implies a persistent 'slow burn' HIV outbreak among PWID in Athens, suggesting that the outbreak was not successfully eliminated. The consistently high proportion of clustered sequences since the onset of the outbreak suggests the persistence of ongoing HIV-1 transmission attributed to injection practices. Our findings underscore the importance of targeted interventions among PWID, considering the ongoing transmission rate and prolonged time from infection to diagnosis.

Distinct SARS-CoV-2 populational immune backgrounds tolerate divergent RBD evolutionary preferences

Immune evasion is a pivotal force shaping the evolution of viruses. Nonetheless, the extent to which virus evolution varies among populations with diverse immune backgrounds remains an unsolved mystery. Prior to the widespread SARS-CoV-2 infections in December 2022 and January 2023, the Chinese population possessed a markedly distinct (less potent) immune background due to its low infection rate, compared to countries experiencing multiple infection waves, presenting an unprecedented opportunity to investigate how the virus has evolved under different immune contexts. We compared the mutation spectrum and functional potential of the newly derived mutations that occurred in BA.5.2.48, BF.7.14 and BA.5.2.49-variants prevalent in China-with their counterparts in other countries. We found that the emerging mutations in the receptor-binding-domain region in these lineages were more widely dispersed and evenly distributed across different epitopes. These mutations led to a higher angiotensin-converting enzyme 2 (ACE2) binding affinity and reduced potential for immune evasion compared to their counterparts in other countries. These findings suggest a milder immune pressure and less evident immune imprinting within the Chinese population. Despite the emergence of numerous immune-evading variants in China, none of them outcompeted the original strain until the arrival of the XBB variant, which had stronger immune evasion and subsequently outcompeted all circulating variants. Our findings demonstrated that the continuously changing immune background led to varying evolutionary pressures on SARS-CoV-2. Thus, in addition to viral genome surveillance, immune background surveillance is also imperative for predicting forthcoming mutations and understanding how these variants spread in the population.

Agaricus subgenus Pseudochitonia in Malakand, Pakistan: An updated phylogeny and description of three new species

Agaricus is a genus with more than 500 species. Most of the new species reported since 2000 are tropical or subtropical. The study area, the Malakand region, located in the north of Pakistan, has a subtropical climate. In this study, nine species, including three new species, of Agaricus subgenus Pseudochitonia, are reported from this region. Description of the new species are based on morphological characteristics and phylogenetic analyses using three DNA regions: nuc ribosomal DNA internal transcribed spacers (ITS), fragments of the large subunit of nuc ribosomal DNA (28S), and the translation elongation factor 1-alpha gene (TEF1). One new species, Agaricus lanosus, with wooly squamules on its cap, forms a lineage within Agaricus sect. Bivelares and cannot be classified with certainty in one of the two subsections (Cupressorum and Hortenses) of this section. Agaricus rhizoideus with rhizoid-like structure at the base of the stipe forms a basal clade in Agaricus sect. Hondenses. Specimens of the third new species, Agaricus malakandensis, form a species-level clade within Agaricus sect. Catenulati and exhibits the morphological characteristics of this section. Due to their similar ITS sequences, two previously unnamed specimens from Thailand (A. sp. LD2012162 and CA799) are considered conspecific with A. malakandensis.

Complete chloroplast genomes of the hemiparasitic genus Cymbaria: Insights into comparative analysis, development of molecular markers, and phylogenetic relationships

The hemiparasitic tribe Cymbarieae (Orobanchaceae) plays a crucial role in elucidating the initial stage of the transition from autotrophism to heterotrophism. However, the complete chloroplast genome of the type genus Cymbaria has yet to be reported. In addition, the traditional Mongolian medicine Cymbaria daurica is frequently subjected to adulteration or substitution because of the minor morphological differences with Cymbaria mongolica. In this study, the complete chloroplast genomes of the two Cymbaria species were assembled and annotated, and those of other published 52 Orobanchaceae species were retrieved for comparative analyses. We found that the Cymbaria chloroplast genomes are characterized by pseudogenization or loss of stress-relevant genes (ndh) and a unique rbcL-matK inversion. Unlike the high variability observed in holoparasites, Cymbaria and other hemiparasites exhibit high similarity to autotrophs in genome size, guanine-cytosine (GC) content, and intact genes. Notably, four pairs of specific DNA barcodes were developed and validated to distinguish the medicinal herb from its adulterants. Phylogenetic analyses revealed that the genus Cymbaria and the Schwalbea-Siphonostegia clade are grouped into the tribe Cymbarieae, which forms a sister clade to the remaining Orobanchaceae parasitic lineages. Moreover, the diversification of monophyletic Cymbaria occurred during the late Miocene (6.72 Mya) in the Mongol-Chinese steppe region. Our findings provide valuable genetic resources for studying the phylogeny of Orobanchaceae and plant parasitism, and genetic tools to validate the authenticity of the traditional Mongolian medicine "Xinba.".

A reproducible workflow for assembling the mitochondrial genome of Acheta domesticus (Orthoptera: Gryllidae)

In this study, we report the assembly and annotation of the mitochondrial genome (mitogenome) of Acheta domesticus from breeding facility, a species commonly known as the house cricket. This species is considered to be an important edible cricket. The mitogenome was assembled using a reproducible protocol implemented on the Galaxy Europe Server, which involved uploading paired-end fastq reads for bioinformatic analysis. The resulting mitogenome is 15,784 base pairs in length and has a GC content of 29.05%. The nucleotide composition of this mitogenome is similar to that of other insect mitogenomes, with A, T, C, and G nucleotides comprising 39.2%, 31.7%, 19.6%, and 9.5% of the mitogenome, respectively. The gene organization of the A. domesticus mitogenome is identical to that of other cricket species. The mitogenome consists of 37 genes, including 13 protein-coding genes, 22 tRNA genes, and two rRNA genes. The congruence between PCA and Bayesian evolutionary tree analysis in clustering the divergent A. domesticus sequences highlights these genomes as candidates for further study to elucidate their distinct features and evolutionary history.

Structure, Evolution, and Mitochondrial Genome Analysis of Mussel Species (Bivalvia, Mytilidae)

Based on the nucleotide sequences of the mitochondrial genome (mitogenome) of specimens taken from two mussel species (Arcuatula senhousia and Mytilus coruscus), an investigation was performed by means of the complex approaches of the genomics, molecular phylogenetics, and evolutionary genetics. The mitogenome structure of studied mussels, like in many other invertebrates, appears to be much more variable than in vertebrates and includes changing gene order, duplications, and deletions, which were most frequent for tRNA genes; the mussel species' mitogenomes also have variable sizes. The results demonstrate some of the very important properties of protein polypeptides, such as hydrophobicity and its determination by the purine and pyrimidine nucleotide ratio. This fact might indirectly indicate the necessity of purifying natural selection for the support of polypeptide functionality. However, in accordance with the widely accepted and logical concept of natural cutoff selection for organisms living in nature, which explains its action against deleterious nucleotide substitutions in the nonsynonymous codons (mutations) and its holding of the active (effective) macromolecules of the polypeptides in a population, we were unable to get unambiguous evidence in favor of this concept in the current paper. Here, the phylogeny and systematics of mussel species from one of the largest taxons of bivalve mollusks are studied, the family known as Mytilidae. The phylogeny for Mytilidae (order Mytilida), which currently has no consensus in terms of systematics, is reconstructed using a data matrix of 26-27 mitogenomes. Initially, a set of 100 sequences from GenBank were downloaded and checked for their gender: whether they were female (F) or male (M) in origin. Our analysis of the new data confirms the known drastic differences between the F/M mitogenome lines in mussels. Phylogenetic reconstructions of the F-lines were performed using the combined set of genetic markers, reconstructing only protein-coding genes (PCGs), only rRNA + tRNA genes, and all genes. Additionally, the analysis includes the usage of nucleotide sequences composed of other data matrices, such as 20-68 mitogenome sequences. The time of divergence from MRCA, estimated via BEAST2, for Mytilidae is close to 293 Mya, suggesting that they originate in the Silurian Period. From all these data, a consensus for the phylogeny of the subfamily of Mytilinae and its systematics is suggested. In particular, the long-debated argument on mussel systematics was resolved as to whether Mytilidae, and the subfamily of Mytilinae, are monophyletic. The topology signal, which was strongly resolved in this paper and in the literature, has refuted the theory regarding the monophyly of Mytilinae.

Evolution of DS-1-like G8P[8] rotavirus A strains from Vietnamese children with acute gastroenteritis (2014-21): Adaptation and loss of animal rotavirus-derived genes during human-to-human spread

Animal rotaviruses A (RVAs) are considered the source of emerging, novel RVA strains that have the potential to cause global spread in humans. A case in point was the emergence of G8 bovine RVA consisting of the P[8] VP4 gene and the DS-1-like backbone genes that appeared to have jumped into humans recently. However, it was not well documented what evolutionary changes occurred on the animal RVA-derived genes during circulation in humans. Rotavirus surveillance in Vietnam found that DS-1-like G8P[8] strains emerged in 2014, circulated in two prevalent waves, and disappeared in 2021. This surveillance provided us with a unique opportunity to investigate the whole process of evolutionary changes, which occurred in an animal RVA that had jumped the host species barrier. Of the 843 G8P[8] samples collected from children with acute diarrhoea in Vietnam between 2014 and 2021, fifty-eight strains were selected based on their distinctive electropherotypes of the genomic RNA identified using polyacrylamide gel electrophoresis. Whole-genome sequence analysis of those fifty-eight strains showed that the strains dominant during the first wave of prevalence (2014-17) carried animal RVA-derived VP1, NSP2, and NSP4 genes. However, the strains from the second wave of prevalence (2018-21) lost these genes, which were replaced with cognate human RVA-derived genes, thus creating strain with G8P[8] on a fully DS-1-like human RVA gene backbone. The G8 VP7 and P[8] VP4 genes underwent some point mutations but the phylogenetic lineages to which they belonged remained unchanged. We, therefore, propose a hypothesis regarding the tendency for the animal RVA-derived genes to be expelled from the backbone genes of the progeny strains after crossing the host species barrier. This study underlines the importance of long-term surveillance of circulating wild-type strains in order to better understand the adaptation process and the fate of newly emerging, animal-derived RVA among the human population. Further studies are warranted to disclose the molecular mechanisms by which spillover animal RVAs become readily transmissible among humans, and the roles played by the expulsion of animal-derived genes and herd immunity formed in the local population.

Annexin gene family in Spirometra mansoni (Cestoda: Diphyllobothriidae) and its phylogenetic pattern among Platyhelminthes of medical interest

The plerocercoid larvae of Spirometra mansoni are etiological agents of human and animal sparganosis. Annexins are proteins with important roles in parasites. However, our knowledge of annexins in S. mansoni is still inadequate. In this study, 18 new members of the Annexin (ANX) family were characterized in S. mansoni. The clustering analysis demonstrated that all the SmANXs were divided into two main classes, consistent with the patterns of conserved motif organization. The 18 SmANXs were detected at all developmental stages (plerocercoid, adult, and egg) and displayed ubiquitous but highly variable expression patterns in all tissues/organs studied. The representative member rSmANX18 was successfully cloned and expressed. The protein was immunolocalized in the tegument and parenchyma of the plerocercoid and in the tegument, parenchyma, uterus and egg shell of adult worms. The recombinant protein can bind phospholipids with high affinity in a Ca2+-dependent manner, shows high anticoagulant activity and combines with FITC to recognize apoptotic cells. Annexin gene polymorphism and conservative core motif permutation were found in both cestodes and trematodes. SmANXs also revealed high genetic diversity among Platyhelminthes of medical interest. Our findings lay a foundation for further studies on the biological functions of ANXs in S. mansoni as well as other taxa in which ANXs occur.

DNA barcoding reveals cryptic species in the sea slater Ligia italica (Crustacea, Isopoda) from Tunisia

Barcoding studies have provided significant insights into phylogenetic relationships among species belonging to the genus Ligia (Crustacea, Isopoda). Herein the diversity of the Italian sea slater Ligia italica from Tunisia is studied for the first time. Samples were collected from 18 localities in Tunisia, and the analysis included previously published sequences from Italy and Greece available in GenBank. Bayesian and Maximum Likelihood phylogenetic analyses were carried out using a fragment of the mitochondrial COI gene. Putative cryptic species were explored using the 'barcode gap' approach in the software ASAP. A genetic landscape shape analysis was carried out using the program Alleles in Space. The analyses revealed highly divergent and well-supported clades of L. italica dispersed across Tunisia (Clades A1 and A2), Greece (Clade B) and Italy (Clades C1 and C2). High genetic dissimilarity among clades suggested that L. italica constitute a cryptic species complex. Divergence among different L. italica lineages (Clades A, B and C) occurred around 7-4.5 Ma. The detected high genetic distances among clades did not result from atypical mitochondrial DNAs or intracellular infection by Wolbachia bacteria. The complex history of the Mediterranean Sea appears to have played a significant role in shaping the phylogeographic pattern of Ligia italica. Additional morphological and molecular studies are needed to confirm the existence of cryptic species in Ligia italica in Mediterranean.

Distinct effects of phyllosphere and rhizosphere microbes on invader Ageratina adenophora during its early life stages

Microbes strongly affect invasive plant growth. However, how phyllosphere and rhizosphere soil microbes distinctively affect seedling mortality and growth of invaders across ontogeny under varying soil nutrient levels remains unclear. In this study, we used the invader Ageratina adenophora to evaluate these effects. We found that higher proportions of potential pathogens were detected in core microbial taxa in leaf litter than rhizosphere soil and thus leaf inoculation had more adverse effects on seed germination and seedling survival than soil inoculation. Microbial inoculation at different growth stages altered the microbial community and functions of seedlings, and earlier inoculation had a more adverse effect on seedling survival and growth. The soil nutrient level did not affect microbe-mediated seedling growth and the relative abundance of the microbial community and functions involved in seedling growth. The effects of some microbial genera on seedling survival are distinct from those on growth. Moreover, the A. adenophora seedling-killing effects of fungal strains isolated from dead seedlings by non-sterile leaf inoculation exhibited significant phylogenetic signals, by which strains of Allophoma and Alternaria generally caused high seedling mortality. Our study stresses the essential role of A. adenophora litter microbes in population establishment by regulating seedling density and growth.

Antimicrobial resistance and population genomics of emerging multidrug-resistant Salmonella 4,[5],12:i:- in Guangdong, China

Salmonella 4,[5],12:i:-, a monophasic variant of Salmonella Typhimurium, has emerged as a global cause of multidrug-resistant salmonellosis and has become endemic in many developing and developed countries, especially in China. Here, we have sequenced 352 clinical isolates in Guangdong, China, during 2009-2019 and performed a large-scale collection of Salmonella 4,[5],12:i:- with whole genome sequencing (WGS) data across the globe, to better understand the population structure, antimicrobial resistance (AMR) genomic characterization, and transmission routes of Salmonella 4,[5],12:i:- across Guangdong. Salmonella 4,[5],12:i:- strains showed broad genetic diversity; Guangdong isolates were found to be widely distributed among the global lineages. Of note, we identified the formation of a novel Guangdong clade (Bayesian analysis of population structure lineage 1 [BAPS1]) genetically diversified from the global isolates and likely emerged around 1990s. BAPS1 exhibits unique genomic features, including large pan-genome, decreased ciprofloxacin susceptibility due to mutation in gyrA and carriage of plasmid-mediated quinolone resistance (PMQR) genes, and the multidrug-resistant IncHI2 plasmid. Furthermore, high genetic similarity was found between strains collected from Guangdong, Europe, and North America, indicating the association with multiple introductions from overseas. These results suggested that global dissemination and local clonal expansion simultaneously occurred in Guangdong, China, and horizontally acquired resistance to first-line and last-line antimicrobials at local level, underlying emergences of extensive drug and pan-drug resistance. Our findings have increased the knowledge of global and local epidemics of Salmonella 4,[5],12:i:- in Guangdong, China, and provided a comprehensive baseline data set essential for future molecular surveillance.IMPORTANCESalmonella 4,[5],12:i:- has been regarded as the predominant pandemic serotype causing diarrheal diseases globally, while multidrug resistance (MDR) constitutes great public health concerns. This study provided a detailed and comprehensive genome-scale analysis of this important Salmonella serovar in the past decade in Guangdong, China. Our results revealed the complexity of two distinct transmission modes, namely global transmission and local expansion, circulating in Guangdong over a decade. Using phylogeography models, the origin of Salmonella 4,[5],12:i:- was predicted from two aspects, year and country, that is, Salmonella 4,[5],12:i:- emerged in 1983, and was introduced from the UK, and subsequently differentiated into the local endemic lineage circa 1991. Additionally, based on the pan-genome analysis, it was found that the gene accumulation rate in local endemic BAPS 1 lineage was higher than in other lineages, and the horizontal transmission of MDR IncHI2 plasmid associated with high resistance played a major role, which showed the potential threat to public health.

Genomes of the Orestias pupfish from the Andean Altiplano shed light on their evolutionary history and phylogenetic relationships within Cyprinodontiformes

BACKGROUND

To unravel the evolutionary history of a complex group, a comprehensive reconstruction of its phylogenetic relationships is crucial. This requires meticulous taxon sampling and careful consideration of multiple characters to ensure a complete and accurate reconstruction. The phylogenetic position of the Orestias genus has been estimated partly on unavailable or incomplete information. As a consequence, it was assigned to the family Cyprindontidae, relating this Andean fish to other geographically distant genera distributed in the Mediterranean, Middle East and North and Central America. In this study, using complete genome sequencing, we aim to clarify the phylogenetic position of Orestias within the Cyprinodontiformes order.

RESULTS

We sequenced the genome of three Orestias species from the Andean Altiplano. Our analysis revealed that the small genome size in this genus (~ 0.7 Gb) was caused by a contraction in transposable element (TE) content, particularly in DNA elements and short interspersed nuclear elements (SINEs). Using predicted gene sequences, we generated a phylogenetic tree of Cyprinodontiformes using 902 orthologs extracted from all 32 available genomes as well as three outgroup species. We complemented this analysis with a phylogenetic reconstruction and time calibration considering 12 molecular markers (eight nuclear and four mitochondrial genes) and a stratified taxon sampling to consider 198 species of nearly all families and genera of this order. Overall, our results show that phylogenetic closeness is directly related to geographical distance. Importantly, we found that Orestias is not part of the Cyprinodontidae family, and that it is more closely related to the South American fish fauna, being the Fluviphylacidae the closest sister group.

CONCLUSIONS

The evolutionary history of the Orestias genus is linked to the South American ichthyofauna and it should no longer be considered a member of the Cyprinodontidae family. Instead, we submit that Orestias belongs to the Orestiidae family, as suggested by Freyhof et al. (2017), and that it is the sister group of the Fluviphylacidae family, distributed in the Amazonian and Orinoco basins. These two groups likely diverged during the Late Eocene concomitant with hydrogeological changes in the South American landscape.

Illuminating the coevolution of photosynthesis and Bacteria

Life harnessing light energy transformed the relationship between biology and Earth-bringing a massive flux of organic carbon and oxidants to Earth's surface that gave way to today's organotrophy- and respiration-dominated biosphere. However, our understanding of how life drove this transition has largely relied on the geological record; much remains unresolved due to the complexity and paucity of the genetic record tied to photosynthesis. Here, through holistic phylogenetic comparison of the bacterial domain and all photosynthetic machinery (totally spanning >10,000 genomes), we identify evolutionary congruence between three independent biological systems-bacteria, (bacterio)chlorophyll-mediated light metabolism (chlorophototrophy), and carbon fixation-and uncover their intertwined history. Our analyses uniformly mapped progenitors of extant light-metabolizing machinery (reaction centers, [bacterio]chlorophyll synthases, and magnesium-chelatases) and enzymes facilitating the Calvin-Benson-Bassham cycle (form I RuBisCO and phosphoribulokinase) to the same ancient Terrabacteria organism near the base of the bacterial domain. These phylogenies consistently showed that extant phototrophs ultimately derived light metabolism from this bacterium, the last phototroph common ancestor (LPCA). LPCA was a non-oxygen-generating (anoxygenic) phototroph that already possessed carbon fixation and two reaction centers, a type I analogous to extant forms and a primitive type II. Analyses also indicate chlorophototrophy originated before LPCA. We further reconstructed evolution of chlorophototrophs/chlorophototrophy post-LPCA, including vertical inheritance in Terrabacteria, the rise of oxygen-generating chlorophototrophy in one descendant branch near the Great Oxidation Event, and subsequent emergence of Cyanobacteria. These collectively unveil a detailed view of the coevolution of light metabolism and Bacteria having clear congruence with the geological record.

New insights into population structure, demographic history, and effective population size of the critically endangered blue shark Prionace glauca in the Mediterranean Sea

The blue shark, Prionace glauca, is the most abundant pelagic shark in the open ocean but its vulnerability remains poorly understood while being one of the most fecund sharks. In the Mediterranean Sea, the blue shark is listed as Critically Endangered (CR) by the International Union for Conservation of Nature. The species is facing a strong decline due to fishing, and scientific data regarding its genetic structure and vulnerability are still lacking. Here, we investigated the genetic diversity, demographic history, and population structure of the blue shark within the Mediterranean Sea, from samples of the Gulf of Lion and Malta, using sequences of the mtDNA control region and 22 microsatellite markers. We also compared our mitochondrial data to previous studies to examine the Atlantic-Mediterranean population structure. We assessed the blue shark's genetic vulnerability in the Mediterranean basin by modelling its effective population size. Our results showed a genetic differentiation between the Atlantic and the Mediterranean basins, with limited gene flow between the two areas, and distinct demographic histories making the Mediterranean population an independent management unit. Within the Mediterranean Sea, no sign of population structure was detected, suggesting a single population across the Western and Central parts of the sea. The estimated effective population size was low and highlighted the high vulnerability of the Mediterranean blue shark population, as the estimated size we calculated might not be sufficient to ensure the long-term persistence of the population. Our data also provide additional evidence that the Gulf of Lion area acts as a nursery for P. glauca, where protection is essential for the conservation strategy of the species in the Mediterranean.

Phylogeographic Structure and Population Dynamics of Baoxing Osmanthus (Osmanthus serrulatus), an Endemic Species from the Southwest Sichuan Basin, China

The mountainous regions of southwest China are recognized as pivotal centers for the origin and evolution of Osmanthus species. Baoxing Osmanthus (Osmanthus serrulatus Rehder), a rare and endemic species known for its spring blooms, is sparsely distributed within the high altitude evergreen broad-leaved forests surrounding the southwestern Sichuan Basin. However, persistent anthropogenic disturbances and habitat fragmentation have precipitated a significant decline in its natural population size, leading to the erosion of genetic resources. To assess the genetic status of O. serrulatus and formulate effective conservation strategies, we conducted sampling across ten wild populations, totaling 148 individuals in their natural habitats. We employed two cpDNA fragments (matK and trnS-trnG) to elucidate the phylogeographic structure and historical population dynamics. The results revealed low species-level genetic diversity, alongside pronounced regional differentiation among populations (FST = 0.812, p < 0.05) and a notable phylogeographic structure (NST = 0.698 > GST = 0.396, p < 0.05). Notably, genetic variation was predominantly observed among populations (81.23%), with no evidence of recent demographic expansion across the O. serrulatus distribution range. Furthermore, divergence dating indicated a timeline of approximately 4.85 Mya, corresponding to the late Miocene to early Pleistocene. This temporal correlation coincided with localized uplift events in the southwestern mountains and heightened Asian monsoons, suggesting pivotal roles for these factors in shaping the current phylogeographic pattern of O. serrulatus. These findings support the effective conservation of O. serrulatus germplasm and offer insights into the impact of Quaternary climate oscillations on companion species within evergreen broad-leaved forests. They also enhance our understanding of the origin and evolution of these forests in the southwestern mountains, aiding biodiversity conservation efforts in the region.

Camelus knoblochi genome reveals the complex evolutionary history of Old World camels

Extant Old World camels (genus Camelus) contributed to the economic and cultural exchanges between the East and West for thousands of years.1,2 Although many remains have been unearthed,3,4,5 we know neither whether the prevalent hybridization observed between extant Camelus species2,6,7 also occurred between extinct lineages and the ancestors of extant Camelus species nor why some populations became extinct while others survived. To investigate these questions, we generated paleogenomic and stable isotope data from an extinct two-humped camel species, Camelus knoblochi. We find that in the mitochondrial phylogeny, all C. knoblochi form a paraphyletic group that nests within the diversity of modern, wild two-humped camels (Camelus ferus). In contrast, they are clearly distinguished from both wild and domesticated (Camelus bactrianus) two-humped camels on the nuclear level. Moreover, the divergence pattern of the three camel species approximates a trifurcation, because the most common topology is only slightly more frequent than the two other possible topologies. This mito-nuclear phylogenetic discordance likely arose due to interspecific gene flow between all three species, suggesting that interspecific hybridization is not exclusive to modern camels but a recurrent phenomenon throughout the evolutionary history of the genus Camelus. These results suggest that the genomic complexity of Old World camels' evolutionary history is underestimated when considering data from only modern species. Finally, we find that C. knoblochi populations began declining prior to the last glacial maximum and, by integrating palaeoecological evidence and stable isotope data, suggest that this was likely due to failure to adapt to a changing environment.

Mitochondrial DNA evidence reflects high genetic divergence of Amynthas aspergillum (Oligochaeta: Megascolecidae) in southern China

Amynthas aspergillum (Perrier, 1872), a natural resource used in traditional Chinese medicine (Guang-dilong) with high economic value, is widely distributed in forests and farmland habitats in the hilly areas of southern China. To investigate the extent of genetic differentiation and diversity in A. aspergillum, a population genetic structure study was performed on 157 samples from 75 locations in southern China using the mitochondrial genes COI, COII, 12S rRNA, 16S rRNA, and NDI. The results indicated that A. aspergillum had a high level of genetic diversity, and variation within populations was the main source of the total variation. Six deeply divergent mitochondrial clades (I-VI) were detected using both phylogenetic tree and haplotype network analyses. This finding was supported by the high Kimura two-parameter genetic distance and the pairwise fixation index value obtained based on the COI gene. No significant phylogeographic structures were observed. The widespread geographic distribution of clades II, IV, and VI suggested a recent demographic expansion based on multiple analysis results. These results include a high level of Hd and low π, star-shaped haplotype network structures with a high number of less frequent haplotypes, significantly negative neutrality test values, and a unimodal mismatch distribution pattern. The divergence time estimates and reconstruction of the ancestral area revealed that A. aspergillum originated in Guangxi Province and underwent initial intraspecific diversification in the early Pliocene to generate clade I. Then, it gradually dispersed eastward and rapidly differentiated into clades II-V during the Pleistocene. The Yunnan-Guizhou Plateau and Nanling and Wuyi Mountains might act as geographical barriers for the spread of A. aspergillum to the west and north.