Punctuated Evolution of Myxoma Virus: Rapid and Disjunct Evolution of a Recent Viral Lineage in Australia.

Myxoma virus (MYXV) has been evolving in a novel host species-European rabbits-in Australia since 1950. Previous studies of viruses sampled from 1950 to 1999 revealed a remarkably clock-like evolutionary process across all Australian lineages of MYXV. Through an analysis of 49 newly generated MYXV genome sequences isolated in Australia between 2008 and 2017, we show that MYXV evolution in Australia can be characterized by three lineages, one of which exhibited a greatly elevated rate of evolutionary change and a dramatic breakdown of temporal structure. Phylogenetic analysis revealed that this apparently punctuated evolutionary event occurred between 1996 and 2012. The branch leading to the rapidly evolving lineage contained a relatively high number of nonsynonymous substitutions, and viruses in this lineage reversed a mutation found in the progenitor standard laboratory strain (SLS) and all previous sequences that disrupts the reading frame of the M005L/R gene. Analysis of genes encoding proteins involved in DNA synthesis or RNA transcription did not reveal any mutations likely to cause rapid evolution. Although there was some evidence for recombination across the MYXV phylogeny, this was not associated with the increase in the evolutionary rate. The period from 1996 to 2012 saw significant declines in wild rabbit numbers, due to the introduction of rabbit hemorrhagic disease and prolonged drought in southeastern Australia, followed by the partial recovery of populations. It is therefore possible that a rapidly changing environment for virus transmission changed the selection pressures faced by MYXV, altering the course and pace of virus evolution.IMPORTANCE The coevolution of myxoma virus (MYXV) and European rabbits in Australia is one of the most important natural experiments in evolutionary biology, providing insights into virus adaptation to new hosts and the evolution of virulence. Previous studies of MYXV evolution have also shown that the virus evolves both relatively rapidly and in a strongly clock-like manner. Using newly acquired MYXV genome sequences from Australia, we show that the virus has experienced a dramatic change in evolutionary behavior over the last 20 years, with a breakdown in clock-like structure, the appearance of a rapidly evolving virus lineage, and the accumulation of multiple nonsynonymous and indel mutations. We suggest that this punctuated evolutionary event may reflect a change in selection pressures as rabbit numbers declined following the introduction of rabbit hemorrhagic disease virus and drought in the geographic regions inhabited by rabbits.

RNA-Seq analysis of chikungunya virus infection and identification of granzyme A as a major promoter of arthritic inflammation.

Chikungunya virus (CHIKV) is an arthritogenic alphavirus causing epidemics of acute and chronic arthritic disease. Herein we describe a comprehensive RNA-Seq analysis of feet and lymph nodes at peak viraemia (day 2 post infection), acute arthritis (day 7) and chronic disease (day 30) in the CHIKV adult wild-type mouse model. Genes previously shown to be up-regulated in CHIKV patients were also up-regulated in the mouse model. CHIKV sequence information was also obtained with up to ≈8% of the reads mapping to the viral genome; however, no adaptive viral genome changes were apparent. Although day 2, 7 and 30 represent distinct stages of infection and disease, there was a pronounced overlap in up-regulated host genes and pathways. Type I interferon response genes (IRGs) represented up to ≈50% of up-regulated genes, even after loss of type I interferon induction on days 7 and 30. Bioinformatic analyses suggested a number of interferon response factors were primarily responsible for maintaining type I IRG induction. A group of genes prominent in the RNA-Seq analysis and hitherto unexplored in viral arthropathies were granzymes A, B and K. Granzyme A-/- and to a lesser extent granzyme K-/-, but not granzyme B-/-, mice showed a pronounced reduction in foot swelling and arthritis, with analysis of granzyme A-/- mice showing no reductions in viral loads but reduced NK and T cell infiltrates post CHIKV infection. Treatment with Serpinb6b, a granzyme A inhibitor, also reduced arthritic inflammation in wild-type mice. In non-human primates circulating granzyme A levels were elevated after CHIKV infection, with the increase correlating with viral load. Elevated granzyme A levels were also seen in a small cohort of human CHIKV patients. Taken together these results suggest granzyme A is an important driver of arthritic inflammation and a potential target for therapy. TRIAL REGISTRATION: ClinicalTrials.gov NCT00281294.

