Sequencing and Analysis of Lumpy Skin Disease Virus Whole Genomes Reveals a New Viral Subgroup in West and Central Africa.

Lumpy skin disease virus (LSDV) is a member of the capripoxvirus (CPPV) genus of the Poxviridae family. LSDV is a rapidly emerging, high-consequence pathogen of cattle, recently spreading from Africa and the Middle East into Europe and Asia. We have sequenced the whole genome of historical LSDV isolates from the Pirbright Institute virus archive, and field isolates from recent disease outbreaks in Sri Lanka, Mongolia, Nigeria and Ethiopia. These genome sequences were compared to published genomes and classified into different subgroups. Two subgroups contained vaccine or vaccine-like samples ("Neethling-like" clade 1.1 and "Kenya-like" subgroup, clade 1.2.2). One subgroup was associated with outbreaks of LSD in the Middle East/Europe (clade 1.2.1) and a previously unreported subgroup originated from cases of LSD in west and central Africa (clade 1.2.3). Isolates were also identified that contained a mix of genes from both wildtype and vaccine samples (vaccine-like recombinants, grouped in clade 2). Whole genome sequencing and analysis of LSDV strains isolated from different regions of Africa, Europe and Asia have provided new knowledge of the drivers of LSDV emergence, and will inform future disease control strategies.

Spatiotemporal genotype replacement of H5N8 avian influenza viruses contributed to H5N1 emergence in 2021/2022 panzootic.

Since 2020, clade 2.3.4.4b highly pathogenic avian influenza H5N8 and H5N1 viruses have swept through continents, posing serious threats to the world. Through comprehensive analyses of epidemiological, genetic, and bird migration data, we found that the dominant genotype replacement of the H5N8 viruses in 2020 contributed to the H5N1 outbreak in the 2021/2022 wave. The 2020 outbreak of the H5N8 G1 genotype instead of the G0 genotype produced reassortment opportunities and led to the emergence of a new H5N1 virus with G1's HA and MP genes. Despite extensive reassortments in the 2021/2022 wave, the H5N1 virus retained the HA and MP genes, causing a significant outbreak in Europe and North America. Furtherly, through the wild bird migration flyways investigation, we found that the temporal-spatial coincidence between the outbreak of the H5N8 G1 virus and the bird autumn migration may have expanded the H5 viral spread, which may be one of the main drivers of the emergence of the 2020-2022 H5 panzootic.IMPORTANCESince 2020, highly pathogenic avian influenza (HPAI) H5 subtype variants of clade 2.3.4.4b have spread across continents, posing unprecedented threats globally. However, the factors promoting the genesis and spread of H5 HPAI viruses remain unclear. Here, we found that the spatiotemporal genotype replacement of H5N8 HPAI viruses contributed to the emergence of the H5N1 variant that caused the 2021/2022 panzootic, and the viral evolution in poultry of Egypt and surrounding area and autumn bird migration from the Russia-Kazakhstan region to Europe are important drivers of the emergence of the 2020-2022 H5 panzootic. These findings provide important targets for early warning and could help control the current and future HPAI epidemics.

West Nile virus spread in Europe: Phylogeographic pattern analysis and key drivers.

