key: cord-0997644-0e2809hr
authors: Barbosa, Carla M.; Durigon, Edison L.; Thomazelli, Luciano M.; Ometto, Tatiana; Marcatti, Roberta; Nardi, Marcello Shiavo; de Aguiar, Daniel M.; Pinho, João Batista; Petry, Maria Virginia; Neto, Isaac Simão; Serafini, Patrícia; Rodrigues, Roberta Costa; de Azevedo Junior, Severino Mendes; Góes, Luiz Gustavo B.; de Araujo, Jansen
title: Divergent coronaviruses detected in wild birds in Brazil, including a central park in São Paulo
date: 2019-03-15
journal: Brazilian Journal of Microbiology
DOI: 10.1007/s42770-019-00065-7
sha: 2851e74df0bf3b3169b499537c05b5e52afb32fd
doc_id: 997644
cord_uid: 0e2809hr

Coronaviruses are single-stranded positive-sense RNA viruses associated with important avian diseases. Their relatively high rates of mutation and recombination frequencies allow them to adapt to new hosts and ecological niches. Although Brazil has 18% of global avian species diversity, studies regarding the presence of avian viral diseases in wild birds in South America are scarce. In this study, we performed a retrospective analysis of the presence of CoVs in 746 wild birds. Oropharyngeal and cloacal swabs were obtained and placed together in vials containing VTM transport medium collected in different regions of Brazil between 2006 and 2013. Screening for viral nucleic acid was performed using conventional RT-PCR and pancoronavirus nested PCR. Positive samples were characterized by partial sequencing of the RNA-dependent RNA polymerase (RdRp) gene, and ensuing phylogenetic analysis was performed to investigate the association between virus epidemiology and bird migration routes. Coronavirus RNA were detected and sequenced from six samples, in which three were related to gammacoronaviruses group and the other three to deltacoronavirus group. Our study documents the presence of CoVs related to avian gamma- and deltacoronaviruses circulating in both urban- and poultry-farm regions of Brazil, implicating wild birds as potential carriers of CoVs which may represent a risk to poultry farms and public health in Brazil. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (10.1007/s42770-019-00065-7) contains supplementary material, which is available to authorized users.

Coronaviruses are RNA single-stranded positive polarity virus associated with important disease in avian poultry [1] .

Infectious bronchitis virus (IBV), known as coronavirus of chickens, causes important economic losses to the poultry industry. Studies about the presence of CoVs in wild birds, especially IBV, in South America are scarce [2] . Coronaviruses

