BEAGLE: an application programming interface and high-performance computing library for statistical phylogenetics

Whole genome sequencing unravels cryptic circulation of divergent dengue virus lineages in the rainforest region of Nigeria

Dengue is often misclassified and underreported in Africa due to inaccurate differential diagnoses of nonspecific febrile illnesses such as malaria, sparsity of diagnostic testing and poor clinical and genomic surveillance. There are limited reports on the seroprevalence and genetic diversity of dengue virus (DENV) in humans and vectors in Nigeria. In this study, we investigated the epidemiology and genetic diversity of dengue in the rainforest region of Nigeria. We screened 515 febrile patients who tested negative for malaria and typhoid fever in three hospitals in Oyo and Ekiti States in southern Nigeria with a combination of anti-dengue IgG/IgM/NS1 rapid test kits and metagenomic sequencing. We found that approximately 28% of screened patients had previous DENV exposure, with the highest prevalence in persons over sixty. Approximately 8% of the patients showed evidence of recent or current infection, and 2.7% had acute infection. Following sequencing of sixty samples, we assembled twenty DENV-1 genomes (3 complete and 17 partial). We found that all assembled genomes belonged to DENV-1 genotype III. Our phylogenetic analyses showed evidence of prolonged cryptic circulation of divergent DENV lineages in Oyo state. We were unable to resolve the source of DENV in Nigeria owing to limited sequencing data from the region. However, our sequences clustered closely with sequences in Tanzania and sequences reported in Chinese with travel history to Tanzania in 2019. This may reflect the wider unsampled bidirectional transmission of DENV-1 in Africa, which strongly emphasizes the importance of genomic surveillance in monitoring ongoing DENV transmission in Africa.

Human Y chromosome haplogroup L1-M22 traces Neolithic expansion in West Asia and supports the Elamite and Dravidian connection

West and South Asian populations profoundly influenced Eurasian genetic and cultural diversity. We investigate the genetic history of the Y chromosome haplogroup L1-M22, which, while prevalent in these regions, lacks in-depth study. Robust Bayesian analyses of 165 high-coverage Y chromosomes favor a West Asian origin for L1-M22 ∼20.6 thousand years ago (kya). Moreover, this haplogroup parallels the genome-wide genetic ancestry of hunter-gatherers from the Iranian Plateau and the Caucasus. We characterized two L1-M22 harboring population groups during the Early Holocene. One expanded with the West Asian Neolithic transition. The other moved to South Asia ∼8-6 kya but showed no expansion. This group likely participated in the spread of Dravidian languages. These South Asian L1-M22 lineages expanded ∼4-3 kya, coinciding with the Steppe ancestry introduction. Our findings advance the current understanding of Eurasian historical dynamics, emphasizing L1-M22's West Asian origin, associated population movements, and possible linguistic impacts.

Phylodynamics unveils invading and diffusing patterns of dengue virus serotype-1 in Guangdong, China from 1990 to 2019 under a global genotyping framework

BACKGROUND

The strong invasiveness and rapid expansion of dengue virus (DENV) pose a great challenge to global public health. However, dengue epidemic patterns and mechanisms at a genetic scale, particularly in term of cross-border transmissions, remain poorly understood. Importation is considered as the primary driver of dengue outbreaks in China, and since 1990 a frequent occurrence of large outbreaks has been triggered by the imported cases and subsequently spread to the western and northern parts of China. Therefore, this study aims to systematically reveal the invasion and diffusion patterns of DENV-1 in Guangdong, China from 1990 to 2019.

METHODS

These analyses were performed on 179 newly assembled genomes from indigenous dengue cases in Guangdong, China and 5152 E gene complete sequences recorded in Chinese mainland. The genetic population structure and epidemic patterns of DENV-1 circulating in Chinese mainland were characterized by phylogenetics, phylogeography, phylodynamics based on DENV-1 E-gene-based globally unified genotyping framework.

RESULTS

Multiple serotypes of DENV were co-circulating in Chinese mainland, particularly in Guangdong and Yunnan provinces. A total of 189 transmission clusters in 38 clades belonging to 22 subgenotypes of genotype I, IV and V of DENV-1 were identified, with 7 Clades of Concern (COCs) responsible for the large outbreaks since 1990. The epidemic periodicity was inferred from the data to be approximately 3 years. Dengue transmission events mainly occurred from Great Mekong Subregion-China (GMS-China), Southeast Asia (SEA), South Asia Subcontinent (SASC), and Oceania (OCE) to coastal and land border cities respectively in southeastern and southwestern China. Specially, Guangzhou was found to be the most dominant receipting hub, where DENV-1 diffused to other cities within the province and even other parts of the country. Genome phylogeny combined with epidemiological investigation demonstrated a clear local consecutive transmission process of a 5C1 transmission cluster (5C1-CN4) of DENV-1 in Guangzhou from 2013 to 2015, while the two provinces of Guangdong and Yunnan played key roles in ongoing transition of dengue epidemic patterns. In contextualizing within Invasion Biology theories, we have proposed a derived three-stage model encompassing the stages of invasion, colonization, and dissemination, which is supposed to enhance our understanding of dengue spreading patterns.

CONCLUSIONS

This study demonstrates the invasion and diffusion process of DENV-1 in Chinese mainland within a global genotyping framework, characterizing the genetic diversities of viral populations, multiple sources of importation, and periodic dynamics of the epidemic. These findings highlight the potential ongoing transition trends from epidemic to endemic status offering a valuable insight into early warning, prevention and control of rapid spreading of dengue both in China and worldwide.

Spatiotemporal transmission of SARS-CoV-2 lineages during 2020-2021 in Pernambuco-Brazil

SARS-CoV-2 virus emerged as a new threat to humans and spread around the world, leaving a large death toll. As of January 2023, Brazil is among the countries with the highest number of registered deaths. Nonpharmacological and pharmacological interventions have been heterogeneously implemented in the country, which, associated with large socioeconomic differences between the country regions, has led to distinct virus spread dynamics. Here, we investigate the spatiotemporal dispersion of SARS-CoV-2 lineages in the Pernambuco state (Northeast Brazil) throughout the distinct epidemiological scenarios that unfolded in the first 2 years of the pandemic. We generated a total of 1,389 new SARS-CoV-2 genomes from June 2020 to August 2021. This sampling captured the arrival, communitary transmission, and the circulation of the B1.1, B.1.1.28, and B.1.1.33 lineages; the emergence of the former variant of interest P.2; and the emergence and fast replacement of all previous variants by the more transmissible variant of concern P.1 (Gamma). Based on the incidence and lineage spread pattern, we observed an East-to-West to inner state pattern of transmission, which is in agreement with the transmission of more populous metropolitan areas to medium- and small-size country-side cities in the state. Such transmission patterns may be partially explained by the main routes of traffic across municipalities in the state. Our results highlight that the fine-grained intrastate analysis of lineages and incidence spread can provide actionable insights for planning future nonpharmacological intervention for air-borne transmissible human pathogens.IMPORTANCEDuring the COVID-19 pandemic, Brazil was one of the most affected countries, mainly due its continental-size, socioeconomic differences among regions, and heterogeneous implementation of intervention methods. In order to investigate SARS-CoV-2 dynamics in the state of Pernambuco, we conducted a spatiotemporal dispersion study, covering the period from June 2020 to August 2021, to comprehend the dynamics of viral transmission during the first 2 years of the pandemic. Throughout this study, we were able to track three significant epidemiological waves of transmission caused by B1.1, B.1.1.28, B.1.1.33, P.2, and P.1 lineages. These analyses provided valuable insights into the evolution of the epidemiological landscape, contributing to a deeper understanding of the dynamics of virus transmission during the early years of the pandemic in the state of Pernambuco.

Extreme genome scrambling in marine planktonic Oikopleura dioica cryptic species

Genome structural variations within species are rare. How selective constraints preserve gene order and chromosome structure is a central question in evolutionary biology that remains unsolved. Our sequencing of several genomes of the appendicularian tunicate Oikopleura dioica around the globe reveals extreme genome scrambling caused by thousands of chromosomal rearrangements, although showing no obvious morphological differences between these animals. The breakpoint accumulation rate is an order of magnitude higher than in ascidian tunicates, nematodes, Drosophila, or mammals. Chromosome arms and sex-specific regions appear to be the primary unit of macrosynteny conservation. At the microsyntenic level, scrambling did not preserve operon structures, suggesting an absence of selective pressure to maintain them. The uncoupling of the genome scrambling with morphological conservation in O. dioica suggests the presence of previously unnoticed cryptic species and provides a new biological system that challenges our previous vision of speciation in which similar animals always share similar genome structures.

Pathophysiology of chikungunya virus infection associated with fatal outcomes

