The evolutionary history of hepaciviruses

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

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

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

Evolutionarily diverse origins of deformed wing viruses in western honey bees

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

Genomic Surveillance of Recent Dengue Outbreaks in Colombo, Sri Lanka

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

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

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

Northern Richness, Southern Dead End-Origin and Dispersal Events of Pseudolycoriella (Sciaridae, Diptera) between New Zealand's Main Islands

Sciaridae (Diptera) is a widespread insect family of which some species can reach high abundances in arboreal habitats. This trait, together with their (passive) mobility, enables them to quickly colonise suitable habitats. To reveal the biogeographic history of the New Zealand members of the sciarid genus Pseudolycoriella, we analysed three molecular markers of selected species and populations in a Bayesian approach. At the intra- and interspecific levels, we detected a pattern of northern richness vs. southern purity, which has probably developed as a result of Pleistocene glacial cycles. Since the late Miocene, we identified 13 dispersal events across the sea strait separating New Zealand's main islands. As nine of these dispersal events were south-directed, North Island can be considered the centre of radiation for this genus. An unequivocal re-colonisation of North Island was only observed once. Based on the inclusion of three undescribed species from Tasmania and on previously published data, three colonisations of New Zealand are likely, all of them assumed to be of Australian origin. One of these most probably took place during the late Miocene, and the other two during the late Pliocene or at the Pliocene-Pleistocene boundary.

Phylogeography of 912 Cherry Accessions Insight into Independent Origins of Fruiting Cherries and Domestication Footprints of Cultivated Chinese Cherry (Prunus pseudocerasus Lindl.)

The subgenus Cerasus (Rosaceae) contain numerous fruit trees and ornamentals with high economic values. The origin and genetic divergence among various types of fruiting cherries always remain a perplexing issue. We employed three plastom fragments and ITS sequence matrices derived from 912 cherry accessions to elucidate the phylogeographic structure and genetic relationship among fruiting cherries, as well as the origin and domestication of cultivated Chinese cherry. The integration of haplotype genealogies, Approximate Bayesian computation (ABC) approach and estimation of genetic differentiation within and between different groups and lineages has facilitated the resolution of several previously unresolved questions. Firstly, distant phylogenetic relationships between Cerasus and Microcerasus accessions, as indicated by both nuclear and chloroplast data, suggested independent origins and evolution for these two taxa. Moreover, two distinct geographic origin centers (Europe and China) have been confirmed, with significant phylogeographic signals and high genetic differentiation observed between cherries from these regions. This may be attributed to long-term geographic isolation caused by Himalaya-Hengduan Mountains. Our phylogeographic analyses and ABC analysis suggested that cherries inhabiting in China may have undergone multiple hybridization events during the glacial refugia of the eastern edge and southern Himalaya-Hengduan Mountains, followed by rapid radiation throughout their current habitats during interglacial period. The discrepancy between nuclear and chloroplast data may be attributed to hybridization events and incomplete lineage sorting. Furthermore, we speculated that the domesticated Chinese cherries were derived from wild accessions in Longmenshan Fault Zones approximately 2600 years ago. We have also traced the domestication processes and dispersal routes of cultivated Chinese cherries.

Retrospective Phylodynamic and Phylogeographic Analysis of the Human Papillomavirus 16 E6 Gene in the Mediterranean Region

Human papillomavirus 16 (HPV16) is considered to be strongly correlated with the development of cervical cancer. Among the 8 HPV16 genes, the E6 constitutes a remarkable marker to follow the evolutionary history and spatial phylodynamics of HPV16 in the Mediterranean basin. Thus, this work aims to decipher the major evolutionary events and crosstalks in the Mediterranean basin with a focus on Tunisian strains regarding the E6 oncogene. In this study, we first extracted the available and annotated Mediterranean strains of HPV16 E6 gene sequences (n = 155) from the NCBI nucleotide database. These sequences were aligned, edited, and used for the downstream phylogenetic analyses. Finally, a Bayesian Markov Chain Monte Carlo approach was applied to reconstruct the evolutionary history of HPV16 migration. Our results showed that the HPV circulating in Tunisia derived from a Croatian ancestor around the year 1987. This starting point spreads to most European countries to reach northern Africa through the Moroccan gateway in 2004.

