The presence of identical deformed wing virus sequence variants in co-occurring Apis species in Northern Thailand may represent a potential epidemiological threat to native honey bees of Southeast Asia.

Widespread native honey bee species in South and East Asia (Apis cerana, Apis dorsata and Apis florea) and the imported western honey bee (Apis mellifera) share habitats and potentially also share pathogens. Chief among the threats facing A. mellifera in Europe and North America is deformed wing virus (DWV), including its two principal genotypes: A and B (DWV-A and DWV-B respectively). Though DWV-A has been recorded in Asia's native Apis species, it is not known if DWV-B, or both DWV-A and DWV-B, are currently widespread in Asia and, if so, whether viral transmission is primarily intraspecific or interspecific. This study aims to fill these knowledge gaps by (i) determining the DWV genotype in four co-occurring Apis host species using qPCR and (ii) inferring viral transmission between them using nucleotide sequences of DWV from Apis host species collected at three independent localities in Northern Thailand. We found DWV-A and -B in all four Apis species, the exotic A. mellifera and the native A. cerana, A. dorsata and A. florea. That DWV-A sequences were identical across Apis species at the same locality, with a similar pattern for DWV-B sequences, suggests that DWV's epidemiology is largely driven by ongoing interspecific transmission (spillover) of DWV across co-occurring native and exotic Apis species. Both genotypes of DWV represent a serious threat to Asia's exotic and native honey bee species.

Transcriptomic Responses Underlying the High Virulence of Black Queen Cell Virus and Sacbrood Virus following a Change in Their Mode of Transmission in Honey Bees (Apis mellifera).

BACKGROUND: Over the last two decades, honey bees (Apis mellifera) have suffered high rates of colony losses that have been attributed to a variety of factors, chief among which are viral pathogens, such as deformed wing virus (DWV), whose virulence has increased because of vector-based transmission by the invasive, ectoparasitic varroa mite (Varroa destructor). A shift in the experimental mode of transmission of the black queen cell virus (BQCV) and sacbrood virus (SBV) from fecal/food-oral (direct horizontal) to vector-mediated (indirect horizontal) transmission also results in high virulence and viral titers in pupal and adult honey bees. Agricultural pesticides represent another factor that acts independently or in interaction with pathogens, and they are also thought to cause colony loss. Understanding the molecular mechanisms underlying the higher virulence following a vector-based mode of transmission provides deeper insight into honey bee colony losses, as does determining whether or not host-pathogen interactions are modulated by exposure to pesticides. METHODS: Through an experimental design with controlled laboratory, we investigated the effects of the modes of transmission of BQCV and SBV (feeding vs. vector-mediated via injection) alone or in combination with chronic exposure to sublethal and field-realistic concentrations of flupyradifurone (FPF), a novel agricultural insecticide, on honey bee survival and transcription responses by using high-throughput RNA sequencing (RNA-seq) analysis. RESULTS: Co-exposure to viruses via feeding (VF) or injection (VI) and FPF insecticide had no statistically significant interactive effect on their survival compared to, respectively, VF or VI treatments alone. Transcriptomic analysis revealed a distinct difference in the gene expression profiles of bees inoculated with viruses via injection (VI) and exposed to FPF insecticide (VI+FPF). The number of differentially expressed genes (DEGs) at log2 (fold-change) > 2.0 in VI bees (136 genes) or/and VI+FPF insecticide (282 genes) was very high compared to that of VF bees (8 genes) or the VF+FPF insecticide treatment (15 genes). Of these DEGs, the expression in VI and VI+FPF bees of some immune-related genes, such as those for antimicrobial peptides, Ago2, and Dicer, was induced. In short, several genes encoding odorant binding proteins, chemosensory proteins, odor receptors, honey bee venom peptides, and vitellogenin were downregulated in VI and VI+FPF bees. CONCLUSIONS: Given the importance of these suppressed genes in honey bees' innate immunity, eicosanoid biosynthesis, and olfactory associative function, their inhibition because of the change in the mode of infection with BQCV and SBV to vector-mediated transmission (injection into haemocoel) could explain the high virulence observed in these viruses when they were experimentally injected into hosts. These changes may help explain why other viruses, such as DWV, represent such a threat to colony survival when transmitted by varroa mites.