Chikungunya virus (CHIKV) is a mosquito-borne alphavirus that causes acute, subacute, and chronic human arthritogenic diseases and, in rare instances, can lead to neurological complications and death. Here, we combined epidemiological, virological, histopathological, cytokine, molecular dynamics, metabolomic, proteomic, and genomic analyses to investigate viral and host factors that contribute to chikungunya-associated (CHIK) death. Our results indicate that CHIK deaths are associated with multi-organ infection, central nervous system damage, and elevated serum levels of pro-inflammatory cytokines and chemokines compared with survivors. The histopathologic, metabolite, and proteomic signatures of CHIK deaths reveal hemodynamic disorders and dysregulated immune responses. The CHIKV East-Central-South-African lineage infecting our study population causes both fatal and survival cases. Additionally, CHIKV infection impairs the integrity of the blood-brain barrier, as evidenced by an increase in permeability and altered tight junction protein expression. Overall, our findings improve the understanding of CHIK pathophysiology and the causes of fatal infections.

Complete Mitochondrial Genome of the Eggplant Fruit and Shoot Borer, Leucinodes orbonalis Guenée (Lepidoptera: Crambidae), and Comparison with Other Pyraloid Moths

The eggplant fruit and shoot borer (EFSB) (Leucinodes orbonalis Guenée) is a devastating lepidopteran pest of eggplant (Solanum melongena L.) in the Philippines. Management of an insect pest like the EFSB requires an understanding of its biology, evolution, and adaptations. Genomic resources provide a starting point for understanding EFSB biology, as the resources can be used for phylogenetics and population structure studies. To date, genomic resources are scarce for EFSB; thus, this study generated its complete mitochondrial genome (mitogenome). The circular mitogenome is 15,244 bp-long. It contains 37 genes, namely 13 protein-coding, 22 tRNA, and 2 rRNA genes, and has conserved noncoding regions, motifs, and gene syntenies characteristic of lepidopteran mitogenomes. Some protein-coding genes start and end with non-canonical codons. The tRNA genes exhibit a conserved cloverleaf structure, with the exception in trnS1. Partitioned phylogenetic analysis using 72 pyraloids generated highly supported maximum likelihood and Bayesian inference trees revealing expected basal splits between Crambidae and Pyralidae, and Spilomelinae and Pyraustinae. Spilomelinae was recovered to be paraphyletic, with the EFSB robustly placed before the split of Spilomelinae and Pyraustinae. Overall, the EFSB mitogenome resource will be useful for delineations within Spilomelinae and population structure analysis.

Evolutionary trajectory of diverse SARS-CoV-2 variants at the beginning of COVID-19 outbreak

Despite extensive scientific efforts directed toward the evolutionary trajectory of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in humans at the beginning of the COVID-19 epidemic, it remains unclear how the virus jumped into and evolved in humans so far. Herein, we recruited almost all adult coronavirus disease 2019 (COVID-19) cases appeared locally or imported from abroad during the first 8 months of the outbreak in Shanghai. From these patients, SARS-CoV-2 genomes occupying the important phylogenetic positions in the virus phylogeny were recovered. Phylogenetic and mutational landscape analyses of viral genomes recovered here and those collected in and outside of China revealed that all known SARS-CoV-2 variants exhibited the evolutionary continuity despite the co-circulation of multiple lineages during the early period of the epidemic. Various mutations have driven the rapid SARS-CoV-2 diversification, and some of them favor its better adaptation and circulation in humans, which may have determined the waxing and waning of various lineages.

Synchrony of Bird Migration with Global Dispersal of Avian Influenza Reveals Exposed Bird Orders

Highly pathogenic avian influenza virus (HPAIV) A H5, particularly clade 2.3.4.4, has caused worldwide outbreaks in domestic poultry, occasional spillover to humans, and increasing deaths of diverse species of wild birds since 2014. Wild bird migration is currently acknowledged as an important ecological process contributing to the global dispersal of HPAIV H5. However, this mechanism has not been quantified using bird movement data from different species, and the timing and location of exposure of different species is unclear. We sought to explore these questions through phylodynamic analyses based on empirical data of bird movement tracking and virus genome sequences of clade 2.3.4.4 and 2.3.2.1. First, we demonstrate that seasonal bird migration can explain salient features of the global dispersal of clade 2.3.4.4. Second, we detect synchrony between the seasonality of bird annual cycle phases and virus lineage movements. We reveal the differing exposed bird orders at geographical origins and destinations of HPAIV H5 clade 2.3.4.4 lineage movements, including relatively under-discussed orders. Our study provides a phylodynamic framework that links the bird movement ecology and genomic epidemiology of avian influenza; it highlights the importance of integrating bird behavior and life history in avian influenza studies.

Whole-genome resequencing of major populations revealed domestication-related genes in yaks

BACKGROUND

The yak is a symbol of the Qinghai-Tibet Plateau and provides important basic resources for human life on the plateau. Domestic yaks have been subjected to strong artificial selection and environmental pressures over the long-term. Understanding the molecular mechanisms of phenotypic differences in yak populations can reveal key functional genes involved in the domestication process and improve genetic breeding.

MATERIAL AND METHOD

Here, we re-sequenced 80 yaks (Maiwa, Yushu, and Huanhu populations) to identify single-nucleotide polymorphisms (SNPs) as genetic variants. After filtering and quality control, remaining SNPs were kept to identify the genome-wide regions of selective sweeps associated with domestic traits. The four methods (π, XPEHH, iHS, and XP-nSL) were used to detect the population genetic separation.

RESULTS

By comparing the differences in the population stratification, linkage disequilibrium decay rate, and characteristic selective sweep signals, we identified 203 putative selective regions of domestic traits, 45 of which were mapped to 27 known genes. They were clustered into 4 major GO biological process terms. All known genes were associated with seven major domestication traits, such as dwarfism (ANKRD28), milk (HECW1, HECW2, and OSBPL2), meat (SPATA5 and GRHL2), fertility (BTBD11 and ARFIP1), adaptation (NCKAP5, ANTXR1, LAMA5, OSBPL2, AOC2, and RYR2), growth (GRHL2, GRID2, SMARCAL1, and EPHB2), and the immune system (INPP5D and ADCYAP1R1).

CONCLUSIONS

We provided there is an obvious genetic different among domestic progress in these three yak populations. Our findings improve the understanding of the major genetic switches and domestic processes among yak populations.

A new perspective on the taxonomy and systematics of Arvicolinae (Gray, 1821) and a new time-calibrated phylogeny for the clade

Background

Arvicoline rodents are one of the most speciose and rapidly evolving mammalian lineages. Fossil arvicolines are also among the most common vertebrate fossils found in sites of Pliocene and Pleistocene age in Eurasia and North America. However, there is no taxonomically robust, well-supported, time-calibrated phylogeny for the group.

Methods

Here we present well-supported hypotheses of arvicoline rodent systematics using maximum likelihood and Bayesian inference of DNA sequences of two mitochondrial genes and three nuclear genes representing 146 (82% coverage) species and 100% of currently recognized arvicoline genera. We elucidate well-supported major clades, reviewed the relationships and taxonomy of many species and genera, and critically compared our resulting molecular phylogenetic hypotheses to previously published hypotheses. We also used five fossil calibrations to generate a time-calibrated phylogeny of Arvicolinae that permitted some reconciliation between paleontological and neontological data.

Results

Our results are largely congruent with previous molecular phylogenies, but we increased the support in many regions of the arvicoline tree that were previously poorly-sampled. Our sampling resulted in a better understanding of relationships within Clethrionomyini, the early-diverging position and close relationship of true lemmings (Lemmus and Myopus) and bog lemmings (Synaptomys), and provided support for recent taxonomic changes within Microtini. Our results indicate an origin of ∼6.4 Ma for crown arvicoline rodents. These results have major implications (e.g., diversification rates, paleobiogeography) for our confidence in the fossil record of arvicolines and their utility as biochronological tools in Eurasia and North America during the Quaternary.

Genomic epidemiology and evolutionary analysis during XBB.1.16-predominant periods of SARS-CoV-2 omicron variant in Bangkok, Thailand: December 2022-August 2023

The growing occurrence of novel recombinants, such as XBB.1.16, has emerged and become predominant, raising concerns about the impact of genomic recombination on the evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). This study investigated the molecular epidemiological trends and evolution of the Omicron XBB.1.16 epidemic in Bangkok between December 2022 and August 2023. Partial spike and complete genome sequencing of SARS-CoV-2 samples collected from collaborating hospitals were performed. The analysis of 491 partial spike sequences identified 15 distinct lineages, with XBB.1.16 dominating the lineages beginning in March 2023. Phylogenetic analysis revealed at least four distinct XBB.1.16 lineages, suggesting multiple independent introductions into Bangkok. The estimated emergence of XBB.1.16 occurred approximately in January 2022, with an evolutionary rate of 0.79 × 10-3 substitutions per site per year. Monitoring the genomic epidemiology and evolution of XBB.1.16 is vital for the early detection of new strains or emerging variants, which may guide vaccine design and the inclusion of new vaccine strains.