BACKGROUND: West Nile virus (WNV) outbreaks in birds, humans, and livestock have occurred in multiple areas in Europe and have had a significant impact on animal and human health. The patterns of emergence and spread of WNV in Europe are very different from those in the US and understanding these are important for guiding preparedness activities. METHODS: We mapped the evolution and spread history of WNV in Europe by incorporating viral genome sequences and epidemiological data into phylodynamic models. Spatially explicit phylogeographic models were developed to explore the possible contribution of different drivers to viral dispersal direction and velocity. A "skygrid-GLM" approach was used to identify how changes in environments would predict viral genetic diversity variations over time. FINDINGS: Among the six lineages found in Europe, WNV-2a (a sub-lineage of WNV-2) has been predominant (accounting for 73% of all sequences obtained in Europe that have been shared in the public domain) and has spread to at least 14 countries. In the past two decades, WNV-2a has evolved into two major co-circulating clusters, both originating from Central Europe, but with distinct dynamic history and transmission patterns. WNV-2a spreads at a high dispersal velocity (88km/yr-215 km/yr) which is correlated to bird movements. Notably, amongst multiple drivers that could affect the spread of WNV, factors related to land use were found to strongly influence the spread of WNV. Specifically, the intensity of agricultural activities (defined by factors related to crops and livestock production, such as coverage of cropland, pasture, cultivated and managed vegetation, livestock density) were positively associated with both spread direction and velocity. In addition, WNV spread direction was associated with high coverage of wetlands and migratory bird flyways. CONCLUSION: Our results suggest that-in addition to ecological conditions favouring bird- and mosquito- presence-agricultural land use may be a significant driver of WNV emergence and spread. Our study also identified significant gaps in data and the need to strengthen virological surveillance in countries of Central Europe from where WNV outbreaks are likely seeded. Enhanced monitoring for early detection of further dispersal could be targeted to areas with high agricultural activities and habitats of migratory birds.

Temporal Dynamics, Discovery, and Emergence of Human-Transmissible RNA Viruses.

Transmissibility, the ability to spread within host populations, is a prerequisite for a pathogen to have epidemic or pandemic potential. Here, we estimate the phylogenies of human infectivity and transmissibility using 1,408 genome sequences from 743 distinct RNA virus species/types in 59 genera. By repeating this analysis using data sets censored by virus discovery date, we explore how temporal changes in the known diversity of RNA viruses-especially recent increases in recognized nonhuman viruses-have altered these phylogenies. Over time, we find significant increases in the proportion of RNA virus genera estimated to have a nonhuman-infective ancestral state, in the fraction of distinct human virus lineages that are purely human-transmissible or strictly zoonotic (compared to mixed lineages), and in the number of human viruses with nearest relatives known not to infect humans. Our results are consistent with viruses that are capable of spreading in human populations commonly emerging from a nonhuman reservoir. This is more likely in lineages that already contain human-transmissible viruses but is rare in lineages that contain only strictly zoonotic viruses.

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.

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

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

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

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

The SARS-CoV-2 Alpha variant was associated with increased clinical severity of COVID-19 in Scotland: A genomics-based retrospective cohort analysis.

OBJECTIVES: The SARS-CoV-2 Alpha variant was associated with increased transmission relative to other variants present at the time of its emergence and several studies have shown an association between Alpha variant infection and increased hospitalisation and 28-day mortality. However, none have addressed the impact on maximum severity of illness in the general population classified by the level of respiratory support required, or death. We aimed to do this. METHODS: In this retrospective multi-centre clinical cohort sub-study of the COG-UK consortium, 1475 samples from Scottish hospitalised and community cases collected between 1st November 2020 and 30th January 2021 were sequenced. We matched sequence data to clinical outcomes as the Alpha variant became dominant in Scotland and modelled the association between Alpha variant infection and severe disease using a 4-point scale of maximum severity by 28 days: 1. no respiratory support, 2. supplemental oxygen, 3. ventilation and 4. death. RESULTS: Our cumulative generalised linear mixed model analyses found evidence (cumulative odds ratio: 1.40, 95% CI: 1.02, 1.93) of a positive association between increased clinical severity and lineage (Alpha variant versus pre-Alpha variants). CONCLUSIONS: The Alpha variant was associated with more severe clinical disease in the Scottish population than co-circulating lineages.

Effect of neonatal immunoglobulin status on the outcomes of spring-born suckler calves.