Virological factors that increase the transmissibility of emerging human viruses.

The early detection of pathogens with epidemic potential is of major importance to public health. Most emerging infections result in dead-end "spillover" events in which a pathogen is transmitted from an animal reservoir to a human but is unable to achieve the sustained human-to-human transmission necessary for a full-blown epidemic. It is therefore critical to determine why only some virus infections are efficiently transmitted among humans whereas others are not. We sought to determine which biological features best characterized those viruses that have achieved sustained human transmission. Accordingly, we compiled a database of 203 RNA and DNA human viruses and used an information theoretic approach to assess which of a set of key biological variables were the best predictors of human-to-human transmission. The variables analyzed were as follows: taxonomic classification; genome length, type, and segmentation; the presence or absence of an outer envelope; recombination frequency; duration of infection; host mortality; and whether or not a virus exhibits vector-borne transmission. This comparative analysis revealed multiple strong associations. In particular, we determined that viruses with low host mortality, that establish long-term chronic infections, and that are nonsegmented, nonenveloped, and, most importantly, not transmitted by vectors were more likely to be transmissible among humans. In contrast, variables including genome length, genome type, and recombination frequency had little predictive power. In sum, we have identified multiple biological features that seemingly determine the likelihood of interhuman viral transmissibility, in turn enabling general predictions of whether viruses of a particular type will successfully emerge in human populations.

Complete genome of Aedes aegypti anphevirus in the Aag2 mosquito cell line.

A novel negative-sense RNA virus, Aedes aegypti anphevirus, was recently identified in wild Aedes aegypti mosquitoes. We show that this virus is also present in the Aag2 Aedes aegypti cell line and characterize its complete genome and evolutionary history. The Aedes aegypti anphevirus genome is estimated to be 12 916 nucleotides in length, contains four genes and has a genome structure similar to that of other anpheviruses. Phylogenetically, Aedes aegypti anphevirus falls within an unclassified group of insect-specific viruses in the order Mononegavirales that form a sister-group to the chuviruses. Notably, the Aag2 cell line used here was also experimentally infected with dengue virus and naturally contained a Phasi Charoen-like virus and cell-fusing agent virus. All four viruses were at relatively high abundance, with 0.5 % of sequence reads assigned to Aedes aegypti anphevirus. The Aag2 cell line is therefore permissive to efficient co-infection with dengue virus and multiple insect-specific viruses.

Dinucleotide Composition in Animal RNA Viruses Is Shaped More by Virus Family than by Host Species.

Viruses use the cellular machinery of their hosts for replication. It has therefore been proposed that the nucleotide and dinucleotide compositions of viruses should match those of their host species. If this is upheld, it may then be possible to use dinucleotide composition to predict the true host species of viruses sampled in metagenomic surveys. However, it is also clear that different taxonomic groups of viruses tend to have distinctive patterns of dinucleotide composition that may be independent of host species. To determine the relative strength of the effect of host versus virus family in shaping dinucleotide composition, we performed a comparative analysis of 20 RNA virus families from 15 host groupings, spanning two animal phyla and more than 900 virus species. In particular, we determined the odds ratios for the 16 possible dinucleotides and performed a discriminant analysis to evaluate the capability of virus dinucleotide composition to predict the correct virus family or host taxon from which it was isolated. Notably, while 81% of the data analyzed here were predicted to the correct virus family, only 62% of these data were predicted to their correct subphylum/class host and a mere 32% to their correct mammalian order. Similarly, dinucleotide composition has a weak predictive power for different hosts within individual virus families. We therefore conclude that dinucleotide composition is generally uniform within a virus family but less well reflects that of its host species. This has obvious implications for attempts to accurately predict host species from virus genome sequences alone.IMPORTANCE Determining the processes that shape virus genomes is central to understanding virus evolution and emergence. One question of particular importance is why nucleotide and dinucleotide frequencies differ so markedly between viruses. In particular, it is currently unclear whether host species or virus family has the biggest impact on dinucleotide frequencies and whether dinucleotide composition can be used to accurately predict host species. Using a comparative analysis, we show that dinucleotide composition has a strong phylogenetic association across different RNA virus families, such that dinucleotide composition can predict the family from which a virus sequence has been isolated. Conversely, dinucleotide composition has a poorer predictive power for the different host species within a virus family and across different virus families, indicating that the host has a relatively small impact on the dinucleotide composition of a virus genome.