The phylogenetic reconstruction of the Neotropical cycad genus Ceratozamia (Zamiaceae) reveals disparate patterns of niche evolution

The cycad genus Ceratozamia comprises 40 species from Mexico, Guatemala, Belize, and Honduras, where cycads occur throughout climatically varied montane habitats. Ceratozamia has the potential to reveal the history and processes of species diversification across diverse Neotropical habitats in this region. However, the species relationships within Ceratozamia and the ecological trends during its evolution remain unclear. Here, we aimed to clarify the phylogenetic relationships, the timing of clade and species divergences, and the niche evolution throughout the phylogenetic history of Ceratozamia. Genome-wide DNA sequences were obtained with MIG-seq, and multiple data-filtering steps were used to optimize the dataset used to construct an ultrametric species tree. Divergence times among branches and ancestral niches were estimated. The niche variation among species was evaluated, summarized into two principal components, and their ancestral states were reconstructed to test whether niche shifts among branches can be explained by random processes, under a Brownian Motion model. Ceratozamia comprises three main clades, and most species relationships within the clades were resolved. Ceratozamia has diversified since the Oligocene, with major branching events occurring during the Miocene. This timing is consistent with fossil evidence, the timing estimated for other Neotropical plant groups, and the major geological events that shaped the topographic and climatic variation in Mexico. Patterns of niche evolution in the genus do not accord with the Brownian Motion model. Rather, non-random evolution with shifts towards more seasonal environments at high latitudes, or shifts towards humid or dry environments at low latitudes explain the diversification of Ceratozamia. We present a comprehensive phylogenetic reconstruction for Ceratozamia and identify for the first time the environmental factors involved in clade and species diversification within the genus. This study alleviates the controversies regarding the species relationships in the genus and provides the first evidence that latitude-associated environmental factors may influence processes of niche evolution in cycads.

Diversity and dissemination of viruses in pathogenic protozoa

Viruses are the most abundant biological entities on Earth and play a significant role in the evolution of many organisms and ecosystems. In pathogenic protozoa, the presence of viruses has been linked to an increased risk of treatment failure and severe clinical outcome. Here, we studied the molecular epidemiology of the zoonotic disease cutaneous leishmaniasis in Peru and Bolivia through a joint evolutionary analysis of Leishmania braziliensis and their dsRNA Leishmania virus 1. We show that parasite populations circulate in tropical rainforests and are associated with single viral lineages that appear in low prevalence. In contrast, groups of hybrid parasites are geographically and ecologically more dispersed and associated with an increased prevalence, diversity and spread of viruses. Our results suggest that parasite gene flow and hybridization increased the frequency of parasite-virus symbioses, a process that may change the epidemiology of leishmaniasis in the region.

Detection of genomic variations and selection signatures in Wagyu using whole-genome sequencing data

Wagyu is recognized for producing marbled beef with high nutritional value and flavor. Reportedly, Wagyu has been widely used to improve the meat quality of local breeds around the world. However, studies on the genetic mechanism of meat quality in Wagyu at the whole-genome level are rarely reported. Here, whole-genome sequencing data of 11 Wagyu and 115 other individuals were used to explore the genomic variations and genes under selection pressure in Wagyu. A total of 31 349 non-synonymous variants and 53 102 synonymous variants were identified in Wagyu. The population structure analysis showed that Wagyu had the closest genetic relationship with Mishima-Ushi cattle and was apparently separated from other cattle breeds. Then, composite likelihood ratio (CLR), integrated haplotype score, fixation index and cross-population composite likelihood ratio (XP-CLR) tests were performed to identify the candidate genes under positive selection in Wagyu. In total, 770 regions containing 312 genes were identified by at least three methods. Among them, 97 regions containing 27 genes were detected by all four methods. We specifically illustrate a list of interesting genes, including LRP2BP, GAA, CACNG6, CXADR, GPCPD1, KLF2, KLF13, SOX5, MYBPC1, SLC25A10, ATP8A1 and MYH15, which are associated with lipid metabolism, fat deposition, muscle development, bone development, feed intake and growth traits in Wagyu. This is the first study to explore the genomic variations and selection signatures of Wagyu at the whole-genome level. These results will provide significant help to beef cattle improvement and breeding.

Leaping through Tree Space: Continuous Phylogenetic Inference for Rooted and Unrooted Trees

Phylogenetics is now fundamental in life sciences, providing insights into the earliest branches of life and the origins and spread of epidemics. However, finding suitable phylogenies from the vast space of possible trees remains challenging. To address this problem, for the first time, we perform both tree exploration and inference in a continuous space where the computation of gradients is possible. This continuous relaxation allows for major leaps across tree space in both rooted and unrooted trees, and is less susceptible to convergence to local minima. Our approach outperforms the current best methods for inference on unrooted trees and, in simulation, accurately infers the tree and root in ultrametric cases. The approach is effective in cases of empirical data with negligible amounts of data, which we demonstrate on the phylogeny of jawed vertebrates. Indeed, only a few genes with an ultrametric signal were generally sufficient for resolving the major lineages of vertebrates. Optimization is possible via automatic differentiation and our method presents an effective way forward for exploring the most difficult, data-deficient phylogenetic questions.

Genetic architecture of ear traits based on association mapping and co-expression networks in maize inbred lines and hybrids

Ear traits are key contributors to grain yield in maize; therefore, exploring their genetic basis facilitates the improvement of grain yield. However, the underlying molecular mechanisms of ear traits remain obscure in both inbred lines and hybrids. Here, two association panels, respectively, comprising 203 inbred lines (IP) and 246 F1 hybrids (HP) were employed to identify candidate genes for six ear traits. The IP showed higher phenotypic variation and lower phenotypic mean than the HP for all traits, except ear tip-barrenness length. By conducting a genome-wide association study (GWAS) across multiple environments, 101 and 228 significant single-nucleotide polymorphisms (SNPs) associated with six ear traits were identified in the IP and HP, respectively. Of these significant SNPs identified in the HP, most showed complete-incomplete dominance and over-dominance effects for each ear trait. Combining a gene co-expression network with GWAS results, 186 and 440 candidate genes were predicted in the IP and HP, respectively, including known ear development genes ids1 and sid1. Of these, nine candidate genes were detected in both populations and expressed in maize ear and spikelet tissues. Furthermore, two key shared genes (GRMZM2G143330 and GRMZM2G171139) in both populations were found to be significantly associated with ear traits in the maize Goodman diversity panel with high-density variations. These findings advance our knowledge of the genetic architecture of ear traits between inbred lines and hybrids and provide a valuable resource for the genetic improvement of ear traits in maize.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01426-9.

Phylodynamic approaches to studying avian influenza virus

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

The episodic resurgence of highly pathogenic avian influenza H5 virus

Highly pathogenic avian influenza (HPAI) H5N1 activity has intensified globally since 2021, increasingly causing mass mortality in wild birds and poultry and incidental infections in mammals1-3. However, the ecological and virological properties that underscore future mitigation strategies still remain unclear. Using epidemiological, spatial and genomic approaches, we demonstrate changes in the origins of resurgent HPAI H5 and reveal significant shifts in virus ecology and evolution. Outbreak data show key resurgent events in 2016-2017 and 2020-2021, contributing to the emergence and panzootic spread of H5N1 in 2021-2022. Genomic analysis reveals that the 2016-2017 epizootics originated in Asia, where HPAI H5 reservoirs are endemic. In 2020-2021, 2.3.4.4b H5N8 viruses emerged in African poultry, featuring mutations altering HA structure and receptor binding. In 2021-2022, a new H5N1 virus evolved through reassortment in wild birds in Europe, undergoing further reassortment with low-pathogenic avian influenza in wild and domestic birds during global dissemination. These results highlight a shift in the HPAI H5 epicentre beyond Asia and indicate that increasing persistence of HPAI H5 in wild birds is facilitating geographic and host range expansion, accelerating dispersion velocity and increasing reassortment potential. As earlier outbreaks of H5N1 and H5N8 were caused by more stable genomic constellations, these recent changes reflect adaptation across the domestic-bird-wild-bird interface. Elimination strategies in domestic birds therefore remain a high priority to limit future epizootics.

SoyDNGP: a web-accessible deep learning framework for genomic prediction in soybean breeding

