Spatiotemporal Phylogenetic Analysis and Molecular Characterisation of Infectious Bursal Disease Viruses Based on the VP2 Hyper-Variable Region

Bayesian phylodynamic analysis reveals the evolutionary history and the dispersal patterns of citrus tristeza virus in China based on the p25 gene

BACKGROUND

Citrus tristeza virus (CTV) is one of the most serious threats to the citrus industry, and is present in both wild and cultivated citrus. The origin and dispersal patterns of CTV is still poorly understood in China.

METHODS

In this study, 524 CTV suspected citrus samples from China were collected, including 354 cultivated citrus samples and 174 wild citrus samples. Finally, 126 CTV coat protein sequences were obtained with time-stamped from 10 citrus origins in China. Bayesian phylodynamic inference were performed for CTV origin and dispersal patterns study in China.

RESULT

We found that CTV was mainly distributed in southern and coastal areas of China. The substitution rate of CTV was 4.70 × 10- 4 subs/site/year (95% credibility interval: 1.10 × 10- 4 subs/site/year ~ 9.10 × 10- 4 subs/site/year), with a slight increasing trend in CTV populations between 1990 and 2006. The CTV isolates in China shared a most common recent ancestor around 1875 (95% credibility interval: 1676.57 ~ 1961.02). The CTV in China was originated from wild citrus in Hunan and Jiangxi, and then spread from the wild citrus to cultivated citrus in the growing regions of Sichuan, Chongqing, Hubei, Fujian, Zhejiang, Guangxi and Guangdong provinces.

CONCLUSIONS

This study has proved that CTV in China was originated from wild citrus in Hunan and Jiangxi. The spatial-temporal distribution and dispersal patterns has uncovered the population and pandemic history of CTV, providing hints toward a better understanding of the spread and origin of CTV in China.

Investigation and sequence analysis of psittacine beak and feather disease virus and avian polyomavirus from companion birds in Windhoek, Namibia

The commercial farming and trading of parrots and ornamental birds as companion animals are important economic activities in many countries. Some of the bird species farmed/traded are captured from the wild or are closely related to wild birds and therefore represent a risk of pathogen exchange/introduction. Beak and feather disease virus (BFDV) and avian poliomavirus (APV) are among the viruses with the biggest impact on companion bird populations and have been detected in different hosts worldwide. Despite their relevance for both domesticated and wild birds, our knowledge of BFDV and APV epidemiology remains limited in several African countries. In the present study, 143 cloacal swabs were collected from companion birds in Windhoek, Namibia, and tested by polymerase chain reaction for BFDV and APV. Of the samples tested, 35/143 (24.48%) tested positive for BFDV; 11/143 (7.69%) were positive for APV; and 6/143 (4.2%) tested positive for both pathogens. Positive amplicons, consisting of segments of the ORF1 and VP1 genes, were sequenced and compared with sequences from viruses identified in other countries. Four Namibian-only clades of BFDV were identified, loosely related to foreign strains, which suggest the occurrence of multiple introduction events in the past, potentially from South Africa, followed by local, independent evolution. In contrast, the Namibian APV sequences were identical to each other and form a single clade. In both instances, no correlation was observed between the sampling host and the viral phylogeny, suggesting the absence of host-specific adaptation and a remarkable, unconstrained viral circulation within Namibian borders. Therefore, while regulations and control measures developed against foreign strain introduction have proven to be effective over time, the spread of BFDV and APV within Namibia's borders appears undeterred. Additional resources should be dedicated to limit strain circulation in commercial farming facilities, markets and small-scale traders.

Analysis of the global origin, evolution and transmission dynamics of the emerging novel variant IBDV (A2dB1b): The accumulation of critical aa-residue mutations and commercial trade contributes to the emergence and transmission of novel variants