Epidemiology of a major honey bee pathogen, deformed wing virus: potential worldwide replacement of genotype A by genotype B.

The western honey bee (Apis mellifera) is of major economic and ecological importance, with elevated rates of colony losses in temperate regions over the last two decades thought to be largely caused by the exotic ectoparasitic mite Varroa destructor and deformed wing virus (DWV), which the mite transmits. DWV currently exists as two main genotypes: the formerly widespread DWV-A and the more recently described and rapidly expanding DWV-B. It is an excellent system to understand viral evolution and the replacement of one viral variant by another. Here we synthesise published results on the distribution and prevalence of DWV-A and -B over the period 2008-2021 and present novel data for Germany, Italy and the UK to suggest that (i) DWV-B has rapidly expanded worldwide since its first description in 2004 and (ii) that it is potentially replacing DWV-A. Both genotypes are also found in wild bee species. Based on a simple mathematical model, we suggest that interference between viral genotypes when co-infecting the same host is key to understanding their epidemiology. We finally discuss the consequences of genotype replacement for beekeeping and for wild pollinator species.

Human population history at the crossroads of East and Southeast Asia since 11,000 years ago.

Past human genetic diversity and migration between southern China and Southeast Asia have not been well characterized, in part due to poor preservation of ancient DNA in hot and humid regions. We sequenced 31 ancient genomes from southern China (Guangxi and Fujian), including two ∼12,000- to 10,000-year-old individuals representing the oldest humans sequenced from southern China. We discovered a deeply diverged East Asian ancestry in the Guangxi region that persisted until at least 6,000 years ago. We found that ∼9,000- to 6,000-year-old Guangxi populations were a mixture of local ancestry, southern ancestry previously sampled in Fujian, and deep Asian ancestry related to Southeast Asian Hòabìnhian hunter-gatherers, showing broad admixture in the region predating the appearance of farming. Historical Guangxi populations dating to ∼1,500 to 500 years ago are closely related to Tai-Kadai and Hmong-Mien speakers. Our results show heavy interactions among three distinct ancestries at the crossroads of East and Southeast Asia.

Statistics and Patterns of Occurrence of Simple Tandem Repeats in SARS-CoV-1 and SARS-CoV-2 Genomic Data.

The data presented in this article is related to the research article entitled "Developing an ultra-efficient microsatellite discoverer to find structural differences between SARS-CoV-1 and Covid-19" [Naghibzadeh et al. 2020]. Simple tandem repeats (microsatellites, STR) are extracted and investigated across all viral families from four main viral realms. An ultra-efficient and reliable software, which is recently developed by the authors and published in the above-mentioned article, is used for extracting STRs. The analysis is done for k-mer tandem repeats where k varies from one to seven. In particular the frequency of trimer STRs is shown to be low in RNA viruses compared with DNA viruses. Special attention is paid to seven zoonotic viruses from family Coronaviridae which caused several severe human crises during last two decades including MERS, SARS 2003 and Covid-19.

Exact results for the probability and stochastic dynamics of fixation in the Wright-Fisher model.

In this work we consider fixation of an allele in a population. Fixation is key to understanding the way long-term evolutionary change occurs at the gene and molecular levels. Two basic aspects of fixation are: (i) the chance it occurs and (ii) the way the gene frequency progresses to fixation. We present exact results for both aspects of fixation for the Wright-Fisher model. We give the exact fixation probability for some different schemes of frequency-dependent selection. We also give the corresponding exact stochastic difference equation that generates frequency trajectories which ultimately fix. Exactness of the results means selection need not be weak. There are possible applications of this work to data analysis, modelling, and tests of approximations. The methodology employed illustrates that knowledge of the fixation probability, for all initial frequencies, fully characterises the dynamics of the Wright-Fisher model. The stochastic equations for fixing trajectories allow insight into the way fixation occurs. They provide the alternative picture that fixation is driven by the injection of one carrier of the fixing allele into the population each generation. The stochastic equations allow explicit calculation of some properties of fixing trajectories and their efficient simulation. The results are illustrated and tested with simulations.