Soybean is a globally significant crop, playing a vital role in human nutrition and agriculture. Its complex genetic structure and wide trait variation, however, pose challenges for breeders and researchers aiming to optimize its yield and quality. Addressing this biological complexity requires innovative and accurate tools for trait prediction. In response to this challenge, we have developed SoyDNGP, a deep learning-based model that offers significant advancements in the field of soybean trait prediction. Compared to existing methods, such as DeepGS and DNNGP, SoyDNGP boasts a distinct advantage due to its minimal increase in parameter volume and superior predictive accuracy. Through rigorous performance comparison, including prediction accuracy and model complexity, SoyDNGP represents improved performance to its counterparts. Furthermore, it effectively predicted complex traits with remarkable precision, demonstrating robust performance across different sample sizes and trait complexities. We also tested the versatility of SoyDNGP across multiple crop species, including cotton, maize, rice and tomato. Our results showed its consistent and comparable performance, emphasizing SoyDNGP's potential as a versatile tool for genomic prediction across a broad range of crops. To enhance its accessibility to users without extensive programming experience, we designed a user-friendly web server, available at http://xtlab.hzau.edu.cn/SoyDNGP. The server provides two features: 'Trait Lookup', offering users the ability to access pre-existing trait predictions for over 500 soybean accessions, and 'Trait Prediction', allowing for the upload of VCF files for trait estimation. By providing a high-performing, accessible tool for trait prediction, SoyDNGP opens up new possibilities in the quest for optimized soybean breeding.

Evolution and distribution of rabies viruses from a panorama view

Rabies kills more than 59,000 people annually, mainly in developing countries. Previous studies on the evolution and distribution of rabies viruses (RABVs) were scattered. Here, we explore the evolution and distribution of this deadly virus from a novel panorama view. Multiple bioinformatic software tools were employed to analyze the phylogenetic diversity, evolution, spatiotemporal, and distribution of RABVs. The analyses were based on 1,202 qualified full-length genomes of RABVs and numerous literatures. Of the 10 distinct phylogenetic clades of RABV that we identified, more frequent intra- and inter-clade recombination occurs in the sequences of Asian-SEA, Arctic, and Cosmopolitan clades isolated from China, while according to existing sequence information, RABV might originate from bats (posterior probability, PP = 0.75, PP = 0.60 inferred from N and L genes, separately) in North America (PP = 0.57, PP = 0.62 inferred from N and L genes, separately). Due to the difference in evolutionary rate of N (2.22 × 10-4 subs/site/year, 95% HPD 1.99-2.47 × 10-4 subs/site/year) and L genes (1.67 × 10-4 subs/site/year, 95% HPD 1.59-1.74 × 10-4 subs/site/year), the root age was 1,406.6 (95% HPD 1,291.2-1,518.2) and 1,122.7 (95% HPD 1,052.4-1,193.9) inferred from N and L genes, separately. Among other findings, Mephitidae plays an important role in the interspecific transmission and communication of RABV, which we found tends to spread to populations genetically proximate to the host. We also identified amino acids under positive selection in different genes of different clades as well as single nucleotide variation sites important for different lineages. IMPORTANCE Rabies virus is widely distributed all over the world, and wild animals are its largest potential reservoir. Our study offers a panorama view about evolution and distribution of rabies viruses and emphasizes the need to monitor the transmission dynamics of animal rabies.

Niche breadth explains the range size of European-centred butterflies, but dispersal ability does not

Aim

The breadth of ecological niches and dispersal abilities have long been discussed as important determinants of species' range sizes. However, studies directly comparing the relative effects of both factors are rare, taxonomically biased and revealed inconsistent results.

Location

Europe.

Time Period

Cenozoic.

Major Taxa

Butterflies, Lepidoptera.

Methods

We relate climate, diet and habitat niche breadth and two indicators of dispersal ability, wingspan and a dispersal tendency index, to the global range size of 369 European-centred butterfly species. The relative effects of these five predictors and their variation across the butterfly phylogeny were assessed by means of phylogenetic generalized least squares models and phylogenetically weighted regressions respectively.

Results

Climate niche breadth was the most important single predictor, followed by habitat and diet niche breadth, while dispersal tendency and wingspan showed no relation to species' range size. All predictors together explained 59% of the variation in butterfly range size. However, the effects of each predictor varied considerably across families and genera.

Main Conclusions

Range sizes of European-centred butterflies are strongly correlated with ecological niche breadth but apparently independent of dispersal ability. The magnitude of range size-niche breadth relationships is not stationary across the phylogeny and is often negatively correlated across the different dimensions of the ecological niche. This variation limits the generalizability of range size-trait relationships across broad taxonomic groups.

Epidemiological and genomic investigation of chikungunya virus in Rio de Janeiro state, Brazil, between 2015 and 2018

Since 2014, Brazil has experienced an unprecedented epidemic caused by chikungunya virus (CHIKV), with several waves of East-Central-South-African (ECSA) lineage transmission reported across the country. In 2018, Rio de Janeiro state, the third most populous state in Brazil, reported 41% of all chikungunya cases in the country. Here we use evolutionary and epidemiological analysis to estimate the timescale of CHIKV-ECSA-American lineage and its epidemiological patterns in Rio de Janeiro. We show that the CHIKV-ECSA outbreak in Rio de Janeiro derived from two distinct clades introduced from the Northeast region in mid-2015 (clade RJ1, n = 63/67 genomes from Rio de Janeiro) and mid-2017 (clade RJ2, n = 4/67). We detected evidence for positive selection in non-structural proteins linked with viral replication in the RJ1 clade (clade-defining: nsP4-A481D) and the RJ2 clade (nsP1-D531G). Finally, we estimate the CHIKV-ECSA's basic reproduction number (R0) to be between 1.2 to 1.6 and show that its instantaneous reproduction number (Rt) displays a strong seasonal pattern with peaks in transmission coinciding with periods of high Aedes aegypti transmission potential. Our results highlight the need for continued genomic and epidemiological surveillance of CHIKV in Brazil, particularly during periods of high ecological suitability, and show that selective pressures underline the emergence and evolution of the large urban CHIKV-ECSA outbreak in Rio de Janeiro.

Dynamics of the rice yellow mottle disease in western Burkina Faso: Epidemic monitoring, spatio-temporal variation of viral diversity, and pathogenicity in a disease hotspot

The rice yellow mottle virus (RYMV) is a model in plant virus molecular epidemiology, with the reconstruction of historical introduction routes at the scale of the African continent. However, information on patterns of viral prevalence and viral diversity over multiple years at a local scale remains scarce, in spite of potential implications for crop protection. Here, we describe a 5-year (2015-9) monitoring of RYMV prevalence in six sites from western Burkina Faso (geographic areas of Bama, Banzon, and Karfiguela). It confirmed one irrigated site as a disease hotspot and also found one rainfed lowland (RL) site with occasional high prevalence levels. Within the studied fields, a pattern of disease aggregation was evidenced at a 5-m distance, as expected for a mechanically transmitted virus. Next, we monitored RYMV genetic diversity in the irrigated disease hotspot site, revealing a high viral diversity, with the current coexistence of various distinct genetic groups at the site scale (ca. 520 ha) and also within various specific fields (25 m side). One genetic lineage, named S1bzn, is the most recently emerged group and increased in frequency over the studied period (from 20 per cent or less in 2015-6 to more than 65 per cent in 2019). Its genome results from a recombination between two other lineages (S1wa and S1ca). Finally, experimental work revealed that three rice varieties commonly cultivated in Burkina Faso were not different in terms of resistance level, and we also found no significant effect of RYMV genetic groups on symptom expression and viral load. We found, however, that infection outcome depended on the specific RYMV isolate, with two isolates from the lineage S1bzn accumulating at the highest level at early infections. Overall, this study documents a case of high viral prevalence, high viral diversity, and co-occurrence of divergent genetic lineages at a small geographic scale. A recently emerged lineage, which comprises viral isolates inducing severe symptoms and high accumulation under controlled conditions, could be recently rising through natural selection. Following up the monitoring of RYMV diversity is required to confirm this trend and further understand the factors driving the local maintenance of viral diversity.

Investigating the Genetic Diversity of H5 Avian Influenza Viruses in the United Kingdom from 2020-2022