Species limits and diversification of the Dendropsophus rubicundulus subgroup (Anura, Hylidae) in Neotropical savannas

Understanding the processes that generate and maintain biodiversity at and below the species level is a central goal of evolutionary biology. Here we explore the spatial and temporal drivers of diversification of the treefrog subgroup Dendropsophus rubicundulus, a subgroup of the D. microcephalus species group, over periods of pronounced geological and climatic changes in the Neotropical savannas that they inhabit. This subgroup currently comprises 11 recognized species distributed across the Brazilian and Bolivian savannas, but the taxonomy has been in a state of flux, necessitating reexamination. Using newly generated single nucleotide polymorphism (SNP) data from restriction-site associated DNA sequencing (RADseq) and mitochondrial 16S sequence data for ∼150 specimens, we inferred phylogenetic relationships, tested species limits using a model-based approach, and estimated divergence times to gain insights into the geographic and climatic events that affected the diversification of this subgroup. Our results recognized at least nine species: D. anataliasiasi, D. araguaya, D. cerradensis, D. elianeae, D. jimi, D. rubicundulus, D. tritaeniatus, D. rozenmani, and D. sanborni. Although we did not collect SNP data for the latter two species, they are likely distinct based on mitochondrial data. In addition, we found genetic structure within the widespread species D. rubicundulus, which comprises three allopatric lineages connected by gene flow upon secondary contact. We also found evidence of population structure and perhaps undescribed diversity in D. elianeae, which warrants further study. The D. rubicundulus subgroup is estimated to have originated in the Late Miocene (∼5.45 million years ago), with diversification continuing through the Pliocene and Early Pleistocene, followed by the most recent divergence of D. rubicundulus lineages in the Middle Pleistocene. The epeirogenic uplift followed by erosion and denudation of the central Brazilian plateau throughout the Pliocene and Pleistocene, in combination with the increasing frequency and amplitude of climatic fluctuations during the Pleistocene, was important for generating and structuring diversity at or below the species level in the D. rubicundulus subgroup.

Molecular and cytogenetic evidence for sibling species in the Chagas disease vector Triatoma maculata

Triatoma maculata (Hemiptera, Reduviidae, Triatominae) occurs across dry-to-semiarid ecoregions of northern South America, where it transmits Trypanosoma cruzi, causative agent of Chagas disease. Using 207 field-caught specimens from throughout the species' range, mitochondrial(mt) DNA sequence data, and cytogenetics, we investigated inter-population genetic diversity and the phylogenetic affinities of T. maculata. Mitochondrial DNA sequence analyses (cytb and nd4) disclosed a monophyletic T. maculata clade encompassing three distinct geographic groups: Roraima formation (Guiana shield), Orinoco basin, and Magdalena basin (trans-Andean). Between-group cytb distances (11.0-12.8%) were larger than the ~7.5% expected for sister Triatoma species; the most recent common ancestor of these T. maculata groups may date back to the late Miocene. C-heterochromatin distribution and the sex-chromosome location of 45S ribosomal DNA clusters both distinguished Roraima bugs from Orinoco and Magdalena specimens. Cytb genealogies reinforced that T. maculata is not sister to Triatoma pseudomaculata and probably represents an early (middle-late Miocene) offshoot of the 'South American Triatomini lineage'. In sum, we report extensive genetic diversity and deep phylogeographic structuring in T. maculata, suggesting that it may consist of a complex of at least three sibling taxa. These findings have implications for the systematics, population biology, and perhaps medical relevance of T. maculata sensu lato.

Ecogeographic Drivers of the Spatial Spread of Highly Pathogenic Avian Influenza Outbreaks in Europe and the United States, 2016-Early 2022

