key: cord-253436-dz84icdc
authors: Wille, Michelle; Muradrasoli, Shaman; Nilsson, Anna; Järhult, Josef D.
title: High Prevalence and Putative Lineage Maintenance of Avian Coronaviruses in Scandinavian Waterfowl
date: 2016-03-03
journal: PLoS One
DOI: 10.1371/journal.pone.0150198
sha: 
doc_id: 253436
cord_uid: dz84icdc

Coronaviruses (CoVs) are found in a wide variety of wild and domestic animals, and constitute a risk for zoonotic and emerging infectious disease. In poultry, the genetic diversity, evolution, distribution and taxonomy of some coronaviruses have been well described, but little is known about the features of CoVs in wild birds. In this study we screened 764 samples from 22 avian species of the orders Anseriformes and Charadriiformes in Sweden collected in 2006/2007 for CoV, with an overall CoV prevalence of 18.7%, which is higher than many other wild bird surveys. The highest prevalence was found in the diving ducks—mainly Greater Scaup (Aythya marila; 51.5%)—and the dabbling duck Mallard (Anas platyrhynchos; 19.2%). Sequences from two of the Greater Scaup CoV fell into an infrequently detected lineage, shared only with a Tufted Duck (Aythya fuligula) CoV. Coronavirus sequences from Mallards in this study were highly similar to CoV sequences from the sample species and location in 2011, suggesting long-term maintenance in this population. A single Black-headed Gull represented the only positive sample from the order Charadriiformes. Globally, Anas species represent the largest fraction of avian CoV sequences, and there seems to be no host species, geographical or temporal structure. To better understand the eitiology, epidemiology and ecology of these viruses more systematic surveillance of wild birds and subsequent sequencing of detected CoV is imperative.

Coronaviruses (CoVs) are found in a wide variety of animals in which they can cause respiratory, enteric, hepatic, and neurological disease of varying severity. In humans, CoVs cause a large fraction of common colds [1] , as well as more rare, but serious disease such as the outbreak of Severe Acute Respiratory Syndrome (SARS, caused by SARS-CoV) in 2003, and the polymerase (RdRp) were generated using methods described in [14] or [15] , and sequences generated in this study have been deposited in GenBank under accession numbers KT882615-28. Resulting sequences were aligned using the MAFFT algorithm [16] implemented in Geneious R7 (Biomatters, New Zealand). The nucleotide substitution model was determined in MEGA 5.2 [17] and maximum likelihood trees were constructed using PhyML [18] with aLRT branch support in SeaView [19] . Trees were projected using FigTree v1.4 [20] .

A total of 764 birds from 11 species of Anserifomes (ducks, geese, swans, n = 691) and 11 species of Charadriiformes (gulls, terns, shorebirds, n = 143) were sampled in this study. Overall, the prevalence of CoV was high, at 18.7%, but species, groups and orders were unevenly represented. Diving ducks had the highest prevalence (39%), driven by high prevalence in Greater Scaup (Aythya marila; 51.5%), however the sample size was small (n = 37). Dabbling ducks of the genus Anas also had high prevalence; Mallard, in particular, had a prevalence of 19.2%. While Anseriformes had the highest prevalence, in the order Charadriiformes, CoV was detected in a Black-headed Gull (Chroicocephalus ridibundus), but absent in tern and wader species sampled (Table 1) . Eleven sequenced Mallard CoV and three sequenced Scaup CoV were gammacoronaviruses.

