Respiratory diseases like infectious bronchitis virus (IBV) and avian metapneumovirus (aMPV) have been held accountable for major losses for poultry production. Nevertheless, scarce information was present dealing with the prevalence and molecular epidemiology of these infections in Greece and the efficacy of currently applied control strategies. To fill this gap, a specific epidemiological study was designed. A total of 106 broiler and layer farms, including 10 backyard and 96 commercial flocks, were sampled between March 2016 and May 2017, and the obtained tracheal swabs were tested for IBV and aMPV using RT-PCR based techniques followed by sequencing. For each farm, data regarding production type, flock features, clinical signs, and vaccination program were also recorded. Different associations between vaccination protocol, production type, animal category, birds density, age, presence of clinical signs, and IBV and/or aMPV infection were tested. Both IBV and aMPV field strain prevalence were proven high, approximately 20 and 30%, respectively, being the GI-19 lineage (14 out of 19; 73.6%) and B subtype (30 out of 30; 100%), the most commonly detected IBV and aMPV genetic types. Infection with IBV field strains was significantly associated with clinical sign presence (odds ratio = 8.55 [95CI = 2.17-42.90]). Remarkably, only the vaccination protocol involving a double vaccination at 1 D of age was proven protective against IBV-induced symptomatology, with the odds of developing disease being 4.14 [95CI = 1.34-14.51] times lower. No association was demonstrated between aMPV infection and clinical outbreaks or between aMPV and IBV detection, suggesting the marginal role of the former pathogen in poultry farming. Globally, the present study provides the first detailed investigation of the epidemiological scenario of 2 viruses traditionally considered of pivotal relevance in poultry farming and demonstrates that remarkable benefits could be obtained with just minor adjustments in vaccination protocols.

Aleutian mink disease virus (AMDV) is one the most relevant pathogens of domestic mink, where it can cause significant economic losses, and wild species, which are considered a threat to mink farms. Despite their relevance, many aspects of the origin, evolution, and geographic and host spreading patterns of AMDV have never been investigated on a global scale using a comprehensive biostatistical approach. The present study, benefitting from a large dataset of sequences collected worldwide and several phylodynamic-based approaches, demonstrates the ancient origin of AMDV and its broad, unconstrained circulation from the initial intercontinental spread to the massive among-country circulation, especially within Europe, combined with local persistence and evolution. Clear expansion of the viral population size occurred over time until more effective control measures started to be applied. The role of frequent changes in epidemiological niches, including different hosts, in driving the high nucleotide and amino acid evolutionary rates was also explored by comparing the strengths of selective pressures acting on different populations. The obtained results suggest that the viral passage among locations and between wild and domesticated animals poses a double threat to farm profitability and animal welfare and health, which is particularly relevant for endangered species. Therefore, further efforts must be made to limit viral circulation and to refine our knowledge of factors enhancing AMDV spread, particularly at the wild-domestic interface.

Similarly to poultry industry, coccidiosis may cause significant economic losses also in the commercial quail industry, an emerging sector undergoing uneven development around the world. Although scant and mostly dated, the available literature reports detailed morphological and morphometric features of both oocysts and sporocysts of the Eimeria species hitherto recognized in Japanese quails, i.e. E. tsunodai, E. uzura, E. bateri, and E. fluminensis. Mixed infections are very common in the field and require an accurate differential diagnosis of diverse species of coccidia, identifying the highly pathogenic ones, in particular E. tsunodai (localized in the caeca), and E. uzura (localized in both caeca and small intestine). This goal is hampered by time-consuming laboratory procedures involving highly qualified staff and facilities, and poorly compatible with routine management practices in farmed quails. A supplemental difficulty is represented by the lack of nucleotide sequences available in GenBank. To overcome these issues, copromicroscopic and molecular analyses (amplifying the 18S rRNA region, and the internal transcribed spacers regions ITS1-5.8rRNA-ITS2) were performed on oocysts populations separately isolated from pools of 12 caecal and 12 cloacal contents collected from 240 naturally infected laying Japanese quails. Data on morphological and morphometric features of 1,000 sporulated oocysts were statistically compared, demonstrating the presence of different Eimeria species colonizing the 2 intestinal tracts. This result was also confirmed by PCR and phylogenetic analysis of the 18S rRNA gene. Overall results allowed to hypothesize the presence of E. uzura in our Japanese quails. Although a certain identification at species level was not obtained, the present study demonstrates that reasonable turnaround times of monitoring procedures performed on Japanese quail farms, shedding light on the in vivo and post-mortem differential diagnosis of coccidiosis can be achieved, and provide obvious benefits in disease understanding and control.