The Chinese IBDV novel variant (nvIBDV), belonging to the genotype A2dB1b, an emerging pathotype that can cause subclinical disease with severe, prolonged immunosuppression, poses a new threat to the poultry industry. The process of the global origin, evolution and transmission dynamics of nvIBDV, however, is poorly understood. In this study, phylogenetic trees, site substitutions of amino acid (aa) and highly accurate protein structure modelling, selection pressure, evolutionary and transmission dynamics of nvIBDV were analysed. Interestingly, nvIBDV was classified into the same genogroup with the early US antigenic variants (avIBDV) but in a new lineage with a markedly different and specific pattern of 17 aa-residual substitutions: 13 in VP2 (77D, 213N, 221K, 222T, 249K, 252I, 253Q, 254N, 284A, 286I, 299S, 318D and 323E) and four in VP1 (141I, 163V, 240E and 508K). Importantly, the aa-residues 299S and 163V may play a key role in cell binding and polymerase activity, respectively. The effective population size of the circulating avIBDV experienced two growth phases, respectively, in the years 1999-2007 (in North America) and 2015-2021 (in Asia), which is consistent with the observed trend of the epidemic outbreaks. The most recent common ancestor (tMRCA) of avIBDV most first originated in the USA and was dated around the 1970s. After its emergence, the ancestor virus of this group probably spread to China around the 1990s and the variants experienced a long-term latent circulation with the accumulation of several critical aa-residue mutations in VP2 until re-emerging in 2016. At present, central China has become the epicentre of nvIBDV spread to other parts of China and Asian countries. Importantly, a strong correlation seems to exist between the transmission patterns of virus and the flow of commercial trade of live poultry and products. These findings provide important insights into the origin, evolution and transmission of the nvIBDV and will assist in the development of programs for control strategies for these emerging viruses.

Assessing the antigenicity of different VP3 regions of infectious bursal disease virus in chickens from South Brazil

Background

Infectious bursal disease (IBD), also known as Gumboro disease, is a viral infection that causes mortality and immunosuppression in chickens (Gallus gallus). VP2 and VP3 are the major structural viral capsid components and are the most immunogenic proteins of IBD virus (IBDV). Reliable diagnostic tests using VP2 and VP3 produced in heterologous systems are important tools to control this infection. One advantage of an IBD diagnostic based on VP3, over those that use VP2, is that VP3 has linear epitopes, enabling its production in bacteria.

Results

We tested the suitability of recombinant VP3 (rVP3) as a diagnostic reagent in an enzyme-linked immunosorbent assay (ELISA). Compared with a commercial test, rVP3 ELISA showed high sensitivity and specificity as a diagnostic tool for vaccinated animals. In addition, rVP3, but not the commercial ELISA, was able to detect antibodies in nonvaccinated chickens, probably developed against circulating IBDV strains. It was possible the assessment of VP3 regions antigenicity using chicken antisera.

Conclusions

The full-length recombinant VP3 can be used to assess post vaccination immunological status of chickens and its production is feasible and inexpensive. The evaluation of VP3 regions as candidates for general use in the diagnosis of IBD in chickens should be conducted with caution. Our work was the first to identify several regions of VP3 recognized by chicken antibodies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12917-021-02956-0.

The Novel Genetic Background of Infectious Bursal Disease Virus Strains Emerging from the Action of Positive Selection

The circulation in Europe of novel reassortant strains of infectious bursal disease virus (IBDV), containing a unique genetic background composition, represents a serious problem for animal health. Since the emergence of this novel IBDV mosaic was first described in Poland, this scenario has become particularly attractive to uncover the evolutionary forces driving the genetic diversity of IBDV populations. This study additionally addressed the phenotypic characterization of these emergent strains, as well as the main features affecting the viral fitness during the competition process of IBDV lineages in the field. Our results showed how different evolutionary mechanisms modulate the genetic diversity of co-existent IBDV lineages, leading to the error catastrophe effect, Muller ratchet effect, or prevalence, depending on their genetic compositions. We also determined that the action of the positive selection pressure, depending on the genomic segment on which it is acting, can drive two main phenotypes for IBDV: immune-escaping strains from the selection on segment A or strains with functional advantages from the selection on segment B. This last group seems to possess an increased fitness landscape in the viral quasispecies composition, presenting better adaptability to dissimilar environmental conditions and likely becoming the dominant population. The reassortant strains also exhibited a lower mortality rate compared with the well-known vvIBDV strains, which can facilitate their spreading.

