key: cord-0693434-gptjqti7
authors: Ball, Christopher; Forrester, Anne; Herrmann, Andreas; Lemiere, Stephane; Ganapathy, Kannan
title: Comparative protective immunity provided by live vaccines of Newcastle disease virus or avian metapneumovirus when co-administered alongside classical and variant strains of infectious bronchitis virus in day-old broiler chicks
date: 2019-12-10
journal: Vaccine
DOI: 10.1016/j.vaccine.2019.09.081
sha: 03fea4a49ce80db75ae5c76eeaa176a87ce668a6
doc_id: 693434
cord_uid: gptjqti7

This study reports on the simultaneous administration of live NDV or aMPV subtype B vaccines alongside two live IBV (Massachusetts-H120 and 793B-CR88) vaccines in day-old maternal-antibody positive commercial broiler chicks. In the first experiment, chicks were divided into four groups; one unvaccinated and three groups vaccinated with live NDV VG/GA-Avinew, live H120 + CR88, or VG/GA-Avinew + H120 + CR88. In the second experiment, live aMPV subtype B vaccine was used in place of NDV. Clinical signs were monitored daily and oropharyngeal swabs were taken at regular intervals for vaccine virus detection. Blood was collected at 21 dpv for serology. 10 chicks from each group were challenged with virulent strains of M41 or QX or aMPV subtype B. For IBV, after 5 days post challenge (dpc), tracheal ciliary protection was assessed. For aMPV, clinical scores were recorded up to 10 dpc. For NDV, haemagglutination inhibition (HI) antibody titres were assayed as an indicator of protective immunity. In both experiments, ciliary protection for IBV vaccinated groups was maintained above 90%. The protection against virulent aMPV challenge was not compromised when aMPV, H120 and CR88 were co-administered. NDV HI mean titres in single and combined NDV-vaccinated groups remained above the protective titre (>3 log(2)). Both experiments demonstrated that simultaneous administration of live NDV VG/GA-Avinew or aMPV subtype B alongside H120 and CR88 vaccines does not interfere with protection conferred against NDV, IBV or aMPV.

Newcastle disease (ND), caused by an avian paramyxovirus serotype 1 (now designated as orthoavulavirus 1) [1] , and infectious bronchitis (IB) caused by an avian coronavirus [2] , are known to give rise to severe diseases in chickens, with disastrous economic losses and welfare concerns [3, 4] . ND and IB have been important and common chicken diseases worldwide since 1926 and 1936 respectively. In the last three decades, publications from both our and other laboratories have highlighted the importance of aMPV [1] , which causes a drop in egg production and quality, and swollen head syndrome in chickens [5, 6] . Infectious bronchitis virus (IBV), avian metapneumovirus (aMPV) and Newcastle Disease virus (NDV) are respiratory RNA viruses that primarily infect the tracheal epithelium of chickens [7] .

A number of countries have introduced the use of live and inactivated NDV, IBV and aMPV vaccines [8] . Simultaneous vaccination with multiple live viral vaccines is carried out in poultry for a number of reasons. These include the avoidance of inducing stress in chicks through excessive handling, and to reduce vaccination costs [9] . Traditionally, NDV and single IBV Massachusetts live vaccines are given to day-old chicks at the hatchery or in the farms, mostly by spray. However, as the disease and losses caused by variant IBVs increase [10] [11] [12] , growing numbers of producers have been including a 793B vaccine in the day-old vaccination programme. All three diseases (IBV, aMPV and NDV) are now routinely controlled in the poultry industry through the use of live attenuated vaccines [13, 14] , with long lasting antibody titres recorded for inactivated NDV vaccinations [15] .

Reports since the 1960 0 s have described how IBV vaccines can impair the efficacy of NDV vaccines [16, 17] , potentially leading to production losses [18] . Gelb vaccine in broiler chickens [7] . With the introduction of aMPV in 1990 0 s, the possibility of heterologous co-vaccination of NDV + aMPV and IBV + aMPV at day-old was investigated [19] [20] [21] [22] [23] . These publications showed that simultaneous administration of aMPV vaccine alongside NDV or IBV vaccine were compatible, with protection against virulent IBV, NDV or aMPV remaining uncompromised. However, in these studies, only a single vaccine strain of IBV (Massachusetts) was used. With the emergence of variant IBVs, inclusion of a variant IBV alongside Massachusetts for dayold vaccination has become a common practice [19, 24, 25] .

Recent work has shown that IBV vaccines can undergo genetic mutations following inoculation into the chicken host [26, 27] . As persistent mutations can alter strain characteristics [28] , it is preferable to monitor potential changes. Furthermore, it will be beneficial to establish if simultaneous vaccine administration has any role in nucleotide mutations occurring in the IBV S1 gene and subsequent amino acid changes. In contrast to IBV vaccine strains, prolonged persistence or genetic mutations in NDV or aMPV vaccine strains have not been reported in chickens.