The genetic variability of Infectious bronchitis virus (IBV) is one of the main challenges for its control, hindering not only the development of effective vaccination strategies but also its classification and, consequently, epidemiology understanding. The 624/I and Q1 genotypes, now recognized to be part of the GI-16 lineage, represent an excellent example of the practical consequences of IBV molecular epidemiology limited knowledge. In fact, being their common origin unrecognized for a long time, independent epidemiological pictures were drawn for the two genotypes. To fix this misinterpretation, the present study reconstructs the history, population dynamics and spreading patterns of GI-16 lineage as a whole using a phylodynamic approach. A collection of worldwide available hypervariable region 1 and 2 (HVR12) and 3 (HVR3) sequences of the S1 protein was analysed together with 258 HVR3 sequences obtained from samples collected in Italy (the country where this genotype was initially identified) since 1963. The results demonstrate that after its emergence at the beginning of the XX century, GI-16 was able to persist until present days in Italy. Approximately in the late 1980s, it migrated to Asia, which became the main nucleus for further spreading to Middle East, Europe and especially South America, likely through multiple introduction events. A remarkable among-country diffusion was also demonstrated in Asia and South America. Interestingly, although most of the recent Italian GI-16 strains originated from ancestral viruses detected in the same country, a couple were closely related to Chinese ones, supporting a backward viral flow from China to Italy. Besides to the specific case-study results, this work highlights the misconceptions that originate from the lack of a unified nomenclature and poor molecular epidemiology data generation and sharing. This shortcoming appears particularly relevant since the described scenario could likely be shared by many other IBV genotypes and pathogens in general.

Feline panleukopenia is a severe disease of cats caused by feline parvovirus (FPV), and marginally canine parvovirus (CPV). Despite being less rapid than CPV, FPV evolution deserves attention, especially since outbreaks of particular severity are currently reported. This apparently different virulence needs monitoring from genetic and clinical points of view. This manuscript explored FPV molecular epidemiology at both Italian and international levels and the possible association between viral phylogeny and disease severity. Sequences from clinical cases of feline panleukopenia in Italy were obtained from 2011 to 2019, and the etiological agent was characterized, distinguishing FPV from CPV. Phylogenetic and phylodynamic analyses were conducted on Italian and international sequences. Moreover, the association between the viral sequence and clinical variables was evaluated on a group of highly characterized patients. After its origin in the 1920s, FPV showed a constant population size until a more recent expansion since 2000. Few long-distance introduction events characterized FPV spreading, however, most of its evolution occurred locally. Although without a strong statistical association, several clinical variables appeared influenced by viral phylogeny, suggesting a differential virulence potentially characterizing FPV strains. These results stress the importance of the continuous study of viral evolution and its repercussions on the disease clinical aspects.

PCV2 has emerged as one of the most devastating viral infections of swine farming, causing a relevant economic impact due to direct losses and control strategies expenses. Epidemiological and experimental studies have evidenced that genetic diversity is potentially affecting the virulence of PVC2. The growing number of PCV2 complete genomes and partial sequences available at GenBank questioned the accepted PCV2 classification.

The swine sector in Africa plays an important role in local economies, contributing to poverty alleviation and community subsistence. In addition, intensive farming is progressively becoming more important in the region. Therefore, any disease affecting swine populations can have detrimental effects on local communities. Porcine Reproductive and Respiratory Syndrome (PRRS) is among the most important infectious diseases affecting swine worldwide, but information on its epidemiology in Africa is extremely limited.