H5Nx highly pathogenic avian influenza (HPAI) viruses of clade 2.3.4.4 have caused outbreaks in Europe among wild and domestic birds since 2016 and were introduced to North America via wild migratory birds in December 2021. We examined the spatiotemporal extent of HPAI viruses across continents and characterized ecological and environmental predictors of virus spread between geographic regions by constructing a Bayesian phylodynamic generalized linear model (phylodynamic-GLM). The findings demonstrate localized epidemics of H5Nx throughout Europe in the first several years of the epizootic, followed by a singular branching point where H5N1 viruses were introduced to North America, likely via stopover locations throughout the North Atlantic. Once in the United States (US), H5Nx viruses spread at a greater rate between US-based regions as compared to prior spread in Europe. We established that geographic proximity is a predictor of virus spread between regions, implying that intercontinental transport across the Atlantic Ocean is relatively rare. An increase in mean ambient temperature over time was predictive of reduced H5Nx virus spread, which may reflect the effect of climate change on declines in host species abundance, decreased persistence of the virus in the environment, or changes in migratory patterns due to ecological alterations. Our data provide new knowledge about the spread and directionality of H5Nx virus dispersal in Europe and the US during an actively evolving intercontinental outbreak, including predictors of virus movement between regions, which will contribute to surveillance and mitigation strategies as the outbreak unfolds, and in future instances of uncontained avian spread of HPAI viruses.

Climate heterogeneity shapes phylogeographic pattern of Hippophae gyantsensis (Elaeagnaceae) in the east Himalaya-Hengduan Mountains

The interaction of recent orographic uplift and climate heterogeneity acted as a key role in the East Himalaya-Hengduan Mountains (EHHM) has been reported in many studies. However, how exactly the interaction promotes clade diversification remains poorly understood. In this study, we both used the chloroplast trnT-trnF region and 11 nuclear microsatellite loci to investigate the phylogeographic structure and population dynamics of Hippophae gyantsensis and estimate what role geological barriers or ecological factors play in the spatial genetic structure. The results showed that this species had a strong east-west phylogeographic structure, with several mixed populations identified from microsatellite data in central location. The intraspecies divergence time was estimated to be about 3.59 Ma, corresponding well with the recent uplift of the Tibetan Plateau. Between the two lineages, there was significant climatic differentiation without geographic barriers. High consistency between lineage divergence, climatic heterogeneity, and Qingzang Movement demonstrated that climatic heterogeneity but not geographic isolation drives the divergence of H. gyantsensis, and the recent regional uplift of the QTP, as the Himalayas, creates heterogeneous climates by affecting the flow of the Indian monsoon. The east group of H. gyantsensis experienced population expansion c. 0.12 Ma, closely associated with the last interglacial interval. Subsequently, a genetic admixture event between east and west groups happened at 26.90 ka, a period corresponding to the warm inter-glaciation again. These findings highlight the importance of the Quaternary climatic fluctuations in the recent evolutionary history of H. gyantsensis. Our study will improve the understanding of the history and mechanisms of biodiversity accumulation in the EHHM region.

Novel Avian Influenza Virus (H5N1) Clade 2.3.4.4b Reassortants in Migratory Birds, China

Two novel reassortant highly pathogenic avian influenza viruses (H5N1) clade 2.3.4.4b.2 were identified in dead migratory birds in China in November 2021. The viruses probably evolved among wild birds through different flyways connecting Europe and Asia. Their low antigenic reaction to vaccine antiserum indicates high risks to poultry and to public health.

Phylogeographic analysis of Siraitia grosvenorii in subtropical China provides insights into the origin of cultivated monk fruit and conservation of genetic resources

Siraitia grosvenorii, an economically important plant species with high medicinal value, is endemic to subtropical China. To determine the population structure and origin of cultivated S. grosvenorii, we examined the variation in three chloroplast DNA regions (trnR-atpA, trnH-psbA, trnL-trnF) and two orthologous nuclear genes (CHS and EDL2) of S. grosvenorii in 130 wild individuals (selected from 13 wild populations across its natural distribution range) and 21 cultivated individuals using a phylogeographic approach. The results showed three distinct chloroplast lineages, which were restricted to different mountain ranges, and strong plastid phylogeographic structure. Our findings suggest that S. grosvenorii likely experienced ancient range expansion and survived in multiple refuges in subtropical China during glacial periods, resulting in population fragmentation in different mountainous areas. Our results also demonstrated that wild populations in Guilin (Guangxi, China) share the same gene pool as cultivated S. grosvenorii, suggesting that current cultivars were collected directly from local wild resources, consistent with the principles of "nearby domestication." The results of this study provide insights into improving the efficiency of S. grosvenorii breeding using a genetic approach and outline measures for the conservation of its genetic resources.