Substitution Model Adequacy and Assessing the Reliability of Estimates of Virus Evolutionary Rates and Time Scales.

Determining the time scale of virus evolution is central to understanding their origins and emergence. The phylogenetic methods commonly used for this purpose can be misleading if the substitution model makes incorrect assumptions about the data. Empirical studies consider a pool of models and select that with the highest statistical fit. However, this does not allow the rejection of all models, even if they poorly describe the data. An alternative is to use model adequacy methods that evaluate the ability of a model to predict hypothetical future observations. This can be done by comparing the empirical data with data generated under the model in question. We conducted simulations to evaluate the sensitivity of such methods with nucleotide, amino acid, and codon data. These effectively detected underparameterized models, but failed to detect mutational saturation and some instances of nonstationary base composition, which can lead to biases in estimates of tree topology and length. To test the applicability of these methods with real data, we analyzed nucleotide and amino acid data sets from the genus Flavivirus of RNA viruses. In most cases these models were inadequate, with the exception of a data set of relatively closely related sequences of Dengue virus, for which the GTR+Γ nucleotide and LG+Γ amino acid substitution models were adequate. Our results partly explain the lack of consensus over estimates of the long-term evolutionary time scale of these viruses, and indicate that assessing the adequacy of substitution models should be routinely used to determine whether estimates are reliable.

Distinct Viral Lineages from Fish and Amphibians Reveal the Complex Evolutionary History of Hepadnaviruses.

UNLABELLED: Hepadnaviruses (hepatitis B viruses [HBVs]) are the only animal viruses that replicate their DNA by reverse transcription of an RNA intermediate. Until recently, the known host range of hepadnaviruses was limited to mammals and birds. We obtained and analyzed the first amphibian HBV genome, as well as several prototype fish HBVs, which allow the first comprehensive comparative genomic analysis of hepadnaviruses from four classes of vertebrates. Bluegill hepadnavirus (BGHBV) was characterized from in-house viral metagenomic sequencing. The African cichlid hepadnavirus (ACHBV) and the Tibetan frog hepadnavirus (TFHBV) were discovered using in silico analyses of the whole-genome shotgun and transcriptome shotgun assembly databases. Residues in the hydrophobic base of the capsid (core) proteins, designated motifs I, II, and III, are highly conserved, suggesting that structural constraints for proper capsid folding are key to capsid protein evolution. Surface proteins in all vertebrate HBVs contain similar predicted membrane topologies, characterized by three transmembrane domains. Most striking was the fact that BGHBV, ACHBV, and the previously described white sucker hepadnavirus did not form a fish-specific monophyletic group in the phylogenetic analysis of all three hepadnaviral genes. Notably, BGHBV was more closely related to the mammalian hepadnaviruses, indicating that cross-species transmission events have played a major role in viral evolution. Evidence of cross-species transmission was also observed with TFHBV. Hence, these data indicate that the evolutionary history of the hepadnaviruses is more complex than previously realized and combines both virus-host codivergence over millions of years and host species jumping. IMPORTANCE: Hepadnaviruses are responsible for significant disease in humans (hepatitis B virus) and have been reported from a diverse range of vertebrates as both exogenous and endogenous viruses. We report the full-length genome of a novel hepadnavirus from a fish and the first hepadnavirus genome from an amphibian. The novel fish hepadnavirus, sampled from bluegills, was more closely related to mammalian hepadnaviruses than to other fish viruses. This phylogenetic pattern reveals that, although hepadnaviruses have likely been associated with vertebrates for hundreds of millions of years, they have also been characterized by species jumping across wide phylogenetic distances.

Fluid Spatial Dynamics of West Nile Virus in the United States: Rapid Spread in a Permissive Host Environment.