BACKGROUND: Successfully rearing beef calves to weaning and beyond determines the economic performance of a beef farm. As such, it is important to understand the factors influencing performance outcomes. METHODS: This study recorded the health events, mortality and growth rates of 674 calves born on 50 commercial beef farms in Great Britain using a postsampling questionnaire. All calves had a known postcolostral serum IgG status. RESULTS: Preweaning mortality in the study population was 1.5% (10/674 calves), while the treatment rate was 6.4% (43/674 calves). Serum IgG, calf sex and dystocia were significant predictors of whether a calf died and/or required treatment. Average daily liveweight gain was calculated for calves where weaning weights were provided (n = 513). Serum IgG and calf sex were consistent predictors of calf growth rates, while birthweight and whether the calf was born to a cow or heifer were predictive in a model where average daily liveweight gain was converted to a binary response variable using the mean average daily liveweight gain on the calf's farm of origin. LIMITATIONS: Morbidity and mortality were lower than comparable studies, potentially due to limitations in the study design. CONCLUSION: Serum IgG and calf sex were significant explanatory variables that affected beef calf average daily liveweight gain. For every 5 g/L increase in serum IgG, the odds ratio of dying and/or requiring treatment decreased by 0.86.

Multiclonal human origin and global expansion of an endemic bacterial pathogen of livestock.

Most new pathogens of humans and animals arise via switching events from distinct host species. However, our understanding of the evolutionary and ecological drivers of successful host adaptation, expansion, and dissemination are limited. Staphylococcus aureus is a major bacterial pathogen of humans and a leading cause of mastitis in dairy cows worldwide. Here we trace the evolutionary history of bovine S. aureus using a global dataset of 10,254 S. aureus genomes including 1,896 bovine isolates from 32 countries in 6 continents. We identified 7 major contemporary endemic clones of S. aureus causing bovine mastitis around the world and traced them back to 4 independent host-jump events from humans that occurred up to 2,500 y ago. Individual clones emerged and underwent clonal expansion from the mid-19th to late 20th century coinciding with the commercialization and industrialization of dairy farming, and older lineages have become globally distributed via established cattle trade links. Importantly, we identified lineage-dependent differences in the frequency of host transmission events between humans and cows in both directions revealing high risk clones threatening veterinary and human health. Finally, pangenome network analysis revealed that some bovine S. aureus lineages contained distinct sets of bovine-associated genes, consistent with multiple trajectories to host adaptation via gene acquisition. Taken together, we have dissected the evolutionary history of a major endemic pathogen of livestock providing a comprehensive temporal, geographic, and gene-level perspective of its remarkable success.

Ecological and evolutionary dynamics of multi-strain RNA viruses.

Potential interactions among co-circulating viral strains in host populations are often overlooked in the study of virus transmission. However, these interactions probably shape transmission dynamics by influencing host immune responses or altering the relative fitness among co-circulating strains. In this Review, we describe multi-strain dynamics from ecological and evolutionary perspectives, outline scales in which multi-strain dynamics occur and summarize important immunological, phylogenetic and mathematical modelling approaches used to quantify interactions among strains. We also discuss how host-pathogen interactions influence the co-circulation of pathogens. Finally, we highlight outstanding questions and knowledge gaps in the current theory and study of ecological and evolutionary dynamics of multi-strain viruses.

A novel strain of lumpy skin disease virus causes clinical disease in cattle in Hong Kong.

Lumpy skin disease virus (LSDV) is an emerging poxviral pathogen of cattle that is currently spreading throughout Asia. The disease situation is of high importance for farmers and policy makers in Asia. In October 2020, feral cattle in Hong Kong developed multi-focal cutaneous nodules consistent with lumpy skin disease (LSD). Gross and histological pathology further supported the diagnosis and samples were sent to the OIE Reference Laboratory at The Pirbright Institute for confirmatory testing. LSDV was detected using quantitative polymerase chain reaction (qPCR) and additional molecular analyses. This is the first report of LSD in Hong Kong. Whole genome sequencing (WGS) of the strain LSDV/Hong Kong/2020 and phylogenetic analysis were carried out in order to identify connections to previous outbreaks of LSD, and better understand the drivers of LSDV emergence. Analysis of the 90 core poxvirus genes revealed LSDV/Hong Kong/2020 was a novel strain most closely related to the live-attenuated Neethling vaccine strains of LSDV and more distantly related to wildtype LSDV isolates from Africa, the Middle East and Europe. Analysis of the more variable regions located towards the termini of the poxvirus genome revealed genes in LSDV/Hong Kong/2020 with different patterns of grouping when compared to previously published wildtype and vaccine strains of LSDV. This work reveals that the LSD outbreak in Hong Kong in 2020 was caused by a different strain of LSDV than the LSD epidemic in the Middle East and Europe in 2015-2018. The use of WGS is highly recommended when investigating LSDV disease outbreaks.

