key: cord-0875764-8il7xiw8
authors: Ayala, Andrea J.; Dimitrov, Kiril M.; Becker, Cassidy R.; Goraichuk, Iryna V.; Arns, Clarice W.; Bolotin, Vitaly I.; Ferreira, Helena L.; Gerilovych, Anton P.; Goujgoulova, Gabriela V.; Martini, Matheus C.; Muzyka, Denys V.; Orsi, Maria A.; Scagion, Guilherme P.; Silva, Renata K.; Solodiankin, Olexii S.; Stegniy, Boris T.; Miller, Patti J.; Afonso, Claudio L.
title: Presence of Vaccine-Derived Newcastle Disease Viruses in Wild Birds
date: 2016-09-14
journal: PLoS One
DOI: 10.1371/journal.pone.0162484
sha: 1a29f97fc037ea6456b41d054e5603fbd4f76e94
doc_id: 875764
cord_uid: 8il7xiw8

Our study demonstrates the repeated isolation of vaccine-derived Newcastle disease viruses from different species of wild birds across four continents from 1997 through 2014. The data indicate that at least 17 species from ten avian orders occupying different habitats excrete vaccine-derived Newcastle disease viruses. The most frequently reported isolates were detected among individuals in the order Columbiformes (n = 23), followed in frequency by the order Anseriformes (n = 13). Samples were isolated from both free-ranging (n = 47) and wild birds kept in captivity (n = 7). The number of recovered vaccine-derived viruses corresponded with the most widely utilized vaccines, LaSota (n = 28) and Hitchner B1 (n = 19). Other detected vaccine-derived viruses resembled the PHY-LMV2 and V4 vaccines, with five and two cases, respectively. These results and the ubiquitous and synanthropic nature of wild pigeons highlight their potential role as indicator species for the presence of Newcastle disease virus of low virulence in the environment. The reverse spillover of live agents from domestic animals to wildlife as a result of the expansion of livestock industries employing massive amounts of live virus vaccines represent an underappreciated and poorly studied effect of human activity on wildlife.

Introduction sequenced NDV wild bird isolates from Ukraine, Bulgaria and Brazil. Third, we performed active NDV surveillance in a wild Rock Pigeon (Columba livia) population within an urban environment (Atlanta, GA, USA) to assess the presence of Newcastle disease viruses. Lastly, we analyzed the compilation of the GenBank and newly obtained wild bird NDV sequences, and compared them to currently used live NDV vaccines.

The fieldwork was conducted under the Georgia Department of Natural Resources Permit Number 29-WJH-16-48, and was conducted according to University of Georgia Institutional Animal Care and Use Committee (IACUC) protocol AUP #: A2012 03-031-Y1-A0. All birds were captured on public land owned and managed by the City of Atlanta. No protected, endangered or threatened species were involved in this research, and Rock Pigeons captured in the summer were all released within thirty minutes of capture, as per permit regulations. Rock Pigeons captured by APHIS followed APHIS Wildlife Services protocols as outlined in the 2011 State of Georgia Environmental Assessment Protocols [49] .

Rock Pigeons (n = 72) were captured at three urban sites inside the city limits of Atlanta, Georgia (Fulton County) between May and October 2012 within a 10 km 2 area. From May to July (summer capture), pigeons were caught at three urban sites, while from August to October 2012 (fall capture), pigeons were captured by the USDA-APHIS Wildlife Services personnel at multiple sites as part of their integrated wildlife damage management program [50] . Pigeons were captured using a combination of hand-nets, drop-nets, mist-nets, and ground traps. During the summer capture period, all individuals were banded with a unique four-color combination to avoid pseudoreplication [51] .

To ensure consistency, assessment of demographic attributes, recording, and analysis of all wild Rock Pigeon physiological characteristics was performed by the same person. Each bird was placed into one of two age classes: hatch-year birds (HY) who were born during the current breeding cycle or after hatch-year birds (AHY) born prior to 2012 [52] . Breeding status was established using reproductive characteristics, such as the presence of a brood patch, cloacal protuberance, or both [53] . The sex of AHY birds was determined by reproductive characteristics (if present) and, when observed, by courtship behaviors [54] . Individual body conditions were assessed using two different criteria: the ratio of bird mass to wing chord length (also known as the wing-loading aspect), and the relative amount of visible fat. The weight-to-wing chord ratio was calculated as the body mass divided by the length of the un-flattened longest primary feather [53, 55] , while the MAPS (Monitoring Avian Productivity and Survival) protocol [56] was used for fat scoring. Only wing-pit fat was assessed due to the density of plumage at the furculum and stomach, and scored on an eight-point scale of 0 to 7. A zero fat score indicated no visible yellow subcutaneous fat at the wing-pit and a score of seven indicated a bulging pocket of fat.