Dolphin morbillivirus (DMV) has been deemed as one of the most relevant threats for fin whales (Balaenoptera physalus) being responsible for a mortality outbreak in the Mediterranean Sea in the last years. Knowledge of the complete viral genome is essential to understand any structural changes that could modify virus pathogenesis and viral tissue tropism. We report the complete DMV sequence of N, P/V/C, M, F and H genes identified from a fin whale and the comparison of primary to quaternary structure of proteins between this fin whale strain and some of those isolated during the 1990-'92 and the 2006-'08 epidemics. Some relevant substitutions were detected, particularly Asn52Ser located on F protein and Ile21Thr on N protein. Comparing mutations found in the fin whale DMV with those occurring in viral strains of other cetacean species, some of them were proven to be the result of diversifying selection, thus allowing to speculate on their role in host adaptation and on the way they could affect the interaction between the viral attachment and fusion with the target host cells.

The discovery of a globally distributed porcine circovirus (Porcine circovirus 3; PCV-3) has led to intense research activity and the production of a large amount of molecular data. Different research groups have proposed several, not always concordant, genotypes for this virus. While such categories could aid an easier interpretation of PCV-3 molecular epidemiology, any classification, to be useful in practical settings, must be univocal and of help in the understanding of underlying biological features and epidemiology. Based on these premises, the possibility of defining PCV-3 genotypes was evaluated on the broadest available dataset of PCV-3 complete genome (n = 357) and open reading frame 2 (ORF2, n = 653) sequences. Genetic distance and phylogenetic clustering were selected as the main objective criteria. Additional factors, including the number of within-cluster sequences, host and geographic clustering, concordance between different genomic regions, and analysis method were also taken in account to generate a classification that could be effectively applied in research and diagnostic settings. A maximum within-genotype genetic distance of 3% at the complete genome and 6% at the ORF2 levels, bootstrap support higher than 90%, and concordance between analysis methods allowed us to clearly define two clades which could be potentially defined as genotypes. Further subdivision was not suggested due to the absence of a meaningful association between PCV-3 and its biological/epidemiological features. Nevertheless, since one of the clades included two strains only, thus far we formally propose the definition of only one PCV-3 genotype (PCV-3a). The established criteria will allow us to automatically recognize other genotypes when more strain sequences are characterized.

Infectious bronchitis virus (IBV) represents one of the poultry industry major threats, particularly in high density producing countries. The emergence and spread of new IBV genotypes have frustrated the various disease control efforts implemented over time. Despite that, few comprehensive and large scale studies have been performed to understand the international and local spreading dynamics of this virus. In the present work, these phenomena were evaluated by implementing a Bayesian phylodynamic approach to reconstruct the epidemiological patterns and population history of the QX genotype (currently renamed GI-19 lineage), the most relevant IBV lineage of the Old-World. Our analysis, based on 807 partial S1 sequences of strains collected from 18 countries between 1993 and 2015, demonstrates that this genotype originated in China well before its first identification. After a prolonged local circulation, it started spreading to other European, Asian and Middle East countries in successive waves, which were mirrored by concomitant fluctuations in viral population size. Interestingly, the within-Europe spread was characterized by a higher estimated migration rate compared with the inter-continental one, potentially reflecting the closer geographic and economic relationships among these countries. Nevertheless, the colonization of new states by the GI-19 lineage appeared to occur mostly by single introduction events in both intra and inter-continental spread, likely because of epidemiological factor and health policy combination which seems to prevent the frequent introduction and mixing of different strains. On the other hand, the within Italy QX circulation reconstruction showed a much more intricate connection network among different locations, evidencing the difficulty in controlling IBV spread especially in highly densely poultry populated areas. The presence of several well supported epidemiological links among distantly related Italian regions testifies that animal transportation and indirect transmission routes rather than local airborne diffusion contribute to the QX success and persistence at local scale. Globally, the spreading dynamics and evolution of the QX genotype were reconstructed from its very origin to nowadays, demonstrating the need of more effective direct control measures, particularly within each country. Unfortunately, the incompleteness of available molecular epidemiology data represents an insurmountable limit which leaves many questions currently unsolved, thus highlighting the compulsoriness of a structured monitoring and data sharing system implementation.