To date, the impact of co-administration of live NDV or aMPV alongside classical and variant live IBV vaccines in IBV-aMPV-NDV maternal antibody positive birds has not been investigated. With the increased use of co-vaccination strategies in commercial broiler chicks, it is important to understand how such vaccinecombinations may impact protection against IBV, aMPV or NDV challenges. This study reports the results of two experiments where live IBV vaccines of Massachusetts and 793B were coadministered alongside a live aMPV subtype B or live Newcastle disease vaccine in commercial broiler chicks.

Chicks: Day-old commercial broiler chicks were obtained from a commercial hatchery, and raised in the University of Liverpool isolation facilities. Chicks were reared on deep litter with water and feed provided ad libitum. No in-feed or in-water antibiotics were used throughout the study. All experimental procedures were conducted according to UK legislation on the use of animals for experiments, and were granted ethical approval by the University of Liverpool ethics committee.

Vaccine strains: Commercially available live IBV Mass type strain H120 (Bioral, Boehringer Ingelheim, Lyon, France) and IBV variant 793B type strain CR88 (Gallivac Ò IB88, Boehringer Ingelheim, Lyon, France), NDV VG/GA-Avinew strain (Avinew Ò NeO, Boehringer Ingelheim, Lyon, France) and aMPV subtype B (Nemovac Ò , Boehringer Ingelheim, Lyon, France) vaccines were administered during this study. All vaccines were prepared prior to administration using sterile distilled water, to allow for 100 ml of vaccine per chick to be administered according to manufacturer's instructions [5, 19] .

Challenge strains: Virulent virus strains of IBV (M41, 10 5.0 Ciliostatic Dose (CD) 50 /ml; QX, 10 6.0 CD 50 /ml) and aMPV (subtype B; 10 4.51 CD 50 /ml) were used for challenge. All challenge strains have been propagated in our laboratory and used in several previous studies [19, 23, 29] . PCR and bacterial culturing confirmed an absence of NDV, aMPV, avian influenza virus (AIV), infectious laryngotracheitis virus (ILTV), infectious bursal disease virus (IBDV), fowl adenovirus (FAdV) and mycoplasmas in all inocula [30] [31] [32] [33] [34] [35] .

Experiment 1: Chicks were randomly divided into four groups (35 chicks per group) and kept in separate isolation units. At day-old, each chick was inoculated via oculo (50 ml) and nasal (50 ml) routes with the following; VG/GA-Avinew (Group A), Mass + 793B (Group B), VG/GA-Avinew + Mass + 793B (Group C) and sterile H 2 O (Group D). Chicks were observed daily for clinical signs [36] ; Mild = Coughing, head shaking and depression of short duration, Severe = gasping, coughing and depression, accompanied by ruffled feathers. Mortality rates and lesions following necropsy were recorded.

At 21 dpv, chicks from each group were separated into two further groups (A1, A2, B1, B2, C1, C2, D1 and D2), with 10 chicks per sub-group. Chicks from A1, B1, C1 and D1 were challenged with 100 ml of virulent M41, and chicks from A2, B2, C2 and D2 were challenged with 100 ml of virulent QX via oculo-nasal route. Chicks were observed daily for clinical signs as previously described, with mortality rates and post mortem lesions recorded.

At 0, 3, 7, 10, 14 and 21 dpv, oropharyngeal (OP) swabs were collected from five chicks in each group. At 7, 14, 21 dpv, five trachea and kidney samples were collected from each group. At 5 dpc, five trachea and kidney samples were collected and processed as previously described [26] . Blood was collected from ten chicks at 0 dpv and eight chicks at 21 dpv for the detection of IBV and NDV antibodies.

Experiment 2: Day-old broiler chicks were randomly divided into four groups and kept in separate isolation units, with 30 chicks per group. At day-old, each chick was inoculated via oculo (50 ml) and nasal (50 ml) routes with the following; aMPV subtype B (Group A), Mass + 793B (Group B), aMPV subtype B + Mass + 793B (Group C) and sterile H 2 O (Group D). Chicks were observed daily for clinical signs as described above.

At 21 days post vaccination (dpv), chicks from each group (A-D) were further split into twelve groups (A1-A3, B1-B3, C-C3, D1-D3), with five chicks per group. Chicks from A1, B1, C1 and D1 were challenged with 100 ml of virulent M41, and chicks in the A2, B2, C2 and D2 groups were challenged with 100 ml of virulent aMPV B via oculo-nasal route. Chicks from A3, B3, C3 and D3 were sham-inoculated with sterile distilled water. All birds were observed daily for clinical signs, with IBV challenge groups scored as above, and aMPV challenge groups scored as follows: 0 = No nasal exudate, 1 = Mild exudate, 2 = Unilateral excessive exudate, 3 = Bilateral excessive exudate [5] .

At 0, 3, 7, 10, 14 and 21 dpv, oropharyngeal (OP) swabs were collected from five chicks in each group. At 7, 14, 21 dpv five trachea and kidney samples were collected. At 5 dpc, five trachea and kidney samples were collected and processed as previously described [26] . Blood was collected from ten chicks at 0 dpv and eight chicks at 21 dpv for the detection of IBV and aMPV antibodies.

