Whole Genome Sequencing of Infectious Bursal Disease Viruses Isolated from a Californian Outbreak Unravels the Underlying Virulence Markers and Highlights Positive Selection Incidence

Outbreaks of the immunosuppressive infectious bursal disease (IBD) are frequently reported worldwide, despite the vaccination regimes. A 2009 Californian IBD outbreak caused by rA and rB isolates was described as very virulent (vv) IBD virus (IBDV); however, molecular factors beyond this virulence were not fully uncovered. Therefore, segments of both isolates were amplified, successfully cloned, whole genome sequenced by Next Generation Sequencing, genotyped, and the leading virulence factors were entirely investigated in terms of phylogenetic and amino acid analysis and protein modeling for positive selection orientation and interaction analysis. rA and rB isolates displayed the highest amino acid identity (97.84-100%) with Genotype 3 strains. Interestingly, rA and rB contained all virulence hallmarks of hypervariable (HVR), including 222A, 242I, 249Q, 256I, 284A, 286T, 294I, 299S, and 318G, as well as the serine-rich heptapeptide sequence. Moreover, we pinpointed the A3B2 genotype of rA and rB, predominant in non-reassortants, and we highlighted the absence of recombination events. Furthermore, gene-wise phylogenetic analysis showed the entire genes of rA and rB clustered with the vvIBDVs and emphasized their share in IBDV virulence. VP5 showed a virulence marker, MLSL (amino acid sequence). VP2 encountered three significant novel mutations apart from the HVR, including G163E in rA and Y173C and V178A in rB, all residing within interacting motifs. VP4 contained 168Y, 173N, 203S, and 239D characteristic for the vv phenotype. A235V mutation was detected at the dsRNA binding domain of VP3. In VP1, the TDN triplet and the mutation (V4I) were detected, characteristic of hypervirulence occurring at the N-terminus responsible for protein priming. Although selection analysis revealed seven sites, codon 222 was the only statistically significant selection site. The VP2 modeling of rA and rB highlighted great structure fitness, with 96.14% Ramachandran favored positioning including the 222A, i.e., not influencing the structure stability. The 222A was found to be non-interface surface residue, associated with no interaction with the attachment-mediated ligand motif. Our findings provide pivotal insights into the evolution and underlying virulence factors and will assist in the development of control strategies via sequence-based continuous monitoring for the early detection of novel vv strains.

Infectious bursal disease virus in Western Europe: the rise of reassortant strains as the dominant field threat

Infectious bursal disease virus (IBDV) is a highly contagious birnavirus causing a burdensome immunosuppressive disease in chickens. IBDV features a remarkable antigenic, pathogenic and genetic heterogeneity, with significant implications on disease manifestation, control measures and diagnostic approaches. The recent proposals of comprehensive phylogenetic classification systems offered the ideal platform for large-scale molecular surveys, which are crucial to gather epidemiological data and inform control efforts. In this study, the IBDV scenario was investigated in most of Western Europe by considering the results of diagnostic activities performed internationally throughout 2021. In total, 470 bursal samples from nine different countries were analysed by RT-PCR targeting the VP2. When a field virus was identified, the VP1 was also characterized. Most of the 132 detected field viruses were highly homologous reassortants featuring a very virulent-like VP2 and a classical-like VP1 (genotype A3B1). Despite emerging recently, these reassortants were already signalled in several countries in North-Western Europe and associated with subclinical infections. Here, we report their further spread in the region, where they currently represent the dominant field threat. Two other IBDV types were found, one in Italy, where all the identified viruses clustered in a clade of the A3B1 genotype previously reported only in Russia and the Middle East, and the other in Portugal, where the recently characterized A9B1 genotype was confirmed to circulate. The obtained data suggest the recent occurrence of a major shift in the Western European epidemiological landscape of IBDV, stressing the importance of steady monitoring and sharing of information among different countries and laboratories.RESEARCH HIGHLIGHTS The IBDV scenario in Western Europe seems to have radically changed in recent years.IBDV reassortants were found to be the dominant field type in the region.Local circulation of two other IBDV types was detected in Italy and Portugal.

Genetic Diversity of Recent Infectious Bursal Disease Viruses Isolated From Vaccinated Poultry Flocks in Malaysia

Vaccination is an essential component in controlling infectious bursal disease (IBD), however, there is a lack of information on the genetic characteristics of a recent infectious bursal disease virus (IBDV) that was isolated from IBD vaccinated commercial flocks in Malaysia. The present study investigated 11 IBDV isolates that were isolated from commercial poultry farms. The isolates were detected using reverse transcription-polymerase chain reaction (RT-PCR) targeting the hypervariable region (HVR) of VP2. Based on the HVR sequences, five isolates (IBS536/2017, IBS624/2017, UPM766/2018, UPM1056/2018, and UPM1432/2019) were selected for whole-genome sequencing using the MiSeq platform. The nucleotide and amino acid (aa) sequences were compared with the previously characterized IBDV strains. Deduced aa sequences of VP2HVR revealed seven isolates with 94–99% aa identity to very virulent strains (genogroup 3), two isolates with 97–100% aa identity to variant strains (genogroup 2), and two strains with 100% identity to the vaccine strain (genogroup 1) of IBDV. The phylogenetic analysis also showed that the isolates formed clusters with the respective genogroups. The characteristic motifs 222T, 249K, 286I, and 318D are typical of the variant strain and were observed for UPM1219/2019 and UPM1432/2019. In comparison, very virulent residues such as 222A, 249Q, 286T, and 318G were found for the vvIBDV, except for the UPM1056/2018 strain with a A222T substitution. In addition, the isolate has aa substitutions such as D213N, G254D, S315T, S317R, and A321E that are not commonly found in previously reported vvIBDV strains. Unlike the other vvIBDVs characterized in this study, UPM766/2018 lacks the MLSL aa residues in VP5. The aa tripeptides 145/146/147 (TDN) of VP1 were conserved for the vvIBDV, while a different motif, NED, was observed for the Malaysian variant strain. The phylogenetic tree showed that the IBDV variant clustered with the American and Chinese variant viruses and are highly comparable to the novel Chinese variants, with 99.9% identity. Based on the sequences and phylogenetic analyses, this is the first identification of an IBDV variant being reported in Malaysia. Further research is required to determine the pathogenicity of the IBDV variant and the protective efficacy of the current IBD vaccines being used against the virus.

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.

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.