Up to 1 mL of blood was taken from the brachial vein of each individual, transferred to a 3 mL serum vacutainer (Becton Dickinson, San Jose, CA, USA), and placed in a cooler with ice packs for transport to the lab. Vacutainers were tilted overnight at room temperature, then centrifuged at 1500 rpm, and the harvested sera were subsequently stored at -20°C [57] . Oral and cloacal swabs were collected from each bird, then transferred into separate 2.0 mL cryovials (Corning Inc., Corning, NY, USA) containing 1.5 mL of Brain-Heart Infusion (Becton Dickinson, San Jose, CA, USA) mixed with Gentamicin (200 μg/mL), Penicillin (2000 units/mL), and Swab medium from each cryovial was inoculated in 9-11 day-old specific-pathogen-free (SPF) embryonated chicken eggs (ECE) in Biosafety Level-2 at the Southeast Poultry Research Laboratory (SEPRL). The recovered allantoic fluid was passaged once more. Protocols for virus isolation and titration followed OIE standards [59] . After a week of incubation, harvested allantoic fluids were tested using the hemagglutination assay (HA). Samples with HA activity by the second passage were considered positive and set aside for further characterization. Viral titers were assessed via 1:10 serial dilutions (10 −5 to 10 −10 ). Five SPF ECE per dilution were inoculated with 100 μL of each dilution, incubated for a week, then evaluated with the HA test. The viral mean embryo infectious dose (EID 50 ) were calculated as described previously [60] .

The virulence of Rock Pigeon isolates was determined from allantoic fluid by ICPI test in 1-day-old SPF chickens [61] . Viruses with ICPI values below 0.7 were considered of low virulence [59] . Lastly, Rock Pigeon sera were individually assessed using hemagglutination inhibition (HI) assay in 96-well plates [59] using antigens for APMV serotypes 1, 2, 3, 4, 6, and 7 [62] . Only samples demonstrating complete HI at the serum dilution 1:16 (2 4 ) or higher were considered positive for presence of NDV antibodies [63] .

Newcastle disease viruses from laboratory repositories in Bulgaria and Ukraine (Table 1) were submitted to SEPRL for further characterization. Each sample was propagated into 9-11 dayold SPF ECE using standard methods [59] . RNA extractions of USA pigeon isolates and viruses from Bulgaria and Ukraine were performed using TRIzol LS (Life Technologies, Carlsbad, CA, USA), as per manufacturer instructions. One-step reverse-transcriptase PCR amplification proceeded using the SuperScript 1 III One-Step RT-PCR System with Platinum 1 Taq DNA Polymerase (Life Technologies, Carlsbad, CA, USA) and previously described primers (4331F/ 5090R, 4911F/5857R, 5669F/6433R, 4961F/5772R) were used for the PCR and sequencing [64] . Amplicons were separated through a 1% agarose gel, with the ensuing DNA bands excised and purified using the QuickClean II Gel Extraction Kit (GenScript, Piscataway, NJ, USA). Nucleotide sequencing and assembly were performed as described previously [39] .

The Brazilian samples were characterized at the University of São Paulo (Brazil). Viral RNA purifications were performed with QIAamp Viral RNA Mini Kit (Qiagen, Hilden, Germany), according to manufacturer's instructions, followed by RT-PCR using SuperScript 1 III One-