Electronic supplementary material The online version of this article (https://doi.org/10.1007/s42770-019-00065-7) contains supplementary material, which is available to authorized users.

(CoVs) possess the largest genomes among RNA viruses and have relatively high mutation rates and recombination frequencies when compared with other single-strand RNA viruses, which may allow them to adapt to new hosts and ecological niches [3] [4] [5] . These characteristics are exemplified by the two emerged highly pathogenic betacoronaviruses, CoV-SARS (severe acute respiratory syndrome) and the CoV-MERS (Middle East respiratory syndrome) [6] .

Even though avian coronaviruses with zoonotic potential have not been reported yet, many species of wild and domestic birds can carry CoVs responsible for substantial economic loss within the poultry industry such as the gammacoronavirus infectious bronchitis virus (IBV). This virus affects the gastrointestinal, respiratory, urinary, and reproductive systems and may lead to mortality rates as high as 30% and sudden drops in egg production [7] . Moreover, many novel CoVs genetically distinct from IBV have been recently identified in different families of wild and domestic birds and mammals but their potential economic losses remain unknown. Recently, Durães-Carvalho et al. (2015) , described the detection of delta-and betacoronaviruses in captive wild birds, from a rehabilitation center, originally captured in São Paulo and Mato Grosso do Sul state with close evolutionary relationships with CoVs isolated from mammals [8] .

Wild migratory birds may facilitate the spread of viruses, as demonstrated by the surveillance of avian influenza virus [9, 10] and West Nile virus [11] [12] [13] , but the circulation of CoV in avian migratory routes in South America, more specifically in Brazil, still need to be elucidated. The potential diversity of CoV in Brazilian bird fauna as consequence of the bird species diversity (18%), together with the large number of endangered species and the high intensity of poultry farms in Brazil, indicates that birds from Brazilian biomes are important subjects for CoV surveillance [2, 14] .

In this study, we performed a retrospective analysis, molecular detection, and phylogenetic analysis of avian coronaviruses using oropharyngeal and cloacal swab samples from wild and domestic birds collected between 2006 and 2013 in different regions of Brazil to investigate the association between virus epidemiology, bird migration routes, and the proximity of urban areas and poultry farms.

In total, 746 samples (Table 1) of birds captured between 2006 and 2013 in five expeditions to north, northeast, midwest, south, and southeast Brazil (Fig. 1) were analyzed. For selection, the migratory routes were investigated through a map analysis and experience of ornithologists [15] . The samples were collected by field team of the Laboratory Clinical and Molecular Virology, University of São Paulo (USP), São Paulo, Brazil, with support of the Wildlife Department of São Paulo (DEPAVE); University of Vale do Rio dos Sinos (UNISINOS); University of Pernambuco (UFRPE); National Center for Research and Conservation of Wild Birds (CEMAVE); Laboratory of Bird Ecology, Institute of Biosciences (UFMT); and Laboratory of Virology and Rickettsioses, Veterinary Hospital (UFMT). Samples were collected from wild birds after capture with mist nets, strategically placed between 30 and 250 cm of the ground, and by hand nets on an open field in public parks. The field works were made according to the permission of the Brazilian Institute of Environment and Renewable Natural Resources (IBAMA) and National Center for Research and Conservation of Wild Birds (CEMAVE) developed and followed protocols approved by the Ethics Committee on Animal Experimentation of the Institute of Biomedical Sciences, University of São Paulo. Sampling from oropharyngeal and cloacal swabs was obtained and placed together in vials containing VTM transport medium [16] . All birds were then released after sampling.

Viral RNA was extracted using the MagMax TM-96 RNA Isolation Kit (Ambion, Austin, TX, USA) according to the manufacturer's instructions. The extracted RNA was eluted in 80 μL of RNase-free water and stored at − 80°C until processed. The cDNA was transcribed using the High-Capacity cDNA Archive Kit (Applied Biosystems, Foster City, CA, USA) following the manufacturer's instructions. The RT reaction was performed by adding, to each reaction, 10 μL of total RNA, 2 μL of RT Buffer 10X, 2 μL of Random Primers (50 ng/μL) (Thermo Fisher Scientific, Waltham, MA, USA), 4.2 μL of DEPC water, and 50 U of MultiScribe™ Reverse Transcriptase (Thermo Fisher Scientific). The RT step conditions were 25°C for 10 min, 37°C for 120 min, and 85°C for 5 min, and cooling to 4°C.

Polymerase chain reaction (PCR) followed by nested PCR were performed according to the protocol for the detection of avian CoVs developed by Chu et al. (2011) [17] , which targets the RNA-dependent RNA polymerase (RdRp) gene in the ORF1 area of the viral genome. The resulting approximately 600-bp PCR fragment and 440-bp nested PCR fragment were separated by agarose gel [1.5%] electrophoresis and viewed under UV light after staining with ethidium bromide (25 μg/100 mL).

Samples were purified using the ExoSap-IT enzyme (GE Healthcare, Chalfont St. Giles, UK) and GeneJET™ Gel Extraction Kit (Thermo Fisher Scientific) according to the Genetic characterization and phylogenetic tree of the RdRp gene Generated chromatograms containing nucleotide sequences were visualized using Bioedit Sequence Alignment Editor (v.7.1.5.0). Nucleotide sequences were aligned using CLC sequencer WorkBench version 7.6 (QIAGEN Aarhus, Denmark), and the nucleotide percentage similarities were calculated using the MegAlign program of the DNASTAR package (DNASTAR, Inc., Madison, WI, USA). The similarities with existing sequences were determined using the Basic Local Alignment Search Tool (BLAST). Phylogenetic analysis was performed using the neighbor-joining method (NJ) implemented in the MEGA 7.0.26 [18] . The nucleotide substitution model used for phylogenetic reconstruction was the p-distance model. The bootstrap values > 55% were obtained in the analysis of 10,000 replicates and are presented at the branching points. Samples related to gammacoronaviruses were obtained from one Calidris alba and one Calidris fuscicollis individual, respectively, related a restricted to sandpiper species (Fig. 2) . Samples related to deltacoronavirus were obtained from a black skimmer (Rynchops niger) carrying a strain that clustered, in phylogenetic analysis, closer to the clade of viruses found in Ibirapuera Park (Fig. 3 ).

The nested PCR assay with degenerate consensus primers was used to many studies. This assay has been used for the identification and differentiation of CoV in the four genera d e s c r i b e d t o d a t e ( a l p h a -, b e t a -, g a m m a -, a n d deltacoronavirus) [17, 19, 20] . Similarly to our report, most of avian CoV diversity studies detected the same coronavirus genera [2, 17, 21, 22] , with the exception of a report of betacoronavirus detection in distinct birds species from a rehabilitation center in Brazil [8] . Further studies are needed to understand the circulation of betacoronavirus in birds. Ibirapuera is the main park in São Paulo, receiving on average 150,000 visitors in a single weekend. Wild birds make up a large portion of the species of animals that inhabit the park (ca. 142 species). The birds are kept free and cohabiting with the human population and other wild species (resident or visitors). There are another 21 species of wild animals, including bats and skunks, and the park also allows pets [23] . Such close contact between so many different species may cause interspecific transmission and recombination events, leading to the appearance of emerging viruses [24] . The three goose samples were collected from geese that live in this same environment in contact with wild birds that use the park as a stopping point. On the other hand, three migratory birds captured in wintering areas at Lagoa do Peixe National Park-RS have shown positive results of CoVs too. Sandpipers are boreal migrants that migrate south from North America in the winter [25] . A CoV related to deltacoronavirus group was obtained from a black skimmer (Rynchops niger), thus highlighting the possibility that migratory birds may be carrying these viruses across migratory routes in Brazil [9] . The Lagoa do Peixe region is an important migratory route [25] , and it is located close (400 miles) from one of the most important region of poultry farm production of Brazil, more specifically region of Porto Alegre in Rio Grande do Sul State [26] . The proximity between the migratory routes of wild birds infected by Fig. 2 Phylogenetic analysis of partial RdRp gene sequence of avian gammacoronaviruses. Phylogenetic tree of partial sequences of the RdRp-encoding region constructed using the neighbor-joining method. In green, CoVs related to gammacoronaviruses sequences detected in samples from Lagoa do Peixe National Park, state of Rio Grande do Sul, southern Brazil. In red, CoVs related to gammacoronaviruses sequences detected in samples from Ibirapuera Park, São Paulo, southeast Brazil. Bootstrap values > 55% were obtained in the analysis of 10,000 replicates and are presented at the branching points gammacoronavirus shows the importance of a continuous surveillance and potential contact control between poultry and migratory birds as measures to avoid any chance of viral transmission between these animals.

The CoVs related to deltacoronaviruses group (DPV 5 and DPV 10) and gammacoronavirus group (DPV 16) from Ibirapuera Park showed the highest similarity to a CoV found by Lau et al. (2012) [27] in Guangdong province (84%), China, and by Lima et al. (2015) [28] (93% nt identity) in a migratory bird (Porphyrula alleni) from Madagascar in 2011, respectively ( Table 2 ). The gammacoronaviruses found at Lagoa do Peixe National Park were more similar to those found by Muradrasoli et al. (2010) [29] in sandpipers (Calidris spp.) from the Bering Strait (99% nt identity). However, the deltacoronavirus was more similar to the strain isolated from a black-headed gull (Chroicocephalus ridibundus) by Hepojoki et al. (2017) in Finland (95% nt identity) (Table 3) . Interestingly, this data were reported about sampling collected during the same years (2010-2013) by the same way, using RT-PCR targeting the RdRp gene to detect the presence of CoV in wild birds from Finland. The prevalence was found around 5.4% (51/939), which corroborates with our findings. Hepojoki et al. (2017) have shown positive samples to diverse gamma-and deltacoronavirus too [22] .

In South America, there is a lack of information about coronavirus in wild birds. These data are interesting, once that, present birds infected in central park in Sao Paulo city and wild birds through migratory routes. In Hong Kong, Chu et al. (2011) detected a high prevalence (12.5%) in aquatic wild birds, reporting gamma and deltacoronavirus, being the first one more prevalent [17] . Similar studies screened before in Sweden (2016), with the number of sampling similar to our study, and reported a high prevalence (18.7%) in wild birds, with emphasis for Anseriformes that had the highest prevalence (51%), all belonging to the gammacoronavirus group [30] . On the 2015) analyzed 303 samples from cloacal swabs, using universal CoV primers, and found only one duck coronavirus, belonging to the gammacoronavirus group too, which suggested a low prevalence in wild aquatic birds in USA, corroborating with our data [21] . As we know, migratory routes in Brazil are linked with flyways coming from North America and can confirm our findings [16] . Although our data have shown low prevalence when compared with studies before in Scandinavian Waterfowl, our data is clear when compared with North American data. In addition, it is the first time that samples related to Gamma and deltacoronavirus group in free wild birds were reported in Brazil.