Since 2020, the United Kingdom and Europe have experienced annual epizootics of high-pathogenicity avian influenza virus (HPAIV). The first epizootic, during the autumn/winter of 2020-2021, involved six H5Nx subtypes, although H5N8 HPAIV dominated in the United Kingdom. While genetic assessments of the H5N8 HPAIVs within the United Kingdom demonstrated relative homogeneity, there was a background of other genotypes circulating at a lower degree with different neuraminidase and internal genes.  Following a small number of detections of H5N1 in wild birds over the summer of 2021, the autumn/winter of 2021-2022 saw another European H5 HPAIV epizootic that dwarfed the prior epizootic. This second epizootic was dominated almost exclusively by H5N1 HPAIV, although six distinct genotypes were defined. We have used genetic analysis to evaluate the emergence of different genotypes and proposed reassortment events that have been observed. The existing data suggest that the H5N1 viruses circulating in Europe during late 2020 continued to circulate in wild birds throughout 2021, with minimal adaptation, but then went on to reassort with AIVs in the wild bird population. We have undertaken an in-depth genetic assessment of H5 HPAIVs detected in the United Kingdom over two winter seasons and demonstrate the utility of in-depth genetic analyses in defining the diversity of H5 HPAIVs circulating in avian species, the potential for zoonotic risk, and whether incidents of lateral spread can be defined over independent incursions of infections from wild birds. This provides key supporting data for mitigation activities. IMPORTANCE High-pathogenicity avian influenza virus (HPAIV) outbreaks devastate avian species across all sectors, having both economic and ecological impacts through mortalities in poultry and wild birds, respectively. These viruses can also represent a significant zoonotic risk. Since 2020, the United Kingdom has experienced two successive outbreaks of H5 HPAIV. While H5N8 HPAIV was predominant during the 2020-2021 outbreak, other H5 subtypes were also detected. The following year, there was a shift in the subtype dominance to H5N1 HPAIV, but multiple H5N1 genotypes were detected. Through the thorough utilization of whole-genome sequencing, it was possible to track and characterize the genetic evolution of these H5 HPAIVs in United Kingdom poultry and wild birds. This enabled us to assess the risk posed by these viruses at the poultry-wild bird and the avian-human interfaces and to investigate the potential lateral spread between infected premises, a key factor in understanding the threat to the commercial sector.

The genomic landscape of swine influenza A viruses in Southeast Asia

Swine are a primary source for the emergence of pandemic influenza A viruses. The intensification of swine production, along with global trade, has amplified the transmission and zoonotic risk of swine influenza A virus (swIAV). Effective surveillance is essential to uncover emerging virus strains; however gaps remain in our understanding of the swIAV genomic landscape in Southeast Asia. More than 4,000 nasal swabs were collected from pigs in Cambodia, yielding 72 IAV-positive samples by RT-qPCR and 45 genomic sequences. We unmasked the cocirculation of multiple lineages of genetically diverse swIAV of pandemic concern. Genomic analyses revealed a novel European avian-like H1N2 swIAV reassortant variant with North American triple reassortant internal genes, that emerged approximately seven years before its first detection in pigs in 2021. Using phylogeographic reconstruction, we identified south central China as the dominant source of swine viruses disseminated to other regions in China and Southeast Asia. We also identified nine distinct swIAV lineages in Cambodia, which diverged from their closest ancestors between two and 15 B.P., indicating significant undetected diversity in the region, including reverse zoonoses of human H1N1/2009 pandemic and H3N2 viruses. A similar period of cryptic circulation of swIAVs occurred in the decades before the H1N1/2009 pandemic. The hidden diversity of swIAV observed here further emphasizes the complex underlying evolutionary processes present in this region, reinforcing the importance of genomic surveillance at the human-swine interface for early warning of disease emergence to avoid future pandemics.

Human-to-swine introductions and onward transmission of 2009 H1N1 pandemic influenza viruses in Brazil

Introduction

Once established in the human population, the 2009 H1N1 pandemic virus (H1N1pdm09) was repeatedly introduced into swine populations globally with subsequent onward transmission among pigs.

Methods

To identify and characterize human-to-swine H1N1pdm09 introductions in Brazil, we conducted a large-scale phylogenetic analysis of 4,141 H1pdm09 hemagglutinin (HA) and 3,227 N1pdm09 neuraminidase (NA) gene sequences isolated globally from humans and swine between 2009 and 2022.

Results

Phylodynamic analysis revealed that during the period between 2009 and 2011, there was a rapid transmission of the H1N1pdm09 virus from humans to swine in Brazil. Multiple introductions of the virus were observed, but most of them resulted in self-limited infections in swine, with limited onward transmission. Only a few sustained transmission clusters were identified during this period. After 2012, there was a reduction in the number of human-to-swine H1N1pdm09 transmissions in Brazil.

Discussion

The virus underwent continuous antigenic drift, and a balance was established between swine-to-swine transmission and extinction, with minimal sustained onward transmission from humans to swine. These results emphasize the dynamic interplay between human-to-swine transmission, antigenic drift, and the establishment of swine-to-swine transmission in shaping the evolution and persistence of H1N1pdm09 in swine populations.

Unprecedented Antimicrobial and Cytotoxic Polyketides from Cultures of Diaporthe africana sp. nov

Four unprecedented polyketides named isoprenylisobenzofuran B (2), isoprenylisobenzofuran C₁/C₂ (3), diaporisoindole F₁/F₂ (4), and isochromophilonol A₁/A₂ (7) were isolated from ethyl acetate extracts of the newly described endophytic fungus Diaporthe africana. Additionally, the previously reported cyclic depsipeptide eucalactam B (1) was also identified, along with the known compounds diaporisoindole A/B (5), tenellone B (6) and beauvericin (8). The taxonomic identification of the fungus was accomplished using a polyphasic approach combining multi-gene phylogenetic analysis and microscopic morphological characters. The structures 1-8 were determined by a detailed analysis of their spectral data, namely high-resolution electrospray ionization mass spectrometry (HR-ESIMS), 1D/2D nuclear magnetic resonance (NMR) spectroscopy, as well as electronic circular dichroism (ECD) spectra. In addition, chemical methods such as Marfey's analysis were also employed to determine the stereochemistry in compound 1. All the compounds obtained were evaluated for antimicrobial and in vitro cytotoxic properties. Compounds 3-8 were active against certain fungi and Gram-positive bacteria with MIC values of 8.3 to 66.6 µg/mL. In addition, 3-5 displayed cytotoxic effects (22.0 ≤ IC₅₀ ≤ 59.2 µM) against KB3.1 and L929 cell lines, whereas compounds 6-8 inhibited the growth of seven mammalian cancer cell lines with IC₅₀ ranging from 17.7 to 49.5 µM (6), 0.9 to 12.9 µM (7) and 1.9 to 4.3 µM (8).

Herbarium specimen sequencing allows precise dating of Xanthomonas citri pv. citri diversification history

Herbarium collections are an important source of dated, identified and preserved DNA, whose use in comparative genomics and phylogeography can shed light on the emergence and evolutionary history of plant pathogens. Here, we reconstruct 13 historical genomes of the bacterial crop pathogen Xanthomonas citri pv. citri (Xci) from infected Citrus herbarium specimens. Following authentication based on ancient DNA damage patterns, we compare them with a large set of modern genomes to estimate their phylogenetic relationships, pathogenicity-associated gene content and several evolutionary parameters. Our results indicate that Xci originated in Southern Asia ~11,500 years ago (perhaps in relation to Neolithic climate change and the development of agriculture) and diversified during the beginning of the 13th century, after Citrus diversification and before spreading to the rest of the world (probably via human-driven expansion of citriculture through early East-West trade and colonization).

Pharmacovariome scanning using whole pharmacogene resequencing coupled with deep computational analysis and machine learning for clinical pharmacogenomics

BACKGROUND

This pilot study aims to identify and functionally assess pharmacovariants in whole exome sequencing data. While detection of known variants has benefited from pharmacogenomic-dedicated bioinformatics tools before, in this paper we have tested novel deep computational analysis in addition to artificial intelligence as possible approaches for functional analysis of unknown markers within less studied drug-related genes.

METHODS

Pharmacovariants from 1800 drug-related genes from 100 WES data files underwent (a) deep computational analysis by eight bioinformatic algorithms (overall containing 23 tools) and (b) random forest (RF) classifier as the machine learning (ML) approach separately. ML model efficiency was calculated by internal and external cross-validation during recursive feature elimination. Protein modelling was also performed for predicted highly damaging variants with lower frequencies. Genotype-phenotype correlations were implemented for top selected variants in terms of highest possibility of being damaging.

RESULTS

Five deleterious pharmacovariants in the RYR1, POLG, ANXA11, CCNH, and CDH23 genes identified in step (a) and subsequent analysis displayed high impact on drug-related phenotypes. Also, the utilization of recursive feature elimination achieved a subset of 175 malfunction pharmacovariants in 135 drug-related genes that were used by the RF model with fivefold internal cross-validation, resulting in an area under the curve of 0.9736842 with an average accuracy of 0.9818 (95% CI: 0.89, 0.99) on predicting whether a carrying individuals will develop adverse drug reactions or not. However, the external cross-validation of the same model indicated a possible false positive result when dealing with a low number of observations, as only 60 important variants in 49 genes were displayed, giving an AUC of 0.5384848 with an average accuracy of 0.9512 (95% CI: 0.83, 0.99).

CONCLUSION

While there are some technologies for functionally assess not-interpreted pharmacovariants, there is still an essential need for the development of tools, methods, and algorithms which are able to provide a functional prediction for every single pharmacovariant in both large-scale datasets and small cohorts. Our approaches may bring new insights for choosing the right computational assessment algorithms out of high throughput DNA sequencing data from small cohorts to be used for personalized drug therapy implementation.

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.