Development of a recombinant VP2 vaccine for the prevention of novel variant strains of infectious bursal disease virus

Since 2017, novel variant strains of infectious bursal disease virus (nvIBDV) have been detected in China, while the current vaccines on the market against very virulent IBDV have limited protection against this subtype virus. In this context, a strain of the virus has been isolated, and sequencing alignment and bird regression experiments showed that the virus was IBDV, belonging to the nvIBDV subtype (and named IBDV FJ-1812). Furthermore, the Escherichia coli expression system was used to successfully express soluble nvIBDV rVP2, which is specifically recognized by an anti-IBDV standard serum and anti-nvIBDV positive serum, and could be assembled into 14 - 17 nm virus-like particles. Based on the purified nvIBDV rVP2, we developed an IBDV FJ-1812 VP2 VLP vaccine at a laboratory scale to evaluate protection by this vaccine; in addition, we also prepared an IBDV JZ 3/02 VP2 subunit vaccine targeting very virulent IBDV and evaluated its cross-protection against nvIBDV. Results of bird experiments showed that the nvIBDV rVP2 vaccine could induce high titres of specific antibodies, completely protect the bursa of Fabricius from viral infection, and provide 100% immune protection to SPF and Ross 308 broiler chickens. Furthermore, the IBDV JZ 3/02 VP2 subunit vaccine targeting very virulent IBDV could provide 60% protection for SPF chickens and 80% protection for Ross 308 broiler chickens. This report provides important technical supports for the prevention and control of nvIBDV in the future.

Unraveling the Global Phylodynamic and Phylogeographic Expansion of Mycoplasma gallisepticum: Understanding the Origin and Expansion of This Pathogen in Ecuador

Mycoplasma gallisepticum (MG) is among the most significant problems in the poultry industry worldwide, representing a serious threat to international trade. Despite the fact that the mgc2 gene has been widely used for diagnostic and molecular characterization purposes, there is a lack of evidence supporting the reliability of this gene as a marker for molecular epidemiology approaches. Therefore, the current study aimed to assess the accuracy of the mgc2 gene for phylogenetic, phylodynamic, and phylogeographic evaluations. Furthermore, the global phylodynamic expansion of MG is described, and the origin and extension of the outbreak caused by MG in Ecuador were tracked and characterized. The results obtained strongly supported the use of the mgc2 gene as a reliable phylogenetic marker and accurate estimator for the temporal and phylogeographic structure reconstruction of MG. The phylodynamic analysis denoted the failures in the current policies to control MG and highlighted the imperative need to implement more sensitive methodologies of diagnosis and more efficient vaccines. Framed in Ecuador, the present study provides the first piece of evidence of the circulation of virulent field MG strains in Ecuadorian commercial poultry. The findings derived from the current study provide novel and significant insights into the origin, diversification, and evolutionary process of MG globally.

Identification and assessment of pathogenicity of a naturally reassorted infectious bursal disease virus from Henan, China