Infectious bronchitis virus GI-23 lineage, although described approximately two decades ago in the Middle East, has recently drawn remarkable attention and is considered an "emerging" lineage due to its current spread to several other regions, including Europe. Despite the relevance, no comprehensive studies are available investigating its epidemiologic and evolutionary pattern. The present phylodynamic study was designed to fill this gap, benefitting from a collection of freely available GI-23 sequences and ad-hoc generated European ones. After a relatively ancient origin in the Middle East, likely in the first half of the previous century, GI-23 circulated largely undetected or underdiagnosed for a long time in this region, likely causing little damage, potentially because of low virulence coupled with limited development of avian industry in the considered years and regions and insufficient diagnostic activity. The following development of the poultry industry and spread to other countries led to a progressive but slow increase of viral population size between the late '90s and 2010. An increase in viral virulence could also be hypothesized. Of note, a big recombinant cluster, likely originating in the Middle East but spreading thereafter, especially to Europe through Turkey, demonstrated a much-marked increase in viral population size compared to previously circulating variants. The extensive available GI-23 sequence datasets allowed to demonstrate several potential epidemiological links among African, Asian, and European countries, not described for other IBV lineages. However, differently from previously investigated IBV lineages, its spread appears to primarily involve neighbouring countries and those with strong economic and political relationships. It could thus be speculated that frequent effective contacts among locations are necessary for efficient strain transmission. Some countries appear to play a major role as a "bridge" among less related locations, being Turkey the most relevant example. The role of vaccination in controlling the viral population was also tentatively evaluated. However, despite some evidence suggesting such an effect, the bias in sequence and data availability and the variability in the applied vaccination protocols prevent robust conclusions and warrant further investigations.

Despite the fact that vaccine resistance has been typically considered a rare phenomenon, some episodes of vaccine failure have been reported with increasing frequency in intensively-raised livestock. Infectious bronchitis virus (IBV) is a widespread avian coronavirus, whose control relies mainly on extensive vaccine administration. Unfortunately, the continuous emergence of new vaccine-immunity escaping variants prompts the development of new vaccines. In the present work, a molecular epidemiology study was performed to evaluate the potential role of homologous vaccination in driving IBV evolution. This was undertaken by assessing IBV viral RNA sequences from the ORF encoding the S1 portion of viral surface glycoprotein (S) before and after the introduction of a new live vaccine on broiler farms in northern-Italy. The results of several biostatistics analyses consistently demonstrate the presence of a higher pressure in the post-vaccination period. Natural selection was detected essentially on sites located on the protein surface, within or nearby domains involved in viral attachment or related functions. This evidence strongly supports the action of vaccine-induced immunity in conditioning viral evolution, potentially leading to the emergence of new vaccine-escape variants. The great plasticity of rapidly-evolving RNA-viruses in response to human intervention, which extends beyond the poultry industry, is demonstrated, claiming further attention due to their relevance for animal and especially human health.

Recombination in Porcine Reproductive and Respiratory Syndrome Virus (PRRSV) is a well-documented phenomenon. A high recombination frequency has been reported in experimental conditions both in vitro and in vivo, and its role in driving viral evolution has been postulated by several authors. However field evidences are rare, mainly obtained from large-scale sampling and typically represented by single sequences rather than by groups of circulating "recombinant progenies". The present work was aimed to investigate the gray area between experimental studies and large-scale epidemiological investigations. The study was performed on ORF5, ORF7 and concatenated sequences obtained in our laboratory or available in GenBank collected between 2009 and 2012 in northern Italy. Six independent recombinant strains out of 66 concatenated sequences (∼9%) were found, demonstrating a high recombination frequency respect to previous field studies but comparable to in vitro experiments. In silico analysis let speculate that this new strain displayed physicochemical features diverse enough to potentially alter its immunological properties. Taken altogether, the results of our study support previous experimental evidences that depict PRRSV to be extremely prone to recombination. The limited temporal and geographical spread of recombinant strains however states in favor of a limited fitness of the recombinant progeny compared to parental strains and the marginal role of this phenomenon in PRRSV evolution.

Variability has been one of the hallmarks of canine parvovirus type 2 (CPV-2) since its discovery, and several lineages and antigenic variants have emerged. Among these, a group of viruses commonly called Asian CPV-2c has recently been reported with increasing frequency in different regions. Currently, its global epidemiology and evolution are essentially unknown. The present work deals with this information gap by evaluating, via sequence, phylodynamic, and phylogeographic analyses, all the complete coding sequences of strains classified as Asian CPV-2c based on a combination of amino acid markers and phylogenetic analysis. After its estimated origin around 2008, this lineage circulated undetected in Asia until approximately 2012, when an expansion in viral population size and geographical distribution occurred, involving Africa, Europe, and North America. Asia was predicted to be the main nucleus of viral dispersal, leading to multiple introduction events in other continents/countries, where infection establishment, persistence, and rapid evolution occurred. Although the dog is the main host, other non-canine species were also involved, demonstrating the host plasticity of this lineage. Finally, although most of the strains showed an amino acid motif considered characteristic of this lineage, several exceptions were observed, potentially due to convergent evolution or reversion phenomena.