Adaptation, spread and transmission of SARS-CoV-2 in farmed minks and associated humans in the Netherlands.

In the first wave of the COVID-19 pandemic (April 2020), SARS-CoV-2 was detected in farmed minks and genomic sequencing was performed on mink farms and farm personnel. Here, we describe the outbreak and use sequence data with Bayesian phylodynamic methods to explore SARS-CoV-2 transmission in minks and humans on farms. High number of farm infections (68/126) in minks and farm workers (>50% of farms) were detected, with limited community spread. Three of five initial introductions of SARS-CoV-2 led to subsequent spread between mink farms until November 2020. Viruses belonging to the largest cluster acquired an amino acid substitution in the receptor binding domain of the Spike protein (position 486), evolved faster and spread longer and more widely. Movement of people and distance between farms were statistically significant predictors of virus dispersal between farms. Our study provides novel insights into SARS-CoV-2 transmission between mink farms and highlights the importance of combining genetic information with epidemiological information when investigating outbreaks at the animal-human interface.

Risk factor analysis for beef calves requiring assisted vaginal delivery in Great Britain.

BACKGROUND: Dystocia has serious consequences for both mother and offspring. This study therefore aimed to identify risk factors for dystocia in vaginally delivered spring born beef calves in Great Britain (GB). METHODS: Information on calving assistance, calf sex, birthweight, breed, twinning, dam parity and body condition score (BCS) was collected from 1131 calves across 84 GB farms. Variables were modelled against calving assistance as a binary response variable. RESULTS: Twins (Odds Ratio [OR] = 5.45), Charolais calves (OR = 3.24), calves from primiparous dams (OR = 5.75) and male calves (OR = 1.75) were at significantly increased risk of requiring calving assistance across all models. Calves born to cows classed as thin (BCS < 2.5/5) were identified in the univariate analysis and in one of the multivariate models (OR = 1.92) as having an increased likelihood of dystocia. CONCLUSIONS: Most beef herds have limited scope to manage cows on the basis of fetal gender and number. However, calf breed, dam body condition and management of primiparous dams can be manipulated to reduce the risk of dystocia and improve supervision. Poor body condition is a novel risk factor for dystocia in beef cows and worthy of further investigation.

Integrating animal movements with phylogeography to model the spread of PRRSV in the USA.

Viral sequence data coupled with phylodynamic models have become instrumental in investigating the outbreaks of human and animal diseases, and the incorporation of the hypothesized drivers of pathogen spread can enhance the interpretation from phylodynamic inference. Integrating animal movement data with phylodynamics allows us to quantify the extent to which the spatial diffusion of a pathogen is influenced by animal movements and contrast the relative importance of different types of movements in shaping pathogen distribution. We combine animal movement, spatial, and environmental data in a Bayesian phylodynamic framework to explain the spatial diffusion and evolutionary trends of a rapidly spreading sub-lineage (denoted L1A) of porcine reproductive and respiratory syndrome virus (PRRSV) Type 2 from 2014 to 2017. PRRSV is the most important endemic pathogen affecting pigs in the USA, and this particular virulent sub-lineage emerged in 2014 and continues to be the dominant lineage in the US swine industry to date. Data included 984 open reading frame 5 (ORF5) PRRSV L1A sequences obtained from two production systems in a swine-dense production region (∼85,000 mi2) in the USA between 2014 and 2017. The study area was divided into sectors for which model covariates were summarized, and animal movement data between each sector were summarized by age class (wean: 3-4 weeks; feeder: 8-25 weeks; breeding: ≥21 weeks). We implemented a discrete-space phylogeographic generalized linear model using Bayesian evolutionary analysis by sampling trees (BEAST) to infer factors associated with variability in between-sector diffusion rates of PRRSV L1A. We found that between-sector spread was enhanced by the movement of feeder pigs, spatial adjacency of sectors, and farm density in the destination sector. The PRRSV L1A strain was introduced in the study area in early 2013, and genetic diversity and effective population size peaked in 2015 before fluctuating seasonally (peaking during the summer months). Our study underscores the importance of animal movements and shows, for the first time, that the movement of feeder pigs (8-25 weeks old) shaped the spatial patterns of PRRSV spread much more strongly than the movements of other age classes of pigs. The inclusion of movement data into phylodynamic models as done in this analysis may enhance our ability to identify crucial pathways of disease spread that can be targeted to mitigate the spatial spread of infectious human and animal pathogens.