Addressing the complex phylogenetic relationship of the Gempylidae fishes using mitogenome data

The Gempylidae (snake mackerels) family, belonging to the order Perciformes, consists of about 24 species described in 16 genera primarily distributed in tropical, subtropical, and temperate seas worldwide. Despite substantial research on this family utilizing morphological and molecular approaches, taxonomy categorization in this group has remained puzzling for decades prompting the need for further investigation into the underlying evolutionary history among the gempylids using molecular tools. In this study, we assembled eight complete novel mitochondrial genomes for five Gempylidae species (Neoepinnula minetomai, Neoepinnula orientalis, Rexea antefurcata, Rexea prometheoides, and Thyrsites atun) using Ion Torrent sequencing to supplement publicly available mitogenome data for gempylids. Using Bayesian inference and maximum-likelihood tree search methods, we investigated the evolutionary relationships of 17 Gempylidae species using mitogenome data. In addition, we estimated divergence times for extant gempylids. We identified two major clades that formed approximately 48.05 (35.89-52.04) million years ago: Gempylidae 1 (Thyrsites atun, Promethichthys prometheus, Nealotus tripes, Diplospinus multistriatus, Paradiplospinus antarcticus, Rexea antefurcata, Rexea nakamurai, Rexea prometheoides, Rexea solandri, Thyrsitoides marleyi, Gempylus serpens, and Nesiarchus nasutus) and Gempylidae 2 (Lepidocybium flavobrunneum, Ruvettus pretiosus, Neoepinnula minetomai, Neoepinnula orientalis, and Epinnula magistralis). The present study demonstrated the superior performance of complete mitogenome data compared with individual genes in phylogenetic reconstruction. By including T. atun individuals from different regions, we demonstrated the potential for the application of mitogenomes in species phylogeography.

Population structure and phylogeography of three closely related tree peonies

Paeonia decomposita, Paeonia rotundiloba, and Paeonia rockii are three closely related species of Sect. Moutan is distributed in the montane area of the Eastern Hengduan Mountain region. Understanding the population history of these three tree peony species could contribute to unraveling the evolutionary patterns of undergrowth species in this hotspot area. We used one nuclear DNA marker (internal transcribed spacer region, ITS) and two chloroplast DNA markers (matK, ycf1) to reconstruct the phylogeographic pattern of the populations. In total, 228 individuals from 17 populations of the three species were analyzed in this study. Three nuclear clades (Clade I - Clade III) and four maternal clades (Clade A - Clade D) were reconstructed. Molecular dating suggested that young lineages diverged during the late Pliocene and early Pleistocene, younger than the uplift of the Hengduan Mountains but older than the last glacial maximum (LGM). Significant population and phylogeographic structures were detected at both markers. Furthermore, the populations of these tree peonies were overall at equilibrium during the climatic oscillations of the Pleistocene. The simulated palaeoranges of the three species during the LGM period mostly overlapped, which could have led to cross-breeding events. We propose an evolutionary scenario in which mountain orogenesis around the Hengduan Mountain area triggered parapatric isolation between maternal lineages of tree peonies. Subsequent climatic fluctuations drove migration and range recontact of these populations along the valleys. This detailed evolutionary history provides new insights into the phylogeographic pattern of species from mountain-valley systems.

Ocean crossers: a tale of disjunctions and speciation in the dwarf-fruticose Lichina (lichenized Ascomycota)

Lichens thrive in rocky coastal areas in temperate and cold regions of both hemispheres. Species of the genus Lichina, which form characteristic black fruiting thalli associated with cyanobacteria, often create distinguishable bands in the intertidal and supralittoral zones. The present study uses a comprehensive specimen dataset and four gene loci to (1) delineate and discuss species boundaries in this genus, (2) assess evolutionary relationships among species, and (3) infer the most likely causes of their current geographic distribution in the Northern and Southern hemispheres. A dated phylogeny describes the time frame in which extant disjunctions of species and populations were established. The results showed that the genus is integrated by four species, with Lichina pygmaea, L. confinis and the newly described L. canariensis from rocky seashores in the Canary Islands, occurring in the Northern Hemisphere, whereas L. intermedia is restricted to the Southern Hemisphere. Lichina intermedia hosted a much higher intraspecific genetic diversity than the other species, with subclades interpreted as species-level lineages by the different species delimitation approaches. However, a conservative taxonomic approach was adopted. This species showed a striking disjunct distribution between Australasia and southern South America. The timing for the observed interspecific and intraspecific divergences and population disjunctions postdated continental plate movements, suggesting that long-distance dispersal across body waters in the two hemispheres played a major role in shaping the current species distributions. Such ocean crossings were, as in L. canariensis, followed by speciation. New substitution rates for the nrITS of the genus Lichina were inferred using a tree spanning the major Ascomycota lineages calibrated using fossils. In conclusion, this work lays the foundation for a better understanding of the evolution through time and space of maritime lichens.