Evolutionary control: Targeted change of allele frequencies in natural populations using externally directed evolution.

Random processes in biology, in particular random genetic drift, often make it difficult to predict the fate of a particular mutation in a population. Using principles of theoretical population genetics, we present a form of biological control that ensures a focal allele's frequency, at a given locus, achieves a prescribed probability distribution at a given time. This control is in the form of an additional evolutionary force that acts on a population. We provide the mathematical framework that determines the additional force. Our analysis indicates that generally the additional force depends on the frequency of the focal allele, and it may also depend on the time. We argue that translating this additional force into an externally controlled process, which has the possibility of being implemented in a number of different ways corresponding to selection, migration, mutation, or a combination of these, may provide a flexible instrument for targeted change of traits of interest in natural populations. This framework may be applied, or used as an informed form of guidance, in a variety of different biological scenarios including: yield and pesticide optimisation in crop production, biofermentation, the local regulation of human-associated natural populations, such as parasitic animals, or bacterial communities in hospitals.

Corrected simulations for one-dimensional diffusion processes with naturally occurring boundaries.

To simulate a diffusion process, a usual approach is to discretize the time in the associated stochastic differential equation. This is the approach used in the Euler method. In the present work we consider a one-dimensional diffusion process where the terms occurring, within the stochastic differential equation, prevent the process entering a region. The outcome is a naturally occurring boundary (which may be absorbing or reflecting). A complication occurs in a simulation of this situation. The term involving a random variable, within the discretized stochastic differential equation, may take a trajectory across the boundary into a "forbidden region." The naive way of dealing with this problem, which we refer to as the "standard" approach, is simply to reset the trajectory to the boundary, based on the argument that crossing the boundary actually signifies achieving the boundary. In this work we show, within the framework of the Euler method, that such resetting introduces a spurious force into the original diffusion process. This force may have a significant influence on trajectories that come close to a boundary. We propose a corrected numerical scheme, for simulating one-dimensional diffusion processes with naturally occurring boundaries. This involves correcting the standard approach, so that an exact property of the diffusion process is precisely respected. As a consequence, the proposed scheme does not introduce a spurious force into the dynamics. We present numerical test cases, based on exactly soluble one-dimensional problems with one or two boundaries, which suggest that, for a given value of the discrete time step, the proposed scheme leads to substantially more accurate results than the standard approach. Alternatively, the standard approach needs considerably more computation time to obtain a comparable level of accuracy to the proposed scheme, because the standard approach requires a significantly smaller time step.

How does reorganization energy change upon protein unfolding? Monitoring the structural perturbations in the heme cavity of cytochrome c.

In several classes of proteins the redox center provides an additional intrinsic biophysical probe that could be used to study the protein structure and function. In present report reorganization energy (lambda, as a parameter describing electron transfer properties) was used to study the protein structural changes around the heme prosthetic group in cytochrome c (cyt c). We attempted to monitor the value of this parameter upon the unfolding process of cyt c by urea, during which it was increased sigmoidally from about 0.52 to 0.82 eV for native and unfold protein, respectively. Results indicate that by structural changes in the heme site, lambda provides a complementary tool for following the unfolding process. Assuming a reversible two-state model for cyt c unfolding, Delta G(H2O), Cm and m values were determined to be 8.32+/-0.7 kcal mol(-1), 1.53+/-0.19 kcalmol(-1)M(-1) and 5.03 M, respectively.