Step RT-PCR System with Platinum 1 Taq DNA Polymerase (Life Technologies, Carlsbad, CA, USA) and previously described primers [64] . Amplicons were visualized in 2% of SYBR Safe (Life Technologies, Carlsbad, CA, USA) low Melting Point Agarose (Life Technologies, Carlsbad, CA, USA). Products were purified using Illustra GFX PCR DNA and Gel Band Purification Kits (GE Health Care and Life Sciences, Buckinghamshire, England). DNA sequencing was performed with BigDye Terminator v3.1 Cycle Sequencing Kit (Life Technologies, Carlsbad, CA, USA) in ABI 3730XL DNA Analyzer (Applied Biosystems, Foster City, CA USA) at the LACTAD Facility at UNICAMP (http://www.lactad.unicamp.br/en).

All available complete fusion-gene (F-gene) sequences of class II NDV were downloaded from GenBank as of July 2015 [48] and aligned using ClustalW [65] , resulting in 1452 complete Table 1 . Collated Isolates from GenBank and SEPRL samples. A total of 54 isolates from the following taxonomic orders are tabulated below: Accipitriformes (n = 1); Anseriformes (n = 13); Charadriiformes (n = 3); Columbiformes (n = 23); Falconiformes (n = 1); Galliformes (n = 4); Passeriformes (n = 2); Pelecaniformes (n = 1); Phoenicopteriformes (n = 1); Psittaciformes (n = 4); Unknown (n = 1). GenBank accession numbers bolded are strains sequenced from this study. fusion protein gene sequences. Initial phylogenetic analyses were performed utilizing the complete F-gene sequences using the Neighbor Joining method based upon 100 bootstrap replicates [66] , as implemented in the MEGA version 6 software [67] (data not shown). To ensure that all viruses used in further analyses were isolated only from wild birds, we subjected the dataset to rigorous selection criteria. Specifically, viruses isolated from domestic non-poultry species such as waterfowl from live bird markets [68] , and sporting, racing, or pet birds where conspecifics may have been vaccinated [69] were excluded from the dataset. Viral sequences of wild bird isolates (n = 54, including the 24 sequenced in this study) ( Table 1 ) that were evolutionarily closely related to reference vaccine strains (n = 5) from genotypes I and II were parsed from the initial compilation. An additional 15 representative sequences from the remaining genotypes (III-XIV and XVI-XVIII) were also included, resulting in a final dataset of 74 sequences. Analysis of the best-fit substitution model was performed and the goodness-of-fit for each model was measured by the corrected Akaike Information Criterion (AICc) and the Bayesian Information Criterion (BIC) [67] . Tamura 3-parameter model with 500 bootstrap replicates was used for constructing the phylogenetic tree [70] . Finally, estimates of the means and pairwise genetic distances were computed using the Maximum Composite Likelihood method as implemented in MEGA6 [71] . The pairwise distances per decade were also calculated as the nucleotide distance per site divided by the number of years separating the isolation of each virus and the respective vaccine strain, multiplied by ten. The rate variation among sites was modeled with a gamma distribution (shape parameter = 4). For each statistical inference, codon positions consisting of 1st+2nd+3rd+Noncoding were retained, while gaps and/or missing data were trimmed.

In the past, both topology and nucleotide distance have been used to demonstrate that the sources (direct or indirect) of isolated viruses were live vaccines used in poultry [38, 46, 72] . It has been previously reported that a nucleotide distance of approximately 1% per decade is the natural rate of NDV evolution [73, 74] . In the present study, more stringent selection criterion to identify vaccine-derived viruses was used. NDV isolates with nucleotide distances lower than 0.1% per decade, when compared to the most closely related vaccine strains, were termed "vaccine-derived" viruses.

All statistical analyses were performed with SAS v. 9.3 [75] . Fisher's Exact Test for small sample sizes and dichotomous variables [76] (a 2 x 2 contingency table) was used to examine whether a significant relationship existed among the ages at which NDV infected Rock Pigeons were found shedding virus. For that purpose, the covariates age and shedding status were each divided into two classes, the hatch-year (HY) and after hatch-year (AHY), and shedding or not shedding, respectively. We next investigated whether a measurable physiological cost was present in the wild Rock Pigeons infected with vaccine-derived NDV. We used two widely accepted and correlated measures of body condition indices in passerines and near passerines to do so: i) weight-to-wing chord ratio, which is measured on a continuous range, and ii) fat scores, which are measured along an ordinal scale [77] . Prior to analyses, the distributions for fat and weightto-wing chord ratio were inspected to confirm that each variable met the assumption for a non-skewed distribution using the Shapiro Wilk's test statistic [78] .

Since the variable fat (S2 Table) did not meet the assumption of a normalized distribution, we utilized the non-parametric, one-way-ANOVA-by-ranks Kruskal-Wallis H test (S3 Table) [79] to determine whether the fat scores of shedding Rock Pigeons were significantly lower than the fat scores of non-shedding Rock Pigeons. While the variable weight-to-wing chord ratio had a normal distribution (S4 Table) , it was necessary to relax assumptions regarding the within-group variances due to the unequal sample sizes of shedding Rock Pigeons versus nonshedding Rock Pigeons [80] . Therefore, we used a one-way generalized linear model (GLM) (S5 Table) to examine whether shedding Rock Pigeons had a measurably smaller weight-towing chord ratio than their non-shedding counterparts.

The complete F-gene sequences (n = 24) of NDV obtained in this study were submitted to Gen-Bank and are available under the accession numbers from KU133351 to KU133365 and KU159667 to KU159675.

Cloacal and oral swabs (n = 72) and blood samples (n = 71) were collected from Rock Pigeons in Atlanta, Georgia. Rock Pigeon body condition indices, age, serological and virus shedding status are provided in S1 Table. Three birds were positive for NDV antibodies, equaling a serological study prevalence of 4.23%. No NDV was isolated from the serologically positive birds. Nine Rock Pigeon samples tested positive via HA for NDV shedding. None of the samples cross-reacted with the additional tested APMV serotypes (data not shown). The detected shedding of NDV from the Rock Pigeon samples suggests active replication (Table 2) . ICPI data of the HA positive samples are presented in Table 2 . The observed low ICPI values (0.03 to 0.34) indicated that the isolated viruses were of low virulence [59] .

Sequencing data from the F-gene analysis showed that the 24 isolates from wild birds in Brazil, Bulgaria, Ukraine and the USA contained fusion protein cleavage sites specific for NDV of low virulence with two basic amino acids between residue positions 113 and 116 and a leucine at position 117 ( 113 RQGR#L 117 ) [59] . As a result of the preliminary phylogenetic analysis including all available sequences from GenBank and the 24 sequences obtained in this study (data not shown), 54 viral sequences that were evolutionarily closely related to reference vaccine strains from genotypes I and II were selected for further analysis. The corresponding data about these selected viruses are presented in Table 1 . Utilizing established criteria, the 24 NDV sequences were classified as members of genotype II of class II [81] . Nine of them clustered together with the vaccine strain chicken/USA/LaSota/1946, while the other 15 isolates grouped with two additional ND vaccine strains of genotype II-chicken/USA/HitchnerB1/1947 and turkey/USA/VG/GA/1989 (Fig 1) . Fig 1 illustrates the distribution of the isolates from wild birds that are genetically closely related to vaccine strains of NDV genotype I and II. Within genotype II, 28 isolates from wild birds originating from Argentina, Bulgaria, China, India, Mexico and Ukraine clustered in a monophyletic branch with the vaccine strain chicken/USA/LaSota/1946. The mean genetic distance between this vaccine strain and the isolates within the branch was considerably low (0.2%). Within the same genotype, another 19 wild bird isolates from Brazil, China and USA grouped together with vaccine strains chicken/USA/HitchnerB1/1947 and turkey/USA/VG/ GA/1989. The genetic analysis showed that the latter isolates of wild-bird origin displayed genetic distance of 0.3-0.4% when compared to these two vaccine strains. Similarly, seven wild-bird NDV isolates of genotype I of class II also showed close phylogenetic relationship with vaccine strains. Five Mexican isolates of wild-bird origin from 2009 clustered with the vaccine strain PHY-LMV42/1966, while two 2009 and 2011 Chinese isolates of duck origin grouped closely together in the phylogram with the vaccine strain chicken/Australia/Queensland/V-4/1966 (Fig 1) . The mean genetic distance between the viruses and the vaccine strains within each of these genotype I branches was 0.2%. Pairwise distance and pairwise distance per decade results among viruses within of the described groups are presented in S7, S8 and S9 Tables. Based on the above results and the pairwise distances per decade (from 0.009% to 0.093%), the 54 viruses from these groups were determined as "vaccine-derived".

Our data indicate that at least 17 species from ten avian orders occupying different habitats excreted vaccine-derived NDV ( Table 1 ). The most frequent isolations occurred in the orders Columbiformes (n = 23) and Anseriformes (n = 13), with occasional isolations from eight other orders. The oldest isolation detected occurred in 1997 and the most recent in 2014 with samples from free-ranging (n = 47) and captive (n = 7) birds. The obtained results demonstrate that vaccine-derived NDV were detected in eight countries on four different continents. Lastly, some of the vaccine-derived viruses were obtained from wild birds representing declining species or species at imminent risk of decline according to international standards and conservation working groups (Table 1 ). These species of special concern include the Highland Guan (Penelopina nigra), Red-lored Amazon (Amazona autumnalis), Yellow-naped Amazon (Amazona auropalliata), Vinaceous-breasted Amazon (Amazona vinacea), Chilean Flamingo (Phoenicopterus chilensis), and an endemic Psittacine from India, a country where all Psittacines are facing declines [82] .

A total of 70 Rock Pigeons were assessed in quantifying whether age served as an explanatory variable for infection with vaccine-derived NDV (as the age of two captured Rock Pigeons could not be determined). Fisher's Exact Test for small sample sizes and dichotomous variables was employed to determine whether the Rock Pigeon "age class" was a significant predictor of NDV infection. The results (Fisher's Exact Test: p = 0.0217) provide evidence that HY birds are statistically more susceptible to infection with the vaccine-derived NDV in this population (S6 Table) .

We hypothesized that Rock Pigeons infected with the vaccine-derived viruses would fall into a poorer body condition class than their non-shedding conspecifics. Only the fat scores and the weight-to-wing chord ratio values of HY birds were used in these analyses as only HY birds were infected with vaccine-derived viruses. HI titers were not added as a covariate to these analyses since all shedding birds were negative for NDV antibodies. We initially proposed that shedding HY birds would have a lower mean ranked fat score than non-shedding HY birds. We determined that there was no statistical difference in visible subcutaneous fat between shedding (n = 9) and non-shedding (n = 36) HY individuals (Kruskal-Wallis, H = 0.1117, p = 0.7382) (S3 Table) . The variable weight-to-wing chord ratio is a well-accepted metric of body condition that is also frequently applied in avian disease ecology field studies (55, 56) . We found no statistical difference in the weight-to-wing chord ratio indices between shedding birds and non-shedding birds (F = 0.17, df = 1, 42, p = 0.6844) (S5 Table) .

The unexpectedly low genetic distances per decade (from 0.009% to 0.093%) between the vaccine strains and those isolates studied here demonstrate the global presence spanning at least 18 years of vaccine-derived ND viruses in wild birds (S7, S8 and S9 Tables). Although our data do not allow for the identification of the direct sources of NDV infection for the studied wild birds, it is reasonable to suggest that these wild bird isolates originated from recent spillovers of live NDV vaccines, instead of representing strains that naturally circulate in these birds for the following reasons: i) if the newly isolated viruses had been naturally circulating and originated from vaccine [70] . The tree with the highest log likelihood (-108983.3717) is shown. A discrete Gamma distribution was used to model evolutionary rate differences among sites (4 categories (+G, parameter = 0.0936). The rate variation model allowed for some sites to be strains (originally isolated in the 1940s and the 1960s) through natural evolution they would have changed significantly for the last 4 to 6 decades and would have presented much higher nucleotide distances than the ones determined in the pairwise distance analysis (S7, S8 and S9 Tables); ii) preservation of live NDV in the environment, unchanged for such a long period of time, is highly unlikely due to the thermal and biological lability of NDV [83] ; and iii) for the majority of the isolates studied here that are almost genetically identical, there is no evidence of a direct epidemiological link, neither geographical nor temporal, between their isolations. Spillover of NDV vaccines into wild birds reflects the most commonly used live NDV vaccines. Four different types of vaccine-derived viruses have been identified in at least 17 wild bird species from 10 different orders. The identification rates of LaSota-and B1-like viruses were substantially higher compared to the rest of the vaccines (28 and 19, respectively, out of a total of 54), corresponding with the most widely utilized vaccines, LaSota and Hitchner B1 [84] [85] [86] . In addition, LaSota, the most pathogenic of the commonly used NDV vaccines [87] , is more likely to be shed in the environment. Spillovers of other NDV vaccines that are more limitedly employed were identified at lower rates (PHY-LMV42 and V4), or were not found at all (NDW, I-2, F, Clone-30, Ulster) [42, 59, 86] .