Canine circovirus (CanineCV) is a recently identified virus affecting both domestic and wild carnivores, including foxes, sometimes in presence of severe clinical signs. Its circulation in wild animals can thus represent a potential threat for endangered species conservation and an infection source for dogs. Nevertheless, no data were available on its circulation in the Alps region of Northern Italy. In the present study, samples collected from 186 foxes in the period 2009-2020 from Valle d'Aosta and Veneto regions were tested using a real-time PCR assay, demonstrating a viral circulation of approximatively 2-5%, depending on the considered regions. Two complete or almost complete genome sequences were obtained, highlighting that the detected strains were part of a so defined "fox only" clade, which suggests that, despite common contact opportunities, Alps foxes are not involved in frequent transmission events to domestic dogs. Such genetic isolation could be at least partially attributed to some sort of independent evolution occurred in the foxes, leading to species barrier. Additionally, CanineCV strains in foxes from Italy were unexpectedly related to those previously identified in foxes from the United Kingdom and Scandinavian area. Combining the history of fox distribution in Europe since the last glacial maximum (LGM) with the viral history allowed us to speculate a long-standing coexistence between European canine circovirus and this host, justifying the peculiar geographic distribution and evolutionary paths of the fox infecting clade.

Avian Metapneumovirus (aMPV) has been recognized as a respiratory pathogen of turkey and chickens for a long time. Recently, a crescent awareness of aMPV, especially subtype B, clinical and economic impact has risen among European researchers and veterinarians. Nevertheless, the knowledge of its epidemiology and evolution is still limited. In the present study, the broadest available collection of partial G gene sequences obtained from European aMPV-B strains was analyzed using different phylodynamic and biostatistical approaches to reconstruct the viral spreading over time and the role of different hosts on its evolution. After aMPV-B introduction, approximatively in 1985 in France, the infection spread was relatively quick, involving the Western and Mediterranean Europe until the end of the 1990s, and then spreading westwards at the beginning of the new millennium, in parallel with an increase of viral population size. In the following period, a wider mixing among aMPV-B strains detected in eastern and western countries could be observed. Most of the within-country genetic heterogeneity was ascribable to single or few introduction events, followed by local circulation. This, combined with the high evolutionary rate herein demonstrated, led to the establishment of genetically and phenotypically different clusters among countries, which could affect the efficacy of natural or vaccine-induced immunity and should be accounted for when planning control measure implementation. On the contrary, while a significant strain exchange was proven among turkey, guinea fowl and chicken, no evidence of differential selective pressures or specific amino-acid mutations was observed, suggesting that no host adaptation is occurring.

Anaplasma phagocytophilum (A. phagocytophilum) is a tick-borne pathogen causing disease in both humans and animals. Human granulocytic anaplasmosis (HGA) is an emerging disease, but despite the remarkable prevalence in European ticks and wild animals, human infection appears underdiagnosed. Several genetic variants are circulating in Europe, including the zoonotic ecotype I. This study investigated A. phagocytophilum occurrence in wild ungulates and their ectoparasites in an area where HGA has been reported. Blood samples from wild ungulates and ectoparasites were screened by biomolecular methods targeting the mps2 gene. The groEL gene was amplified and sequenced to perform genetic characterization and phylogenetic analysis. A total of 188 blood samples were collected from different wild ungulates species showing an overall prevalence of 63.8% (88.7% in wild ruminants and 3.6% in wild boars). The prevalence of A. phagocytophilum DNA in ticks (manly Ixodes ricinus), and keds collected from wild ruminants was high, reflecting the high infection rates obtained in their hosts. Among ticks collected from wild boars (Hyalomma marginatum and Dermacentor marginatus) no DNA was detected. Phylogenetic analysis demonstrated the presence of ecotype I and II. To date, this is the first Italian report of ecotype I in alpine chamois, mouflon, and wild boar species. These findings suggest their role in HGA epidemiology, and the high prevalence detected in this study highlights that this human tick-borne disease deserves further attention.