UNLABELLED: The introduction of West Nile virus (WNV) into North America in 1999 is a classic example of viral emergence in a new environment, with its subsequent dispersion across the continent having a major impact on local bird populations. Despite the importance of this epizootic, the pattern, dynamics, and determinants of WNV spread in its natural hosts remain uncertain. In particular, it is unclear whether the virus encountered major barriers to transmission, or spread in an unconstrained manner, and if specific viral lineages were favored over others indicative of intrinsic differences in fitness. To address these key questions in WNV evolution and ecology, we sequenced the complete genomes of approximately 300 avian isolates sampled across the United States between 2001 and 2012. Phylogenetic analysis revealed a relatively star-like tree structure, indicative of explosive viral spread in the United States, although with some replacement of viral genotypes through time. These data are striking in that viral sequences exhibit relatively limited clustering according to geographic region, particularly for those viruses sampled from birds, and no strong phylogenetic association with well-sampled avian species. The genome sequence data analyzed here also contain relatively little evidence for adaptive evolution, particularly of structural proteins, suggesting that most viral lineages are of similar fitness and that WNV is well adapted to the ecology of mosquito vectors and diverse avian hosts in the United States. In sum, the molecular evolution of WNV in North America depicts a largely unfettered expansion within a permissive host and geographic population with little evidence of major adaptive barriers. IMPORTANCE: How viruses spread in new host and geographic environments is central to understanding the emergence and evolution of novel infectious diseases and for predicting their likely impact. The emergence of the vector-borne West Nile virus (WNV) in North America in 1999 represents a classic example of this process. Using approximately 300 new viral genomes sampled from wild birds, we show that WNV experienced an explosive spread with little geographical or host constraints within birds and relatively low levels of adaptive evolution. From its introduction into the state of New York, WNV spread across the United States, reaching California and Florida within 4 years, a migration that is clearly reflected in our genomic sequence data, and with a general absence of distinct geographical clusters of bird viruses. However, some geographically distinct viral lineages were found to circulate in mosquitoes, likely reflecting their limited long-distance movement compared to avian species.

Cross-validation to select Bayesian hierarchical models in phylogenetics.

BACKGROUND: Recent developments in Bayesian phylogenetic models have increased the range of inferences that can be drawn from molecular sequence data. Accordingly, model selection has become an important component of phylogenetic analysis. Methods of model selection generally consider the likelihood of the data under the model in question. In the context of Bayesian phylogenetics, the most common approach involves estimating the marginal likelihood, which is typically done by integrating the likelihood across model parameters, weighted by the prior. Although this method is accurate, it is sensitive to the presence of improper priors. We explored an alternative approach based on cross-validation that is widely used in evolutionary analysis. This involves comparing models according to their predictive performance. RESULTS: We analysed simulated data and a range of viral and bacterial data sets using a cross-validation approach to compare a variety of molecular clock and demographic models. Our results show that cross-validation can be effective in distinguishing between strict- and relaxed-clock models and in identifying demographic models that allow growth in population size over time. In most of our empirical data analyses, the model selected using cross-validation was able to match that selected using marginal-likelihood estimation. The accuracy of cross-validation appears to improve with longer sequence data, particularly when distinguishing between relaxed-clock models. CONCLUSIONS: Cross-validation is a useful method for Bayesian phylogenetic model selection. This method can be readily implemented even when considering complex models where selecting an appropriate prior for all parameters may be difficult.

Low-frequency drug-resistant HIV-1 and risk of virological failure to first-line NNRTI-based ART: a multicohort European case-control study using centralized ultrasensitive 454 pyrosequencing.

OBJECTIVES: It is still debated if pre-existing minority drug-resistant HIV-1 variants (MVs) affect the virological outcomes of first-line NNRTI-containing ART. METHODS: This Europe-wide case-control study included ART-naive subjects infected with drug-susceptible HIV-1 as revealed by population sequencing, who achieved virological suppression on first-line ART including one NNRTI. Cases experienced virological failure and controls were subjects from the same cohort whose viraemia remained suppressed at a matched time since initiation of ART. Blinded, centralized 454 pyrosequencing with parallel bioinformatic analysis in two laboratories was used to identify MVs in the 1%-25% frequency range. ORs of virological failure according to MV detection were estimated by logistic regression. RESULTS: Two hundred and sixty samples (76 cases and 184 controls), mostly subtype B (73.5%), were used for the analysis. Identical MVs were detected in the two laboratories. 31.6% of cases and 16.8% of controls harboured pre-existing MVs. Detection of at least one MV versus no MVs was associated with an increased risk of virological failure (OR = 2.75, 95% CI = 1.35-5.60, P = 0.005); similar associations were observed for at least one MV versus no NRTI MVs (OR = 2.27, 95% CI = 0.76-6.77, P = 0.140) and at least one MV versus no NNRTI MVs (OR = 2.41, 95% CI = 1.12-5.18, P = 0.024). A dose-effect relationship between virological failure and mutational load was found. CONCLUSIONS: Pre-existing MVs more than double the risk of virological failure to first-line NNRTI-based ART.