Infectious bursal disease virus (IBDV) is still a vital etiological agent in poultry farms. IBDV outbreaks occasionally occur due to the presence of very virulent, reassortment or variant strains. Vaccine immunization has played crucial roles in IBD control for decades. However, survival pressure of IBDV from the vaccine immunization also increases the reassortments of circulating viruses. In this study, an IBDV strain was isolated from several broiler farms in Henan Province, central part of China, and named IBDV HN strain. Based on the results of RT-PCR, sequencing and phylogenic analyses of VP1 and VP2 genes, the IBDV HN strain is a novel reassortment strain in the Henan region. Segment A of this strain appears to originate from the very virulent IBDV strain, while segment B comes from the other field reassortment strains. This may be the result of natural reassortant of virus circulating in the field. About 60% (6/10) of experimentally infected specific pathogen-free chickens died after 3 to 5 d post-infection with typical symptom and pathological lesions. The IBDV HN strain was prone to horizontal transmission, which poses a serious threat to the chicken industry. Further investigation on the prevalence, virulence, and evolution of HN strain IBDV will provide a foundation for the prevention and control of the disease in this region.

Emergence and expansion of novel pathogenic reassortant strains of infectious bursal disease virus causing acute outbreaks of the disease in Europe

Infectious bursal disease virus (IBDV) is the aetiological agent of a highly contagious chicken immunodeficiency disorder known as Gumboro disease, which cause severe economic loses to the poultry worldwide. The emergence of very virulent IBDV strains (vvIBDV) during the late 80s resulted in drastic changes to the epidemiology of IBDV with a dramatic increase in the mortality of the animals affected. Molecular studies determined that the emergence of the vvIBDV was a consequence of a genomic reorganization of IBDV known as reassortant event by which the virus combined two emergent genetic background vvIBDV for segment A and vvIBDV for segment B. In the current study, a retrospective analysis was conducted, and samples collected during acute outbreaks of Gumboro disease in Poland during 1992-2015 were submitted to sequencing and further molecular and phylogenetic analyses. The results obtained not only revealed a high genetic diversity for Polish IBDV strains but a new population of IBDV was identified. These novel reassortant strains with a unique genetic background contain the segment A from very virulent strains and segment B from an unidentified source, phylogenetically segregated and classified as 'transitional lineage'. The results obtained also showed the presence of this new lineage in Finland, evidencing the expansion of this new genomic reorganized viral strain in Europe representing an additional threat to the global situation of IBDV.

Origin and global spreading of an ancestral lineage of the infectious bursal disease virus

Infectious bursal disease virus (IBDV) is an economically relevant and widespread pathogen that produces immunosuppression in young chickens. IBDV is genetically classified into seven genogroups (G1-G7), where the traditional classic, variant and very virulent strains correspond to G1, G2 and G3, respectively. The G4 strains, also known as 'distinct' (dIBDV), have recently acquired increased relevance because of their prevalence and notorious impair to the poultry industry in South America. Here, worldwide dIBDV strains were studied using phylogenetic and phylodynamic approaches. The phylogenetic analyses performed using partial and complete sequences of both viral segments (A and B) consistently clustered the dIBDV strains in a monophyletic group. The analyses of the VP5, polyprotein and VP1 coding regions identified amino acid residues that act as markers for the identification of the entire dIBDV group or different sub-populations. The phylodynamic analyses performed using the hypervariable region of VP2 indicated that the dIBDV strains emerged in the early 1930s in Eastern Europe, shortly after the emergence of classic strains (1927) and before variant (1949) and very virulent strains (1967). The analysis of the migration routes indicated that after its emergence, the dIBDV strains spread to Eastern Asia around 1959, to Brazil around 1963, and to Argentina around 1990. These inter-continental migrations resulted in three sub-populations that are currently represented by strains from (a) Brazil, (b) Eastern Asia and Canada, and (c) Eastern Europe, Argentina and Uruguay. Taken together, our results highlight the complex evolutionary history of IBDV and the importance of new phylodynamic data to unravel and nearly follow the different evolutionary pathways taken by this important poultry pathogen.

Phylodynamic analyses of Brazilian antigenic variants of infectious bursal disease virus