Outbreaks of inclusion body hepatitis (IBH) and adenoviral gizzard erosion have been anecdotally reported in Greece since approximately 2011. However, a relevant increase in clinical outbreaks compatible with IBH has been described since 2014. Unfortunately, with limited exceptions, only serological assays were performed, and involved strains were not properly characterized. In the present study, 35 outbreaks were investigated in the period between July 2017 and February 2018 in Greece. In addition to clinical and histopathological diagnosis, fowl adenovirus (FAdV) presence was investigated by PCR and sequencing. Thirty-four out of 35 samples tested FAdV positive. Twenty-nine (85.29%) and 5 (14.71%) strains were classified as FAdV-E and FAdV-D, respectively. Fowl adenovirus-E strains were genetically homogeneous and formed an independent cluster of Greek-only sequences, including the sole previously available sequence, suggesting the prolonged circulation of this species in Greece. On the contrary, FAdV-D strains were more heterogeneous and closely related to strains sampled in other European countries, testifying the occurrence of multiple introduction events. The evaluation of phylogenetic relationships, geographic clustering, age of infection, and origin of the broiler breeder flocks suggests that both vertical and horizontal transmission are important in FAdV epidemiology in Greece and highlights the limited efficacy of currently implemented control measures. Of note, a significantly higher mortality was observed in precociously infected flocks, likely because of the higher susceptibility of younger animals. This evidence stresses the need of preventing vertical and/or early infection to limit the economic impact of adenovirus-induced diseases.

Bovine coronavirus (BoCV) is an important pathogen of cattle, causing severe enteric disease and playing a role in the bovine respiratory disease complex. Similar to other coronaviruses, a remarkable variability characterizes both its genome and biology. Despite their potential relevance, different aspects of the evolution of BoCV remain elusive. The present study reconstructs the history and evolution of BoCV using a phylodynamic approach based on complete genome and spike protein sequences. The results demonstrate high mutation and recombination rates affecting different parts of the viral genome. In the spike gene, this variability undergoes significant selective pressures-particularly episodic pressure-located mainly on the protein surface, suggesting an immune-induced selective pressure. The occurrence of compensatory mutations was also identified. On the contrary, no strong evidence in favor of host and/or tissue tropism affecting viral evolution has been proven. The well-known plasticity is thus ascribable to the innate broad viral tropism rather than mid- or long-term adaptation. The evaluation of the geographic spreading pattern clearly evidenced two clusters: a European cluster and an American-Asian cluster. While a relatively dense and quick migration network was identified in the former, the latter was dominated by the primary role of the United States (US) as a viral exportation source. Since the viral spreading pattern strongly mirrored the cattle trade, the need for more intense monitoring and preventive measures cannot be underestimated as well as the need to enforce the vaccination of young animals before international trade, to reduce not only the clinical impact but also the transferal and mixing of BoCV strains.

The identification of a new circovirus (Porcine circovirus 3, PCV-3) has raised a remarkable concern because of some analogies with Porcine circovirus 2 (PCV-2). Preliminary results suggest an extremely recent PCV-3 emergence and high mutation rate. Retrospective studies prove its circulation at least since the early 1990s, revealing that PCV-3 could have been infecting pigs for an even longer period. Therefore, a new evaluation, based on an updated collection of PCV-3 sequences spanning more than 20 years, is performed using a phylodynamic approach. The obtained results overrule the previous PCV-3 history concept, indicating an ancient origin. These evidences are associated with an evolutionary rate far lower (10-5-10-6 substitution/site/year) than the PCV-2 one. Accordingly, the action of selective pressures on PCV-3 open reading frames (ORFs) seems to be remarkably lower compared to those acting on PCV-2, suggesting either a reduced PCV-3 plasticity or a less efficient host-induced natural selection. A complex and not-directional viral flow network is evidenced through phylogeographic analysis, indicating a long lasting circulation rather than a recent emergence followed by spreading. Being recent emergence has been ruled out, efforts should be devoted to understand whether its recent discovery is simply due to improved detection capabilities or to the breaking of a previous equilibrium.