Molecular epidemiology to aid virtual elimination of HIV transmission in Australia.

The Global UNAIDS 95/95/95 targets aim to increase the percentage of persons who know their HIV status, receive antiretroviral therapy, and have achieved viral suppression. Achieving these targets requires efforts to improve the public health response to increase access to care for those living with HIV, identify those yet undiagnosed with HIV early, and increase access to prevention for those most at risk of HIV acquisition. HIV infections in Australia are among the lowest globally having recorded significant declines in new diagnoses in the last decade. However, the HIV epidemic has changed with an increasing proportion of newly diagnosed infections among those born outside Australia observed in the last five years. Thus, the current prevention efforts are not enough to achieve the UNAIDS targets and virtual elimination across all population groups. We believe both are possible by including molecular epidemiology in the public health response. Molecular epidemiology methods have been crucial in the field of HIV prevention, particularly in demonstrating the efficacy of treatment as prevention. Cluster detection using molecular epidemiology can provide opportunities for the real-time detection of new outbreaks before they grow, and cluster detection programs are now part of the public health response in the USA and Canada. Here, we review what molecular epidemiology has taught us about HIV evolution and spread. We summarize how we can use this knowledge to improve public health measures by presenting case studies from the USA and Canada. We discuss the successes and challenges of current public health programs in Australia, and how we could use cluster detection as an add-on to identify gaps in current prevention measures easier and respond quicker to growing clusters. Lastly, we raise important ethical and legal challenges that need to be addressed when HIV genotypic data is used in combination with personal data.

Novel siRNA therapeutics demonstrate multi-variant efficacy against SARS-CoV-2.

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) is a respiratory virus that causes COVID-19 disease, with an estimated global mortality of approximately 2%. While global response strategies, which are predominantly reliant on regular vaccinations, have shifted from zero COVID to living with COVID, there is a distinct lack of broad-spectrum direct acting antiviral therapies that maintain efficacy across evolving SARS-CoV-2 variants of concern. This is of most concern for immunocompromised and immunosuppressed individuals who lack robust immune responses following vaccination, and others at risk for severe COVID and long-COVID. RNA interference (RNAi) therapeutics induced by short interfering RNAs (siRNAs) offer a promising antiviral treatment option, with broad-spectrum antiviral capabilities unparalleled by current antiviral therapeutics and a high genetic barrier to antiviral escape. Here we describe novel siRNAs, targeting highly conserved regions of the SARS-CoV-1 and 2 genome of both human and animal species, with multi-variant antiviral potency against eight SARS-CoV-2 lineages - Ancestral VIC01, Alpha, Beta, Gamma, Delta, Zeta, Kappa and Omicron. Treatment with our siRNA resulted in significant protection against virus-mediated cell death in vitro, with >97% cell survival (P < 0.0001), and corresponding reductions of viral nucleocapsid RNA of up to 99.9% (P < 0.0001). When compared to antivirals; Sotrovimab and Remdesivir, the siRNAs demonstrated a more potent antiviral effect and similarly, when multiplexing siRNAs to target different viral regions simultaneously, an increased antiviral effect was observed compared to individual siRNA treatments (P < 0.0001). These results demonstrate the potential for a highly effective broad-spectrum direct acting antiviral against multiple SARS-CoV-2 variants, including variants resistant to antivirals and vaccine generated neutralizing antibodies.

EH3 (ABHD9): the first member of a new epoxide hydrolase family with high activity for fatty acid epoxides.

