Production of infectious inoculum. The ABBV1 isolate utilized in this experiment was collected from the brain of a naturally infected Canada goose (GenBank number MK966418) and passaged in immortalized duck embryo fibroblasts (American Type Culture Collection CCL-141)37.
The infectious inoculum was produced as described previously15-17,24. Briefly, CCL-141 cells persistently infected with ABBV1 were grown to confluency in groups of 100 mm dishes, and the virus was released via osmotic shock (using distilled water) followed by one cycle of freeze / thaw in an ultra-cold freezer (-80 °C). The lysate was clarified by centrifugation at 3000 x g for 10 min, and the supernatant concentrated by precipitation with 30 % w/v polyethylene glycol (PEG)-8000 overnight at 4 °C, followed by centrifugation at 3,200 x g for 15 min. The pellet was resuspended in 1 x phosphate-buffered saline (PBS) with 1 % fetal bovine serum (FBS) to obtain an original supernatant-to-final resuspended volume ratio of approximately 100:1, which was aliquoted and kept at -80 °C until inoculation. The amount of infectious virus present in the stock was quantified by preparing a series of ten-fold dilutions in 2 % FBS/ Dulbecco’s Modified Eagle Medium (DMEM; 200 µL/well, Corning, New York, United States) which were added to CCL-141 cells in an array of wells using a 96-well-plate format (100 µL per well), and incubated at 37 °C for 5 to 7 days. Due to the lack of cytopathic effects, the infection status of each well was determined by immunofluorescence assay, as previously described37. A 50 % tissue culture infectious dose (TCID50) was calculated according to the Spearman-Karber method38 and converted into focus forming units (FFUs).
Animal experiment and sample collection
Experimental design. The layout of the experimental plan is summarized in Supplementary Figure 1. A total of 205, 1-day-old domestic turkeys (Meleagris gallopavo domesticus) were purchased from a commercial hatchery in Ontario (Hendrix Genetics, Kitchener, Ontario, Canada) and delivered to the University of Guelph Production Animal Isolation Facility (Guelph, ON, Canada) in 2 groups, 6 weeks apart. The first cohort (old) consisted of 92, 1-day-old birds (80 birds plus 15 % to account for attrition) that were randomly assigned into 2 rooms (virus-inoculated (VI) and control (CO), each n = 46. At 56 days of age, because of their larger size, CO birds were moved to a larger room, and those in the VI group were split between 2 rooms. The second cohort (young) consisted of 113, 1-day-old birds (80 birds plus 41 % attrition), which were randomly assigned to an additional 2 rooms (VI and CO, respectively, n = 57 and 56). Increased poult numbers for attrition in the young cohort was based on increased unexpected mortality in the previous cohort. As these birds were kept for up to 12 weeks of age, additional room space was not needed.
All poults were housed in negative pressure rooms. Upon being received, poults were neck tagged, placed on the floor with wood shaving, and provided with a standardized photoperiod, as well as ad libitum access to feed and water, according to hatchery guidelines39. For the first two weeks, poults were kept in brooding rings and provided supplemental heat using heat lamps approximately half a meter from the floor. Feeders were checked daily and replenished as needed. Waterers were equipped with a trough, flushed daily, and refilled automatically.
In order to synchronize the time of inoculation between the two cohorts (Supplementary Figure 1), the old cohort was kept for 6 weeks before inoculation (i.e., poults were inoculated at 6 weeks of age), while the young cohort was received from the producer 24 hours before inoculation (i.e., poults were inoculated between approximately 24 - 48 hrs of age). A total of 40 old and 50 young turkeys were inoculated in the left gastrocnemius muscle with 6.9 x 105 focus forming units (FFUs) of ABBV1 in 100 µL volume. The additional 7 young turkeys were inoculated with the same amount of virus into 250 µL (due to a less concentrated batch of the stock), delivered in 100 and 150 µL aliquots in the left and right gastrocnemius muscle, respectively. The control birds (44 old and 56 young) were inoculated with 100 µL of carrier only (PBS with 1 % FBS). Note the number of inoculated old turkeys is lower compared to the number initially purchased due to initial attrition before inoculation and sampling of bird sera for ELISA optimization (pre-inoculation mortalities are not reported).
Turkeys were monitored daily for the development of clinical signs for the entire duration of the study (18 weeks). Emergency euthanasia (i.e., humane euthanasia outside of the scheduled sampling time) was elected in cases of severe lameness, moderate lameness intractable to medical management, lethargy, moderate to marked decrease in body condition score, neurologic signs (any combination of head tremors, torticollis, paresis/paralysis), severe feather pecking, or ascites. These culls were classified as emergency birds. Animals with minor clinical signs, such as ruffled feathers, huddling with conspecifics, mild decrease in body condition score (based on breast muscle examination), feather pecking, mild to moderate lameness (but able to maintain access to feed and water) were monitored closely. Where medically appropriate (e.g., for lame animals), a small portion of the trial room was partitioned to create a separate pen where convalescent animals were placed and given ad libitum access to food and water. Lame animals were given daily meloxicam injections as needed (1mg / kg, IM, alternating right and left pectoralis muscle). Pine tar was applied to animals experiencing mild to moderate feather pecking. During periods of feather pecking, photoperiod was temporarily reduced to 7 hours of light per day with a light intensity of 10 lux.
Euthanasia. At 1, 4, 8, and 12 weeks post-inoculation (wpi), n = 10 turkeys from each treatment group were selected and euthanized for sample collection. At each sampling time, birds demonstrating mild to moderate signs of illness were intentionally selected, while the remaining turkeys were collected without specific criteria. Randomized selection based on neck tag number was not used to avoid the stress of handling and the occurrence of sudden deaths due to underlying cardiovascular issues (see incidental lesions). Poults that were less than or equal to 4-week-old were anesthetized with isoflurane in a 7 L vented induction chamber (VetEquip©, Livermore, California, USA) and then euthanized by carbon dioxide (CO2) inhalation. Older turkeys were weighed, sedated (dexmedetomidine [0.03 mg / kg; Dexdomitor, Zoetis, Parsippany-Troy Hills, New Jersey, USA] and ketamine [30 mg/kg; Narketan, Vetoquinol, Lavaltrie, Québec, Canada] delivered intramuscularly), and euthanized by pentobarbital overdose (100 mg / kg; Euthasol, Virbac, Cambridge, Ontario, Canada) delivered intravenously in the right or left basilic vein.
Postmortem Examination, Swabs, and Tissue Collection. Immediately after euthanasia, 2 – 6 mL of blood was collected from the jugular vein using a 5 - 10 ml syringe, and injected into a red-top tube (for serum separation). Next, the choanal slit and cloaca were swabbed using sterile cotton-tipped applicators and separately placed in cryovials containing 0.7 ml of RNA preserving solution (20 mM ethylenediaminetetraacetic acid [EDTA], 25 mM sodium citrate, and 70 % [weight per volume] ammonium sulfate with a pH of 5.2). For each turkey, a routine postmortem examination was conducted, gross lesions were recorded, and the following 5 tissues were sampled for both RNA extraction and histology (partial sampling, Supplementary Figure 1): brain, lumbar spinal cord, kidneys, proventriculus, and gonads. Separate samples were collected either in 10 % neutral buffered formalin or in 1 ml sterile screwcap tubes containing 0.7 mL of RNA preserving solution. Samples in preserving solution were kept at -80 °C until needed for RNA extraction. For the brain, small portions of the cerebrum, cerebellum and brainstem were collected from one half of the brain for RNA extraction, while the remaining tissue was fixed in formalin. For more complete histopathological assessment, several additional tissues were collected from 3 turkeys in each cohort per sampling point (full sampling). These tissues were directly fixed in formalin and included: globe, vagus nerve, brachial plexus, sciatic nerve, cervical spinal cord, thoracic spinal cord, skeletal muscle, bone marrow, spleen, skin, uropygial gland, trachea, lung, esophagus, crop, ventriculus, duodenum, jejunum, ileum, ceca, colon, cloaca, liver, bursa, thymus, Meckel’s diverticulum, thyroid gland, adrenal gland, pancreas, and heart.
Unexpected mortalities or turkeys euthanized outside of sampling schedule (emergency birds) underwent partial sampling, unless excessively autolyzed. For these birds, test results were added to those of the subsequent sampling time point for group data analysis.
Animal Use Approval. The Animal Care and Use Committee of the University of Guelph granted approval for all experimental procedures (Animal Utilization Protocol #3978). The conducted activities adhered to the applicable regulations established by the Canadian Council on Animal Care.
Sample processing and Data Analysis
Histological Assessment and Scoring. Tissues for histology were fixed in 10 % neutral buffered formalin, and routinely processed for staining with hematoxylin and eosin (HE) by the Animal Health Laboratory histotechnology laboratory. Over the course of the trial, tissues were collected as follows: 34 birds prior to planned endpoint (20 emergency and 14 found dead birds) underwent partial sampling, 107 underwent partial sampling at the planned endpoint, and 48 underwent full sampling at the planned endpoint (Supplementary Figure 1). Tissues were processed for histology from all the emergency birds, and all the turkeys euthanized in-time point that underwent full sampling. Of the birds euthanized in-time point that underwent partial sampling, all those euthanized at 12 wpi (n = 28) were processed, while tissues from 1, 4 and 8 wpi were not processed due to lack of CNS inflammation as seen in the full sampling birds from the same time point.
A member of the investigative team (LG) identified and tallied the microscopic lesions across all tissues and scored inflammation in the brain and spinal cord (CNS) in a blind fashion. Scoring was performed using the same rubric devised previously for chickens17. Briefly, 4 regions of the brain (cerebrum, cerebellum, optic lobe, and brainstem) and all available segments of the spinal cord (lumbar +/- cervical and thoracic) were scored separately by assessing the average thickness of the perivascular cuffs (PVC; intensity of inflammation) and the average number of vessels affected (distribution of inflammation) in 10 randomly selected 100x fields (one 100x field with F22 lens aperture = 3.80 mm2). Each region received a score (regional score, RS) that was the sum of the distribution (no PVCs, score 0; 1-10 PVCs, score 1; 11-20 PVCs, score 2; > 21 PVCs, score 3) and intensity (no inflammation, score 0; 1-2 layers wide, score 1; 3-4 layers wide, score 2; > 5 layers wide, score 3) scores, ranging between 0 - 6. The score for the brain and spinal cord (organ scores, OS) were calculated respectively as the simple average between all assessed brain regions (i.e., average of RS for cerebrum, cerebellum, optic lobe, and brainstem), or all assessed spinal cord segments (i.e., average RS for lumbar +/- thoracic and cervical). Lesions outside the CNS were not scored, and their presence was arranged into contingency tables.
Immunohistochemistry. Immunohistochemistry (IHC) for ABBV1 nucleoprotein (N) was performed by the Animal Health Laboratory (Guelph, Ontario, Canada), and conducted in all of the tissues collected from the VI birds that underwent full sampling at 12 wpi, as well as the five young VI emergency birds collected between 8 and 12 wpi (Table 5). Immunohistochemistry for CD3 (cluster of differentiation 3; T cell marker) and PAX5 (paired box 5; B cell marker) was conducted on a single full sampling turkey sampled at 12 wpi from the young and old VI groups. Tissues from a single young and old control bird (T251 and T141, respectively) were used as negative controls. Positive controls, employed in each run, consisted of the brain of an ABBV1-infected Canada goose (for the N protein) and a canine lymph node (for the CD3 and PAX5 antigens). Information about the antibodies is summarized in Supplementary Table 6. Briefly, sections were mounted on positively charged slides (Fisherbrand™ Superfrost™ Plus Microscope Slides, Thermo Fisher Scientific, Waltham, Massachusetts, United States), de-paraffinized and rehydrated, and antigens retrieved using either a citrate-based solution (pH 6) at 120 °C for 2 min, tris(hydroxymethyl)amino-methane-ethylenediaminetetraacetic acid (TRIS-EDTA) (pH 9) at 98 °C for 25 min, or Proteinase K (Dako, Agilent Technologies, Santa Clara, California, United States) at room temperature for 12 min. Endogenous peroxidases were quenched with 3 % hydrogen peroxide for 10 min, and non-specific binding was blocked with serum-free protein solution (Dako, Agilent Technologies, Santa Clara, California, United States). Primary antibodies were incubated at room temperature in a humid chamber, and after washes with Tris wash buffer (Dako, Agilent Technologies, Santa Clara, California, United States), sections were incubated with secondary antibodies conjugated to horseradish peroxidase (HRP) (Dako, Agilent Technologies, Santa Clara, California, United States). The reaction was revealed by a chromogenic substrate (NovaRED Horseradish Peroxidase substrate; Vector Laboratories, Newark, California, United States), and sections counterstained with hematoxylin and mounted using Surgipath Micromount (Leica Biosystems, Wetzlar, Germany).
3.3.3 Viral RNA Quantification in Tissues and Swabs. Presence of ABBV1 RNA in tissues and swabs was determined by RT-qPCR. RNA was extracted from approximately 300 µg of tissue sample and 150 µL of swab solution using the E.Z.N.A® RNA Kit II (Omega Bio-Tek, Norcross, Georgia, United States). The target segment of the ABBV1 N gene was amplified using a set of primers coupled with a fluorescent probe, using the Luna® Universal Probe One-Step RT- qPCR kit (New England Biolabs, Whitby, Ontario, Canada) as previously described15-17. Each sample was run in triplicate using a 96-well plate on a Roche LightCycler® System (Rotkreuz, Switzerland), and the cycle threshold (Ct) values were interpolated as copy numbers using a standard curve made of serial dilutions of a 500 basepair (bp) gBlock gene cassette (integrated DNA technologies) containing the target sequence of the ABBV1 N gene. A known positive sample and nuclease-free water were used as positive and negative controls, respectively. The RT-qPCR output consisted of copy number per 150 ng of total RNA (for tissues), or copy number per 150 µL of swab fluid. Samples were considered positive for Ct < 35.
Serology. The blood collected from the jugular vein in the red-top tube was kept at room temperature and allowed to coagulate (up to 8 h post collection). The serum was separated from the clot by centrifugation at 2000 x g for 10 min, pipetted into sterile 1 mL screw top vials and preserved at -80 °C until used. As in previous trials, seroconversion was assessed by indirect ELISA across all groups and 1, 8, and 12 wpi time points15-17. A 96-well microtitre plate was coated for 2 hours with recombinant full-length ABBV1 N protein (Biomatik, Kitchener, Ontario, Canada) at 25 ng / well at 37 °C and blocked overnight at 4 °C with 200 μL / well of 5 % fish gelatin in PBS. Turkey serum (50 μL; diluted 1:200 in PBS-T with 1 % fish gelatin) was added to the wells of the microtitre plate and incubated for 2 hours at 37 °C. After washes (3 times with PBS-T and 2 times with water), a secondary goat anti-Turkey immunoglobulin Y (IgY) (H+L)-HRP conjugated antibody (SouthernBiotech, Birmingham, Alabama, United States) was applied at 1:4000 dilution for 1 hour at 37 °C. After a last series of washes (3 times with PBS-T and 2 times with water) the chromogenic substrate (Pierce 1-Step Ultra TMB-ELISA Substrate Solution, Thermo Fisher Scientific, Waltham, Massachusetts, United States) was added to the plate (50 μL / well) to detect bound serum antibodies. Chromogenic substrate reaction was stopped by addition of 0.2 N sulfuric acid solution (50 μL / well). The intensity of the signal (optical density, OD) was measured using an EnSpire multimode plate reader (PerkinElmer, Waltham, Massachusetts, United States) at a wavelength of 450 nm. All samples were run in triplicate, and a single sample exhibiting moderate to high reactivity was run in triplicate for each plate, serving as a calibrator. To minimize variation across plates, the OD value of each well on a plate was divided by the average OD value of the calibrator wells on that plate, and then multiplied by the mean OD value of all calibrators across plates, to obtain normalized values. The threshold for classifying sera as positive, regardless of age cohort, was established as the average OD value calculated from all control sera of the young and old cohorts, plus three times the standard deviation (SD). All analyses were based on the averages from triplicate wells.
Statistical Analysis.
Week mortality rates were calculated using the total number of culled or spontaneously dead birds over the total number of at-risk bird-weeks40. Proportions (contingency tables) were compared by Fisher’s exact test. Differences in the average pathology scores (RS and OS) or normalized (i.e., calibrated) OD values between groups and time points were tested using the non-parametric Kruskal-Wallis test, followed by the post-hoc Dunn’s test for multiple comparisons. Group differences in the concentration of ABBV1 RNA copies (log-10 transformed) in tissues were tested by a two-way analysis of variance (ANOVA) fitting a full model (variables: age group and time post-inoculation), followed by Tukey’s test for multiple comparisons. Analyses were carried out using GraphPad Prism (version 9.5.1 for Windows, GraphPad Software, Boston, Massachusetts USA, www.graphpad.com), with significance for all tests was set at p < 0.05.
Multivariable regression analysis. Multivariable logistic regression analysis, implemented using Firth’s penalized likelihood method41, was conducted to evaluate correlation between infectious status (Model 1) of the turkeys or inflammation (Model 2) in the spinal cord (outcome categorical variables; presence / absence), with a series of independent variables. These included age cohort (young / old, categorical variable), virus RNA amount in the bran / spinal cord (continuous), brain and spinal cord pathology scores (continuous), magnitude of seroconversion (OD values), clinical signs (categorical), and presence of incidental lesions (gross and microscopic; categorical). The magnitude of seroconversion was quantified using OD values multiplied by 100. This scaling was implemented because utilizing the raw OD values in the analysis yielded excessively large odds ratios and broad 95% confidence intervals, potentially compromising the reliability of statistical interpretations40. To avoid having categories with < 5 % frequency, clinical signs were grouped as a single category (presence / absence); for gross lesions, all cardiovascular and musculoskeletal lesions were merged, while dilated cardiomyopathy (DCM) and angular limb deformity (ALD) were kept separately. Microscopic incidental findings were classified in the following 2 categories: interstitial inflammation of the kidney, and interstitial inflammation of the gonads.
A relaxed p-value (p ≤ 0.2) was initially implemented to seek out explanatory variables, which were incorporated into the multivariable models if they were shown to be significantly associated with the dependent variable using univariable logistic models. Explanatory variables were assessed for collinearity by implementing pairwise Spearman’s correlation tests as in previous studies. In brief, the multivariable models were constructed with the variables with the smallest p-value when two variables were highly correlated (rho > 0.70; p < 0.05). Virus RNA concentration values were log10-converted, and all zero values (Ct value ≥ 35) were recorded as 1.0 x 10-5. Analyses were carried out using R (version 4.3.1, Vienna, Austria, http://www.r-project.org).