Evaluation of tracheal protection following M41 challenge: At 5 days post challenge (dpc), five chicks from each group were humanely killed by wing vein injection of sodium pentobarbital. Tracheas were removed from each chick and processed for ciliary and percentage protection assessment as previously described [14, 19] . Briefly, 10 rings from each bird were examined, with a maximum possible ciliary score of 40 indicating no protection (no cilia beating). The mean ciliary score for each bird was calculated and percentage protection for each group was calculated.

Measurement of antibody levels: Commercial ELISA kits [IBV (IDEXX, Maine, United States) and aMPV (BioChek, ER Reeuwijk, Netherlands)] were used according to manufacturer's instructions, with a titre result higher than 396 (IBV) or 1655 (aMPV) indicating a positive protection titre. Haemagluttination inhibition (HI) was utilised to determine IBV M41 and NDV protection status as previously described [37] .

IBV, NDV and aMPV detection by RT-PCR: Swabs were pooled and dipped as a single sample per group, per day in a sterile bijou containing 1.5 ml of Eagles serum-free minimum essential medium with glutamine, streptomycin [50 mg/ml] and penicillin [50 IU/ml], which was stored at À70°C until required.

RNA was extracted from the swab and tissue samples using the phenol chloroform method [38] and subjected to RT-PCR for IBV, NDV and aMPV [32, 39, 40] . Positive amplicons were purified and sent for bi-directional commercial Sanger sequencing (Source BioScience, Nottingham, UK).

IBV and aMPV viral load in trachea and kidney: Total RNA was extracted from each group's five trachea and kidney samples using the RNeasy Plus Mini Kit (Qiagen, Manchester, UK). Quantification of viral RNA was determined by quantitative RT-PCR (qRT-PCR) for both IBV and aMPV [41, 42] . The IBV qRT-PCR was performed using the Qiagen OneStep RT-PCR kit and 40 ng of total RNA per reaction. The aMPV assay was performed using the Qiagen Quantitect 2xRT-PCR no ROX master mix and 40 ng of total RNA per reaction. Ct values were converted to log 2 of Relative Equivalent Units (REU) of viral RNA, using previously determined standard curves for IBV and aMPV [25] .