Epoxide hydrolases are a small superfamily of enzymes important for the detoxification of chemically reactive xenobiotic epoxides and for the processing of endogenous epoxides that act as signaling molecules. Here, we report the identification of two human epoxide hydrolases: EH3 and EH4. They share 45% sequence identity, thus representing a new family of mammalian epoxide hydrolases. Quantitative RT-PCR from mouse tissue indicates strongest EH3 expression in lung, skin, and upper gastrointestinal tract. The recombinant enzyme shows a high turnover number with 8,9-, 11,12-, and 14,15-epoxyeicosatrienoic acid (EET), as well as 9,10-epoxyoctadec-11-enoic acid (leukotoxin). It is inhibited by a subclass of N,N'-disubstituted urea derivatives, including 12-(3-adamantan-1-yl-ureido)-dodecanoic acid, 1-cyclohexyl-3-dodecylurea, and 1-(1-acetylpiperidin-4-yl)-3-(4-(trifluoromethoxy)phenyl)urea, compounds so far believed to be selective inhibitors of mammalian soluble epoxide hydrolase (sEH). Its sensitivity to this subset of sEH inhibitors may have implications on the pharmacologic profile of these compounds. This is particularly relevant because sEH is a potential drug target, and clinical trials are under way exploring the value of sEH inhibitors in the treatment of hypertension and diabetes type II.

Tailored enrichment strategy detects low abundant small noncoding RNAs in HIV-1 infected cells.

BACKGROUND: The various classes of small noncoding RNAs (sncRNAs) are important regulators of gene expression across divergent types of organisms. While a rapidly increasing number of sncRNAs has been identified over recent years, the isolation of sncRNAs of low abundance remains challenging. Virally encoded sncRNAs, particularly those of RNA viruses, can be expressed at very low levels. This is best illustrated by HIV-1 where virus encoded sncRNAs represent approximately 0.1-1.0% of all sncRNAs in HIV-1 infected cells or were found to be undetected. Thus, we applied a novel, sequence targeted enrichment strategy to capture HIV-1 derived sncRNAs in HIV-1 infected primary CD4+ T-lymphocytes and macrophages that allows a greater than 100-fold enrichment of low abundant sncRNAs. RESULTS: Eight hundred and ninety-two individual HIV-1 sncRNAs were cloned and sequenced from nine different sncRNA libraries derived from five independent experiments. These clones represent up to 90% of all sncRNA clones in the generated libraries. Two hundred and sixteen HIV-1 sncRNAs were distinguishable as unique clones. They are spread throughout the HIV-1 genome, however, forming certain clusters, and almost 10% show an antisense orientation. The length of HIV-1 sncRNAs varies between 16 and 89 nucleotides with an unexpected peak at 31 to 50 nucleotides, thus, longer than cellular microRNAs or short-interfering RNAs (siRNAs). Exemplary HIV-1 sncRNAs were also generated in cells infected with different primary HIV-1 isolates and can inhibit HIV-1 replication. CONCLUSIONS: HIV-1 infected cells generate virally encoded sncRNAs, which might play a role in the HIV-1 life cycle. Furthermore, the enormous capacity to enrich low abundance sncRNAs in a sequence specific manner highly recommends our selection strategy for any type of investigation where origin or target sequences of the sought-after sncRNAs are known.

Reemergence of an atypical bluetongue virus strain in goats, Sardinia, Italy.

Bluetongue virus (BTV) is the etiologic agent of bluetongue, a WOAH (founded as Office International des Épizooties, OIE)-notifiable economically important disease of ruminants. BTV is transmitted by Culicoides biting midges and 24 different "classical" serotypes have been reported to date. In recent years, several putative novel BTV serotypes, often referred to as "atypical" BTVs, have been documented. These are characterized by unusual biological characteristics, most notably avirulence and vector-independent transmission. Here, we describe the recurrence of such an atypical virus strain BTV-X ITL2021 detected in goats six years after its first discovery in Sardinia, Italy. Combined serological and genome analysis results clearly suggest that the two strains belong to the same BTV serotype. However, unlike the 2015 strain, BTV-X ITL2021 was successfully isolated in BSR cell-culture allowing further serological characterization. Lastly, seropositivity for BTV-X ITL2021 was detected by virus-neutralization in approximately 74% of animals tested, suggesting that this atypical BTV serotype has been circulating undetected in asymptomatic animals for years.