Mapping the Dynamics of Contemporary PRRSV-2 Evolution and Its Emergence and Spreading Hotspots in the U.S. Using Phylogeography

The repeated emergence of new genetic variants of PRRSV-2, the virus that causes porcine reproductive and respiratory syndrome (PRRS), reflects its rapid evolution and the failure of previous control efforts. Understanding spatiotemporal heterogeneity in variant emergence and spread is critical for future outbreak prevention. Here, we investigate how the pace of evolution varies across time and space, identify the origins of sub-lineage emergence, and map the patterns of the inter-regional spread of PRRSV-2 Lineage 1 (L1)-the current dominant lineage in the U.S. We performed comparative phylogeographic analyses on subsets of 19,395 viral ORF5 sequences collected across the U.S. and Canada between 1991 and 2021. The discrete trait analysis of multiple spatiotemporally stratified sampled sets (n = 500 each) was used to infer the ancestral geographic region and dispersion of each sub-lineage. The robustness of the results was compared to that of other modeling methods and subsampling strategies. Generally, the spatial spread and population dynamics varied across sub-lineages, time, and space. The Upper Midwest was a main spreading hotspot for multiple sub-lineages, e.g., L1C and L1F, though one of the most recent emergence events (L1A(2)) spread outwards from the east. An understanding of historical patterns of emergence and spread can be used to strategize disease control and the containment of emerging variants.

Unraveling the epidemiology of Mycobacterium bovis using whole-genome sequencing combined with environmental and demographic data

When studying the dynamics of a pathogen in a host population, one crucial question is whether it transitioned from an epidemic (i.e., the pathogen population and the number of infected hosts are increasing) to an endemic stable state (i.e., the pathogen population reached an equilibrium). For slow-growing and slow-evolving clonal pathogens such as Mycobacterium bovis, the causative agent of bovine (or animal) and zoonotic tuberculosis, it can be challenging to discriminate between these two states. This is a result of the combination of suboptimal detection tests so that the actual extent of the pathogen prevalence is often unknown, as well as of the low genetic diversity, which can hide the temporal signal provided by the accumulation of mutations in the bacterial DNA. In recent years, the increased availability, efficiency, and reliability of genomic reading techniques, such as whole-genome sequencing (WGS), have significantly increased the amount of information we can use to study infectious diseases, and therefore, it has improved the precision of epidemiological inferences for pathogens such as M. bovis. In this study, we use WGS to gain insights into the epidemiology of M. bovis in Cameroon, a developing country where the pathogen has been reported for decades. A total of 91 high-quality sequences were obtained from tissue samples collected in four abattoirs, 64 of which were with complete metadata. We combined these with environmental, demographic, ecological, and cattle movement data to generate inferences using phylodynamic models. Our findings suggest M. bovis in Cameroon is slowly expanding its epidemiological range over time; therefore, endemic stability is unlikely. This suggests that animal movement plays an important role in transmission. The simultaneous prevalence of M. bovis in co-located cattle and humans highlights the risk of such transmission being zoonotic. Therefore, using genomic tools as part of surveillance would vastly improve our understanding of disease ecology and control strategies.

Spatial genetic patterns and distribution dynamics of Begonia grandis (Begoniaceae), a widespread herbaceous species in China

Introduction

Begonia L., one of the 10 largest plant genera, contains over 2,100 species, most of which have a very limited distribution range. Understanding the spatial genetic structure and distribution dynamics of a widespread species in this genus will contribute to clarifying the mechanism responsible for Begonia speciation.