The overall prevalence of CoV found in this study was 0.8% (6/746 birds), with 3.7% (3/81) from Ibirapuera Park and 1.6% (3/192) from migratory birds captured in Lagoa do Peixe National Park. Other studies report an overall prevalence of CoV in aquatic and migratory birds between 0.33 and 18.5% [8, 17, 21, 30, 31] . This data reinforces that CoVs are circulating across major migratory routes worldwide.

Our report points to the circulation of coronaviruses among migratory and resident birds sampled from distinct regions of Brazil. Coronavirus RNA were detected in birds species of Anatidae family (Chinese goose) in Ibirapuera Park, São Paulo, southeast Brazil. Anatidae birds were already described as hosts of coronavirus all over the globe. However, interestingly, we report for the first time the detection in two sandpiper species (Calidris alba and Calidris fuscicollis) from Lagoa do Peixe National Park, a wintering area in southern Brazil. Although our data indicate a low prevalence of coronavirus in birds sampled, the genetic similarity between coronavirus detected at birds from distinct continents may indicate that wild birds are possibly carrying coronavirus throughout migratory routes as observed for avian influenza virus. We report the detection of samples related to delta-and gammacoronaviruses in different species of birds, some sampled in a park located inside the largest city in South America, indicating the circulation of coronavirus in urban areas (Table 3 ). In addition, all birds tested in this study were in apparently healthy populations. Further studies are required to determine the species and full genome similarity with other avian CoVs. However, our findings indicate that CoVs may be circulating among one of the major avian migratory routes in the world.