Putative maintenance, despite no global spatial, temporal, host species patterns A global phylogeny of all sequenced avian gammacoronaviruses (not IBV-like) RdRp show no spatial, temporal or host species differentiation between clades. That is, all clades contain sequences from multiple geographic locations, different years of sampling, and different host species (Fig 1A, S1 and S2 Figs) . Indeed, viruses sequenced in this study were placed in a clade containing sequences from Hong Kong, China, Beringia, the United States of America and Sweden (Fig 1B and S2 Table) . However, despite limited global patterns of clade differentiation, 8/10 sequences from Mallards in 2007 were most closely related (within 98-99% identity) to CoV sequences from Mallards at Ottenby in 2011 (Fig 1B and S1 Table) , suggesting local circulation of these RdRp. Mallard CoV sequence 69998 had the highest pairwise identity to viruses from Hong Kong, but is nested in the broad and highly similar clade containing the largest number of sequences from Ottenby in 2011 ( Fig 1B) . A single Mallard CoV sequence (69740) was in a differentiated clade, most similar to sequences from Hong Kong and more distantly China (Fig 1B and S2 Table and S2 Fig) . Sequences from Scaup CoV were not similar to Mallard CoV; Scaup 67699 was similar to sequences from Beringia and China (within 98% identity; S1 Table) . Two Scaup CoV sequences (67693 and 67703) were highly similar to only one sequence in Genbank, which was collected from a Tufted Duck (Aythya fuligula) in Hong Kong (JN788847). All other sequences were less than 92% similar, indicating a differentiated lineage that is poorly sampled (Fig 1B and S1 

In this study we aimed to further contribute to the natural history of avian coronaviruses by screening an array of waterbirds for these viruses. We found a prevalence of 18.7% CoV, which is higher than the 0-15% reported previously in wild bird studies [11, 14, 15, [21] [22] [23] [24] [25] [26] [27] . This may be due to the high proportion of dabbling ducks, particularly Mallard, in our study, and the temporal and spatial features of the dabbling duck sampling. More specifically, sampling of dabbling ducks occurred at a migratory stopover site, and thus high prevalence could be driven by migrants which not only import viruses through infected individuals, but also replenish the pool of susceptible individuals across the migratory period, resulting in high prevalence at these locations [28, 29] . Directly comparing prevalence estimates to previous studies, however, is challenging due to the array of different methods utilized, ranging from conventional PCR to multiplex real-time PCR, in addition to species distributions and seasonal variations. Regardless, this study corroborates previous studies suggesting the importance of dabbling ducks (Anas sp) as hosts of avian coronaviruses [11] [12] [13] [14] [15] [21] [22] [23] [24] [25] , but also highlighting the importance of diving ducks. The high diving duck prevalence in this study was driven by high prevalence in Greater Scaup, which could be due to congregation in very large flocks facilitating transmission, but could also be due to other factors such as a local variation in disease prevalence or reflect outbreak dynamics wherein the virus rapidly expands across the susceptible population. Overall, diving ducks have only been sampled in this study and by Chu et al. (2011) , so it is unclear whether this pattern of high prevalence is present globally. We found a low prevalence in waders, which is in contrast to [14] and [21] , which found a very high prevalence (>20%) in waders, but corroborates findings in [22, 25, 26] . This may be due to differences in congregation and feeding patterns of waders between the locations, but has to be interpreted with caution, as numbers of waders sampled in our study, and others, were small. Finally, while assessment of gulls is limited, this study corroborates findings with Muradrasoli et al (2010) indicating that this group requires further scrutiny. In order to better assess CoV dynamics, resampling the same site across time is imperative, and this is the first study to do so. Coronaviruses from Mallards were previously assessed using samples collected at Ottenby in 2011 [15] and we detected high similarity between viruses from 2007 and those from 2011. Indeed, 8/10 CoV were most closely related to sequences from 2011, which strongly suggests maintenance and/or location transmission of these RdRp lineages at this location. This is unlikely to be a coincidence as Swedish Mallard CoV make up less than 1/3 of those available sequences, and all current RdRp clades are comprised of sequences from Asia. The waterfowl host appears important in these samples, however assessing host species biases is challenging due to the overwhelming number of sequences from Anas species. However, coronaviruses from diving ducks were rather different from existing diversity as demonstrated by the similarity between two of three Scaup CoV sequences from this study and Tufted Duck CoV sequence from Hong Kong. The limited number of sequences in this clade could be due to sampling shortcomings, or that currently developed primers do not amplify this clade effectively.

Despite few studies, small samples sizes and differences in prevalence, what is clear, is that in the Northern Hemisphere waterfowl species, especially dabbling and diving ducks are important in the epidemiology of avian CoVs. It is interesting to note that these patterns are very similar to those found in low pathogenic influenza A viruses: high prevalence in waterfowl and gulls in the Northern Hemisphere [30] , and little host species and temporal structuring within waterfowl derived viruses in the conserved polymerase genes (such as PB2, PB1) [31] . Further, detection in fecal or cloacal samples suggest these viruses may also replicate in the gastrointestinal tract, which is the true for Turkey Coronavirus [9] . As demonstrated by years of influenza sampling [32] , to better understand the eitiology, epidemiology and ecology of these viruses more systematic surveillance needs to be undertaken and more viruses need to be sequenced, particularly full genomes. Given the importance of coronaviruses as human pathogens, as exemplified by SARS and MERS-CoV and the potential for wild birds as reservoirs, spreaders and mixing vessels, further studies of coronaviruses in wild birds are warranted. 

Human coronaviruses: what do they cause? Antivir Ther

Coronaviruses: important emerging human pathogens

Bats are natural reservoirs of SARS-like coronaviruses

Middle East respiratory syndrome (MERS): a new zoonotic viral pneumonia. Virulence

Virus Taxonomy: classification and nomenclture of viruses: Ninth Report of the Internaltional Committee on Taxonomy of Viruses

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

Discovery of seven novel mammalian and avian coronaviruses in the genus Deltacoronavirus supports bat coronaviruses as the gene source of alphacoronavirus and betacoronavius and avian coronavirus as the gene source of gammacoronavirus and deltacoronavirus

Evolutionary insights into the ecology of coronaviruses

Turkey coronavirus is more closely related to avian infectious bronchitis virus than to mammalian coronaviruses: a review

Pathogenicity of turkey coronavirus in turkeys and chickens

Avian coronavirus in wild aquatic birds

Genomic Analysis and Surveillance of the Coronavirus Dominant in Ducks in China

Broadly targeted multiprobe QPCR for detection of coronaviruses: coronavirus is common among mallard ducks (Anas platyrhynchos)

Prevalence and phylogeny of coronaviruses in wild birds from the Bering Strait area (Beringia)

Temporal dynamics, diversity, and interplay in three components of the viriodiversity of a Mallard population: Influenza A virus, avian paramyxovirus and avian coronavirus

Multiple Alignment of DNA Sequences with MAFFT

Molecular Evolutionary Genetics Analysis using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods

Estimating maximum likelihood phylogenies with PhyML

SeaView Version 4: A Multiplatform Graphical User Interface for Sequence Alignment and Phylogenetic Tree Building

FigTree v1.1.1: Tree figure drawing tool

Diverse gammacoronaviruses detected in wild birds from Madagascar

Detection and molecular characterization of infectious bronchitis-like viruses in wild bird populations

Genetically diverse coronaviruses in wild bird populations of northern England

Molecular identification and characterization of novel coronaviruses infecting graylag geese (Anser anser), feral pigeons (Columbia livia) and mallards (Anas platyrhynchos)

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

Surveillance of avian coronaviruses in wild bird populations of Korea

Absence of coronaviruses, paramyxoviruses, and influenza A viruses in seabirds in the southwestern Indian Ocean

Animal migration and infectious disease risk

Juveniles and migrants as drivers for seasonal epizootics of avian influenza virus

Global patterns of influenza a virus in wild birds

The evolutionary genetics and emergence of avian influenza A viruses in wild birds

Spatial, temporal, and species variation in prevalence of influenza A viruses in wild migratory birds

We wish to thank and the duck trappers at Ottenby Bird Observatory and Jonas Waldenström for collecting and providing samples used in this study, Jonas Blomberg for kindly providing sequence, Mallard CoV sequences generated in this study are indicated with a filled circle and Scaup CoV sequences with an asterisk. Scale bar represents number of substitutions per site. Traditional projection of panel A with support values is presented in S2 Fig.