Re-emergence of H5N8 highly pathogenic avian influenza virus in wild birds, China.

In mid-November 2020, deaths of whooper swan were reported in the Yellow River Reservoir Area, China. In the present study, we describe the genetic characterizations and phylogenetic relationships of four clade 2.3.4.4b H5N8 highly avian influenza viruses (HPAIVs) identified from a sick whooper swan and environmental samples collected in the Yellow River Reservoir Area in late November 2020. They were closely related to recent H5Nx HPAIVs causing outbreaks in Eurasia in the 2020-2021 influenza season, suggesting these isolates might be imported into China via migratory birds. The newly identified H5N8 HPAIVs possessed Q226 and G228 (H3 numbering), indicating that they prefer to avian-like receptors. However, they had three mutations falling within known antigenic regions, including T144A in antigenic region A, T192I in antigenic region B, and N240D in antigenic region D. Our study highlights the risk of the rapid global spread of H5N8 HPAIVs and the necessity for continuous monitoring of avian influenza viruses in wild birds.

Phylogenetic Structure and Sequential Dominance of Sub-Lineages of PRRSV Type-2 Lineage 1 in the United States.

The genetic diversity and frequent emergence of novel genetic variants of porcine reproductive and respiratory syndrome virus type-2 (PRRSV) hinders control efforts, yet drivers of macro-evolutionary patterns of PRRSV remain poorly documented. Utilizing a comprehensive database of >20,000 orf5 sequences, our objective was to classify variants according to the phylogenetic structure of PRRSV co-circulating in the U.S., quantify evolutionary dynamics of sub-lineage emergence, and describe potential antigenic differences among sub-lineages. We subdivided the most prevalent lineage (Lineage 1, accounting for approximately 60% of available sequences) into eight sub-lineages. Bayesian coalescent SkyGrid models were used to estimate each sub-lineage's effective population size over time. We show that a new sub-lineage emerged every 1 to 4 years and that the time between emergence and peak population size was 4.5 years on average (range: 2-8 years). A pattern of sequential dominance of different sub-lineages was identified, with a new dominant sub-lineage replacing its predecessor approximately every 3 years. Consensus amino acid sequences for each sub-lineage differed in key GP5 sites related to host immunity, suggesting that sub-lineage turnover may be linked to immune-mediated competition. This has important implications for understanding drivers of genetic diversity and emergence of new PRRSV variants in the U.S.

Monitoring SARS-CoV-2 Circulation and Diversity through Community Wastewater Sequencing, the Netherlands and Belgium.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has rapidly become a major global health problem, and public health surveillance is crucial to monitor and prevent virus spread. Wastewater-based epidemiology has been proposed as an addition to disease-based surveillance because virus is shed in the feces of ≈40% of infected persons. We used next-generation sequencing of sewage samples to evaluate the diversity of SARS-CoV-2 at the community level in the Netherlands and Belgium. Phylogenetic analysis revealed the presence of the most prevalent clades (19A, 20A, and 20B) and clustering of sewage samples with clinical samples from the same region. We distinguished multiple clades within a single sewage sample by using low-frequency variant analysis. In addition, several novel mutations in the SARS-CoV-2 genome were detected. Our results illustrate how wastewater can be used to investigate the diversity of SARS-CoV-2 viruses circulating in a community and identify new outbreaks.