First report of Colletotrichum fioriniae infections in brown marmorated stink bugs, Halyomorpha halys

An epizootic caused by fungal pathogens occurred among Halyomorpha halys, brown marmorated stink bugs, while they were overwintering, with infections also occurring after overwintering. We report that one of the two pathogens responsible was Collectotrichum fioriniae (Marcelino & Gouli), Pennycook; a species well known as a plant pathogen and endophyte and which has only previously been reported naturally infecting elongate hemlock scales, Fiorinia externa. To prove pathogenicity, H. halys adults challenged with conidia died from infections and the fungus subsequently produced conidia externally on cadavers.

The evolutionary dynamics history of canine distemper virus through analysis of the hemagglutinin gene during 1930-2020

Canine distemper virus (CDV) is a lethal viral disease of carnivores which is considered to be a serious threat to domestic and wild species. Despite the widespread use of vaccines, CDV still occurs in vaccinated animals and current vaccines does not guarantee complete protection. In this study, a total of 286 hemagglutinin (H) gene sequences of the virus isolated in 25 countries during 90 years (1930-2020) were analyzed by Bayesian maximum likelihood analysis to estimate the population dynamics. We identified the most recent common ancestor (TMRCA) of the virus in 1868 in the USA which arrived in continental Europe in 1948, and from there, the virus spread rapidly to other continents. The Canidae family was identified as the original host as well as a source of the subsequent spread. We identified 11 lineages of geographic co-circulating strains globally. The effective population size experienced a two-phase-exponential growth between 2000-2005 and 2010-2012. Our findings provide a novel insight into the epidemic history of canine distemper virus which may facilitate more effective disease management. This study uses a large set of sequencing data on the H gene of CDV to identify distinct lineages of the virus, track its geographic spread over time, analyze its likelihood of transmission within and between animal families, and provide suggestions for improved strategies to combat the virus.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10344-023-01685-z.

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.

Emergence of a novel genotype of class II New Castle Disease virus in North Eastern States of India

Outbreaks of New Castle Disease from three north eastern states of India were confirmed by clinico-pathological examination followed by reverse transcription-PCR detection of F gene of ND Virus (NDV). Irrespective of vaccination, the outbreaks resulted 90-100% mortality in the affected flocks. The analysis of fusion protein sequences from ten field isolates revealed them as the velogenic or highly virulent strain. Phylogenetic analyses based on the complete F gene nucleotide sequences of the isolates have characterized only one of the isolate (OK149201) in the genotype XIII.2.2. The rest of the nine isolates are depicted in a distinct monophyletic group with average nucleotide distances from the other 20 genotypes ranged from 10.90 - 20.70. The nine isolates were further divided into two sub branches with the bootstrap support value of 100% at the nodes that define the two subgroups with an average evolutionary nucleotide distance of 6.00between the isolates in the two subgroups. As per the recommendation put forth in recently updated unified phylogenetic classification system for NDV, our findings clearly indicates emergence of a novel genotype of class II NDV in the biodiversity hot spot region of NER, India. The isolates in the newly identified genotype is designated with next available Roman numerals XXII. Further, the two subgroups within the genotype are designated as XXII.1 and XXII.2.

Redisposition of apiosporous genera Induratia and Muscodor in the Xylariales, following the discovery of an authentic strain of Induratia apiospora

BACKGROUND

The genus Induratia is based on Induratia apiospora, a xylarialean pyrenomycete from New Zealand with clypeate uniperitheciate stromata, hyaline apiospores and a nodulisporium-like anamorph. However, because of the lack of DNA data from the generic type, its phylogenetic affinities have remained unresolved. Recently, two fungal species with teleomorphs strikingly similar to Induratia were discovered in Thailand. However, they did not produce an anamorph and were found to be phylogenetically close to the species classified within the hyphomycete genus Muscodor, which was described after Induratia. Therefore, in 2020 the species of Muscodor were transferred to Induratia, and a new family Induratiaceae was proposed.

RESULTS

We have encountered an unpublished ex-holotype strain of Induratia apiospora among the holdings of the ATCC collection, enabling detailed morphological and molecular phylogenetic investigations. We observed the characteristic nodulisporium-like anamorph described in the original publication. Phylogenetic analyses of multigene sequence data revealed a close relationship of Induratia apiospora to the Barrmaeliaceae, while a close relationship to the Induratia species formerly classified within Muscodor was rejected.

CONCLUSIONS

We here classify Induratia apiospora within the Barrmaeliaceae and consider Induratiaceae to be synonymous with the former. As the holotype specimen of Induratia apiospora is apparently lost, an isotype specimen from WSP is selected as lectotype. We also propose that the genus Muscodor is resurrected within the Xylariaceae, and formally transfer several Induratia species to Muscodor.

Estimates of genomic inbreeding and identification of candidate regions in Beijing-You chicken populations

Runs of homozygosity (ROHs) has become an effective method for analysing inbreeding in livestock populations. Moreover, ROHs is well-suited to detect signatures of selection via ROH islands. This study aimed to investigate the occurrence and distribution of ROHs, compare the genomic inbreeding coefficients and identify the genomic regions with high ROH frequencies in different Beijing-You chicken (BY) populations, including a random conservation population (BY_R), a pedigree conservation population (BY_P), and a commercial population obtained from the market (BY_C). Among them, BY_R in 2010 and 2019 were BY_R1 and BY_R2 respectively. A total of 27 916 ROHs were identified. The average number of ROHs per individual across the three BY populations ranged from 213 (BY_P) to 161 (BY_C), and the average length of ROHs ranged from 0.432 Mb (BY_R2) to 0.451 Mb (BY_P). The highest inbreeding coefficient calculated based on ROHs (F_ROH ) was 0.1 in BY_P, whereas the lowest F_ROH was 0.0743 in BY_C. In addition, the inbreeding coefficient of BY_R2 (F_ROH  = 0.0798) was higher than that of BY_R1 (F_ROH  = 0.0579). Furthermore, the highest proportion of long ROH fragments (>4 Mb) was observed in BY_P and BY_C. This study showed the top 10 ROH islands of each population, and these ROH islands harboured 53 genes, some of which were related to limb development, body size and immune response. These findings contribute to the understanding of genetic diversity and population demography, and might help improve breeding and conservation strategies for BY populations.

Comparative pathogenicity of a genotype XXI.1.2 pigeon Newcastle disease virus isolate in pigeons and chickens

Here, we performed molecular and pathogenic characterization of a Newcastle disease virus (NDV) isolate from pigeons in Bangladesh. Molecular phylogenetic analysis based on the complete fusion gene sequences classified the three study isolates into genotype XXI (sub-genotype XXI.1.2) together with recent NDV isolates obtained from pigeons in Pakistan (2014-2018). The Bayesian Markov Chain Monte Carlo analysis revealed that the ancestor of Bangladeshi pigeon NDVs and the viruses from sub-genotype XXI.1.2 existed in the late 1990s. Pathogenicity testing using mean embryo death time pathotyped the viruses as mesogenic, while all isolates carried multiple basic amino acid residues at the fusion protein cleavage site. Experimental infection of chickens and pigeons revealed no or minimum clinical signs in chickens, while a relatively high morbidity (70%) and mortality (60%) were observed in pigeons. The infected pigeons showed extensive and systemic lesions including hemorrhagic and/or vascular changes in the conjunctiva, respiratory and digestive system and brain, and atrophy in the spleen, while only mild congestion in the lungs was noticed in the inoculated chickens. Histologically, consolidation in the lungs with collapsed alveoli and edema around the blood vessels, hemorrhages in the trachea, severe hemorrhages and congestion, focal aggregation of mononuclear cells, and single hepatocellular necrosis in the liver, severe congestion, multifocal tubular degeneration, and necrosis, as well as mononuclear cell infiltration in the renal parenchyma, encephalomalacia with severe neuronal necrosis with neuronophagia were noticed in the brain in infected pigeons. In contrast, only slight congestion was found in lungs of the infected chickens. qRT-PCR revealed the replication of the virus in both pigeons and chickens; however, higher viral RNA loads were observed in oropharyngeal and cloacal swabs, respiratory tissues, and spleen of infected pigeons than the chickens. In conclusion, genotype XXI.1.2 NDVs are circulating in the pigeon population of Bangladesh since 1990s, produce high mortality in pigeons with pneumonia, hepatocellular necrosis, renal tubular degeneration, and neuronal necrosis in pigeons, and may infect chickens without overt signs of clinical disease and are likely to shed viruses via the oral or cloacal routes.