Sequencing analysis: IBV partial-S1 gene sequences were trimmed to 320 bp, analysed in ChromasPRO v1.7.3 (http://technelysium.com.au/), and BLAST searches were conducted to confirm the isolate identity. Alignments were carried out in MEGA6 [43] , using Clustal W [44] . Single nucleotide polymorphism (SNP) and insertion and deletion (indel) detection was carried out in MEGA6, following alignment to sequenced vaccine strains. SNPs were characterised as non-synonymous (ns) if it led to an amino acid change, and synonymous (s) if it led to no amino acid changes. Insertions and deletions (indels) were defined as the insertion or deletion of a nucleotide that altered the sequence length. The d s /d ns ratio, to identify positive or purifying selection pressure for each group was calculated per day. This was done using the Nei-Gojobori method [45] where <1.0 indicated the recovered strains to be under positive selection pressure, and >1.0 indicated strains to be under purifying selection [26, 46] . Positive or purifying selection was considered significant at P < 0.05. Translated amino acid variations were also identified and variations that resulted in a change in hydrophobicity were noted according to the Kyte and Doolittle scale [47] .

Statistical analysis: Statistical analysis was conducted using one-way analysis of variance (ANOVA), followed by Turkey's test to examine the differences between pairs of means. Differences were considered to be significant when p < 0.05. Analysis was carried out using SPSS Statistics 22.

Following day-old vaccination, in both IBV-vaccinated groups, mild respiratory signs were recorded from 2 to 14 dpv. After M41 or QX challenge, groups that received no IBV vaccination had mild signs throughout the experiment. The IBV-vaccinated groups, which were M41-challenged or QX-challenged showed mild signs for a total of two days, at 3-4 dpc.

Ciliostasis results showed that all IBV-vaccinated groups (B1, B2, C1, C2) have high protection against both IBV challenge strains, ranging from 94 to 99% (C1 -94%; C2 -96.76%; B1 -99%; B2 -98.56%).

The NDV HI titres in groups that had not received an NDV vaccine were below 2 log 2 (Table 1) , whereas the groups that received NDV vaccine singly or in combination with IBV were 6.63 and 5.63 respectively. A significant increase (p < 0.05) in titre was seen for Group A (NDV vaccinated) when compared to Group C (combined vaccination) at 21 dpv.

When the same sera were assayed against IBV HI antibodies, IBV-vaccinated groups showed an increase in HI log 2 titre at 21 dpv when compared to 0 dpv (Table 1) . Furthermore, there was a significant increase (p < 0.05) in titre for both Groups A and B when compared to Group C for the IBV (4/91) antigen. There was a stronger response towards the 4/91 compared to the M41 HI antigen for all IBV vaccinated groups (Group B -6.75 and 4.00; Group C -5.88 and 5.63). The mean IBV ELISA antibody titre at day-old was 838.38 (±87), indicating the presence of maternally derived antibodies (MDA). By 21 dpv, antibody levels declined to below the detectable titre (3 9 6) in all groups (Group A -23.3, Group B -89, Group C -71.8 and Group D -27.7). Despite being below detectable levels, the two groups receiving IBV vaccines (B and C) had higher antibody titres compared to the IBVunvaccinated groups (A and D).

IBV RNA was detected at all sampling time points (3, 7, 10, 14 and 21 dpv) in Groups B and C, and was absent from the IBVunvaccinated groups (A and D). The Mass-type vaccine was detected up until 10 dpv in Group B and 14 dpv in Group C. After this time, only the 793B-type vaccine was identified in all subsequent samples in both IBV-vaccinated groups. All Mass-type identified strains had partial-S1 nucleotide similarity to the original vaccine, ranging from 92.97 to 99.74%, whereas the 793B-type had 93.72-99.47% relatedness to the CR88 vaccine strain. NDV was not detected at any time.

Total SNP counts remained low, with the exception of samples from Groups B and C at 14 dpv, which contained greater total SNP counts (n = 21) ( Table 2 ). This translated to 14 amino acid (AA) changes in Group B and 16 in Group C. All changes in 793Blike samples were from hydrophobic to hydrophilic properties, and the majority (75%) of changes in the Mass-like samples were hydrophilic to hydrophobic. The average d s /d ns ratio was 1.96 for 793B strains and 2.81 for Mass strains, indicating that the Mass genotype was significantly under purifying selection pressure (p < 0.05). No indels were present for any of the sequences analysed.

3.1.6. Molecular detection of IBV and NDV from trachea and kidney tissues Up to 21 dpv: IBV vaccinal strains were recovered from both trachea and kidney tissue in Group B at 7 and 14 dpv, whereas all sampling days were positive in Group C (Table 3) . No trachea samples were IBV-positive at 21 dpv and no Mass-type strains were recovered in the kidney past 7 dpv. We did not detect NDV in any samples. Sequence similarity to the IBV vaccine strains were between 98 and 100%, with a single 793B exception of 96% at 21 dpv. The majority of strains recovered from the kidney and trachea (84.62%) had a d s /d ns ratio of over 1.00, demonstrating that the majority of partial S1 nucleotide variations had no effect on the translated amino acid composition. The exceptions were the two 793B-like strains from Group C kidney samples at 21 dpv. This group had an average of five non-synonymous changes per strain, which translated to three hydrophobicity changes within all 13 recovered strains.

IBV M41 challenge at 5 dpc: IBV was detected in the trachea of Group D1 (Non-vaccinated), and in the trachea and kidney of Group A1 (NDV vaccinated) ( Table 4 ). IBV vaccine strains were detected in Groups B1 (Mass vaccinated À 99% similarity) and C1 (793B vaccinated À 96% similarity). The majority of variations led to changes in the translated amino acids (n = 14; average d s /d ns ratio = 0.85), however only one change in hydrophobic properties (hydrophilic to hydrophobic) was identified from a sample identified as highly similar to the virulent strain (Group A1; trachea). IBV QX challenge at 5 dpc: Virulent QX strains were detected in the trachea and kidney of Group D2 (Non-vaccinated), with 98-99% nucleotide similarity to the virulent strain. However, only tracheal samples were IBV-positive for Group A2 (NDV-vaccinated) ( Table 4 ). Vaccinal strains were detected in Groups B2 (793B; 98-99% similarity) and C2 (793B and Mass; 99% similarity). The minority of nucleotide variations caused an amino acid change (n = 29; average d s /d ns ratio = 0.97), which resulted in seven hydrophobicity changes (hydrophilic to hydrophobic = 4; hydrophobic to hydrophilic = 3). Vaccine strains were recovered from kidney samples in Groups B and C, with the majority of changes in 793B-like samples being non-synonymous. From a total of 18 amino acid changes, only six caused a change in hydrophobicity.

In the IBV vaccinated tracheal samples, viral load reduced from 7 dpv to 14 dpv (Group B = 1.49 to 0 log REU; Group C = 0.97 to 0.29 log REU), however IBV presence were significantly higher (p < 0.05) at 21 dpv for both groups (Group B = 2.92; Group C = 1.55) when compared to 7 dpv (Fig. 1) . At both 7 and 21 dpv, there was a significantly higher tracheal viral load in the IBVvaccination group when compared to the NDV + IBV vaccine group. Post challenge, high viral loads were present in the trachea (Fig. 2 For the kidney samples at 14 and 21 dpv, Groups B (IBV vaccinated) and C (NDV + IBV vaccinated) increased from 14 to 21 dpv (0.81 to 3.13 log REU and 1.15 to 1.42 log REU respectively) (Fig. 1) . Post challenge (Fig. 2) , significantly lower viral loads were present in non-IBV vaccinated groups (A1 -0.44 log REU; A2 -0.25 log REU; C1 -0.1 log REU; D1 -0.18 log REU) compared to IBV vaccinated groups (B1 -2.15 log REU; B2 -1.22 log REU; C2 -2.62 log REU). Following M41 challenge, viral load did not change in Group B (IBV vaccinated). The virulent QX strain was detected at a greater level in the control group (Group D) when compared to the M41 challenge control group. However, the viral load in both NDV vaccinated groups (A1 and A2) remained low (<0.45 log REU) following challenge.

Up to 21 dpv: All IBV vaccinated groups demonstrated mild clinical signs (prominent snicking and sneezing), starting at 6 and 7 dpv (Group B -Mass + 793B; Group C -aMPV + Mass + 793B respectively), and continuing up to 15 dpv. Group A (aMPV vaccinated) and Group D (Control) showed no clinical signs. No groups presented with moderate clinical signs at any point.

Up to 5 dpc: Following IBV M41 challenge, Groups A1 (aMPV vaccinated) and D1 (Control) showed signs from 1 dpc to 5 and 6 dpc respectively. Mild signs in Groups B1 and C1 subsided after 2 dpc. For aMPV challenge, greater clinical signs were observed in non-vaccinated birds (Group D2), or birds only receiving the IBV vaccines (Group B2), compared to the group receiving the combined vaccination (Group C2). All signs were cleared from aMPV challenge groups by 7 dpc. Groups receiving no challenge virus (A3, B3, C3 and D3) were absent of clinical signs.

All IBV vaccinated groups showed high protection against the M41 challenge strains (91.12-97.84% protection score). The combined groups (C1 and C3) demonstrated a similar level of protection percentage (91.12 and 97.84 respectively) when compared to the single vaccination groups (B1 and B3, both 97.84%).

The mean serum aMPV antibody titre at 0 dpv was 7217.54 (±572.83), indicating the presence of maternally derived antibodies (MDA) against aMPV (Fig. 3) . By 21 dpv, antibody levels declined to below the detectable titre for the IBV vaccinated group (Group B). However, the aMPV vaccinated and combined groups remained at positive titres (3074 and 1668.78 respectively). By 10 dpc, only Groups B1, B3, D1 and D3 showed antibody levels below the detectable limit. Following aMPV challenge, Groups A2 (aMPV vaccinated) and B2 (IBV vaccinated) presented with higher titres (6367.4 and 4543.57 respectively) when compared with Group C2, which received the combined aMPV + IBV vaccination (3509.67). The aMPV vaccinated-aMPV challenged group (A2), and the combined vaccinated-aMPV challenged group (C2) had significantly higher (p < 0.05) titres compared to the non-challenged sub groups A3 and C3. 

H120 (1) 

We were able to detect both IBV and aMPV as early as 3 dpv. However, while IBV was identified until 21 dpv, we did not detect aMPV beyond 7 dpv in either single or combined vaccine groups.

Partial S1 sequencing of IBV positive samples demonstrated a high similarity to the original inoculum, with identities ranging from 97.22 to 99.42% in Group B to 98.21-99.54% in Group C. No Mass-types were identified after 10 dpv in either group. On aver- Fig. 1 . Experiment 1: Quantification of IBV viral load in the trachea and kidney at 21 dpv following vaccination with either (Group B) H120 + CR88 or (Group C) VG/GA-Avinew + H120 + CR88. Data is presented ± standard error margins (SEM). Significance testing carried out within sampling point of the same tissue type, with differences between groups (p < 0.05) labelled with different letters. Data is presented ± standard error margins (SEM). Significance testing carried out between groups of the same tissue type, with differences between groups (p < 0.05) labelled with different letters. The d s /d ns ratio was an average of 3.38 for the Mass-like strains and 1.46 for the 793B samples ( Table 2 ). The two 793B-positive samples with greater than average SNP counts also showed an increase in amino acid changes, and were the only samples with hydrophobicity changes, with the majority (n = 4/6) being a change from hydrophilic to hydrophobic.

Up to 21 dpv: IBV vaccinal strains were recovered from Groups B and C at all sampling points, with the exception of the kidney at 21 dpv (Table 3 ). All detected IBVs had a high partial-S1 nucleotide similarity with the initial vaccine strain (over 99.1%), with only a single non-synonymous change at all time points for Groups B and C. The tracheal sample from Group C at 21 dpv had 20 SNPs, which translated into nine amino acid variations (containing a single hydrophobic to hydrophilic change). The aMPV vaccine was detected up until 21 dpv in tracheal samples, with the kidney remaining negative throughout.

IBV M41 challenge at 5 dpc: IBV was detected in the kidney of Group B1 (IBV-vaccinated) and Group C1 (aMPV + IBVvaccinated), with samples having a high similarity with the vaccine strain (over 99%) ( Table 5 ). All samples from Groups A1 (aMPV vaccinated) and D1 (Non-vaccinated) were positive for IBV, with 100% identity to the challenge strain. Groups not receiving an M41 challenge were IBV negative and no aMPV was detected in birds challenged with M41.

aMPV subtype B challenge at 5 dpc: No aMPV was detected in Group A2 (aMPV vaccinated) or Group C2 (aMPV + IBV vaccinated). The virus was detected in the trachea of birds in Group B2 (IBV vaccinated) and both trachea and kidneys of Group D2 (nonvaccinated).

Post vaccination in the IBV vaccinated tracheal samples, viral load increased from 7 dpv to 14 dpv (Group B = 2.48 to 2.85 log REU; Group C = 1.53 to 1.99 log REU) (Fig. 4) . However, IBV presence was significantly lower (p < 0.05) at 21 dpv for Group B when compared to both 7 and 14 dpv. The IBV vaccination group (Group B) had a significantly higher tracheal viral load when compared to the combined (Group C) vaccine group at both 7 and 21 dpv. While both groups were negative at 7 dpv, viral load in the kidney was significantly higher at 14 and 21 dpv (p < 0.05) in Group B compared to the combined Group C.

Post challenge with IBV M41, groups vaccinated with IBV were showed to be negative by qRT-PCR, whereas Group A (aMPV vaccinated) and Group D (non-vaccinated) showed readings of 2.55 and 2.42 log REU respectively (Fig. 5) .

Post vaccination, kidney samples were negative for IBV from all groups at 7 dpv. However, viral load increased from 14 to 21 dpv in all IBV vaccinated groups (Group B -2.24 to 2.99 log REU; Group C -1.50 to 2.08 log REU). Post challenge with IBV M41, groups not receiving an IBV vaccine showed virus presence (Group A1 -2.10 log REU; Group D1 -2.36 log REU), with vaccinated groups having low detection levels (Group B1 -0.25 log REU; Group C1 -0.06 log REU).

Co-vaccination of day-old broiler chicks with NDV or aMPV alongside two different live IBV (Mass and 793B) vaccines does not impair protection conferred against either IBV or aMPV challenge. In both experiments, following IBV challenge, mild respiratory signs were observed for the first two days and none thereafter, demonstrating clinical protection conferred by IBV Mass H120 + 793B CR88, with or without NDV or aMPV. Ciliostasis analysis was also used to confirm IBV protection [14] , where a higher protection percentage suggests better protection against a virulent IBV challenge. Ciliostasis results in the current study highlighted excellent IBV protection against M41 and QX challenges in groups receiving H120 + 793B, or the same IBV vaccines given along with aMPV or NDV vaccines. Previous work has shown that the combination of Mass + 793B conferred excellent ciliary protection against Asian (QX and Q1) and Middle East (IS/1494/06 and IS/885/00) IBV strains [25, [48] [49] [50] [51] . When IBV H120 and CR88 vaccine strains are combined with the NDV or aMPV vaccine, protection levels remained at or above 94%. This reconfirms that such combined vaccination did not compromise the protection conferred against virulent IBVs. In addition to the cilia-stopping test, detection of virus and viral load in the trachea and kidney of challenged groups was attempted to demonstrate protection against IBV M41 and QX. No vaccine viruses were detected in the trachea at 5 dpc. However, both virulent challenge viruses were present, suggesting host clearance of IBV vaccines by at least 25 dpv [52, 53] . The co-administration of either NDV or aMPV alongside the IBV vaccine caused a lower IBV vaccine virus load in both the trachea and kidney at 21 dpv. However, post-challenge with virulent IBV, the combined IBV with NDV or aMPV vaccination had no effect on viral load when compared with IBV vaccination alone. For both experiments, despite relatively higher tracheal viral load, a low REU value was shown in kidney tissue from non-IBV vaccinated groups following M41 and QX challenge. These results, along with the cilia-stopping test, demonstrated that comparable levels of protection are found when IBV H120 and CR88 vaccine strains are given either alone or alongside an NDV or aMPV vaccine.

Control of Newcastle disease is of paramount importance in endemic countries. Control is normally achieved by the use of live and inactivated vaccines [54] . Haemagluttination inhibition titres were traditionally used to confirm protection against NDV challenge [55, 56] , and the mean NDV HI titre in both single and combined vaccinated groups described here were >5.63 log 2 , which is above the protective titre against virulent NDVs [19, 23, 57, 58] . Current work showed that simultaneous vaccination of NDV + H120 + CR88 in commercial broiler chicks does not interfere with the induction of protective immunity against NDV. In this study, following challenge with virulent aMPV, no clinical signs were found in the single (aMPV) or combined (aMPV + H120 + CR88) vaccinated groups, confirming clinical protection against the challenge virus [19, 23, 59] . Furthermore, following challenge with the virulent strain, no antigen was detected by RT-PCR from tissue samples in either of the vaccinated groups. The aMPV ELISA antibody titres in the triple-vaccinated group was lower than the single vaccinated group, echoing findings from previous work [23] . However, it has been reported that levels of humoral antibodies have no association with clinical protection against aMPV [60, 61] . It is evident that experimental vaccination of commercial broilers (with aMPV and IBV MDAs) with aMPV subtype B and two IBV strains (H120 and CR88) had no adverse effects on clinical protection conferred against a virulent challenge.

In both experiments, the persistence of vaccinal strains following individual IBV, NDV, aMPV and combined (H120 + CR88 + NDV; H120 + CR88 + aMPV) vaccination were monitored at intervals. For NDV and aMPV, classical RT-PCRs [32, 40] were applied to RNA extracted from oropharyngeal swabs, with no detection of the NDV vaccine at any sampling points. Previous work has shown detection of NDV vaccine virus in SPF chickens for up to four days post vaccination [19, 62] . In this experiment, first sampling was carried out at 7 dpv, and the failure to detect NDV may have been due to rapid clearance of the vaccine [19] . For aMPV, it was possible to detect the virus up to 7 dpv using dry oropharyngeal swabs and until 21 dpv from tracheal tissue. While identified the aMPV vaccine at 21 dpv in a combined IBV + aMPV vaccinated group [19] , there was no difference between single and combined vaccinated groups in the current study, suggesting that combined vaccination does not impair a host's ability to rapidly clear the aMPV vaccine strain. IBV vaccinal strains have previously been shown to persist in chickens for a longer period [53, 63] . To further investigate this, the growth kinetics and genetic characteristics of live IBV vaccines used in this study were cross compared with previous reports [26, 64] . Vaccinal strains were detected in OP swabs throughout the current experiment, with the variant 793B strain dominating at the later ages of the birds. With both experiments, Sanger sequencing showed that the genotype of strains recovered using OP swabs altered from Mass to 793B after 14 dpv. As reported previously [26, 27] , while the Mass genotype made up the majority of detections during early sampling days (3-14 dpv), the 793B strain became the dominant strain by 21 dpv. The reason for the shifting of genotypes from Mass to 793B is not known, however, previous work has highlighted that different IBV strains can differentially induce host expression pathways, in particular TLR7 [65] and IFN-b [66] . Such differences may have played a role in a host's ability to clear the Mass vaccine virus at an earlier time point when compared to the 793B vaccine. Further work is required to understand potential underlying mechanisms.

At no point did the nucleotide similarity for recovered Mass vaccine strains drop below 98% when compared to the original vaccine sequence. This suggests that under experimental conditions, the vaccine strain underwent minimal alteration following inoculation [26] . Similarly, the 793B vaccine strain remained high for all but one sampling point (99-100% with the exception of 96% at 21 dpv in Experiment 1). Limited genetic changes were noted in the tissues samples, with 3-4 and 4-7.5 SNPs per sequence in the Mass and 793B vaccine strains respectively. A higher number of polymorphisms were seen in the 793B sequences, compared to Massachusetts, indicating that the persistence of this particular genotype may contribute to greater genetic variations within local virus populations [64, 67, 68] . However, the d s /d ns ratios suggest purifying selection for both genotypes [69] , and therefore evolutionary constraint [70] . Previous work implied 793B amino acid variations do not follow a consistent pattern over time [71] , which was also witnessed in the current study. No persistent alteration was witnessed from any tissues in any of the vaccination programs.

In conclusion, this study indicates that the combined vaccination consisting of NDV (VG/GA-Avinew) or aMPV (subtype B) alongside two strains IBV (Massachusetts -H120 and 793B -CR88) does not impair protective immunity against the globally important IBV, NDV or aMPV strains that were included in this study.

Taxonomy of the order Mononegavirales: update 2019

Coronavirus avian infectious bronchitis virus

Newcastle disease

Infectious Bronchitis. 13th ed. New Jersey: Whiley-Blackwell

Vaccination of commercial broiler chicks against avian metapneumovirus infection: a comparison of drinking-water, spray and oculo-oral delivery methods

Avian Metapneumovirus circulation in Italian broiler farms

Evaluating viral interference between infectious bronchitis virus and Newcastle disease virus vaccine strains using quantitative reverse transcription-polymerase chain reaction

Viral respiratory diseases (ILT, aMPV infections, IB): are they ever under control?

Effect of associated vaccines on the interference between Newcastle disease virus and infectious bronchitis virus in broilers

Isolation and Molecular Characterization of Novel Infectious Bronchitis Virus Variants from Vaccinated Broiler Flocks in Egypt

Emergence of a genotype I variant of avian infectious bronchitis virus from Northern part of India

Seroprevalence and genotyping of avian infectious bronchitis virus detected from Iranian unvaccinated backyard chickens

The use of vaccination in poultry production

Breadth of protection of the respiratory tract provided by different live-attenuated infectious bronchitis vaccines against challenge with infectious bronchitis viruses of heterologous serotypes

Review of Newcastle disease virus with particular references to immunity and vaccination. World's

Infectious bronchitis virus interference with growth of Newcastle disease virus. II. Interference in chickens

Antibody responses to a live Newcastle disease vaccine when combined with a live infectious bronchitis vaccine

Impact of mild Newcastle disease vaccines on control of IBV

Immune responses and interactions following simultaneous application of live Newcastle disease, infectious bronchitis and avian metapneumovirus vaccines in specific-pathogen-free chicks

Infectious bronchitis virus vaccine interferes with the replication of avian pneumovirus vaccine in domestic fowl

Interaction between live avian pneumovirus and Newcastle disease virus vaccines in specific pathogen free chickens

Protection Conferred by a Live Avian Metapneumovirus Vaccine when Co-Administered with Live La Sota Newcastle Disease Vaccine in Chicks

Protection in specific pathogen free chickens with live avian metapneumovirus and Newcastle disease virus vaccines applied singly or in combination

Increased level of protection of respiratory tract and kidney by combining different infectious bronchitis virus vaccines against challenge with nephropathogenic Brazilian genotype subcluster 4 strains

Cellular, and Humoral Immune Responses Induced by Different Live Infectious Bronchitis Virus Vaccination Regimes and Protection Conferred against Infectious Bronchitis Virus Q1 Strain

Infectious bronchitis vaccine virus detection and part-S1 genetic variation following single or dual inoculation in broiler chicks

Genetic mutations in live infectious bronchitis vaccine viruses following single or dual in vitro infection of tracheal organ cultures

A 15-year analysis of molecular epidemiology of avian infectious bronchitis coronavirus in China

Immunopathogenesis of infection in SPF chicks and commercial broiler chickens of a variant infectious bronchitis virus of economic importance

The use of FTA Ò filter papers for diagnosis of avian influenza virus

Detection and differentiation of Newcastle disease virus (avian paramyxovirus type 1)

Longitudinal field studies of infectious bronchitis virus and avian pneumovirus in broilers using typespecific polymerase chain reactions

Diagnosis of a naturally occurring dual infection of layer chickens with fowlpox virus and gallid herpesvirus 1 (infectious laryngotracheitis virus)

Evaluation of FTA Ò paper and phenol for storage, extraction and molecular characterization of infectious bursal disease virus

Hexon based PCRs combined with restriction enzyme analysis for rapid detection and differentiation of fowl adenoviruses and egg drop syndrome virus

Pathogenicity of infectious bronchitis virus isolates from Ontario chickens

Procedures for the haemagglutination and the haemagglutination inhibition tests for avian infectious bronchitis virus

The single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction: twenty-something years on

Genotypes of infectious bronchitis viruses circulating in the Middle East between

A molecular epidemiological study of avian paramyxovirus type 1 (Newcastle disease virus) isolates by phylogenetic analysis of a partial nucleotide sequence of the fusion protein gene

Development and validation of RT-PCR tests for the detection and S1 genotyping of infectious bronchitis virus and other closely related gammacoronaviruses within clinical samples

Laboratory evaluation of a quantitative real-time reverse transcription PCR assay for the detection and identification of the four subgroups of avian metapneumovirus

MEGA6: Molecular Evolutionary Genetics Analysis version 6.0

CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice

Molecular Evolution and Phylogenetics

Positive Darwinian evolution in human influenza A viruses

A simple method for displaying the hydropathic character of a protein

Characterization of three infectious bronchitis virus isolates from China associated with proventriculus in vaccinated chickens

Isolation and identification of glandular stomach type IBV (QX IBV) in chickens

Protection conferred by live infectious bronchitis vaccine viruses against variant Middle East IS/885/ 00-like and IS/1494/06-like isolates in commercial broiler chicks

Heterologous live infectious bronchitis virus vaccination in day-old commercial broiler chicks: clinical signs, ciliary health, immune responses and protection against variant infectious bronchitis viruses

Comparative dynamic distribution of avian infectious bronchitis virus M41, H120, and SAIBK strains by quantitative real-time RT-PCR in SPF chickens

Establishment of persistent avian infectious bronchitis virus infection in antibody-free and antibody-positive chickens

Control of Newcastle disease by vaccination

Immune response of chicks to oral vaccination against Newcastle disease and fowl pox

Protection of the reproductive tract of young chicks by Newcastle disease virus-induced haemagglutinationinhibition antibodies

Immunization with a thermostable newcastle disease virus K148/08 strain originated from wild mallard duck confers protection against lethal viscerotropic velogenic newcastle disease virus infection in chickens

Serology-based potency test for inactivated Newcastle disease vaccines

Demonstration of loss of attenuation and extended field persistence of a live avian metapneumovirus vaccine

Investigations on the protective role of passively transferred antibodies against avian metapneumovirus infection in turkeys

A live attenuated turkey rhinotracheitis virus vaccine. 2. The use of the attenuated strain as an experimental vaccine

Comparison of viral shedding following vaccination with inactivated and live Newcastle disease vaccines formulated with wild-type and recombinant viruses

A reverse transcriptase-polymerase chain reaction survey of infectious bronchitis virus genotypes in Western Europe from

Avian coronavirus infectious bronchitis attenuated live vaccines undergo selection of subpopulations and mutations following vaccination

Early immune responses and development of pathogenesis of avian infectious bronchitis viruses with different virulence profiles

Differential innate immune responses induced by classical and variant infectious bronchitis viruses in specific pathogen free chicks

Infectious bronchitis virus: dominance of ArkDPI-type strains in the United States broiler industry during the last decade. Revista Brasileira de Ciência Avícola

Rapid selection in chickens of subpopulations within ArkDPI-derived infectious bronchitis virus vaccines

Likelihood models for detecting positively selected amino acid sites and applications to the HIV-1 envelope gene

Genetic diversity and purifying selection in West Nile virus populations are maintained during host switching

Variation in the spike protein of the 793/B type of infectious bronchitis virus, in the field and during alternate passage in chickens and embryonated eggs

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.