Infectious bursal disease virus (IBDV) is a very important pathogen to poultry production and it is classified into three main groups: classical virulent (cvIBDV), very virulent (vvIBDV) and antigenic variants (avIBDV). This last group is composed by five different genetic lineages (recently classified in genogroups G2, G4, G5, G6, and G7) distributed in specific regions around the world. Brazil is one of the biggest poultry producers in the world and the present study aimed to investigate the evolutionary history of avIBDVs of the genogroup G4 in Brazil. A total of 5331 IBDV positive bursa samples, from different Brazilian poultry flocks, were genotyped in a period of ten years (2005 to 2014) and 1888 (35.42%) were identified as local avIBDVs. The highly variable region of the viral protein 2 (hvvp2) gene of 28 avIBDVs was sequenced and used in phylogenetic analyses and evaluation of local amino acid signatures. In addition, all complete and partial IBDV vp2 gene sequences, with local and year of collection information available on GenBank, were retrieved. Phylogenetic analyses were carried out based on a maximum likelihood method for the classification of genogroups occurring in Brazil. Based on a Maximum Likelihood (ML) phylogenetic tree, all Brazilian avIBDVs grouped into the genogroup 4. Bayesian phylodynamics analysis demonstrated the ancestor virus of this group was probably introduced in South America in 1968 (1960 to 1974, 95% HPD) and in Brazil in 1974 (1968 to 1977, 95% HPD) and the most likely source was East Europe (Hungary or Poland). All Brazilian avIBDV sequences, as well as the other genogroup 4 sequences, showed a specific pattern of amino acid: S222, T272, P289, I290, and F296. This report brings new insights about the IBDV epidemiology in Brazil and South America.

Phylogeny and evolution of infectious bursal disease virus circulating in Turkish broiler flocks

The emergence of new infectious bursal disease virus (IBDV) variants can threaten poultry health and production all over the world causing significant economic losses. Therefore, this study was performed to determine IBDV molecular epidemilogy, VP2 gene variation, and corresponding pathological lesions in IBDV infected chickens in Turkey. For this, 1855 bursa of Fabricius samples were collected from 371 vaccinated broiler flocks. Atrophia and haemorrhages were seen in the bursa Fabricius of very virulent IBDV (vvIBDV) infected chickens. Partial VP2 gene was sequenced and phylogenetic, recombination, and evolutionary analyses were performed. 1548 (83.5%) out of 1855 of bursa of Fabricius samples were IBDV positive and 1525 of those could be sequenced. The recombination analysis did not detect occurrence of any recombination event among the Turkish strains. Among 1525 sequenced samples, 1380 of them were found to be classical strains. Among 1380 classical strains, 1317 were similar to IBDV 2512, 11 to Faragher 52/70, 40 to 228 E, and 12 to Lukert strain. Out of 1525 reverse transcriptase ploymerase chain reaction positive samples, 144 of them were found to be similar to vvIBDV-VP2 gene reported to GenBank previously. The phylogenetic tree performed on a broad sequence dataset demonstrated grouping of vvIBDV Turkish strains in three different clusters, including sequences collected also from Iraq and Kuwait (Cluster 1), Indian (Cluster 2), and a distinct Turkish-only cluster (Cluster 3). The evolutionary rate estimation on branches/clades including Turkish strain mirrored the expected one for RNA viruses and no significant differences were found among different considered branches. In conclusion, results of this study indicate that vvIBDV strains similar to those circulating in various countries in the Middle East are present and undergoing evolution in chickens from Turkish broiler flocks. This point needs to be taken into account in planning adequate control strategies.

Phylogenetic analysis of Infectious Bursal Disease viruses according to newly proposed model of classification into geno-groups

BACKGROUND

Infectious bursal disease virus (IBDV) is the causative agent of Infectious Bursal Disease (IBD), the disease causes immunosuppression which leads to secondary infections among rearing poultry flocks. Characterization of the virus is important for its control and eradication. The circulating IBDVs are classified on the basis of their antigenic and pathogenic properties. The virus is categorised as classical, variant and very virulent IBDV (vvIBDV). IBDV is a non-envelop, icosahedral double stranded virus. Viral protein 2 (VP2) is the major structural protein of capsid that determines the host-pathogen relationship. The aim of this study was to characterise the IBD virus of Pak-Asian region.