Positive Control-Avian Infectious bronchitis virus strain Ma5

PNLP115/Rynchops_niger/coronavirus/Brasil

DPV_5/Anser_cygnoides/coronavirus/Brasil

JQ065049.1/Common-moorhen/coronavirus/HKU21

JQ065046.1/Magpie-robin/coronavirus/HKU18

DPV_10/Anser_cygnoides/coronavirus/Brasil

Global distributions and strain diversity of avian infectious bronchitis virus: a review

Phylodynamic analysis of avian infectious bronchitis virus in South America

The molecular biology of coronaviruses

Coronavirus diversity, phylogeny and interspecies jumping

Epidemiology, genetic recombination, and pathogenesis of coronaviruses

Interspecies transmission and emergence of novel viruses: lessons from bats and birds

Infectious Bronchitis. In: Saif YM (ed) Dis. Poult, 12th edn

Coronaviruses detected in Brazilian wild birds reveal close evolutionary relationships with Beta-and Deltacoronaviruses isolated from mammals

Avian influenza virus (H11N9) in migratory shorebirds wintering in the Amazon region

North Atlantic migratory bird flyways provide routes for intercontinental movement of avian influenza viruses

West Nile virus surveillance, Brazil

Prevalence of West Nile virus in migratory birds during spring and fall migration

Birds, migration and emerging zoonoses: West Nile virus, Lyme disease, influenza A and enteropathogens

Molecular surveillance of the Newcastle disease virus in domestic and wild birds on the north eastern coast and Amazon biome of Brazil

Relatório Anual de rotas e áreas de concentraçào de aves migratórias no Brasil

Migratory birds in southern Brazil are a source of multiple avian influenza virus subtypes

Avian coronavirus in wild aquatic birds

MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets

Biogenesis and dynamics of the coronavirus replicative structures

Genetic diversity of bats coronaviruses in the Atlantic Forest hotspot biome, Brazil

Identification of avian coronavirus in wild aquatic birds of the central and eastern USA

Prevalence and genetic diversity of coronaviruses in wild birds

Guia dos parques municipais de São Paulo

Coronaviruses post-SARS: update on replication and pathogenesis

Seasonal variation in shorebird abundance in the state of Rio Grande Do Sul, southern Brazil

Embrapa Suinos e Aves

Discovery of seven novel mammalian and avian coronaviruses in the genus Deltacoronavirus supports bat coronaviruses as the gene source of Alphacoronavirus and Betacoronavirus and avian coronaviruses as the gene source of Gammacoronavirus and Deltacoronavi

Diverse gammacoronaviruses detected in wild birds from Madagascar

Prevalence and phylogeny of coronaviruses in wild birds from the bering strait area (Beringia)

High prevalence and putative lineage maintenance of avian coronaviruses in Scandinavian waterfowl

Novel reassortant influenza A(H5N8) viruses, South Korea

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations

Conflict of interest The authors declared that they have no conflict of interest.Ethical approval All applicable international, national, and/or institutional guidelines for the care and use of animals were followed. The procedures involving wild birds were approved by the Brazilian Society