Molecular and epidemiological characterization of avian influenza viruses from gulls and dabbling ducks in Norway
Virology Journal volume 10, Article number: 112 (2013)
Wild aquatic birds constitute the natural reservoir for avian influenza viruses (AIVs). Separate Eurasian and American AIV gene pools exist. Here, the prevalence and diversity of AIVs in gulls and dabbling ducks in Norway were described. The influence of host species and temporal changes on AIV prevalence was examined. Five AIVs from Norway, including three from common gull (Larus canus), were analyzed along with 10 available AIV genomes from gulls in Eurasia to search for evidence of intracontinental and intercontinental reassortment of gene segments encoding the internal viral proteins.
Swabs collected from 2417 dabbling ducks and gulls in the south-west of Norway during five ordinary hunting seasons (August-December) in the period 2005–2010 were analyzed for presence of AIV. Multivariate linear regression was used to identify associations between AIV prevalence, host species and sampling time. Five AIVs from mallard (Anas platyrhynchos) (H3N8, H9N2) and common gull (H6N8, H13N2, H16N3) were full-length characterized and phylogenetically analyzed together with GenBank reference sequences.
Low pathogenic AIVs were detected in 15.5% (CI: 14.1–17.0) of the samples. The overall AIV prevalence was lower in December compared to that found in August to November (p = 0.003). AIV was detected in 18.7% (CI: 16.8–20.6) of the dabbling ducks. A high AIV prevalence of 7.8% (CI; 5.9–10.0) was found in gulls. A similar temporal pattern in AIV prevalence was found in both bird groups. Thirteen hemagglutinin and eight neuraminidase subtypes were detected. No evidence of intercontinental reassortment was found. Eurasian avian (non H13 and H16) PB2 or PA genes were identified in five reference Eurasian gull (H13 and H16) AIV genomes from GenBank. The NA gene from the Norwegian H13N2 gull isolate was of Eurasian avian origin.
The similar temporal pattern in AIV prevalence found in dabbling ducks and gulls, the relatively high virus prevalence detected in gulls and the evidence of intracontinental reassortment in AIVs from gulls indicate that gulls that interact with dabbling ducks are likely to be mixing vessels for AIVs from waterfowl and gulls. Our results support that intercontinental reassortment is rare in AIVs from gulls in Eurasia.
The natural reservoir for 16 hemagglutinin (HA) (H1-H16) and 9 neuraminidase (NA) (N1-N9) subtypes of influenza A viruses is wild aquatic birds, primarily ducks, gulls and shorebirds [1, 2]. Influenza A viruses that cause epidemics and pandemics in humans, mammals and domestic poultry evolve directly or indirectly from the natural reservoir of avian influenza viruses (AIVs) .
Most AIVs are adapted to wild waterfowl and do not cause apparent clinical disease in their hosts. AIVs replicate in the intestinal and respiratory tract, but their tropism varies between bird species . In ducks, virus transmission via the fecal-oral route through intake of virus-contaminated surface water is believed to be important , and infectious virus can persist in water and lake sediments for months [5, 6]. In the northern hemisphere, the AIV prevalence in dabbling ducks is high in the fall, whereas a lower prevalence is found in the spring . In gulls, the AIV prevalence is generally low [7–9], although higher virus prevalence has occasionally been detected in gull colonies [10, 11].
Due to limited contact between birds from the different continents, AIVs have evolved into Eurasian and American phylogenetic lineages. However, this distinction is not absolute. In general, intercontinental reassortment of virus genes occurs relatively infrequently [12–20], but is more commonly detected in AIVs originating from areas of the world where migratory flyways overlap [14, 15, 18]. AIVs can also move with wild birds that cross the Atlantic Ocean [19, 20]. Transmission of influenza virus genes between continents can lead to introduction and establishment of new virus lineages [19, 21]. So far, intercontinental reassortment has involved exchange of single or multiple gene segments . Transfer of complete AIV genomes has not been reported . Within continents, exchange of virus genes is far more common due to frequent interspecies AIV transmissions that result in the formation of influenza viruses with transient genome constellations .
AIVs of the H13 and H16 subtypes are primarily restricted to gulls and shorebirds [2, 22] and have originally been described to form separate American and Eurasian phylogenetic “gull-specific” lineages . It was recently proposed to combine these two lineages into a “Charadriiformes specific” lineage due to the apparent high frequency of genetic reassortment occurring between viruses from these lineages [19, 23]. Eurasian gene segments are commonly present in influenza A viruses isolated from gulls in North America, especially in the H13 and H16 subtypes, and it has therefore been suggested that long-distance migrating gull species play an important role in the transfer of AIV genes from Eurasia into America [18, 19]. On the contrary, American gene segments have only been found on a few occasions in AIVs isolated from gulls in Eurasia [13, 23]. However, the total number of fully sequenced AIV genomes from gulls in Eurasia is still low and more sequence information, particularly from coastal regions in Europe, has been requested .
Norway has the longest coastline in Europe and is located along the East-Atlantic flyway for migratory birds . As a consequence of the emergence of highly pathogenic AIV subtype H5N1 in Southeast-Asia and its spread towards Europe, a surveillance program for avian influenza in wild birds was implemented in Norway in 2005 . Except for a one-year gap in 2008, the program continued to 2010. In the present study, epidemiological analysis was performed on data obtained through five years of AIV surveillance. In addition, five virus isolates, including three from common gull (Larus canus) were sequenced. The aims of this study were to: i) describe the prevalence and subtype diversity of AIVs in dabbling ducks and gulls in the south-west (SW) of Norway, ii) test whether the virus prevalence was influenced by host species and sampling time, iii) determine the genetic similarity of the Norwegian AIVs to previously characterized AIVs, and to iv) search for evidence of intracontinental and intercontinental reassortment of gene segments encoding the internal viral proteins in AIVs from gulls in Eurasia.
Collection of samples
Swabs were collected by hunters from hunter-harvested dabbling ducks (n = 1709) and gulls (n = 708), in Rogaland County on the SW coast of Norway, during five licensed hunting seasons (2005–2007 and 2009–2010) from August to December. This is the most important area both for wintering wild aquatic birds and for poultry farming in Norway. The gulls were hunted regardless of the surveillance program. A part of the gulls were shot due to safety reasons to prevent bird strikes at the airport in Rogaland County, since the airport is situated in the middle of a wetland area. Others were shot since they were considered to be a problem for agriculture or in industrial and commercial areas. From each bird, cloacal and tracheal swabs were collected. Both swabs were placed in a single tube of virus transport medium (VTM) that was sent by postal mail to the Norwegian Veterinary Institute for analysis. The surveillance program for AIV was approved by the Norwegian Food Safety Authority (http://www.mattilsynet.no). An overview of the number of samples collected from August to December each year is presented in Figure 1.
RNA extraction and rRT-PCR
Swabs collected from 2005 to 2007 were previously analyzed [24, 25]. For samples collected in 2009 and 2010, RNA was isolated  and screened for AIV using real-time RT-PCR (rRT-PCR) targeting the matrix (M) gene, with modifications as described in . If found AIV positive, an H5-specific rRT-PCR was performed . Samples that tested positive for AIV, but negative for H5, were subtyped by sequencing of the HA2 part of the HA gene  and the full-length NA gene .
Virus isolation and hemagglutination assay
Virus isolation was performed on a selection of AIV positive swabs. The samples, each containing one ml VTM, were thawed, vortexed, briefly centrifuged and 12.5–25 μl Gentamicin (50 mg/ml, Gibco, Paisley, UK) was added to the supernatant. After one hour incubation at room temperature, 100–200 μl medium was inoculated into the allantoic cavity of nine to 11 days old embryonated chicken eggs in three parallels and incubated for four days at 37°C. A minimum of two passages were made. The hemagglutination assay was performed using 1% chicken erythrocytes in a PBS de Boer solution . The M gene specific rRT-PCR was used to confirm that the observed hemagglutination was caused by AIV . Subtyping was performed as described above [27, 28] and AIV positive allantoic fluid was stored at -80°C.
One H13, one H16 and one non-H13/H16 isolate from common gull were chosen for further characterization: A/common gull/Norway/1602/2009(H6N8), A/common gull/Norway/1313/2009(H13N2) and A/common gull/Norway/1617/2006(H16N3). In addition, two mallard (Anas platyrhyncos) isolates i.e. A/mallard/Norway/779/2009(H3N8) and A/mallard/Norway/1537/2009(H9N2) were randomly selected as references of Eurasian AIVs. Allantoic fluid from first (H6, H13), second (H3, H16) and third (H9) egg passages was used. The five virus isolates represent AIVs of subtypes frequently circulating in dabbling ducks and gulls in Norway.
Sequencing of AIV genomes
For full-length amplification of the viral segments, a two-step RT-PCR was performed, using SuperScript® III (Invitrogen, Carlsbad, CA, USA) and the Advantage 2 PCR kit (Clontech, Mountain View, CA, USA). Primers  and cycling conditions  were used as previously described. As full-length amplification of all polymerase basic protein 1 (PB1) and polymerase basic protein 2 (PB2) genes was not successful, a new set of primers were designed for this purpose. The PCR amplicons were purified, using Zymoclean™ Gel DNA recovery kit (Zymo Research, Orange, CA), and cloned into pCR2.1 TOPO vectors according to the instructions of the manufacturer (Invitrogen). Sequencing was performed by using a commercial sequencing service (GATC Biotech, Constance, Germany). Primers for internal sequencing of the largest amplicons, PB2 and PB1, were designed. For each viral segment, three to six clones were sequenced. The sequences have been submitted to the EMBL Nucleotide Sequence Database (EMBL) [EMBL: HE802704-HE802743]. All primer sequences designed in this study are available upon request.
Similarity searches and phylogenetic analysis
For the AIV genomes sequenced in this study, the sequences were assembled using Vector NTI Advance ver.11 (Invitrogen). Similarity searches were performed using the Basic Local Alignment Search Tool (BLAST) .
A search in the Influenza Virus Resource database at NCBI  for complete AIV genomes from gulls in Europe and Asia identified 13 genomes (as of February 2012) that could be included in the phylogenetic analysis along with the five virus genomes sequenced in the present study (Additional file 1: Table S1). Three virus isolates from the NCBI search were excluded from the analysis due to subtype (HPAI H5N1), suboptimal sequence quality (A/slaty-backed gull/Japan/6KS0191/2006(H4N8) or because it clustered in a separate genetic lineage very different from all the other sequences (A/slaty-backed gull/Japan/6KS0185/2006(H4N8)) .
To search for evidence of reassortment of the internal viral gene segments (PB2, PB1, PA, NP, M and NS), MEGA 4 was used to align the full-length nucleotide sequences by Clustal W and to create phylogenetic trees, using the Neighbour-Joining method with the Maximum composite likelihood model . Bootstrap values of 1000 replicates were used to assess the nodal support. Reference sequences of American avian [GenBank: CY004904, CY004906-10, CY016144], American gull [GenBank: CY005070-71, CY005390, CY004562-63, CY003895, CY01498], Eurasian avian [GenBank: CY060363-65, CY060367, AB274986-87, CY060370] and Eurasian gull [GenBank: AY684874, AY684878, AY684883, M30753, EU580589, EU580578] AIVs were included in the analysis to evaluate the phylogeographical affiliation of the viral segments . Eurasian and American gull AIVs were defined as AIVs of the H13 and H16 subtypes, whereas Eurasian and American AIVs were defined as all AIV subtypes, except H13 and H16 .
Data (bird species and time of sampling) from AI surveillance in dabbling ducks and gulls in Rogaland from 2005 to 2007 [24, 25] and from 2009 to 2010 were included in the statistical analysis. The overall AIV prevalence was estimated for each bird species and for the two bird groups; dabbling ducks and gulls. For these bird groups, the virus prevalence was also studied according to year and month. For each month, the data from 2005 to 2010 were pooled. Annual variations in AIV prevalence during fall (August to December) were studied in detail in three species of dabbling ducks i.e. mallard, common teal (Anas crecca) and Eurasian wigeon (A. penelope), and also in the two most sampled gull species, i.e. herring gull (Larus argentatus) and common gull. Descriptive statistics were performed using Excel 2007 (Microsoft, Redmond, WA). Furthermore, multivariate logistic regression was used to describe the associations between AIV prevalence and host species or bird group and year, month or sampling date and their intercepts. R  was used for this analysis and p-values larger than 0.05 were considered as significant. Models were compared using Akaike information criteria (AIC). Binomial 95% confidence intervals (CIs) were calculated for the AIV prevalence estimates.
AIV prevalence and subtypes
From 2005 to 2010 a total of 15.5% (CI: 14.1–17.0) of the 2417 birds tested positive for AIV by rRT-PCR (Table 1). Comparison of AIC values for different multivariate logistic regression models showed that bird group (dabbling ducks or gulls), year of sampling and December as sampling month or not were the variables that best described the AIV status of each bird. No interaction of the variables improved the AIC value, hence both bird groups showed similar pattern in AIV prevalence based on annual and monthly variation. Sampling date did not explain the data, probably due to clustered sampling periods. The prevalence of AIV was lower in December (p = 0.003) and varied among years with less number of AIV positive birds in 2006 (p < 0.001) and 2007 (p = 0.049) as compared to other years (Figures 2A,B). The overall AIV prevalence was higher in dabbling ducks than in gulls (p < 0.001) with 18.7% (CI: 16.8–20.6) and 7.8% (CI: 5.9–10.0), respectively.
A multivariate logistic regression was also fitted to the dabbling ducks only. Duck species, December or not as sampling month and year gave the best model based on AIC value. Less birds were AIV positive in December (p = 0.008). Moreover, the AIV prevalence was lower in dabbling ducks in 2006 (p < 0.001) and 2007 (p = 0.027) than the other years. AIV was most frequently detected in common teal, closely followed by mallard, with no significant differences between the two species (Figure 2C, Table 1). The prevalence of AIV was significantly lower (p < 0.001) in Eurasian wigeon compared to that in both mallard and common teal (Table 1). The annual fluctuations in virus prevalence were similar for all the three species of dabbling ducks (Figure 2C). In gulls, the highest prevalence was detected in black-headed gull (Chroicocephalus ridibundus) (Table 1), however, the number of samples collected from this species was low. Most samples from gulls were collected from herring gull and common gull (Table 1). In herring gull and common gull, no significant difference in annual virus prevalence was observed (Figure 2D). A logistic regression with only these two gull species showed no significant difference between species, sampling date, month or year.
Of the 374 AIV positive samples, 223 were successfully subtyped. Of these, 105 samples were both HA and NA subtyped, whereas 109 and 9 samples were only HA or NA subtyped, respectively. Thirteen HA subtypes were detected. H6 was most frequently found, followed by H5, in 40.2% (CI: 33.6–47.1) and 13.1% (CI: 8.9–18.4) of the HA characterized samples, respectively. H3 and H9 were both found in 9.8% (CI: 6.2–14.6) of the samples. In dabbling ducks, H1-H9, H11 and H12 were detected, whereas in gulls, H4, H6, H13 and H16 were found (Figure 3). The H4 and H6 subtypes were found in both dabbling ducks and gulls (Figure 3). Phylogenetic analysis of the H4 and H6 genes showed no significant differences between the viruses from ducks and gulls (data not shown). The highest subtype diversity was found in mallard and common teal (data not shown). The H5 subtype was detected in 7.8% (CI: 5.3–10.9) (29/374) of the AIV positive samples and was found in mallard (n = 22) and common teal (n = 7). All H5 viruses had PQRETR/GLF or PQKETR/GLF HA cleavage sites and were therefore considered low pathogenic. H7 was only detected in one mallard, but its cleavage site was not determined. Eight NA subtypes (N1-N6, N8-N9) were detected. All of these were found in ducks, but only N1-N3, N6 and N8 were found in gulls. The N2 and N8 subtypes were most common and occurred in 66.7% (CI: 57.2–75.2) and 13.2% (CI: 7.6–20.1) of the NA subtyped samples, respectively.
The gene segments from the H3N8, H6N8 and H9N2 isolates were similar to Eurasian AIV gene segments (Additional file 2: Table S2). For the H13N2 and H16N3 viruses from common gull, the majority of genes were highly similar to Eurasian gull AIVs. The NA segment of the H13N2 virus was of Eurasian avian origin. For the PB2, NP and NA genes of the H16N3 virus, the results were more complicated to interpret, since highest identities were shared with the PB2, NP and NA genes of H13 and H16 AIVs isolated from American shorebirds containing gene segments previously proposed to have originated from Eurasian gull AIVs [18, 19, 35].
Phylogenetic analysis of the five virus isolates sequenced in this study showed that all internal gene segments were of Eurasian origin, and neither these isolates nor the reference AIVs from Eurasian gulls contained evidence of genetic reassortment between Eurasian and American AIVs (Figure 4). Evidence of reassortment between Eurasian avian and Eurasian gull AIVs was observed in two of the PB2 reference sequences [GenBank: GQ907309 and GQ907301], and in three of the PA reference sequences [GenBank: GQ907315, GQ907323 and JF775475]. No linkage between the reassortment of the PB2 and PA genes of these viruses was found. The Eurasian gull PB1, PA and M sequences were distinct from the American gull sequences and more closely related to the Eurasian avian sequences. In contrast, the evolutionary distance for American and Eurasian gull AIVs PB2, NP and NS was shorter, particularly for the NS gene (Figure 4). The NS genes of the five Norwegian virus isolates all belonged to allele A.
AIV prevalence and subtypes
In this study, the AIV prevalence was concomitantly monitored in both dabbling ducks and gulls over several months and years in the SW of Norway. The overall AIV prevalence was 15.5%. Statistical analysis showed that the AIV prevalence was significantly lower in December than during the other months of sampling (August to November). It is not clear whether ducks are a source of AIV for gulls (or the opposite) or whether two separate transmission cycles occur in ducks and gulls, with only sporadic virus transmissions between these two groups of birds. Based on the multivariate analysis of our complete data set, no interactions i.e. between bird group and year or month, respectively, were found to be statistically significant. These results therefore suggest that the AIV prevalence follows the same annual and monthly patterns from August to December in both dabbling ducks and gulls. This further indicates that a certain exchange of AIVs between these two bird groups probably occurs during this period. The high virus prevalence observed during August to November each year is most likely a result of extensive interspecies mixing and virus transmission, combined with gathering of a large number of immunological naïve juvenile birds . The relatively cold climate in Norway may also increase the environmental persistence of these viruses . It is unknown if the individuals tested were resident, overwintering or migrating birds, and information about the age of the birds included in this study was not available.
Similar virus prevalences were found in common teal (23.8%) and mallard (19.2%), whereas a significantly lower AIV prevalence was found in Eurasian wigeon (8.7%). The largest subtype diversity was found in mallard and common teal, which is in accordance with previous reports . In dabbling ducks in general, there were no significant differences in virus prevalence within the period August to November. Other studies have reported on a peak in AIV prevalence during October-November [7, 25, 36]. In dabbling ducks, host traits like surface-feeding and a high density of filtering plates on the internal surface of the beak increase the exposure to AIV . Mallards, common teals and Eurasian wigeons are known to gather in large flocks during fall and this could facilitate interspecies transmission of virus. The Eurasian wigeon is more terrestrial compared to the mallard and the common teal, and this behavior could lower the risk of getting infected via virus-contaminated water for this species.
In this study, a relatively large number of gulls (n = 708) were tested for AIV, and the overall prevalence of AIV in gulls was high (7.8%) compared to prevalence values reported in other studies [7–9]. The AIV prevalence in gulls has previously been reported to peak in late summer and early fall . One possible explanation for the high prevalence of AIV in gulls may be virus spill-over from ducks to gulls due to partially overlapping habitat use during fall. As mentioned above, our complete data set indicated that the AIV prevalence in dabbling ducks and gulls follows the same annual and monthly patterns from August to December. There was a tendency for a decrease in virus prevalence in December in gulls, however no statistically significant differences in AIV prevalence between months of sampling were found when analyzing the data set from gulls alone. The high virus prevalence in gulls could be caused by spill-over of H4 and H6 AIVs from ducks, as these subtypes were detected in both dabbling ducks and gulls. No significant genetic differences between the H4 and H6 viruses from ducks and gulls were found. The H6 subtype is known to have a broad host-range  and was the most frequently detected HA subtype in this study.
Thirteen HA subtypes, i.e. all except H10, H14, H15, and eight NA subtypes, i.e. all except N7, were identified. During the same period, H10  and N7 were detected in other parts of Norway, whereas H14 and H15 have never been found (personal observations). In the present study, H5, H6, H9 and H3 were the most frequently detected HA subtypes, whereas N2 and N8 were the most commonly detected NA subtypes. Several attempts were made to subtype the samples from 2005 to 2007. In 2009 and 2010, only one attempt per sample was performed. In this study, a H5 rRT-PCR specific for detection of Eurasian H5 viruses was used  which means that it is uncertain whether this test could detect American H5 AIVs potentially introduced by migratory birds into Eurasia. However, the HA gene from positive samples that tested negative in the H5 specific rRT-PCR were amplified and sequenced using generic primers targeting HA2  allowing for possible identification of H5 from American lineages. In addition, non-specific NA primers that amplify the N1, N2, N4, N5 and N8 genes were often the first choice used for NA subtyping instead of NA subtype specific primers that detect N3, N6, N7 and N9 . This may have introduced a bias in the relative detection rate of the different HA and NA subtypes. However, the occurrence of the different HA and NA subtypes is in accordance with previous results from Northern Europe [7, 13].
The similarity searches showed that the majority of gene segments from the five complete AIV genomes in our study were of Eurasian origin and that they were most similar to North European AIV sequences. The PB2, NP and NA genes of the Norwegian H16N3 isolate were most similar to AIV genes from American H13 and H16 viruses isolated from shorebirds and gulls. However, these American virus isolates have previously been described as reassortant viruses containing Eurasian gene segments [19, 20, 35]. The NA and NP genes from the Norwegian H16N3 isolate were also highly similar to the N3 and NP genes of A/glaucous gull/Alaska/44198-027/2006(H16N3), for which Ramey et al. (2010) were not able to determine certain continental affiliations . The NA gene of the Norwegian H16N3 virus was also similar to A/glaucous-winged gull/Southcentral Alaska/9JR0783R0/2009(H16N3). Wille et al. (2011) reported that the N3 genes of the Alaskan gull viruses mentioned above were more similar to sequences from Eurasian viruses than to those from other gull viruses in America. However, they were still described as distinct from the most closely related Eurasian lineage, and this was suggested to be caused by few available H16N3 virus sequences or by the presence of a unique lineage of the N3 gene in western North America .
The phylogenetic analysis showed that for the PB1, PA and M genes, the Eurasian gull sequences were distinct from the American gull sequences and more closely related to the Eurasian avian sequences. In contrast, for PB2, NP and NS, the evolutionary distance between American and Eurasian gull AIVs was shorter, particularly for the NS gene. This confirms the findings reported by Fouchier et al. (2005) who found that influenza A viruses isolated from gulls could be distinguished from other AIVs based on the PB2, NP, and NS genes but not necessarily based on the PB1, PA, and M genes .
No evidence of intercontinental gene reassortment between Eurasian and American AIVs was found in the internal genes of the Eurasian gull viruses analyzed in the present study. These results confirm the findings by Wille et al. (2011), who found evidence of transcontinental gene transfer in American gull AIVs, but not in European gull AIVs . The apparent one-way transfer of viral segments from Eurasia to America should be further studied in Eurasian gull AIV isolates, particularly from pelagic, far migrating gull species, since they are considered to be most important for intercontinental movement of AIV genes . Recently, H6N1 and H5N2 AIVs that contained American avian M genes were isolated from herring gulls in Belgium . The three virus isolates from gulls sequenced in this study were all from common gull. This species is not known to migrate across the Atlantic Ocean , and therefore, it might be a less likely long-distance AIV vector, compared to pelagic gull species like for instance the black-legged kittiwake (Rissa tridactyla) and the great black-backed gull (Larus marinus) [20, 39]. Some long-distance migratory shorebird species within Scolopacidae, for example Calidris sp. and ruddy turnstone (Arenaria interpres), may also contribute to the intercontinental gene reassortment between American and Eurasian AIVs.
The similarity search showed that the N2 gene of the Norwegian H13N2 gull virus was most similar to Eurasian avian viruses, which indicates that intracontinental reassortment i.e. reassortment between a Eurasian avian and a Eurasian gull H13 AIV has occurred. A similar finding was reported in a H13N2 AIV, isolated from great black-backed gull in Newfoundland, where the NA gene was most similar to American avian viruses .
Eurasian avian PA and PB2 genes were observed in half of the Eurasian gull genomes obtained from GenBank. All reassortant viruses were previously described as Mongolian H13 and H16 viruses . Recently, Eurasian avian PA segments were also observed in two H13N6 viruses isolated in Belgium from yellow-legged gull (Larus michaellis) and herring gull, respectively . Together this indicates that multiple infections of gulls with both Eurasian avian and Eurasian gull viruses are common. Gulls that frequently interact with waterfowl probably act as mixing vessels for Eurasian avian and Eurasian gull AIVs. Alternatively, but maybe less likely, viruses mix in other species and spill-over back to gulls. Especially segment 3, encoding the polymerase protein PA, seems to be frequently transferred from Eurasian avian to Eurasian gull AIVs [23, 40], which could indicate that H13 and H16 AIVs are not dependent of gull adapted PA genes to replicate.
To summarize, a high AIV prevalence was found in both dabbling ducks and gulls in the SW of Norway. Multiple AIV subtypes were present within a relatively small geographical area. Only low pathogenic viruses were found and most viral segments were highly similar to the genes of contemporary Eurasian viruses. Intracontinental reassortment events were identified in H13 and H16 AIVs from gulls. A similar temporal pattern in AIV prevalence was found in dabbling ducks and gulls. This together with the high virus prevalence in gulls and the intracontinental reassortment suggests that gull species frequently interacting with dabbling ducks are likely to be mixing vessels for AIVs from waterfowl and gulls. Evidence of intercontinental reassortment was not found, which supports that these events are rare in Eurasian gull AIVs. However, this apparent one-way transfer of virus segments from Eurasia to America should be further studied.
Webster RG, Bean WJ, Gorman OT, Chambers TM, Kawaoka Y: Evolution and ecology of influenza A viruses. Microbiol Rev 1992, 56: 152-179.
Fouchier RAM, Munster V, Wallensten A, Bestebroer TM, Herfst S, Smith D, Rimmelzwaan GF, Olsen B, Osterhaus ADME: Characterization of a novel influenza A virus hemagglutinin subtype (H16) obtained from black-headed gulls. J Virol 2005, 79: 2814-2822. 10.1128/JVI.79.5.2814-2822.2005
Jourdain E, van RD, Munster VJ, Kuiken T, Waldenstrom J, Olsen B, Ellstrom P: The pattern of influenza virus attachment varies among wild bird species. PLoS One 2011, 6: e24155. 10.1371/journal.pone.0024155
Webster RG, Yakhno M, Hinshaw VS, Bean WJ, Murti KG: Intestinal influenza: replication and characterization of influenza viruses in ducks. Virology 1978, 84: 268-278. 10.1016/0042-6822(78)90247-7
Brown JD, Goekjian G, Poulson R, Valeika S, Stallknecht DE: Avian influenza virus in water: infectivity is dependent on pH, salinity and temperature. Vet Microbiol 2009, 136: 20-26. 10.1016/j.vetmic.2008.10.027
Nazir J, Haumacher R, Ike AC, Marschang RE: Persistence of avian influenza viruses in lake sediment, duck feces, and duck meat. Appl Environ Microbiol 2011, 77: 4981-4985. 10.1128/AEM.00415-11
Munster VJ, Baas C, Lexmond P, Waldenstrom J, Wallensten A, Fransson T, Rimmelzwaan GF, Beyer WE, Schutten M, Olsen B: Spatial, temporal, and species variation in prevalence of influenza A viruses in wild migratory birds. PLoS Pathog 2007, 3: e61. 10.1371/journal.ppat.0030061
Olsen B, Munster VJ, Wallensten A, Waldenstrom J, Osterhaus AD, Fouchier RA: Global patterns of influenza a virus in wild birds. Science 2006, 312: 384-388. 10.1126/science.1122438
Gaidet N, Ould El Mamy AB, Cappelle J, Caron A, Cumming GS, Grosbois V, Gil P, Hammoumi S, de Almeida RS, Fereidouni SR: PLoS One. 2012, 7: e46049. 10.1371/journal.pone.0046049
Velarde R, Calvin SE, Ojkic D, Barker IK, Nagy E: Avian influenza virus H13 circulating in ring-billed gulls (Larus delawarensis) in southern Ontario, Canada. Avian Dis 2010, 54: 411-419. 10.1637/8808-040109-Reg.1
Tønnessen R, Valheim M, Rimstad E, Jonassen CM, Germundsson A: Experimental incoulation of chickens with gull-derived low pathogenic avian influenza virus subtype H16N3 causes limited infection. Avian Dis 2011, 55: 680-685. 10.1637/9701-030411-ResNote.1
Wallensten A, Munster VJ, Elmberg J, Osterhaus AD, Fouchier RA, Olsen B: Multiple gene segment reassortment between Eurasian and American lineages of influenza A virus (H6N2) in Guillemot (Uria aalge). Arch Virol 2005, 150: 1685-1692. 10.1007/s00705-005-0543-8
Lebarbenchon C, Chang CM, Gauthier-Clerc M, Thomas F, Renaud F, van der Werf S: H9N2 avian influenza virus in a Mediterranean gull. J Mol Genet Med 2008, 3: 121-123.
Krauss S, Obert CA, Franks J, Walker D, Jones K, Seiler P, Niles L, Pryor SP, Obenauer JC, Naeve CW: Influenza in migratory birds and evidence of limited intercontinental virus exchange. PLoS Pathog 2007, 3: e167. 10.1371/journal.ppat.0030167
Koehler AV, Pearce JM, Flint PL, Franson JC, Ip HS: Genetic evidence of intercontinental movement of avian influenza in a migratory bird: the northern pintail (Anas acuta). Mol Ecol 2008, 17: 4754-4762. 10.1111/j.1365-294X.2008.03953.x
Dugan VG, Chen R, Spiro DJ, Sengamalay N, Zaborsky J, Ghedin E, Nolting J, Swayne DE, Runstadler JA, Happ GM: The evolutionary genetics and emergence of avian influenza viruses in wild birds. PLoS Pathog 2008, 4: e1000076. 10.1371/journal.ppat.1000076
Manzoor R, Sakoda Y, Mweene A, Tsuda Y, Kishida N, Bai GR, Kameyama K, Isoda N, Soda K, Naito M: Phylogenic analysis of the M genes of influenza viruses isolated from free-flying water birds from their Northern Territory to Hokkaido, Japan. Virus Genes 2008, 37: 144-152. 10.1007/s11262-008-0248-7
Ramey AM, Pearce JM, Ely CR, Sheffield Guy LM, Irons DB, Derksen DV, Ip HS: Transmission and reassortment of avian influenza viruses at the Asian-North American interface. Virology 2010, 406: 352-359. 10.1016/j.virol.2010.07.031
Wille M, Robertson GJ, Whitney H, Bishop MA, Runstadler JA, Lang AS: Extensive geographic mosaicism in avian influenza viruses from gulls in the northern hemisphere. PLoS One 2011, 6: e20664. 10.1371/journal.pone.0020664
Wille M, Robertson GJ, Whitney H, Ojkic D, Lang AS: Reassortment of American and Eurasian genes in an influenza A virus isolated from a great black-backed gull (Larus marinus), a species demonstrated to move between these regions. Arch Virol 2011, 156: 107-115. 10.1007/s00705-010-0839-1
zu Dohna H, Li J, Cardona CJ, Miller J, Carpenter TE: Invasions by Eurasian avian influenza virus H6 genes and replacement of the virus’ North American clade. Emerg Infect Dis 2009, 15: 1040-1045. 10.3201/eid1507.090245
Hinshaw VS, Air GM, Gibbs AJ, Graves L, Prescott B, Karunakaran D: Antigenic and genetic characterization of a novel hemagglutinin subtype of influenza A viruses from gulls. J Virol 1982, 42: 865-872.
Van Borm S, Rosseel T, Vangeluwe D, Vandenbussche F, van den Berg T, Lambrecht B: Phylogeographic analysis of avian influenza viruses isolated from Charadriiformes in Belgium confirms intercontinental reassortment in gulls. Arch Virol 2012, 157: 1509-1522. 10.1007/s00705-012-1323-x
Jonassen CM, Handeland K: Avian influenza virus screening in wild waterfowl in Norway, 2005. Avian Dis 2007, 51: 425-428. 10.1637/7555-033106R1.1
Germundsson A, Madslien KI, Hjortaas MJ, Handeland K, Jonassen CM: Prevalence and subtypes of influenza A viruses in wild waterfowl in Norway 2006-2007. Acta Vet Scand 2010, 52: 28. 10.1186/1751-0147-52-28
Slomka MJ, Pavlidis T, Banks J, Shell W, McNally A, Essen S, Brown IH: Validated H5 Eurasian real-time reverse transcriptase-polymerase chain reaction and its application in H5N1 outbreaks in 2005-2006. Avian Dis 2007, 51: 373-377. 10.1637/7664-060906R1.1
Phipps LP, Essen SC, Brown IH: Genetic subtyping of influenza A viruses using RT-PCR with a single set of primers based on conserved sequences within the HA2 coding region. J Virol Methods 2004, 122: 119-122. 10.1016/j.jviromet.2004.08.008
Hoffmann E, Stech J, Guan Y, Webster RG, Perez DR: Universal primer set for the full-length amplification of all influenza A viruses. Arch Virol 2001, 146: 2275-2289. 10.1007/s007050170002
WHO Manual of Animal Influenza Diagnosis and Surveillance. http://www.who.int/vaccine_research/diseases/influenza/WHO_manual_on_animal-diagnosis_and_surveillance_2002_5.pdf
BLAST: Basic local alignment search tool. http://blast.ncbi.nlm.nih.gov/
Influenza virus resource. http://www.ncbi.nlm.nih.gov/genomes/FLU/Database/nph-select.cgi?go=genomeset
Bui VN, Ogawa H, Xininigen , Karibe K, Matsuo K, Awad SS, Minoungou GL, Yoden S, Haneda H, Ngo LH: H4N8 subtype avian influenza virus isolated from shorebirds contains a unique PB1 gene and causes severe respiratory disease in mice. Virology 2012, 423: 77-88. 10.1016/j.virol.2011.11.019
Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol 2007, 24: 1596-1599. 10.1093/molbev/msm092
R Core Team: R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2012.
Pearce JM, Ramey AM, Ip HS, Gill RE Jr: Limited evidence of trans-hemispheric movement of avian influenza viruses among contemporary North American shorebird isolates. Virus Res 2010, 148: 44-50. 10.1016/j.virusres.2009.12.002
Wallensten A, Munster VJ, Latorre-Margalef N, Brytting M, Elmberg J, Fouchier RA, Fransson T, Haemig PD, Karlsson M, Lundkvist A: Surveillance of influenza A virus in migratory waterfowl in northern Europe. Emerg Infect Dis 2007, 13: 404-411. 10.3201/eid1303.061130
Hill NJ, Takekawa JY, Cardona CJ, Ackerman JT, Schultz AK, Spragens KA, Boyce WM: Waterfowl ecology and avian influenza in California: do host traits inform us about viral occurrence? Avian Dis 2010, 54: 426-432. 10.1637/8912-043009-Reg.1
Bakken V, Runde O, Tjørve E: Norwegian bird ringing atlas. Divers - Auks. Stavanger Museum: Stavanger; 2003.
Toennessen R, Germundsson A, Jonassen CM, Haugen I, Berg K, Barrett RT, Rimstad E: Virological and serological surveillance for type A influenza in the black-legged kittiwake (Rissa tridactyla). Virol J 2011, 8: 21. 10.1186/1743-422X-8-21
Spackman E, Swayne DE, Gilbert M, Joly DO, Karesh WB, Suarez DL, Sodnomdarjaa R, Dulam P, Cardona C: Characterization of low pathogenicity avian influenza viruses isolated from wild birds in Mongolia 2005 through 2007. Virol J 2009, 6: 190. 10.1186/1743-422X-6-190
We thank Inger-Lise Larsen, Lone T. Engerdahl, Irene Haugen and Wanda Gajowniczek at the Norwegian Veterinary Institute for excellent technical assistance. We are also grateful to Knut Madslien for his role in organizing the sampling and to all hunters for collecting and submitting the samples. We also acknowledge the Norwegian Research Council (project number 178275), Animalia Meat and Poultry Research Centre and the Norwegian Food Safety Authority for financial support.
The authors declare that they have no competing interests.
RT, AGH, CMJ and ER were responsible for conception and design of the study. MJH, AGH and RT performed the initial PCR analysis and subtyping of the samples. RT and AGH were responsible for virus isolation and hemagglutination assays. RT, EFH and AGH cloned and sequenced the virus isolates. RT performed the similarity searches and the phylogenetic analysis. ABK and RT performed the statistical analysis. The article was drafted by RT, AGH, ER, CMJ, MJH, EFH and ABK. All authors read, revised and approved the final manuscript.
Electronic supplementary material
Additional file 2: Table S2: Results from the BLAST similarity searches. Results from the BLAST similarity searches performed on the full-length nucleotide sequences from the five virus isolates (H3N8, H6N8, H9N2, H13N2 and H16N3) sequenced in this study. Sequence with maximum identities is shown in brackets. The squares are colored green for Eurasian avian-like and blue for Eurasian gull-like segments. (PDF 29 KB)
About this article
Cite this article
Tønnessen, R., Kristoffersen, A.B., Jonassen, C.M. et al. Molecular and epidemiological characterization of avian influenza viruses from gulls and dabbling ducks in Norway. Virol J 10, 112 (2013). https://doi.org/10.1186/1743-422X-10-112