METHODOLOGY

IBDV suspected flocks were examined in Punjab, Pakistan from 2014-2018. Two hundred and fifty samples were collected with complete history of the disease. The suspected samples were collected from broiler, layer and rural poultry farms. RNA was extracted and hyper-variable region of VP2 gene was amplified using specific primers. Nucleotide sequence of the VP2 gene was determined and its Amino Acid sequence was deduced. Moreover, phylogenetic analysis was also performed.

RESULTS

The current classifications based in a hyper-variable region of the capsid protein VP2 (hvVP2), classification of IBDVs is split into newly proposed geno-groups according to Jackwood group. Among these prevailing, some IBDVs are limited geographically whereas, others are reported cosmopolitan. Genetic alterations are continuously playing role in evolution of new strains of the virus.

CONCLUSION

During this study it was found that isolates of IBDV fall in first three geno-groups. Most of the geno-groups are prevalent around the world, whereas the mutated and re-assorted ones are confined in particular areas of the globe.

Low evolutionary rate of infectious pancreatic necrosis virus (IPNV) in Italy is associated with reduced virulence in trout

Infectious pancreatic necrosis virus (IPNV) is a naked double-stranded RNA virus with a bi-segmented genome that is classified within the family Birnaviridae, genus Aquabirnavirus. IPNV was first detected in Italian trout farms in the late 1970s and ultimately became endemic. To characterize the evolution of IPNV circulating in Italy, particularly whether there is a link between evolutionary rate and virulence, we obtained and analyzed the VP1 (polymerase) and the pVP2 (major capsid protein precursor) sequences from 75 IPNV strains sampled between 1978 and 2017. These data revealed that the Italian IPNV exhibit relatively little genetic variation over the sampling period, falling into four genetic clusters within a single genogroup (group 2 for VP1 and genogroup V for pVP2) and contained one example of inter-segment reassortment. The mean evolutionary rates for VP1 and pVP2 were estimated to be 1.70 and 1.45 × 10⁻⁴ nucleotide substitutions per site, per year, respectively, and hence significantly lower than those seen in other Birnaviruses. Similarly, the relatively low ratios of non-synonymous (d_N) to synonymous (d_S) nucleotide substitutions per site in both genes indicated that IPNV was subject to strong selective constraints, again in contrast to other RNA viruses infecting salmonids that co-circulate in the same area during the same time period. Notably, all the Italian IPNV harbored a proline at position 217 (P₂₁₇) and a threonine at position 221 (T₂₂₁) in pVP2, both of which are associated with a low virulence phenotype. We therefore suggest the lower virulence of IPNV may have resulted in reduced rates of virus replication and hence lower rates of evolutionary change. The data generated here will be of importance in understanding the factors that shape the evolution of Aquabirnaviruses in nature.

Impact of RNA Degradation on Viral Diagnosis: An Understated but Essential Step for the Successful Establishment of a Diagnosis Network

The current global conditions, which include intensive globalization, climate changes, and viral evolution among other factors, have led to an increased emergence of viruses and new viral diseases; RNA viruses are key drivers of this evolution. Laboratory networks that are linked to central reference laboratories are required to conduct both active and passive environmental surveillance of this complicated global viral environment. These tasks require a continuous exchange of strains or field samples between different diagnostic laboratories. The shipment of these samples on dry ice represents both a biological hazard and a general health risk. Moreover, the requirement to ship on dry ice could be hampered by high costs, particularly in underdeveloped countries or regions located far from each other. To solve these issues, the shipment of RNA isolated from viral suspensions or directly from field samples could be a useful way to share viral genetic material. However, extracted RNA stored in aqueous solutions, even at -70 °C, is highly prone to degradation. The current study evaluated different RNA storage conditions for safety and feasibility for future use in molecular diagnostics. The in vitro RNA-transcripts obtained from an inactivated highly pathogenic avian influenza (HPAI) H5N1 virus was used as a model. The role of secondary structures in the protection of the RNA was also explored. Of the conditions evaluated, the dry pellet matrix was best able to protect viral RNA under extreme storage conditions. This method is safe, cost-effective and assures the integrity of RNA samples for reliable molecular diagnosis. This study aligns with the globally significant "Global One Health" paradigm, especially with respect to the diagnosis of emerging diseases that require confirmation by reference laboratories.

