Complete genome analysis and time scale evolution of very virulent infectious bursal disease viruses isolated from recent outbreaks in Morocco

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.

Fulllength coding sequence analysis of genome segments A and B of infectious bursal disease virus in Thailand: identification of Chinese-like and recombinant virus in the field

RESEARCH HIGHLIGHTS

For the first time, this work demonstrated a recombinant IBDV strain in Thailand.Two genogroups of IBDV were found in Thailand: including HLJ-504-like and recombinant virus.Analysis of the full coding sequence is essential for monitoring emerging variant IBDV.

Phylogenotyping of infectious bursal disease virus in Vietnam according to the newly unified genotypic classification scheme

Since 1987, infectious bursal disease virus (IBDV) has circulated and evolved in Vietnam, but little is known about the genotypes present. IBDV samples were collected in 1987, 2001-2006, 2008, 2011, 2015-2019, and 2021 in 18 provinces. We conducted phylogenotyping analysis based on an alignment of 143 VP2-HVR (hypervariable region) sequences from 64 Vietnamese isolates (26 previous and 38 additional sequences and two vaccines, and alignment of 82 VP1 B-marker sequences, including one vaccine and four Vietnamese field strains. The analysis identified three A-genotypes, A1, A3, and A7, and two B-genotypes, B1 and B3, among the Vietnamese IBDV isolates. The lowest average evolutionary distance (8.6%) was seen between the A1 and A3 genotypes, and the highest (21.7%) was between A5 and A7, while there was a distance of 14% between B1 and B3 and 17% between B3 and B2. Unique signature residues were observed for genotypes A2, A3, A5, A6, and A8, which could be used for genotypic discrimination. A timeline statistical summary revealed that the A3-genotype predominated (79.8% presence) in Vietnam from 1987 to 2021 and that it remained the dominant IBDV genotype over the last five years (2016-2021). The current study contributes to a better understanding of the circulating genotypes and evolution of IBDV in Vietnam and worldwide.

Development of a One-Step Real-Time TaqMan Reverse Transcription Polymerase Chain Reaction (RT-PCR) Assay for the Detection of the Novel Variant Infectious Bursal Disease Virus (nVarIBDV) Circulating in China

The novel variant IBDV (nVarIBDV, especially genotype A2dB1) mainly affects broilers in China. It causes an infection characterized by the atrophy of the bursa, a decrease in the level of lymphocytes, proliferation of fibrous tissue around the follicle, and severe atrophy of the follicle in the bursa. Poultry vaccinated with live IBDV vaccines do not have the challenge present with bursa atrophy, which is misdiagnosed for nVarIBDV because of the lack of other gross clinical symptoms. The present study sought to explore the potential and reliability of the real-time TaqMan analysis method for the detection and discrimination of the nVarIBDV genotype from that of the non-nVarIBDV, especially in live vaccine strains. This method will help monitor vaccinated poultry to control and manage infection with the nVarIBDV IBDVs. The nucleotide polymorphism in the 5'-UTR region and the vp5/vp2 overlapping region of the segment A sequences of IBDV were used to establish a one-step real-time TaqMan reverse transcription polymerase chain reaction (RT-PCR) method in this study. The results showed that the method accurately distinguished the nVarIBDV and non-nVarIBDV strains (especially live vaccine strains), and there were no cross-reactions with the infectious bronchitis virus (IBV), Newcastle disease virus (NDV), avian influenza virus (AIV), infectious laryngotracheitis virus (ILTV), fowlpox virus (FPV), Mycoplasma gallisepticum (M. gallisepticum), Mycoplasma synoviae (M. synoviae), and IBDV-negative field samples. The method showed a linear dynamic range between 10² and 10⁷ DNA copies/reaction, with an average R² of 0.99 and an efficiency of 93% for nVarIBDV and an average R² of 1.00 and an efficiency of 94% for non-nVarIBDV. The method was also used for the detection of 84 clinical bursae of chickens vaccinated with the live vaccine. The results showed that this method accurately distinguished the nVarIBDV and non-nVarIBDV strains (vaccine strains), compared with a strategy based on the sequence analysis of HVRs at the vp2 gene or the reverse transcription PCR (RT-PCR) for the vp5 gene. These findings showed that this one-step real-time TaqMan RT-PCR method provides a rapid, sensitive, specific, and simple approach for detection of infections caused by nVarIBDV and is a useful clinical diagnostic tool for identifying and distinguishing nVarIBDV from non-nVarIBDV, especially live vaccine strains.

Genomic Analysis of Infectious Bursal Disease Virus in Nigeria: Identification of Unique Mutations of Yet Unknown Biological Functions in Both Segments A and B

Infectious bursal disease (IBD) is a viral poultry disease known worldwide for impacting the economy and food security. The disease is endemic in Nigeria, with reported outbreaks in vaccinated poultry flocks. To gain insight into the dynamics of infectious bursal disease virus (IBDV) evolution in Nigeria, near-complete genomes of four IBDVs were evaluated. Amino acid sequences in the hypervariable region of the VP2 revealed conserved markers (222A, 242I, 256I, 294I and 299S) associated with very virulent (vv) IBDV, including the serine-rich heptapeptide motif (SWSASGS). Based on the newly proposed classification for segments A and B, the IBDVs clustered in the A3B5 group (where A3 are IBDVs with vvIBDV-like segment A, and where B5 are from non-vvIBDV-like segment B) form a monophyletic subcluster. Unique amino acid mutations with yet-to-be-determined biological functions have been observed in both segments. Amino acid sequences of the Nigerian IBDVs showed that they are reassortant viruses. Circulation of reassortant IBDVs may be responsible for the vaccination failures observed in the Nigerian poultry population. Close monitoring of changes in the IBDV genome is recommended to nip deleterious changes in the bud through the identification and introduction of the most appropriate vaccine candidates and advocacy/extension programs for properly implementing disease control.

Epidemiology of Infectious Bursal Disease Virus in Poland during 2016-2022

Infectious bursal disease virus is an immunosuppressive ubiquitous pathogen that causes serious economic losses in poultry production. The virus is prone to genetic changes through mutations and reassortment, which drive the emergence of new variants and lead to a change in the epidemiological situation in a field. Such a situation is currently being reported due to a large wave of IBDV A3B1 reassortant infections in northwestern Europe. On the other hand, in Poland, which is the largest producer of chicken meat in the EU, the IBDVs of genotypes A3B2 and A3B4 were circulating just before the emergence of A3B1 reassortants. The purpose of the presented study was to update the IBDV epidemiological situation. The performed molecular survey based on the sequence of both genome segments showed the presence of very virulent strains (A3B2) and reassortants of genotypes A3B4 and A3B1; moreover, two of these genotypes are newly introduced IBDV lineages. In addition, a number of amino acid substitutions were demonstrated, including within antigenic epitopes and virulence determinants. In conclusion, the results obtained indicated a dynamic epidemiological situation in Poland, which highlights the need for further monitoring studies in the region and verification of protection conferred by the vaccines used against infection with detected IBDV.

Genetics and Pathogenicity of Natural Reassortant of Infectious Bursal Disease Virus Emerging in Latvia

Infectious bursal disease virus is an immunosuppressive pathogen that, despite applied vaccination, is affecting the poultry industry worldwide. This report presents the genetic and pathotypic characterization of a natural reassortant emerging in Europe (Latvia). Genetic characterization showed that strain 25/11/Latvia/2011 represents genotype A3B1, whose segment A is derived from very virulent strains, while segment B is from the classic-like genogroup. Phylogenetic maximum likelihood inference of the B-segment sequence clustered the reassortant strain together with the US antigenic variant E strain. However, the obtained full-length sequence of 25/11/Latvia/2011 revealed that not only reassortment but also dozens of mutations shaped the unique genetic makeup. Phenotypic characterization showed no mortality and no clinical signs of disease but a severe bursa of Fabricius atrophy and splenomegaly in the convalescent birds at 10 days post infection. The results obtained indicate that the acquired genetic constellation contributed to a decrease in virulence; nevertheless, the infection causes severe damage to lymphoid organs, which can lead to impaired immune responses.

Full-length genome sequencing of a very virulent infectious bursal disease virus isolated in Tunisia

Infectious bursal disease (IBD), an acute, highly contagious, and immunosuppressive avian disease, is caused by infectious bursal disease virus (IBDV) and constitutes one of the main threats to the poultry industry, worldwide. This study was performed to isolate and characterize IBDV isolates circulating in Tunisia. Eleven collected bird samples were identified using an SYBR Green–based one-step real-time reverse transcriptase polymerase chain reaction. The full-length genome sequencing of 7 of the 11 IBDV isolates has been realized. VP2 gene data showed limited sequence variations for all the 7 tested samples. The few nucleotide changes were silent and the deduced amino acid sequences were identical with the exception of a unique and characteristic nonsilent mutation (C₁₂₀₃) detected for the TN37/19 isolate, with a change of amino acid (L) to (F) at position 401. In addition, the serine-rich heptapeptide SWSASGS, characteristic of virulent IBDV, as well the amino acid residues, conserved in most very virulent IBDV (vvIBDV) strains, were detected in all the Tunisian tested isolates. Nucleotide sequences of VP5 gene revealed the presence of 5 substitutions leading to changes in the amino acid sequences of the virus. Two of these mutations were unique and characteristic of the Tunisian isolates. Besides, the alternative AUG start codon, characteristic of vvIBDV, was observed in all obtained VP5 gene sequences. The Tunisian protein sequences of VP1 showed E242 and the TDN triplet at positions 145, 146, and 147, a motif specific of vvIBDV. Phylogenetic analyses of the 5 genes confirmed the sequence alignment results and showed that the Tunisian strains are closely related to the very virulent Algerian IBDV strains.