Cardoso J, Gonçalves-dos-Santos ME, Oliveira RS, Aquino-Teixeira SM, Campos AAS, Almeida MAB, Simonini-Teixeira D, da P. Sevá A, Temponi AOD, Magalhães FM, da Silva Menezes AS, Lopes BT, Almeida HP, Pedroso AL, Gonçalves GP, Chaves DCC, de Menezes GG, Bernal-Valle S, Müller NFD, Janssen L, dos Santos E, Mares-Guia MA, Albuquerque GR, Romano APM, Franco AC, Ribeiro BM, Roehe PM, Lourenço-de-Oliveira R, de Abreu FVS. Yellow Fever Virus Maintained by Sabethes Mosquitoes during the Dry Season in Cerrado, a Semiarid Region of Brazil, in 2021

In recent decades, waves of yellow fever virus (YFV) from the Amazon Rainforest have spread and caused outbreaks in other regions of Brazil, including the Cerrado, a savannah-like biome through which YFV usually moves before arriving at the Atlantic Forest. To identify the vectors involved in the maintenance of the virus in semiarid environments, an entomological survey was conducted after confirmation of yellow fever (YF) epizootics at the peak of the dry season in the Cerrado areas of the state of Minas Gerais. In total, 917 mosquitoes from 13 taxa were collected and tested for the presence of YFV. Interestingly, mosquitoes of the Sabethes genus represented 95% of the diurnal captured specimens, displaying a peak of biting activity never previously recorded, between 4:30 and 5:30 p.m. Molecular analysis identified three YFV-positive pools, two from Sabethes chloropterus-from which near-complete genomes were generated-and one from Sa. albiprivus, whose low viral load prevented sequencing. Sa. chloropterus was considered the primary vector due to the high number of copies of YFV RNA and the high relative abundance detected. Its bionomic characteristics allow its survival in dry places and dry time periods. For the first time in Brazil, Sa. albiprivus was found to be naturally infected with YFV and may have played a role as a secondary vector. Despite its high relative abundance, fewer copies of viral RNA were found, as well as a lower Minimum Infection Rate (MIR). Genomic and phylogeographic analysis showed that the virus clustered in the sub-lineage YFVPA-MG, which circulated in Pará in 2017 and then spread into other regions of the country. The results reported here contribute to the understanding of the epidemiology and mechanisms of YFV dispersion and maintenance, especially in adverse weather conditions. The intense viral circulation, even outside the seasonal period, increases the importance of surveillance and YFV vaccination to protect human populations in affected areas.

A comprehensive phylogeny and revised taxonomy illuminate the origin and diversification of the global radiation of Papilio (Lepidoptera: Papilionidae)

The swallowtail genus Papilio (Lepidoptera: Papilionidae) is species rich, distributed worldwide, and has broad morphological habits and ecological niches. Because of its elevated species richness, it has been historically difficult to reconstruct a densely sampled phylogeny for this clade. Here we provide a taxonomic working list for the genus, resulting in 235 Papilio species, and assemble a molecular dataset of seven gene fragments representing ca. 80% of the currently described diversity. Phylogenetic analyses reconstructed a robust tree with highly supported relationships within subgenera, although a few nodes in the early history of the Old World Papilio remain unresolved. Contrasting with previous results, we found that Papilio alexanor is sister to all Old World Papilio and that the subgenus Eleppone is no longer monotypic. The latter includes the recently described Fijian Papilio natewa with the Australian Papilio anactus and is sister to subgenus Araminta (formerly included in subgenus Menelaides) occurring in Southeast Asia. Our phylogeny also includes rarely studied (P. antimachus, P. benguetana) or endangered species (P. buddha, P. chikae). Taxonomic changes resulting from this study are elucidated. Molecular dating and biogeographic analyses indicate that Papilio originated ca. 30 million years ago (Oligocene), in a northern region centered on Beringia. A rapid early Miocene radiation in the Paleotropics is revealed within Old World Papilio, potentially explaining their low early branch support. Most subgenera originated in the early to middle Miocene followed by synchronous southward biogeographic dispersals and repeated local extirpations in northern latitudes. This study provides a comprehensive phylogenetic framework for Papilio with clarification of subgeneric systematics and species taxonomic changes enumerated, which will facilitate further studies to address questions on their ecology and evolutionary biology using this model clade.

Emergence and spread of two SARS-CoV-2 variants of interest in Nigeria

Identifying the dissemination patterns and impacts of a virus of economic or health importance during a pandemic is crucial, as it informs the public on policies for containment in order to reduce the spread of the virus. In this study, we integrated genomic and travel data to investigate the emergence and spread of the SARS-CoV-2 B.1.1.318 and B.1.525 (Eta) variants of interest in Nigeria and the wider Africa region. By integrating travel data and phylogeographic reconstructions, we find that these two variants that arose during the second wave in Nigeria emerged from within Africa, with the B.1.525 from Nigeria, and then spread to other parts of the world. Data from this study show how regional connectivity of Nigeria drove the spread of these variants of interest to surrounding countries and those connected by air-traffic. Our findings demonstrate the power of genomic analysis when combined with mobility and epidemiological data to identify the drivers of transmission, as bidirectional transmission within and between African nations are grossly underestimated as seen in our import risk index estimates.

Impact and mitigation of sampling bias to determine viral spread: Evaluating discrete phylogeography through CTMC modeling and structured coalescent model approximations

Bayesian phylogeographic inference is a powerful tool in molecular epidemiological studies, which enables reconstruction of the origin and subsequent geographic spread of pathogens. Such inference is, however, potentially affected by geographic sampling bias. Here, we investigated the impact of sampling bias on the spatiotemporal reconstruction of viral epidemics using Bayesian discrete phylogeographic models and explored different operational strategies to mitigate this impact. We considered the continuous-time Markov chain (CTMC) model and two structured coalescent approximations (Bayesian structured coalescent approximation [BASTA] and marginal approximation of the structured coalescent [MASCOT]). For each approach, we compared the estimated and simulated spatiotemporal histories in biased and unbiased conditions based on the simulated epidemics of rabies virus (RABV) in dogs in Morocco. While the reconstructed spatiotemporal histories were impacted by sampling bias for the three approaches, BASTA and MASCOT reconstructions were also biased when employing unbiased samples. Increasing the number of analyzed genomes led to more robust estimates at low sampling bias for the CTMC model. Alternative sampling strategies that maximize the spatiotemporal coverage greatly improved the inference at intermediate sampling bias for the CTMC model, and to a lesser extent, for BASTA and MASCOT. In contrast, allowing for time-varying population sizes in MASCOT resulted in robust inference. We further applied these approaches to two empirical datasets: a RABV dataset from the Philippines and a SARS-CoV-2 dataset describing its early spread across the world. In conclusion, sampling biases are ubiquitous in phylogeographic analyses but may be accommodated by increasing the sample size, balancing spatial and temporal composition in the samples, and informing structured coalescent models with reliable case count data.

Ecological, Genetic, and Phylogenetic Aspects of YFV 2017-2019 Spread in Rio de Janeiro State

In Brazil, a yellow fever (YF) outbreak was reported in areas considered YF-free for decades. The low vaccination coverage and the increasing forest fragmentation, with the wide distribution of vector mosquitoes, have been related to yellow fever virus (YFV) transmission beyond endemic areas since 2016. Aiming to elucidate the molecular and phylogenetic aspects of YFV spread on a local scale, we generated 43 new YFV genomes sampled from humans, non-human primates (NHP), and primarily, mosquitoes from highly heterogenic areas in 15 localities from Rio de Janeiro (RJ) state during the YFV 2016-2019 outbreak in southeast Brazil. Our analysis revealed that the genetic diversity and spatial distribution of the sylvatic transmission of YFV in RJ originated from at least two introductions and followed two chains of dissemination, here named the YFV RJ-I and YFV RJ-II clades. They moved with similar dispersal speeds from the north to the south of the RJ state in parallel directions, separated by the Serra do Mar Mountain chain, with YFV RJ-I invading the north coast of São Paulo state. The YFV RJ-I clade showed a more significant heterogeneity across the entire polyprotein. The YFV RJ-II clade, with only two amino acid polymorphisms, mapped at NS1 (I1086V), present only in mosquitoes at the same locality and NS4A (I2176V), shared by all YFV clade RJ-II, suggests a recent clustering of YFV isolates collected from different hosts. Our analyses strengthen the role of surveillance, genomic analyses of YVF isolated from other hosts, and environmental studies into the strategies to forecast, control, and prevent yellow fever outbreaks.

Bidirectional Movement of Emerging H5N8 Avian Influenza Viruses Between Europe and Asia via Migratory Birds Since Early 2020