Tracking the Origin and Deciphering the Phylogenetic Relationship of Porcine Epidemic Diarrhea Virus in Ecuador

In 2010, new Chinese strains of porcine epidemic diarrhea virus (PEDV), clinically more severe than the classical strains, emerged. These strains were spread to United States in 2013 through an intercontinental transmission from China with further spreading across the world, evidencing the emergent nature of these strains. In the present study, an analysis of PEDV field sequences from Ecuador was conducted by comparing all the PEDV S gene sequences available in the GenBank database. Phylogenetic comparisons and Bayesian phylogeographic inference based on complete S gene sequences were also conducted to track the origin and putative route of PEDV. The sequence from the PED-outbreak in Ecuador was grouped into the clade II of PEDV genogroup 2a together with other sequences of isolates from Mexico, Canada, and United States. The phylogeographic study revealed the emergence of the Chinese PEDV strains, followed by spreading to US in 2013, from US to Korea, and later the introduction of PEDV to Canada, Mexico, and Ecuador directly from the US. The sources of imports of live swine in Ecuador in 2014 were mainly from Chile and US. Thus, this movement of pigs is suggested as the main way for introducing PEDV to Ecuador.

Effect of different vaccination strategies on IBV QX population dynamics and clinical outbreaks

•

Two large scale molecular epidemiology studies performed in Italy and Spain.

•

The population dynamics trends differed between the two countries.

•

Relevant association between changes in vaccination strategies and viral population.

•

Strong association between viral population dynamics and outbreak frequency.

•

Hill’s criteria for causation were accomplished.

The extreme variability and rapid evolution of Infectious bronchitis virus (IBV) has always represented the key challenge for its control because of the limited cross-protection among different strains. Several experimental trials have proven a broadening of the protection spectrum when animals are vaccinated with multiple genotypes. Nevertheless, the conditions of vaccine administration in field are so different that the generalization of experimental results is, at least, questionable. In the present study a large scale epidemiological-phylodynamic approach was used to reconstruct the demographic history of the major field genotype (i.e. the QX one) circulating in Italy and Spain. These two countries were selected because, even if they share a comparable epidemiological scenario, the implemented vaccination protocols did not vary in Spain while changed dramatically in Italy over the time period considered. One hundred and ninety-five Italian and 98 Spanish non-recombinant sequences of the hyper-variable region of the S1 gene obtained between 2012 and 2016 were analyzed using a serial coalescent-based approach to reconstruct viral population history over time. While the IBV QX population dynamics remained constant in Spain, a much more complex pattern was evidenced in Italy; both in terms of viral population size and clinical outbreak frequency. Remarkably, a strong association with changes in vaccination strategies was recognized. This allowed demonstrating, by accomplishing all Hill’s criteria for causation, the cause-effect relationship between the vaccine administration/withdrawal and the variation in viral population dynamics and, above all, IBV related outbreaks. Thus, a robust confirmation about the efficacy of IBV vaccination in field conditions was provided. Additionally, the history herein reported testifies the primary importance of rigorously planning not only the intervention strategies but also their monitoring and evaluation.

Evaluation of a Phylogenetic Marker Based on Genomic Segment B of Infectious Bursal Disease Virus: Facilitating a Feasible Incorporation of this Segment to the Molecular Epidemiology Studies for this Viral Agent

Background