Methods

In this study, we used three chloroplast DNA markers (ndhF-rpl32, atpI-atpH, and ndhA intron), coupled with species distribution modeling (SDM), to investigate the population genetic structure and distribution dynamics of Begonia grandis Dryand., the species of Begonia with the widest distribution in China.

Results

Thirty-five haplotypes from 44 populations clustered into two groups, and haplotype divergence began in the Pleistocene (1.75 Mya). High genetic diversity (H _d = 0.894, H _T = 0.910), strong genetic differentiation (F _ST = 0.835), and significant phylogeographical structure (G _ST/N _ST = 0.848/0.917, P < 0.05) were observed. The distribution range of B. grandis migrated northwards after the last glacial maximum, but its core distribution area remained stable.

Discussion

Combined, the observed spatial genetic patterns and SDM results identified the Yunnan-Guizhou Plateau, the Three Gorges region, and the Daba Mountains as potential refugia of B. grandis. BEAST-derived chronogram and haplotype network analysis do not support the Flora Reipublicae Popularis Sinicae and Flora of China for subspecies classification based on morphological characteristics. Our results support the hypothesis that population-level allopatric differentiation may be an important speciation process for the Begonia genus and a key contributor to its rich diversity.

Genetic Diversity and Dispersal of DENGUE Virus among Three Main Island Groups of the Philippines during 2015-2017

Dengue has been one of the major public health concerns in the Philippines for more than a century. The annual dengue case burden has been increasing in recent years, exceeding 200,000 in 2015 and 2019. However, there is limited information on the molecular epidemiology of dengue in the Philippines. We, therefore, conducted a study to understand the genetic composition and dispersal of DENV in the Philippines from 2015 to 2017 under UNITEDengue. Our analyses included 377 envelope (E) gene sequences of all 4 serotypes obtained from infections in 3 main island groups (Luzon, Visayas, and Mindanao) of the Philippines. The findings showed that the overall diversity of DENV was generally low. DENV-1 was relatively more diverse than the other serotypes. Virus dispersal was evident among the three main island groups, but each island group demonstrated a distinct genotype composition. These observations suggested that the intensity of virus dispersal was not substantive enough to maintain a uniform heterogeneity among island groups so that each island group behaved as an independent epidemiological unit. The analyses suggested Luzon as one of the major sources of DENV emergence and CAR, Calabarzon, and CARAGA as important hubs of virus dispersal in the Philippines. Our findings highlight the importance of virus surveillance and molecular epidemiological analyses to gain deep insights into virus diversity, lineage dominance, and dispersal patterns that could assist in understanding the epidemiology and transmission risk of dengue in endemic regions.

Genomic, spatial and morphometric data for discrimination of four species in the Mediterranean Tamus clade of yams (Dioscorea, Dioscoreaceae)

BACKGROUND AND AIMS

Among the numerous pantropical species of the yam genus, Dioscorea, only a small group occurs in the Mediterranean basin, including two narrow Pyrenean endemics (Borderea clade) and two Mediterranean-wide species (D. communis and D. orientalis, Tamus clade). However, several currently unrecognized species and infraspecific taxa have been described in the Tamus clade due to significant morphological variation associated with D. communis. Our overarching aim was to investigate taxon delimitation in the Tamus clade using an integrative approach combining phylogenomic, spatial and morphological data.

METHODS

We analysed 76 herbarium samples using Hyb-Seq genomic capture to sequence 260 low-copy nuclear genes and plastomes, together with morphometric and environmental modelling approaches.

KEY RESULTS

Phylogenomic reconstructions confirmed that the two previously accepted species of the Tamus clade, D. communis and D. orientalis, are monophyletic and form sister clades. Three subclades showing distinctive geographic patterns were identified within D. communis. These subclades were also identifiable from morphometric and climatic data, and introgression patterns were inferred between subclades in the eastern part of the distribution of D. communis.

CONCLUSIONS