Migratory birds play a critical role in the rapid spread of highly pathogenic avian influenza (HPAI) H5N8 virus clade 2.3.4.4 across Eurasia. Elucidating the timing and pattern of virus transmission is essential therefore for understanding the spatial dissemination of these viruses. In this study, we surveyed >27,000 wild birds in China, tracked the year-round migration patterns of 20 bird species across China since 2006, and generated new HPAI H5N8 virus genomic data. Using this new data set, we investigated the seasonal transmission dynamics of HPAI H5N8 viruses across Eurasia. We found that introductions of HPAI H5N8 viruses to different Eurasian regions were associated with the seasonal migration of wild birds. Moreover, we report a backflow of HPAI H5N8 virus lineages from Europe to Asia, suggesting that Europe acts as both a source and a sink in the global HPAI virus transmission network.

Distribution of Bemisia tabaci in different agro-ecological regions in Uganda and the threat of vector-borne pandemics into new cassava growing areas

Previous studies in sub-Saharan Africa have showed the spread of cassava mosaic disease (CMD) and cassava brown streak disease (CBSD) pandemics into different cassava growing regions by high Bemisia tabaci populations. Studies did indicate that there were stark differences in some whitefly species, yet they have not looked extensively across agroecologies. Members of B. tabaci species complex termed sub-Saharan Africa 1 (SSA1) and SSA2 have been linked to the spread of CMD and CBSD viruses. During the period of a severe CMD pandemic in the 1990s, SSA2 was the most predominant until the resurgence of SSA1, particularly SSA1-subgroup1 (SSA1-SG1) from the early 2000s to date. Cassava being a drought resilient crop has become an important food security crop and has been introduced into new areas and regions. Considering the role B. tabaci in the spread of cassava virus pandemics into neighboring regions, we investigated the genetic diversity and distribution of B. tabaci in nine different agro-ecological regions of Uganda in 2017. Adult whiteflies were collected from cassava and 33 other host plants from cassava-growing areas, those with limited cassava and areas with no cassava, where it is being introduced as a food security crop. The partial sequences of the mitochondrial cytochrome oxidase 1 (mtCO1) gene (657 bp) were used to determine the phylogenetic relationships between the sampled B. tabaci. Cassava B. tabaci SSA1 (-SG1, -SG2, -Hoslundia (previously called SSA1-SG1/2), -SG3), SSA2 and SSA3; non-cassava B. tabaci SSA6, SSA10, SSA11, SSA12, SSA13, MED-ASL, MED-Q1, MEAM1, Indian Ocean; and other Bemisia species, Bemisia afer and Bemisia Uganda1 were identified in the study. SSA3, one of the key B. tabaci species that occurs on cassava in West Africa, was identified for the first time in Uganda. The SSA1-SG1 was widely distributed, predominated on cassava and was found on 17 other host-plants. The ability of SSA1-SG1 to exist in environments with limited or no cassava growing poses the risk of continued spread of virus pandemics. Therefore, measures must be put in place to prevent the introduction of diseased materials into new areas, since the vectors exist.

Independent recruitment of FRUITFULL-like transcription factors in the convergent origins of vernalization-responsive grass flowering

Flowering in response to low temperatures (vernalization) has evolved multiple times independently across angiosperms as an adaptation to match reproductive development with the short growing season of temperate habitats. Despite the context of a generally conserved flowering time network, evidence suggests that the genes underlying vernalization responsiveness are distinct across major plant clades. Whether different or similar mechanisms underlie vernalization-induced flowering at narrower (e.g., family-level) phylogenetic scales is not well understood. To test the hypothesis that vernalization responsiveness has evolved convergently in temperate species of the grass family (Poaceae), we carried out flowering time experiments with and without vernalization in several representative species from different subfamilies. We then determined the likelihood that vernalization responsiveness evolved through parallel mechanisms by quantifying the response of Pooideae vernalization pathway FRUITFULL (FUL)-like genes to extended periods of cold. Our results demonstrate that vernalization-induced flowering has evolved multiple times independently in at least five grass subfamilies, and that different combinations of FUL-like genes have been recruited to this pathway on several occasions.

New TNT routines for parallel computing with MPI

Phylogenetic inference, which involves time-consuming calculations, is a field where parallelization can speed up the resolution of many problems. TNT (a widely used program for phylogenetic analysis under parsimony) allows parallelization under the PVM system (Parallel Virtual Machine). However, as the basic aspects of the implementation remain unpublished, few studies have taken advantage of the parallelization routines of TNT. In addition, the PVM system is deprecated by many system administrators. One of the most common standards for high performance computing is now MPI (Message Passing Interface). To facilitate the use of the parallel analyses offered by TNT, this paper describes the basic aspects of the implementation, as well as a port of the parallelization interface of TNT into MPI. The use of the new routines is illustrated by reanalysis of seven significant datasets, either recent phylogenomic datasets with many characters (up to 2,509,064 characters) or datasets with large numbers of taxa (up to 13,921 taxa). Versions of TNT including the MPI functionality are available at: http://www.lillo.org.ar/phylogeny/tnt/.

RNA structure-altering mutations underlying positive selection on Spike protein reveal novel putative signatures to trace crossing host-species barriers in Betacoronavirus

Similar to other RNA viruses, the emergence of Betacoronavirus relies on cross-species viral transmission, which requires careful health surveillance monitoring of protein-coding information as well as genome-wide analysis. Although the evolutionary jump from natural reservoirs to humans may be mainly traced-back by studying the effect that hotspot mutations have on viral proteins, it is largely unexplored if other impacts might emerge on the structured RNA genome of Betacoronavirus. In this survey, the protein-coding and viral genome architecture were simultaneously studied to uncover novel insights into cross-species horizontal transmission events. We analysed 1,252,952 viral genomes of SARS-CoV, MERS-CoV, and SARS-CoV-2 distributed across the world in bats, intermediate animals, and humans to build a new landscape of changes in the RNA viral genome. Phylogenetic analyses suggest that bat viruses are the most closely related to the time of most recent common ancestor of Betacoronavirus, and missense mutations in viral proteins, mainly in the S protein S1 subunit: SARS-CoV (G > T; A577S); MERS-CoV (C > T; S746R and C > T; N762A); and SARS-CoV-2 (A > G; D614G) appear to have driven viral diversification. We also found that codon sites under positive selection on S protein overlap with non-compensatory mutations that disrupt secondary RNA structures in the RNA genome complement. These findings provide pivotal factors that might be underlying the eventual jumping the species barrier from bats to intermediate hosts. Lastly, we discovered that nearly half of the Betacoronavirus genomes carry highly conserved RNA structures, and more than 90% of these RNA structures show negative selection signals, suggesting essential functions in the biology of Betacoronavirus that have not been investigated to date. Further research is needed on negatively selected RNA structures to scan for emerging functions like the potential of coding virus-derived small RNAs and to develop new candidate antiviral therapeutic strategies.

Ancestral reconstruction reveals catalytic inactivation of activation-induced cytidine deaminase concomitant with cold water adaption in the Gadiformes bony fish

BACKGROUND

Antibody affinity maturation in vertebrates requires the enzyme activation-induced cytidine deaminase (AID) which initiates secondary antibody diversification by mutating the immunoglobulin loci. AID-driven antibody diversification is conserved across jawed vertebrates since bony and cartilaginous fish. Two exceptions have recently been reported, the Pipefish and Anglerfish, in which the AID-encoding aicda gene has been lost. Both cases are associated with unusual reproductive behavior, including male pregnancy and sexual parasitism. Several cold water fish in the Atlantic cod (Gadinae) family carry an aicda gene that encodes for a full-length enzyme but lack affinity-matured antibodies and rely on antibodies of broad antigenic specificity. Hence, we examined the functionality of their AID.

RESULTS

By combining genomics, transcriptomics, immune responsiveness, and functional enzymology of AID from 36 extant species, we demonstrate that AID of that Atlantic cod and related fish have extremely lethargic or no catalytic activity. Through ancestral reconstruction and functional enzymology of 71 AID enzymes, we show that this enzymatic inactivation likely took place relatively recently at the emergence of the true cod family (Gadidae) from their ancestral Gadiformes order. We show that this AID inactivation is not only concordant with the previously shown loss of key adaptive immune genes and expansion of innate and cell-based immune genes in the Gadiformes but is further reflected in the genomes of these fish in the form of loss of AID-favored sequence motifs in their immunoglobulin variable region genes.

CONCLUSIONS

Recent demonstrations of the loss of the aicda gene in two fish species challenge the paradigm that AID-driven secondary antibody diversification is absolutely conserved in jawed vertebrates. These species have unusual reproductive behaviors forming an evolutionary pressure for a certain loss of immunity to avoid tissue rejection. We report here an instance of catalytic inactivation and functional loss of AID rather than gene loss in a conventionally reproducing vertebrate. Our data suggest that an expanded innate immunity, in addition to lower pathogenic pressures in a cold environment relieved the pressure to maintain robust secondary antibody diversification. We suggest that in this unique scenario, the AID-mediated collateral genome-wide damage would form an evolutionary pressure to lose AID function.