The presence of vaccine-derived NDV in non-target species is likely to be underestimated in surveillance and passive diagnostic studies. For instance, the recovery of vaccine-derived viruses from samples submitted to diagnostic laboratories is usually neglected since only the identification of pathogenic or virulent isolates are a priority for animal health. Newcastle disease vaccines, although often ignored, may behave similarly to more widely studied vaccines. There are multiple live vaccines that are widely used in human and animal medicine that are documented as being shed from inoculated patients (varicella, vaccinia, polio, distemper) into the environment [88] [89] [90] [91] . Considering that NDV has the capacity to infect at least 250 species of birds and all avian species are considered susceptible [42, 92] , these vaccines have the potential to be easily transmitted across species undetected.

Although our results are preliminary, the high number of cases and the biological and behavioral characteristics of the order Columbiformes (comprising primarily Rock Pigeons) [93, 94] suggest that these birds have a high degree of association with vaccine-derived viruses, and that they may be used as sentinels for spillover of live NDV vaccines. Furthermore, the synanthropic nature of many birds from the order Columbiformes, as well as their worldwide distribution [95] [96] [97] and presence in anthropogenic affected areas, increase the probability for cross-host transmission of NDV between poultry and wild birds. The hatch-year cohort of Rock Pigeons (Atlanta, GA, USA) represented 100% of all incidences of identification of vaccine-derived viruses within that population. This may be associated with the higher vulnerability to disease of juveniles due to their increased movement requirements for foraging, in concert with naïve immune systems [98, 99] . In the studied dataset, the order Anseriformes also had a high incidence of shedding of vaccine-derived NDV. The high rate of detection of vaccine-derived NDV from this order in Asia may reflect a higher degree of contact between wild birds and domestic members of this order as a result of the larger populations of ducks and geese used as poultry in this continent. evolutionarily invariable ([+I], 39.7777% sites). The tree is drawn to scale with branch lengths measured in the number of substitutions per site and the percentage of trees in which the associated taxa clustered together are shown below the branches. The analysis involved 81 nucleotide sequences with a total of 1662 positions in the final dataset. Isolates studied in this work are designated in front of the taxa name as follows: USA-•; Ukraine-; Brazil-□, Bulgaria-■. Evolutionary analyses were conducted in MEGA6 [67] . The Roman numerals presented in the taxa names in the phylogenetic trees represent the respective genotype for each isolate, followed by the GenBank identification number, host name (if available), country of isolation, strain designation and country of isolation. Although the analyses of fat scores and weight-to-wing chord ratios among the studied wild Rock Pigeons did not demonstrate a negative physiological impact derived from infection with vaccine-derived NDV, the ecological effects have yet to be addressed in other species. Newcastle disease virus vaccines are substantially beneficial for the control of the disease, but they have also been documented to have some mild to moderate negative side effects in poultry [33, [100] [101] [102] . The lack of recognition of adverse events or insufficient sampling efforts in wild avifauna [18, 19] may be the reasons why reports of such side effects are absent in the literature. An additional outcome of this study is that species reported to shed vaccine-derived NDV are listed by the International Union for Conservation of Nature as taxa of conservation concern. Some of the vaccine-derived viruses were obtained from wild birds representing declining species or species at imminent risk of decline and the ecological implications of this finding are currently unknown (Table 1) .

Our analyses did not utilize a random dataset, and the possibility of sampling bias cannot be discarded. It is noteworthy to recognize that some degree of reporting bias may have influenced our results, as we are limited to analyzing self-submitted data [103] . The nature of epidemiological source data implies that convenience samples are often all that are available, and can be valuable in cross-sectional, observational studies [104] [105] [106] . Molecular epidemiology metaanalyses with access to convenience datasets comprised of GenBank genomic-associated data are widespread and useful tools. For example, a limited dataset of self-reported isolates from GenBank (n = 97) was used to pinpoint the source and direction of spread for Methicillin-resistant Staphylococcus aureus within hospitals of Florida [107] . Sequences from convenience samples of Leishmania infantum from Brazil (n = 45) submitted to GenBank were evaluated to determine that no statistical association existed between the host species (canids or humans) and genetic variability within the causative agent [108] . The emergence of severe acute respiratory syndrome (SARS) within densely populated Hong Kong prompted a study in which the only available dataset was a series of convenience samples of GenBank isolates from recently acquired clinical cases from China (n = 168). The authors not only identified two co-circulating SARS viral clusters, but also identified the host travel method for the viral isolate detected in North America [109] .

It is difficult to quantify the magnitude of reporting bias since it is unknown how many birds have historically been positive for vaccine-derived viruses, yet were not reported to Gen-Bank. The analysis of 54 isolates is a small sample size from which to draw inferences about the potential impact of spillovers of NDV vaccines; however, the referenced studies did not indicate that vaccine viruses were specifically targeted for sampling as opposed to wild-type NDV, and it was concluded that each vaccine-derived virus was an incidental research discovery. We utilized a convenience sample to make a causal inference. In public health and epidemiology, causal inference and risk management often use the best available data to identify when intervention is feasible and necessary [110] .

Further studies are necessary to evaluate if transmission of these and other vaccines or infectious agents from poultry operations into free-ranging avifauna would have significant ecological impact.

Supporting Information S1 

Diseases at the livestock-wildlife interface: Status, challenges, and opportunities in the United States

One Health in the shrinking world: Experiences with tuberculosis at the human-livestock-wildlife interface

Emerging infectious diseases of wildlife-Threats to biodiversity and human health

Bridge hosts for avian influenza viruses at the wildlife/domestic interface: an eco-epidemiological framework implemented in southern Africa

Agricultural intensification, priming for persistence and the emergence of Nipah virus: a lethal bat-borne zoonosis

Veterinary aspects of ecological monitoring: the natural history of emerging infectious diseases of humans, domestic animals and wildlife. Tropical animal health and production

Unhealthy landscapes: Policy recommendations on land use change and infectious disease emergence. Environmental Health Perspectives

Infectious animal diseases: the wildlife/livestock interface. Revue Scientifique et Technique-Office International des Epizooties

Assessing disease risk at the wildlife-livestock interface: A study of Sierra Nevada bighorn sheep. Biological Conservation

Wildlife and Emerging Zoonotic Diseases: The Biology, Circumstances and Consequences of Cross-Species Transmission. Current Topics in Microbiology and Immunology

Global trends in emerging infectious diseases

Zoonosis emergence linked to agricultural intensification and environmental change

Cross-Species Virus Transmission and the Emergence of New Epidemic Diseases

Diseases of humans and their domestic mammals: pathogen characteristics, host range and the risk of emergence

Pathogen spillover in disease epidemics

Success factors for avian influenza vaccine use in poultry and potential impact at the wild bird-agricultural interface

Ecology of avian influenza viruses in a changing world

Comparison of pathogens in broiler and backyard chickens on the Galapagos Islands: Implications for transmission to wildlife

Disease emergence in birds: challenges for the twenty-first century

Animal Migration and Infectious Disease Risk

Poultry Plants Lowering Production Costs and Increasing Variety

Globalisation: a study of the poultry-meat supply chain

Impacts of poultry vaccination on viruses of wild bird. Current Opinion in Virology

Current status of vaccines against infectious bursal disease

The long view: 40 years of infectious bronchitis research

Molecular epidemiology of infectious laryngotracheitis: a review

Antigenic differences among Newcastle disease virus strains of different genotypes used in vaccine formulation affect viral shedding after a virulent challenge

Challenges and recent advancements in infectious laryngotracheitis virus vaccines

Advancement in vaccination against Newcastle disease: recombinant HVT NDV provides high clinical protection and reduces challenge virus shedding with the absence of vaccine reactions

Attenuated vaccines can recombine to form virulent field viruses

Imperfect Vaccination Can Enhance the Transmission of Highly Virulent Pathogens

Epidemiology and prevention of vaccine-preventative diseases

The use of vaccination in poultry production

Phylogenetic Relationships among Virulent Newcastle Disease Virus Isolates from the 2002-2003 Outbreak in California and Other Recent Outbreaks in North America

Pathogenesis of Newcastle Disease in chickens experimentally infected with viruses of different virulence

A summary of taxonomic changes recently approved by ICTV

A molecular epidemiological investigation of isolates of the variant avian paramyxovirus type 1 virus (PPMV-1) responsible for the 1978 to present panzootic in pigeons

Identification of new sub-genotypes of virulent Newcastle disease virus with potential panzootic features

Food and Agriculture Organization of the United Nations (FAO)

Newcastle Disease. Developments in Veterinary Virology

Diseases of Poultry

Newcastle disease in cormorants. The Canadian Veterinary Journal

Phylogenetic Diversity among Low-Virulence Newcastle Disease Viruses from Waterfowl and Shorebirds and Comparison of Genotype Distributions to Those of Poultry-Origin Isolates

Newcastle disease virus in penguins from King George Island on the Antarctic region

Characterization of Newcastle Disease Viruses in Wild and Domestic Birds in Luxembourg from

Molecular epidemiology of Newcastle disease in Mexico and the potential spillover of viruses from poultry into wild bird species

United States Department of Agriculture

Decision and finding of no significant impact: Environmental assessment-reducing pigeon, starling, and sparrow damage through an integrated wildlife damage managment program in the state of Georgia. USDA/APHIS/WS, School of Forestry and Natural Resources

Pseudoreplication and the Design of Ecological Field Experiments. Ecological Monographs

The Costs, Effects, and Host Response to Parasitism in a Model Cosmopolitan Avian Host Species: Rock Pigeon (Columba livia

University of Georgia

Identification guide to North American birds. Part I

Geographical Variation in Blood Parasites in Feral Pigeons: The Role of Vectors

Variation of West Nile Virus Antibody Prevalence in Migrating and Wintering Hawks in Central California

Protocol: Instructions for the Establishment and Operation of Constant-Effort Bird Banding Stations as Part of the Monitoring Avian Productivity and Survivorship (MAPS) Program. Point Reyes Station, CA: The Institute for Bird Populations

Evaluation of a commercial enzyme-linked immunosorbent assay for detection of antibodies against the H5 subtype of Influenza A virus in waterfowl. Influenza and Other Respiratory Viruses

Role of Poultry in the Spread of Novel H7N9 Influenza Virus in China

Manual of Diagnostic Tests and Vaccines for Terrestrial Animals: Mammals, Birds and Bees. 1, Part 2

Newcastle disease virus (strain Herts 33/56) in tissues and organs of chickens infected experimentally

A Laboratory Manual for the Isolation and Identification of Avian Pathogens

Evidence for a new avian paramyxovirus serotype 10 detected in rockhopper penguins from the Falkland Islands

Effects of Newcastle disease virus vaccine antibodies on the shedding and transmission of challenge viruses

International Biological Engagement Programs Facilitate Newcastle Disease Epidemiological Studies

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

The neighbor-joining method: a new method for reconstructing phylogenetic trees

MEGA6: Molecular Evolutionary Genetics Analysis Version 6.0

Exchange of Newcastle disease viruses in Korea: The relatedness of isolates between wild birds, live bird markets, poultry farms and neighboring countries. Infection, Genetics and Evolution

Newcastle Disease in Captive Falcons in the Middle East: A Review of Clinical and Pathologic Findings

Estimation of the number of nucleotide substitutions when there are strong transition-transversion and G+C-content biases

Prospects for Inferring Very Large Phylogenies by Using the Neighbor-Joining Method

Development of an improved vaccine evaluation protocol to compare the efficacy of Newcastle disease vaccines

A longitudinal study of velogenic Newcastle disease virus genotypes isolated in Italy between

The occurrence of five major Newcastle disease virus genotypes

Fisher's Exact Test

An Evaluation of Condition Indices and Predictive Models for Noninvasive Estimates of Lipid Mass of Migrating Common Yellowthroats, Ovenbirds, and Swainson's Thrushes

An Analysis of Variance Test for Normality (Complete Samples)

The Kruskal-Wallis Test and Stochastic Homogeneity

Generalized Linear Models with Varying Dispersion

Genetic diversity of avian paramyxovirus type 1: Proposal for a unified nomenclature and classification system of Newcastle disease virus genotypes. Infection, Genetics and Evolution

Updates in the IUCN Red List of Threatened Birds of India

Newcastle Disease Viruses Causing Recent Outbreaks Worldwide Show Unexpectedly High Genetic Similarity to Historical Virulent Isolates from the 1940s

Guest editorial: Historical note on the origin of the LaSota strain of Newcastle disease virus

Guest editorial: Serendipity in science: Discovery of the B-1 Strain of Newcastle disease virus

Control of Newcastle Disease by Vaccination

Comparison of the immunogenicity of Newcastle disease virus strains V4, B1 and LaSota in chickens

Spread of Vaccine-Derived Poliovirus from a Paralytic Case in an Immunodeficient Child: an Insight into the Natural Evolution of Oral Polio Vaccine

How Contagious Is Vaccinia?

Transmission of Vaccine Strain Varicella-Zoster Virus from a Healthy Adult with Vaccine-Associated Rash to Susceptible Household Contacts

Distemper encephalitis in pups after vaccination of the dam

Newcastle disease and other avian paramyxoviruses. Revue Scientifique et Technique-Office International des Epizooties

Artificial selection, naturalization, and fitness: Darwin's pigeons revisited

Involvement of free-flying wild birds in the spread of the viruses of avian influenza, Newcastle disease and infectious bursal disease from poultry products to commercial poultry. World's

The role of wildlife in transboundary animal diseases

Diseases of poultry and endemic birds in Galapagos: implications for the reintroduction of native species

Stability of flock composition in urban pigeons

Parasite mediated mortality and host immune response explain age-related differences in blood parasitism in birds

Competition for food in urban pigeons: the cost of being juvenile. The Condor

An Ulster 2C strain-derived Newcastle disease vaccine: vaccinal reaction in comparison with other lentogenic Newcastle disease vaccines

Increased Mortality in Turkeys Selected for Increased Body Weight Following Vaccination with a Live Newcastle Disease Virus Vaccine

Revue scientifique et technique (International Office of Epizootics)

Sampling uncertainty in medical research

Comparison of Convenience Sampling and Purposive Sampling

Essentials of Infectious Disease Epidemiology

Sampling Considerations for Disease Surveillance in Wildlife Populations. The Journal of Wildlife Management

Molecular epidemiology of community-associated methicillin-resistant Staphylococcus aureus in the genomic era: a cross-sectional study

Genetic homogeneity among Leishmania (Leishmania) infantum isolates from dog and human samples in Belo Horizonte Metropolitan Area (BHMA), Minas Gerais, Brazil. Parasite and Vectors

Molecular epidemiology of the novel coronavirus that causes severe acute respiratory syndrome

Causal inference in public health

Successful establishment and global dispersal of genotype VI avian paramyxovirus serotype 1 after cross species transmission. Infection

Molecular Characterization and Phylogenetic Analysis of Newcastle Disease Virus Isolates from Healthy Wild Birds

Pathotypical Characterization and Molecular Epidemiology of Newcastle Disease Virus Isolates from Different Hosts in China from 1996 to

Phylogenetic and pathogenic analysis of Newcastle disease virus isolated from house sparrow (Passer domesticus) living around poultry farm in southern China

Surveillance and molecular characterization of Newcastle disease virus in seafowl from coastal areas of China in 2011

The authors greatly appreciate the contribution of Dr. Vanessa Ezenwa from the University of Georgia for the surveillance and sampling of Rock Pigeons in Atlanta, GA. The authors also acknowledge Tim Olivier of SEPRL, and Matthew Breithaupt, Shreyas Vangalas and Allison Bradwell from the Odum School of Ecology for their superb technical support.The mention of trade names or commercial products in this publication is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the U.S. Department of Agriculture. The USDA is an equal opportunity provider and employer.