First influenza D virus full-genome sequence retrieved from livestock in Namibia, Africa.

Influenza D virus (IDV) was first isolated in 2011 in the USA and has since been shown to circulate in cattle, pigs, sheep, wild boar, and camels. In Africa, there is limited data on the epidemiology of IDV and, so, we investigated the presence of IDV among domestic ruminants and wild animals in Namibia by screening nasal swabs using an IDV-specific molecular assay. IDV RNA was detected in bovines (n=2), giraffes (n=2) and wildebeest (n=1). The hemagglutinin-esterase-fusion (HEF) gene from one of the bovine and the wildebeest samples was successfully sequenced as well as the full genome for the second bovine sample. Phylogenetic analysis of the HEF gene positioned the African virus variants within the D/OK lineage but with a long branch. The African variant had an amino acid diversity of 2.41% and most likely represents a distinct genotype within the lineage. Notably, the polymerase acidic protein gene (PA) was more closely related to a different lineage (D/660), indicative of potential reassortment. This is the first genetic characterization of IDV in Africa and it adds important data to our understanding of the global IDV distribution.

Epidemiological and genomic findings of the first documented Italian outbreak of SARS-CoV-2 Alpha variant of concern.

From 24 December 2020 to 8 February 2021, 163 cases of SARS-CoV-2 Alpha variant of concern (VOC) were identified in Chieti province, Abruzzo region. Epidemiological data allowed the identification of 14 epi-clusters. With one exception, all the epi-clusters were linked to the town of Guardiagrele: 149 contacts formed the network, two-thirds of which were referred to the family/friends context. Real data were then used to estimate transmission parameters. According to our method, the calculated Re(t) was higher than 2 before the 12 December 2020. Similar values were obtained from other studies considering Alpha VOC. Italian sequence data were combined with a random subset of sequences obtained from the GISAID database. Genomic analysis showed close identity between the sequences from Guardiagrele, forming one distinct clade. This would suggest one or limited unspecified viral introductions from outside to Abruzzo region in early December 2020, which led to the diffusion of Alpha VOC in Guardiagrele and in neighbouring municipalities, with very limited inter-regional mixing.

Full Genome Characterization of Respiratory Syncytial Virus Causing a Fatal Infection in an Immunocompromised Patient in Tunisia.

Human orthopneumovirus (HRSV) is a virus belonging to the Pneumovirus genus that causes lower respiratory tract infections (LRTI) in infants worldwide. In Tunisia, thousands of infants hospitalized for LRTI are found to be positive for HRSV but no whole genome sequences of HRSV strains circulating in this country are available thus far. In this study, five nasal swab samples collected at different time points from a three-month-old female baby with severe immunodeficiency that was hospitalized for acute bronchiolitis were investigated by next generation sequencing. The Tunisian sequences from this study originated from samples collected in 2021, belong to the ON1 genotype of HRSV-A, and are clustered with European sequences from 2019 and not from 2020 or 2021. This is most likely related to local region-specific transmission of different HRSV-A variants due to the COVID-19 related travel restrictions. Overall, this is the first report describing the whole genome sequence of HRSV from Tunisia. However, more sequence data is needed to better understand the genetic diversity and transmission dynamic of HRSV.

Epizootic Haemorrhagic Disease Virus Serotype 8 in Tunisia, 2021.

Epizootic haemorrhagic disease (EHD) is a Culicoides-borne viral disease caused by the epizootic haemorrhagic disease virus (EHDV) associated with clinical manifestations in domestic and wild ruminants, primarily white-tailed deer (Odocoileus virginianus) and cattle (Bos taurus). In late September 2021, EHDV was reported in cattle farms in central/western Tunisia. It rapidly spread throughout the country with more than 200 confirmed outbreaks. We applied a combination of classical and molecular techniques to characterize the causative virus as a member of the serotype EHDV-8. This is the first evidence of EHDV- 8 circulation since 1982 when the prototype EHDV-8 strain was isolated in Australia. This work highlights the urgent need for vaccines for a range of EHDV serotypes.