Infectious bursal disease (IBD) is a highly contagious and acute viral disease, which has caused high mortality rates in birds and considerable economic losses in different parts of the world for more than two decades and it still represents a considerable threat to poultry. The current study was designed to rigorously measure the reliability of a phylogenetic marker included into segment B. This marker can facilitate molecular epidemiology studies, incorporating this segment of the viral genome, to better explain the links between emergence, spreading and maintenance of the very virulent IBD virus (vvIBDV) strains worldwide.

Methodology/Principal Findings

Sequences of the segment B gene from IBDV strains isolated from diverse geographic locations were obtained from the GenBank Database; Cuban sequences were obtained in the current work. A phylogenetic marker named B-marker was assessed by different phylogenetic principles such as saturation of substitution, phylogenetic noise and high consistency. This last parameter is based on the ability of B-marker to reconstruct the same topology as the complete segment B of the viral genome. From the results obtained from B-marker, demographic history for both main lineages of IBDV regarding segment B was performed by Bayesian skyline plot analysis. Phylogenetic analysis for both segments of IBDV genome was also performed, revealing the presence of a natural reassortant strain with segment A from vvIBDV strains and segment B from non-vvIBDV strains within Cuban IBDV population.

Conclusions/Significance

This study contributes to a better understanding of the emergence of vvIBDV strains, describing molecular epidemiology of IBDV using the state-of-the-art methodology concerning phylogenetic reconstruction. This study also revealed the presence of a novel natural reassorted strain as possible manifest of change in the genetic structure and stability of the vvIBDV strains. Therefore, it highlights the need to obtain information about both genome segments of IBDV for molecular epidemiology studies.

Molecular epidemiology study of swine influenza virus revealing a reassorted virus H1N1 in swine farms in Cuba

In this report, we describe the emergence of reassorted H1N1 swine influenza virus, originated from a reassortment event between the H1N1 pandemic influenza virus (H1N1p/2009) and endemic swine influenza virus in Cuban swine population. In November 2010, a clinical respiratory outbreak was reported on a pig fattening farm in Cuba. Phylogenetic analysis showed that all the genes of one of the isolate obtained, with the exception of neuraminidase, belonged to the H1N1p/2009 cluster. This finding suggests that H1N1pdm has been established in swine and has become a reservoir of reassortment that may produce new viruses with both animal and public health risks.

Phylodynamic analysis of porcine reproductive and respiratory syndrome virus (PRRSV) in Italy: action of selective pressures and interactions between different clades

Porcine reproductive and respiratory syndrome (PRRS) is the most relevant and challenging infectious disease to affect swine breeding. Despite this, several aspects of the virus' evolution and virus-host interaction are still poorly understood and largely based on knowledge obtained through in vitro or in vivo experimental infections. Due to peculiar experimental conditions, our understanding is often contradictory and difficult to infer with respect to actual field conditions. Our phylodynamic study, based on ORF5 sequences of 141 samples collected in Italy from 1993 to 2012, explores different aspects of PRRSV epidemiology, evolution, and virus-host interaction. Two major clades, belonging to Type 1 subtype 1, were demonstrated to co-circulate while harboring a relevant intra- and inter-clade genetic diversity. Most Recent Common Ancestor (MRCA), evolution rates, and population dynamics were estimated using a serial coalescent-based approach, and different demographic histories were reconstructed for the two clades. Analysis of selective pressure revealed that sites subjected to diversifying selection were mainly located in the region of glycoprotein 5 (GP5) exposed to the host environment. Similarly, the vast majority of strains were highly glycosylated, confirming the proposed protective role of the glycan shield against the humoral immune response. Overall, our study reports both interactions among the viral populations as well as between virus and host, and their relevance in shaping viral evolution: different population dynamics over time seem to reflect a competition between clades. Some evidence argues in favor of the role of immune pressure in affecting GP5 evolution, including frequent changes in the region exposed to the host immune response, and preserving glycosylation profiles that can hamper humoral immunity.