We propose a taxonomy that maintains D. orientalis, endemic to the eastern Mediterranean region, and splits D. communis sensu lato into three species: D. edulis, endemic to Macaronesia (Canary Islands and Madeira); D. cretica, endemic to the eastern Mediterranean region; and D. communis sensu stricto, widespread across western and central Europe. Introgression inferred between D. communis s.s. and D. cretica is likely to be explained by their relatively recent speciation at the end of the Miocene, disjunct isolation in eastern and western Mediterranean glacial refugia and a subsequent westward recolonization of D. communis s.s. Our study shows that the use of integrated genomic, spatial and morphological approaches allows a more robust definition of species boundaries and the identification of species that previous systematic studies failed to uncover.

Recent Advances in Genetics and Genomics of Snub-Nosed Monkeys (Rhinopithecus) and Their Implications for Phylogeny, Conservation, and Adaptation

The snub-nosed monkey genus Rhinopithecus (Colobinae) comprises five species (Rhinopithecus roxellana, Rhinopithecus brelichi, Rhinopithecus bieti, Rhinopithecus strykeri, and Rhinopithecus avunculus). They are range-restricted species occurring only in small areas in China, Vietnam, and Myanmar. All extant species are listed as endangered or critically endangered by the International Union for Conservation of Nature (IUCN) Red List, all with decreasing populations. With the development of molecular genetics and the improvement and cost reduction in whole-genome sequencing, knowledge about evolutionary processes has improved largely in recent years. Here, we review recent major advances in snub-nosed monkey genetics and genomics and their impact on our understanding of the phylogeny, phylogeography, population genetic structure, landscape genetics, demographic history, and molecular mechanisms of adaptation to folivory and high altitudes in this primate genus. We further discuss future directions in this research field, in particular how genomic information can contribute to the conservation of snub-nosed monkeys.

Population Genetic Structure and Diversity of Cryptic Species of the Plant Genus Macrocarpaea (Gentianaceae) from the Tropical Andes

The Pleistocene climatic oscillations (PCO) that provoked several cycles of glacial-interglacial periods are thought to have profoundly affected species distribution, richness and diversity around the world. While the effect of the PCO on population dynamics at temperate latitudes is well known, considerable questions remain about its impact on the biodiversity of neotropical mountains. Here, we use amplified fragment length polymorphism molecular markers (AFLPs) to investigate the phylogeography and genetic structure of 13 plant species belonging to the gentian genus Macrocarpaea (Gentianaceae) in the tropical Andes. These woody herbs, shrubs or small trees show complex and potentially reticulated relationships, including cryptic species. We show that populations of M. xerantifulva in the dry system of the Rio Marañón in northern Peru have lower levels of genetic diversity compared to other sampled species. We suggest that this is due to a recent demographic bottleneck resulting from the contraction of the montane wet forests into refugia because of the expansion of the dry system into the valley during the glacial cycles of the PCO. This may imply that the ecosystems of different valleys of the Andes might have responded differently to the PCO.

Genome-wide phylogeography reveals cryptic speciation in the circumglobal planktonic calcifier Limacina bulimoides

Little is known about when and how planktonic species arise and persist in the open ocean without apparent dispersal barriers. Pteropods are planktonic snails with thin shells susceptible to dissolution that are used as bio-indicators of ocean acidification. However, distinct evolutionary units respond to acidification differently and defining species boundaries is therefore crucial for predicting the impact of changing ocean conditions. In this global population genomic study of the shelled pteropod Limacina bulimoides, we combined genetic (759,000 single nucleotide polymorphisms) and morphometric data from 161 individuals, revealing three major genetic lineages (F_ST = 0.29 to 0.41): an 'Atlantic lineage' sampled across the Atlantic, an 'Indo-Pacific lineage' sampled in the North Pacific and Indian Ocean, and a 'Pacific lineage' sampled in the North and South Pacific. A time-calibrated phylogeny suggests that the lineages diverged about one million years ago, with estimated effective population size remaining high (~10 million) throughout Pleistocene glacial cycles. We do not observe any signatures of recent hybridisation, even in areas of sympatry in the North Pacific. While the lineages are reproductively isolated, they are morphologically cryptic, with overlapping shell shape and shell colour distributions. Despite showing that the circumglobal L. bulimoides consists of multiple species with smaller ranges than initially thought, we found that these pteropods still possess high levels of genetic variability. Our study adds to the growing evidence that speciation is often overlooked in the open ocean, and suggests the presence of distinct biological species within many other currently defined circumglobal planktonic species.