Skip to main content
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Research and News from ASM Journals
    • mSphere of Influence: Commentaries from Early Career Microbiologists
    • Archive
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mSphere
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
mSphere
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Research and News from ASM Journals
    • mSphere of Influence: Commentaries from Early Career Microbiologists
    • Archive
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mSphere
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
Research Article | Host-Microbe Biology

The Cloacal Microbiome of Five Wild Duck Species Varies by Species and Influenza A Virus Infection Status

Sarah M. Hird, Holly Ganz, Jonathan A. Eisen, Walter M. Boyce
Julia Oh, Editor
Sarah M. Hird
aDepartment of Molecular and Cell Biology, University of Connecticut, Storrs, Connecticut, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Sarah M. Hird
Holly Ganz
bAnimalBiome, Oakland, California, USA
cGenome Center, University of California, Davis, Davis, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jonathan A. Eisen
cGenome Center, University of California, Davis, Davis, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Walter M. Boyce
dPathology, Microbiology and Immunology, University of California, Davis, Davis, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Julia Oh
The Jackson Laboratory for Genomic Medicine
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/mSphere.00382-18
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

Waterfowl, especially ducks of the genus Anas, are natural reservoir species for influenza A virus (IAV). Duck populations contain nearly all the known diversity of IAVs, and the birds are asymptomatic to infection. Previous work established that IAV infection status is correlated with changes in the cloacal microbiome in juvenile mallards. Here, we analyze five Anas species to determine whether these duck species have similar IAV+ and IAV− cloacal microbiomes, or if the relationships among a host, influenza virus, and the microbiome are species specific. We assessed taxonomic composition of the microbiome, alpha diversity, and beta diversity and found very few patterns related to microbiome and infection status across species, while detecting strong differences within species. A host species-specific signal was stronger in IAV− ducks than IAV+ ducks, and the effect size of host species on the microbiome was three times higher in IAV− birds than IAV+ birds. The mallards and the northern shovelers, the species with highest sample sizes but also with differing feeding ecology, showed especially contrasting patterns in microbiome composition, alpha diversity, and beta diversity. Our results indicate that the microbiome may have a unique relationship with influenza virus infection at the species level.

IMPORTANCE Waterfowl are natural reservoir species for influenza A virus (IAV). Thus, they maintain high levels of pathogen diversity, are asymptomatic to the infection, and also contribute to the risk of a global influenza pandemic. An individual’s microbiome is a critical part in how a vertebrate manages pathogens and illness. Here, we describe the cloacal microbiome of 300 wild ducks, from five species (four with previously undescribed microbiomes), including both IAV-negative and IAV-positive individuals. We demonstrate that there is not one consistent “flu-like” microbiome or response to flu across species. Individual duck species appear to have unique relationships between their microbiomes and IAV, and IAV-negative birds have a stronger tie to host species than the IAV-positive birds. In a broad context, understanding the role of the microbiome in IAV reservoir species may have future implications for avian disease management.

INTRODUCTION

A vertebrate’s microbiome consists of millions to trillions of microorganisms that exist symbiotically with the host’s immune system, yet all vertebrates must also recognize and combat microbial pathogens. How the microbiome interacts with pathogens is of increasing interest due to recently established causal relationships. For example, a host’s microbiome provides a barrier to incoming pathogens and disruption of this barrier increases colonization by pathogens (1). Pathogens can influence the microbiome through systemic inflammation (such as in HIV infection, reviewed in reference 2) or the altering of microenvironments in local tissue. Changes in the microbiome to a state of dysbiosis can simultaneously be the result of viral infection and the cause of increased susceptibility to pathogens (3).

One globally important virus is influenza A virus (here IAV, family Orthomyxoviridae), which causes significant annual illness and mortality in humans and animals, especially poultry. The IAV genome consists of eight single-stranded RNA segments that code for up to 17 proteins (4). Two glycoproteins play a key role in the ability of the virus to infect a cell. Hemagglutinin (HA) allows the virus to get into a cell, and neuraminidase (NA) allows progeny virus to get out of the cell. The HA and NA subtypes of a virus designate the HA-NA subtype (e.g., H1N1). Birds, specifically waterfowl and shorebirds, are the natural reservoir species for IAV, as nearly all known HA (H1 to H16) and NA (N1 to N9) diversity can be found in their populations, including most subtypes (5–7). Additionally, infection with low-pathogenic avian influenza, or LPAI, causes little to no pathology in waterfowl (7, 8), although LPAI infection may be associated with lower body mass in mallards (9) and altered host behaviors, including feeding and timing of migration, in Bewick’s swans (10). Highly pathogenic avian influenza (HPAI) viruses can evolve in poultry infected with LPAI of H5 or H7 subtypes (11). HPAI infections have more aggressive pathology and higher mortality rates and infect respiratory tissues and other organs, instead of just the intestinal epithelium, as is the case in LPAI (12). HPAI viruses are highly infectious for poultry and can be passed to wild waterfowl and humans (13–16). HPAI can cause illness and mortality in some wild waterfowl species (12), which is concerning because birds migrate and represent a globally connected network of IAV hosts (6, 17). Species common to major flyways (e.g., American and Eurasian), such as Anas acuta (northern pintail), may provide an “intercontinental bridge” for IAVs (6), and their tendency to cross continents during migration is correlated with the amount of mixture between viruses of Eurasian origin and American origin (18, 19). IAVs can cross species and even class boundaries within the Vertebrata, infecting many birds and mammals, including humans (20, 21).

Correlative studies, like this one, can be a first step toward establishing causation between dysbiosis and disease. Given a correlation between viral infection and microbiome dysbiosis, one causal hypothesis is that infection by a pathogen can directly cause dysbiosis in the microbiome, the “IAV→ΔMB” hypothesis (Fig. 1A and B). This hypothesis has been demonstrated between IAV and the intestinal microbiota of chickens using time series data from experimentally infected hosts (22). Notably, chickens experience IAV infection more strongly than wild ducks (with some HPAI subtypes having mortality rates of up to 100% in chickens [11, 23, 24]), and the viral infection dynamics are not identical in the two systems. A second hypothesis is that hosts with an altered or dysbiotic microbiome are more prone to pathogen infection (e.g., see reference 25), the “ΔMB→IAV” hypothesis (Fig. 1C and D). Perhaps existing in mixed flocks can alter a wild duck’s microbiome, thus also altering susceptibility to infection. Third, it is possible that both mechanisms are occurring: dysbiotic states can be caused by infection while simultaneously increasing the susceptibility to pathogen invasion, the “ΔMB↔IAV” hypothesis (Fig. 1E and F). This hypothesis has been demonstrated with IAV in a mouse model (3), and overall, it may be the most plausible hypothesis as it acknowledges the various feedbacks and cross talk between a host’s immune system and the microbiome (26).

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

Conceptual model of the relationship between IAV infection and the microbiome. (A) The IAV→ΔMB hypothesis, depicted using a hypothetical ordination showing a central group of “normal” duck microbiomes (−) and IAV+ (+) samples surrounding the normal group. The “−” with the arrow shows that as the infection is cleared, the dysbiotic samples will move closer to the normal group as the microbiome recovers. Similarly, the “+” with the arrow shows that a newly infected duck will move further from “normal” as the infection progresses. (B) Over time, the viral load of IAV will increase, and shortly thereafter, the microbiome will also change. As the infection begins to clear, the microbiome will return to the “normal” state. (C) The ΔMB→IAV hypothesis. A dysbiotic microbiome (outlined in red circles) may make ducks more susceptible to infection, but infection itself will not alter the microbiome. (D) As IAV is contracted, proliferates, and is cleared by the host immune system, the microbiome is unaffected. (E) The ΔMB↔IAV hypothesis. Here, the landscape of pathogen susceptibility varies, and species-specific microbiomes have the lowest susceptibility. (F) Dysbiosis makes the host more susceptible to secondary infection, and additional infections further alter the microbiome (i.e., maintain dysbiosis).

The relationship between IAV and the microbiome is complex in wild ducks because many extrinsic factors are known to influence both the microbiome and virus dynamics. Many factors contribute to the dynamics of IAV infection in wild birds, including biotic and abiotic properties of a locality, time, space, and the host (27). Genetic similarity of LPAI viruses correlates with geographic space across duck species (28). The virus is likely spread through water, by virions that originated in the feces of an infected duck (29), or through virus shedding from the oropharynx (30). Viruses can persist in water for up to a month, and perhaps longer, depending on the starting viral load, temperature, pH, and salinity of the water (29, 31, 32). Infectious IAV can be isolated from feathers after an even longer time (33). IAV can also persist in sediments and lake ice and may be able to overwinter in aquatic environments (29). Aquatic plants, invertebrates, and filter feeders are also involved in the persistence of IAV in aquatic environments, in some cases sequestering virus in the environment (34, 35).

Additionally, there are influences of space, time, and the host’s immune system on viral system dynamics. Infection in ducks is temporally cyclic, with prevalence of infection, as well as particular subtypes, showing annual patterns that can vary by over 25% between years (36). The cyclical nature of subtypes may be due to the interaction between IAV subtypes and the avian immune system. Infection by a subtype can cause homosubtypic as well as heterosubtypic immunity (37), making reinfection by that subtype in migrating ducks returning to a locality less likely in the following year (although individual immunity varies [38]). The host can also influence virus transmission: age of infection contributes to the amount of virus shed during infection (39), and juveniles who have not been exposed previously typically have the highest prevalence. Host species is also important, as different ducks vary substantially in infection prevalence (6) and viral shedding (40).

Recent work on the microbiome has revealed its important role in health and disease for many vertebrates. How bacteria and viruses interact within the microbiome and how those interactions relate to pathogenicity are active areas of research (41, 42). Many of the same factors that influence flu dynamics are important for wild avian microbiome composition, including age (43), species (44), and sampling locality (45). Migratory status may also influence host health and has been associated with the microbiome as well (46). Understanding how these variables interact across complex, wild systems is important for evolutionary biology, infectious disease studies, and possibly global health.

Previously, we demonstrated that there is a strong correlation between the cloacal microbiome and influenza virus infection status in juvenile mallards from a single collecting locality (47). Ganz et al. (47) were able to detect robust associations between IAV and the microbiome because the study controlled for age, sampling locality, and host species. Here, we investigate the dynamics of IAV infection in five duck species (including the 122 mallards from reference 47) across geographic space and time, with the intention of elucidating the roles of these variables in the cloacal microbiome of IAV+ and IAV− ducks. We were particularly interested in whether the cloacal microbiomes of different species bear similar signatures of influenza virus infection status.

Despite the number of potentially confounded variables in our data set, we had several expectations for signal in the data. First, we expected to see decreases in alpha diversity for the IAV+ ducks within each species and across the data set, due to influenza virus infection being a state of dysbiosis. Second, we expected there to be compositional differences between the microbiomes of IAV− and IAV+ individuals within a species. Changes shared across species were expected to indicate flu-mediated changes. Finally, we expected to detect some taxa (OTUs) that were associated with influenza virus infection, within and across host species.

RESULTS

Sequencing.A total of 10,158,296 reads were produced from 300 birds, belonging to five species (Table 1 and also Table S1 in the supplemental material). After removing 54,019 sequences for being identified as chimeras, 66,072 sequences that were identified as being from chloroplasts, and 37 sequences identified as being from mitochondria, the final data set contained 10,038,168 reads. One hundred twenty-three ducks tested negative for IAV infection (i.e., “IAV−”), and 177 samples screened positive for IAV (i.e., “IAV+”). Of the IAV+ samples, 36 unique subtypes were represented, including 12 HAs and 9 NAs (Fig. 2). Four individuals were infected with two HAs, six individuals were infected with two NAs, and one additional individual was infected with two HAs and two NAs.

View this table:
  • View inline
  • View popup
TABLE 1

Duck species and influenza infection summary

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

Raw counts of influenza A virus subtypes from the different duck species. (Top) Hemagglutinin (HA) subtype frequencies. (Bottom) Neuraminidase (NA) subtype frequencies.

TABLE S1

Relative abundances for bacterial phyla in all samples. Download Table S1, XLSX file, 0.1 MB.
Copyright © 2018 Hird et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Taxonomic diversity.Six bacterial phyla were found in all 124 IAV− birds: Firmicutes (average relative abundance, 33.7%), Proteobacteria (32.7%), Bacteroidetes (13.8%), Fusobacteria (11.6%), Actinobacteria (3.1%), and Tenericutes (2.2%) (Fig. 3A). Sequences unable to be assigned to a bacterial phylum comprised an average of 1.6% of each sample. Another 24 phyla were found in at least two of the IAV− birds (Table S1).

FIG 3
  • Open in new tab
  • Download powerpoint
FIG 3

(A) Relative abundance of the six dominant bacterial phyla identified by the RDP classifier in 16S rRNA gene sequences per bird species, broken down by influenza virus status. Note the general lack of pattern between IAV+ and IAV− birds across species. The exception is that green-winged teals and mallards show same directionality for all phyla and northern shovelers show the opposite directionality from the green-winged teals and mallards. (B) Bar and whisker plot showing one measure of alpha diversity (observed OTUs) between IAV− and IAV+ individuals in each species. Note that there is no consistent pattern across species.

We were interested in whether all birds showed the same relationship to IAV infection. This was easily viewed at the phylum level (Fig. 3A); within each host species, each IAV infection status had dramatically different relative abundances of bacterial phyla. However, there was no consistent pattern to the direction of the differences between IAV− and IAV+ birds across species.

Alpha diversity.We used three measures of alpha diversity (Chao1, observed OTUs, and phylogenetic diversity [PD]; Table S2) to examine whether IAV+ individuals had a consistent difference in alpha diversity. Generally, all three measures were consistent within a species at P < 0.05 (exception: Anas carolinensis). Only Anas platyrhynchos and A. carolinensis showed a significant decrease in alpha diversity in the IAV+ individuals. The other species were not significantly different (and Anas clypeata even showed a nonsignificant trend of increased alpha diversity in IAV+ individuals, Fig. 3B).

TABLE S2

Alpha diversity data. Download Table S2, XLSX file, 0.03 MB.
Copyright © 2018 Hird et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Beta diversity.To examine how IAV infection correlated to the microbial communities (as a whole), we performed beta diversity analyses using unweighted UniFrac distance matrices (Fig. 4), which were similar in signal to beta diversity calculated with Bray-Curtis dissimilarities and weighted UniFrac distances (Fig. S1). Microbiomes did not strongly group by host species, and although PERMANOVA tests (via Adonis function in vegan package in R) for statistical differences in group centroids were significant for the combined data set (R2 = 0.029, P < 0.001), the dispersions of the groups were also significantly different (unweighted UniFrac, P = 0.003; weighted UniFrac, P = 0.002), indicating that the significance of the centroid differences may be inflated due to the different dispersions of the groups.

FIG 4
  • Open in new tab
  • Download powerpoint
FIG 4

Principal-coordinate analysis of unweighted UniFrac distances (a phylogenetic distance metric) between all samples. Largest panel shows all ducks (species denoted by shape); IAV+ (red) versus IAV− (black) samples. Smaller panels highlight each species and are still colored by influenza virus status. Adonis results indicate effect size (R2) and significance of influenza virus as an explanatory variable for each species; Adonis results for the categories “species” and “influenza status” are shown for all ducks.

FIG S1

NMDS of entire data set using Bray-Curtis and weighted UniFrac distance metrics. Download FIG S1, PDF file, 0.1 MB.
Copyright © 2018 Hird et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Influenza virus was always a significant factor (P < 0.02) with low effect size (R2 < 0.1). Here, dispersion tests show insignificant differences between the groups, indicating that differences between IAV− and IAV+ individuals are not driven by difference in dispersion of the groups. Age of the bird, sampling locality, and sample collection date were all associated with microbiome, as well as the hemagglutinin subtype (HA), neuraminidase subtype (NA), and HA-NA subtype (HANA, Table 2). Once the mallards were removed from the data set, however, only IAV infection status, species, location, and sample collection date remained significant. Removing mallards also caused species to lose significance in the data set that included only IAV+ individuals (Table 2).

View this table:
  • View inline
  • View popup
TABLE 2

Effect size and significance of PERMANOVA tests (via the Adonis function in vegan in R) of the metadata in our data seta

A RandomForest classifier was used to learn microbiome types from known samples and then assign unknown samples to different categories. We assigned unknown individuals within species to infection status with the following out-of-box (OOB) error rates: Anas acuta (northern pintail) = 26.6%, Anas americana (American wigeon) = 28%, A. carolinensis (green-winged teal) = 21.0%, A. clypeata (northern shoveler) = 10.7%, A. platyrhynchos (mallard) = 11.0%.

Common microbiome.A common microbiome of nine OTUs was defined as the OTUs that were present in 90% of individuals in each species (Table 3). The OTUs present in 90% of the IAV+ birds and the OTUs present in 90% of the IAV− birds were both subsets of the OTUs found in 90% of all species.

View this table:
  • View inline
  • View popup
TABLE 3

The OTUs that were common to 90% of individuals within each species; taxonomic assignments were made using RDP classifier (see text for more details)

To further examine the phylogenetic placement of these OTUs, we searched representative sequences from the OTUs using blastn against the NCBI nr and 16S rRNA gene databases (which include additional sequences not included in the database used for the taxonomic annotation shown in Table 3) and also constructed phylogenetic trees based on the blast matches. Four of these OTUs had 100% identical matches to sequences in the NCBI nr database. New.ReferenceOTU451 was found to be 100% identical to a single clone isolated from coastal surface waters and fell within the Campylobacter canadensis. OTU 4439398 was 100% identical to an uncultured bacterium from a human ileum (found in the nr database) and in phylogenetic analysis of sequence hits from the 16S rRNA gene database was found to be sister to two sequences belonging to Fusobacterium mortiferum (48). New.ReferenceOTU712 was 100% identical to a sequence found in a study of Canada Goose feces (49). OTU1122504 was 100% identical to Flavobacterium ceti, as well as many uncultured Flavobacterium clones. The remaining OTUs did not have any 100% matches to any sequences in either database, but phylogenetic analysis of the matches provides some additional insight into their placement. New.ReferenceOTU552 was sister to a clade containing Pasteurellaceae bacteria and uncultured bacteria (see Fig. S2A in the supplemental material). OTU 284123 fell within a clade of Campylobacter and was sister to Campylobacter cuniculorum (Fig. S2B). New.ReferenceOTU81 fell within a clade containing Streptobacillus hongkongensis, Leptotrichia goodfellowii, and Sebaldella termitidis and was sister to a clade containing Sneathia sanguinegens, several Streptobacillus species, and Oceanivirga salmonicida (Fig. S2C). New.ReferenceOTU568 was sister to a clade containing several Veillonella species as well as Megasphaera micronuciformis (Fig. S2D). This clade belonged to a larger clade containing more Veillonella and Dialister species. New.ReferenceOTU701 was sister to a large clade within Firmicutes containing many species (Fig. S2E).

FIG S2

NCBI Blast to Tree information for the poorly identified common OTUs. Download FIG S2, PDF file, 0.8 MB.
Copyright © 2018 Hird et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Differentially abundant OTUs.Within species, we detected many OTUs that were differentially abundant between IAV+ and IAV− birds (Fig. 5A). Compared across the five species, IAV+ northern shovelers and IAV− mallards contained many of the same significant taxa/OTUs. Many OTUs were detected that were significantly associated with infection status in more than one species. One pattern worth noting is that northern shovelers and mallards appear to share many OTUs that are significantly associated with flu status but in opposite directions: e.g., box 1 on Fig. 5 depicts a Rothia mucilaginosa OTU that is significantly more abundant in IAV+ northern shovelers and IAV− mallards.

FIG 5
  • Open in new tab
  • Download powerpoint
FIG 5

Northern shovelers may have unique microbiome differences with respect to IAV infection, relative to the other duck species. (A) Differentially abundant OTUs within each species. Bars indicate individual OTUs belonging to the taxa indicated by the bar at top of the figure, where each column going up indicates an OTU significantly higher in IAV+ birds within a species and columns going down indicate OTUs significantly higher in IAV− birds. Color indicates bird species. Gray bars highlight where an identical OTU was found to be significant in more than one species, with the highest taxonomic resolution of these OTUs numbered and labeled across the bottom. (B) NMDS ordination of Bray-Curtis distances on northern shoveler and mallard samples only. Left panel highlights the IAV− ducks and includes the multivariate t-distributed ellipses; right panel highlights the IAV+ ducks. (C) Venn diagrams of the OTUs found in 90% of the IAV+ and IAV− individuals for mallards and two other species. (Upper) Mallard OTUs compared to northern pintail OTUs, showing that there are some OTUs that overlap in the IAV+/IAV+ region and in the IAV−/IAV− region with no OTUs in the IAV+/IAV− regions. These data meet the expectation that there is a similar infection-based pattern to the OTUs in these two species. (Lower) Mallard OTUs compared to northern shoveler OTUs, showing that none of the previous expectations are met and the IAV+/IAV− regions contain several OTUs.

Mallard and northern shoveler comparisons.We focused next on just the mallards and northern shovelers, as they had the highest sample sizes and to minimize noise in the data. A RandomForest classifier was used to learn microbiome types from known samples and then assign unknown samples to different categories. The classifier had an out-of-box (OOB) error rate of 2.61% for assigning IAV− birds to species and an OOB error rate of 16.5% for IAV+ birds. An ordination of just the mallards and northern shovelers shows that IAV− individuals comprise distinct clusters, whereas the IAV+ birds are largely overlapping (Fig. 5B). To further look at the differences between northern shovelers and the other species, we created Venn diagrams of the OTUs present in 90% of the IAV+ and IAV− individuals within each species (Fig. 5C).

DISCUSSION

The interactions among a host organism and its microbiome and viral pathogens are complex and likely of evolutionary relevance. Here, we have described the IAV− and IAV+ cloacal microbiomes of five Anas ducks and found major differences between individuals of different infection statuses within a species, yet little commonality to those differences across species. The communities were largely comprised of six dominant bacterial phyla: Proteobacteria and Firmicutes together constituted greater than 50% of the bacteria in all bird species and infection statuses. These two phyla have previously been associated with avian microbiomes, including avian intestinal samples (e.g., references 45 and 50) and cloacal swabs specifically (47, 51). Fusobacteria, Actinobacteria, Tenericutes, and Bacteroidetes were the other dominant phyla.

The IAV− and IAV+ birds within each species showed dramatically different microbial community composition (Fig. 3A). Within each species, infection status explained between 4 and 7% of the microbiome community variation (P values < 0.017). Analysis of all birds together showed that infection status explained 1.29% of the variation (P = 0.001); however, when mallards were removed due to suspicion that their large sample size might be affecting the significance of the tests, infection status continued to have a small (R2 = 1.66%) but significant (P = 0.001) association. There were no within-species differences between IAV− and IAV+ individuals that are shared across all five species. For example, alpha diversity is significantly lower in the IAV+ individuals for mallards and green-winged teals, whereas no significant difference was found for the other three species. In fact, there was a nonsignificant increase in alpha diversity in the IAV+ northern shovelers.

Northern shovelers, who are the only filter feeders in the data set, appeared to show an opposite signal relative to the other ducks in a variety of metrics. For example, Fusobacteria were higher in the IAV+ individuals within each species except the northern shovelers and many of the OTUs that were significantly higher in one infection status in northern shovelers showed the opposite infection status in the other species. When we compared just the mallards and the northern shovelers, the IAV− individuals from both species were distinct, whereas the IAV+ individuals were largely overlapping. This was corroborated by the random forest analysis that revealed a higher rate of correct assignment to species in the IAV− birds (2.6% error) than the IAV+ birds (16.5% error). In the data sets including all five species, host species was significant (P = 0.001) in both the IAV− (R2 = 10.2%) and IAV+ data sets (R2 = 3.7%); however, removal of the mallards caused the effect size in IAV− birds to decrease to 8.5% (P = 0.001) and in the IAV+ data set to increase (R2 = 8.9%) and become nonsignificant (P = 0.048). We take this to mean that mallards may have been driving the species signal in the IAV+ data set and that there is a much stronger species-level association in the microbiome of IAV− individuals.

In addition to the infection status and host species being significantly associated with the metadata, sampling locality (R2 = 5.1%) and time (season and year, R2 = 7.5%) were all significantly associated with the microbiomes (P = 0.001) and remained significant even after removal of the mallards (although effect sizes decreased by ∼1%, Table 2). This is unsurprising, given the strong a priori correlations between IAV subtype, season and year sampled, and host species collected at the various localities, as well as the known interactions between these variables, e.g., references 6, 11, 27, and 38), including specifically in mallards (52). Sampling locality can be more important thant ecology or phylogeny for shaping the avian microbiome (45), and we may be seeing this in the ducks as well. When we analyzed three species with both infection statuses from a single locality (Fig. 6), species had both the highest effect size (6.89%) and the lowest P value (0.001). Season/year that the samples were collected had the next highest effect size (5.75%, P = 0.12); IAV infection status had the lowest effect size (2.93%) and a low P value (P = 0.006).

FIG 6
  • Open in new tab
  • Download powerpoint
FIG 6

A single NMDS plot of Bray-Curtis dissimilarity of 16S rRNA gene sequences from cloacal swab of ducks from a single locality. Panels highlight ducks of different species (A), IAV infection status (B), and time collected (C).

The significance of IAV subtype on the microbiome remains unknown; our data set contained 36 IAV subtypes, with generally low sample sizes for any single HA or NA subtype. Additionally, the subtypes were correlated to sampling site, time of collection, and species. However, HA subtype, NA subtype, and HA-NA subtype were all significantly associated with the microbiome in the full data set, with subtype (HANA) showing an effect size of 25% in the IAV+ birds (P = 0.001). The significance of all three rose to greater than 0.13, though, when the mallards were removed. Subtypes were not significantly correlated with the microbiome within any of the species except the mallard (HA, NA, and HANA, P = 0.001) and the association between the northern shoveler and the HA subtype (P = 0.003). This indicates that HA, NA, and HA-NA subtype may be correlated with the microbiome; however, how significant that association is cannot be established.

Sample size may be essential for determining which metadata are associated with the microbiome across infection statuses—species with larger sample sizes (mallards and northern shovelers) had more categorical data significantly associated with them than the lower-sample-size species (Table 2). This is unsurprising, as the complex interactions between these likely confounded variables make for a challenging statistical problem. We hypothesize that increasing the sample size for green-winged teals, northern pintails, and American wigeons would likely uncover more statistically robust associations between birds and their microbiomes.

Do these data bear on the causative hypotheses? All the data are congruent with the “IAV→ΔMB” hypothesis, particularly the fact that the strength of the host species association is higher in IAV− birds than in IAV+ birds. This implies that there is a species-level microbiome for these duck species that deteriorates during IAV infection. It also implies that over time, ducks may return to their “species state,” which is predicted by the “IAV→ΔMB” hypothesis. Our results, however, are also compatible with the “ΔMB→IAV” hypothesis: including all dysbiotic states as one would cause the dispersion of the IAV+ groups to be large and the IAV− groups to be smaller. It may also be true that whom a duck associates with (or is simply near) may influence its likelihood of IAV infection, including the subtype. Mixed species groups may have increased exposure to flu but also have a more similar microbiome. Time series data in these wild reservoir species would be ideal to determine a causative role of IAV on dysbiosis or dysbiosis on IAV susceptibility.

MATERIALS AND METHODS

Sampling.Samples were collected from both living ducks (all mallard samples) and those that were recently deceased (<2 h) (for all other species). Samples were collected via swabs of cloaca. Total number of cloacal swabs for each species were as follows: 31 northern pintails (Anas acuta), 25 American wigeons (Anas americana), 19 green-winged teals (Anas carolinensis), 57 northern shovelers (Anas clypeata), and 168 mallards (Anas platyrhynchos), including the 122 juvenile individuals from the work of Ganz et al. (47). Detailed information about the samples is found in Table 1 and Table S3. Swabs were placed in separate vials containing 2 ml of ice-cold virus transport medium (VTM; medium 199 with Earle’s salts, l-glutamine, and sodium bicarbonate, plus 2 mU/liter penicillin G, 200 mg/liter streptomycin, 2 mU/liter polymyxin B, 250 mg/liter gentamicin, 0.5 mU/liter nystatin, 60 mg/liter ofloxacin, 200 mg/liter sulfamethoxazole, and 0.5% bovine serum albumin V). The samples were transported on ice to the laboratory where they were stored at −80°C. Samples underwent one or more freeze-thaw cycles as aliquots were removed for influenza virus testing prior to microbiome analysis.

TABLE S3

Sample information. Download Table S3, XLSX file, 0.1 MB.
Copyright © 2018 Hird et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Ethics statement.All animal research was conducted with the approval of the University of California Davis Animal Care and Use Committee, approval number 19760. The protocols adhere to federal law, specifically the Animal Welfare Act and the Health Research Extension Act of 1985.

Influenza virus screening and typing.Samples were screened for the influenza virus matrix gene by RT-PCR, viruses were isolated from matrix-positive samples by egg inoculation, and full genome sequences were generated as described in reference 53. Samples were classified as “negative” if they (i) did not yield virus on egg inoculation and (ii) were matrix RT-PCR negative.

DNA extraction and 16S rRNA gene sequencing.DNA extraction and sequencing of the V4 region of 16S rRNA genes were performed in the same manner as in reference 47. DNA was extracted from cloacal swabs using the Mo Bio PowerSoil 96-well soil DNA isolation kit (Carlsbad, CA). Samples in 96-well plates were incubated at 65°C for ten minutes after addition of C1 solution. Each 96-well plate was then vortexed for 3 min using a plate shaker on high (2,600 rpm), and then the standard kit protocol was followed. After elution, DNA was quantified using Qubit fluorometric quantitation (Invitrogen, South San Francisco, CA). We used nested PCR to amplify and sequence the 16S rRNA gene. DNA was characterized for bacterial diversity based on the V4 region of the16S rRNA gene following the methods of reference 54. When possible, 1.0 to 5.0 μl of template DNA was used for PCR. However, due to low DNA concentrations for some samples, bacterial DNA was amplified by a two-step PCR enrichment using the primers 27F−YM+3 and 1492R to target V1 to V4 of the 16S rRNA gene (55). The 7-fold-degenerate primer 27f−YM+3 is composed of four parts 27f−YM (AGAGTTTGATYMTGGCTCAG), plus one part each of primers specific for the amplification of Bifidobacteriaceae, Borrelia, and Chlamydiales sequences (55). For the second PCR, primers used were the bacteria/archaeal primers 515F/806R (54) modified by addition of Illumina adaptor and an in-house barcode system (described in reference 56). After amplification. magnetic beads (Agencourt AMPure XP, Beckman Coulter, Indianapolis, IN) were used to clean PCR mixtures. Amplicons were quantified and characterized using Qubit fluorometric quantitation, qPCR, and the Bioanalyzer (Agilent Technologies, Santa Clara, CA) prior to sequencing. Libraries were sequenced using an Illumina MiSeq system, generating 250-bp paired-end amplicon reads. The amplicon data were multiplexed using dual barcode combinations for each sample.

Data analysis.We used a custom script (available in a GitHub repository, https://github.com/gjospin/scripts/blob/master/Demul_trim_prep.pl) to assign each pair of reads to their respective samples when parsing the raw data. This script allows for 1-bp difference per barcode. The paired reads were then aligned, and a consensus was computed using FLASH (57) with maximum overlap of 120 and a minimum overlap of 70 (other parameters were left as default). The custom script automatically demultiplexes the data into fastq files, executes FLASH, and parses its results to reformat the sequences with appropriate naming conventions for Quantitative Insights in Microbial Ecology (QIIME v.1.9.1 [58]) in fasta format; QIIME and the R package phyloseq (59) were subsequently used for most analyses.

First, representative operational taxonomic units (OTUs) were chosen by de novo clustering of sequences at the >97% identity level, the taxonomic composition of each cluster was assigned with RDP classifier, and an OTU table was constructed for each sample using QIIME’s pick_otus_through_otu_table.py script. Alpha diversity (observed OTUs, phylogenetic diversity, Chao1) was estimated using the core_diversity_analyses.py script. Pairwise distance matrices were constructed for unweighted UniFrac distance (60), weighted UniFrac distance, and Bray-Curtis dissimilarity for beta diversity analyses.

Differentially abundant OTUs were inferred using the program DESeq2 (61) through the program phyloseq. Each species was analyzed separately for OTUs that were significantly higher (adjusted P value < 0.001) in either IAV+ or IAV− individuals. Only OTUs with a log fold increase greater than 2 were retained. The R package “randomForest” (62) was used to perform machine learning classifications within species and infection statuses.

Accession number(s).Sequence data have been uploaded to the NCBI SRA under BioProject ID PRJNA464410.

ACKNOWLEDGMENTS

We acknowledge Magdalena Plancarte, Winston Vickers, Matt Rolston and Guillaume Jospin for their contributions to this project. We also thank the United States Geological Survey, the California Waterfowl Association, and the California Department of Fish and Wildlife for their cooperation and assistance in sample collection.

This study was supported in part by funding from NIH (CEIRS—HHSN272201400008C and HHSN266200700010C) (W.M.B.) and by a grant to J.A.E. from the Alfred P. Sloan Foundation as part of their “Microbiology of the Built Environment” program. S.M.H. was funded through the UC Davis Chancellor’s Postdoctoral Fellowship and startup funds from the University of Connecticut.

We declare a conflict of interest. We declare that H.G. is CEO of a commercial company, AnimalBiome, and that Jonathan Eisen is a scientific advisor to the companies AnimalBiome, Zymo Research, Phylagen, and Indigo Agriculture. This has not altered the motivation for the research, the full and objective presentation of the research, or our adherence to all ASM policies. The research was not funded by any of these companies.

FOOTNOTES

    • Received July 11, 2018.
    • Accepted September 26, 2018.
  • Copyright © 2018 Hird et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

REFERENCES

  1. 1.↵
    1. Sekirov I,
    2. Russell SL,
    3. Antunes LCM,
    4. Finlay BB
    . 2010. Gut microbiota in health and disease. Physiol Rev 90:859–904. doi:10.1152/physrev.00045.2009.
    OpenUrlCrossRefPubMed
  2. 2.↵
    1. Zilberman-Schapira G,
    2. Zmora N,
    3. Itav S,
    4. Bashiardes S,
    5. Elinav H,
    6. Elinav E
    . 2016. The gut microbiome in human immunodeficiency virus infection. BMC Med 14:83. doi:10.1186/s12916-016-0625-3.
    OpenUrlCrossRef
  3. 3.↵
    1. Yildiz S,
    2. Mazel-Sanchez B,
    3. Kandasamy M,
    4. Manicassamy B,
    5. Schmolke M
    . 2018. Influenza A virus infection impacts systemic microbiota dynamics and causes quantitative enteric dysbiosis. Microbiome 6:9. doi:10.1186/s40168-017-0386-z.
    OpenUrlCrossRef
  4. 4.↵
    1. Dubois J,
    2. Terrier O,
    3. Rosa-Calatrava M
    . 2014. Influenza viruses and mRNA splicing: doing more with less. mBio 5:e00070-14. doi:10.1128/mBio.00070-14.
    OpenUrlAbstract/FREE Full Text
  5. 5.↵
    1. Fouchier RAM,
    2. Munster V,
    3. Wallensten A,
    4. Bestebroer TM,
    5. Herfst S,
    6. Smith D,
    7. Rimmelzwaan GF,
    8. Olsen B,
    9. Osterhaus ADME
    . 2005. Characterization of a novel influenza A virus hemagglutinin subtype (H16) obtained from black-headed gulls. J Virol 79:2814–2822. doi:10.1128/JVI.79.5.2814-2822.2005.
    OpenUrlAbstract/FREE Full Text
  6. 6.↵
    1. Olsen B,
    2. Munster VJ,
    3. Wallensten A,
    4. Waldenström J,
    5. Osterhaus ADME,
    6. Fouchier RAM
    . 2006. Global patterns of influenza A virus in wild birds. Science 312:384–388. doi:10.1126/science.1122438.
    OpenUrlAbstract/FREE Full Text
  7. 7.↵
    1. Webster RG,
    2. Bean WJ,
    3. Gorman OT,
    4. Chambers TM,
    5. Kawaoka Y
    . 1992. Evolution and ecology of influenza A viruses. Microbiol Rev 56:152–179.
    OpenUrlAbstract/FREE Full Text
  8. 8.↵
    1. Kida H,
    2. Yanagawa R,
    3. Matsuoka Y
    . 1980. Duck influenza lacking evidence of disease signs and immune response. Infect Immun 30:547–553.
    OpenUrlAbstract/FREE Full Text
  9. 9.↵
    1. Latorre-Margalef N,
    2. Gunnarsson G,
    3. Munster VJ,
    4. Fouchier RAM,
    5. Osterhaus ADME,
    6. Elmberg J,
    7. Olsen B,
    8. Wallensten A,
    9. Haemig PD,
    10. Fransson T,
    11. Brudin L,
    12. Waldenström J
    . 2009. Effects of influenza A virus infection on migrating mallard ducks. Proc Biol Sci 276:1029–1036. doi:10.1098/rspb.2008.1501.
    OpenUrlAbstract/FREE Full Text
  10. 10.↵
    1. van Gils JA,
    2. Munster VJ,
    3. Radersma R,
    4. Liefhebber D,
    5. Fouchier RAM,
    6. Klaassen M
    . 2007. Hampered foraging and migratory performance in swans infected with low-pathogenic avian influenza A virus. PLoS One 2:e184. doi:10.1371/journal.pone.0000184.
    OpenUrlCrossRefPubMed
  11. 11.↵
    1. Alexander DJ
    . 2000. A review of avian influenza in different bird species. Vet Microbiol 74:3–13. doi:10.1016/S0378-1135(00)00160-7.
    OpenUrlCrossRefPubMedWeb of Science
  12. 12.↵
    1. Short KR,
    2. Veldhuis Kroeze EJB,
    3. Reperant LA,
    4. Richard M,
    5. Kuiken T
    . 2014. Influenza virus and endothelial cells: a species specific relationship. Front Microbiol 5:653. doi:10.3389/fmicb.2014.00653.
    OpenUrlCrossRefPubMed
  13. 13.↵
    1. Liu J,
    2. Xiao H,
    3. Lei F,
    4. Zhu Q,
    5. Qin K,
    6. Zhang X-W,
    7. Zhang X-L,
    8. Zhao D,
    9. Wang G,
    10. Feng Y,
    11. Ma J,
    12. Liu W,
    13. Wang J,
    14. Gao GF
    . 2005. Highly pathogenic H5N1 influenza virus infection in migratory birds. Science 309:1206. doi:10.1126/science.1115273.
    OpenUrlAbstract/FREE Full Text
  14. 14.↵
    1. Subbarao K,
    2. Klimov A,
    3. Katz J,
    4. Regnery H,
    5. Lim W,
    6. Hall H,
    7. Perdue M,
    8. Swayne D,
    9. Bender C,
    10. Huang J,
    11. Hemphill M,
    12. Rowe T,
    13. Shaw M,
    14. Xu X,
    15. Fukuda K,
    16. Cox N
    . 1998. Characterization of an avian influenza A (H5N1) virus isolated from a child with a fatal respiratory illness. Science 279:393–396. doi:10.1126/science.279.5349.393.
    OpenUrlAbstract/FREE Full Text
  15. 15.↵
    1. Li KS,
    2. Guan Y,
    3. Wang J,
    4. Smith GJD,
    5. Xu KM,
    6. Duan L,
    7. Rahardjo AP,
    8. Puthavathana P,
    9. Buranathai C,
    10. Nguyen TD,
    11. Estoepangestie ATS,
    12. Chaisingh A,
    13. Auewarakul P,
    14. Long HT,
    15. Hanh NTH,
    16. Webby RJ,
    17. Poon LLM,
    18. Chen H,
    19. Shortridge KF,
    20. Yuen KY,
    21. Webster RG,
    22. Peiris JSM
    . 2004. Genesis of a highly pathogenic and potentially pandemic H5N1 influenza virus in eastern Asia. Nature 430:209. doi:10.1038/nature02746.
    OpenUrlCrossRefPubMedWeb of Science
  16. 16.↵
    1. Chen H,
    2. Smith GJD,
    3. Zhang SY,
    4. Qin K,
    5. Wang J,
    6. Li KS,
    7. Webster RG,
    8. Peiris JSM,
    9. Guan Y
    . 2005. Avian flu: H5N1 virus outbreak in migratory waterfowl. Nature 436:191–192. doi:10.1038/nature03974.
    OpenUrlCrossRefPubMedWeb of Science
  17. 17.↵
    1. Wallensten A,
    2. Munster VJ,
    3. Elmberg J,
    4. Osterhaus ADME,
    5. Fouchier RAM,
    6. Olsen B
    . 2005. Multiple gene segment reassortment between Eurasian and American lineages of influenza A virus (H6N2) in guillemot (Uria aalge). Arch Virol 150:1685–1692. doi:10.1007/s00705-005-0543-8.
    OpenUrlCrossRefPubMedWeb of Science
  18. 18.↵
    1. Pearce JM,
    2. Reeves AB,
    3. Ramey AM,
    4. Hupp JW,
    5. Ip HS,
    6. Bertram M,
    7. Petrula MJ,
    8. Scotton BD,
    9. Trust KA,
    10. Meixell BW,
    11. Runstadler JA
    . 2011. Interspecific exchange of avian influenza virus genes in Alaska: the influence of trans-hemispheric migratory tendency and breeding ground sympatry. Mol Ecol 20:1015–1025. doi:10.1111/j.1365-294X.2010.04908.x.
    OpenUrlCrossRefPubMed
  19. 19.↵
    1. Ramey AM,
    2. Pearce JM,
    3. Flint PL,
    4. Ip HS,
    5. Derksen DV,
    6. Franson JC,
    7. Petrula MJ,
    8. Scotton BD,
    9. Sowl KM,
    10. Wege ML,
    11. Trust KA
    . 2010. Intercontinental reassortment and genomic variation of low pathogenic avian influenza viruses isolated from northern pintails (Anas acuta) in Alaska: examining the evidence through space and time. Virology 401:179–189. doi:10.1016/j.virol.2010.02.006.
    OpenUrlCrossRefPubMed
  20. 20.↵
    1. Root JJ,
    2. Shriner SA,
    3. Ellis JW,
    4. VanDalen KK,
    5. Sullivan HJ,
    6. Franklin AB
    . 2015. When fur and feather occur together: interclass transmission of avian influenza A virus from mammals to birds through common resources. Sci Rep 5:14354. doi:10.1038/srep14354.
    OpenUrlCrossRef
  21. 21.↵
    1. Short KR,
    2. Richard M,
    3. Verhagen JH,
    4. van Riel D,
    5. Schrauwen EJA,
    6. van den Brand JMA,
    7. Mänz B,
    8. Bodewes R,
    9. Herfst S
    . 2015. One health, multiple challenges: the inter-species transmission of influenza A virus. One Health 1:1–13. doi:10.1016/j.onehlt.2015.03.001.
    OpenUrlCrossRefPubMed
  22. 22.↵
    1. Yitbarek A,
    2. Weese JS,
    3. Alkie TN,
    4. Parkinson J,
    5. Sharif S
    . 2018. Influenza A virus subtype H9N2 infection disrupts the composition of intestinal microbiota of chickens. FEMS Microbiol Ecol 94. doi:10.1093/femsec/fix165.
    OpenUrlCrossRef
  23. 23.↵
    1. Naeem K
    . 2003. The avian influenza H7 N3 outbreak in south central Asia. Avian Dis 47:31–35.
    OpenUrl
  24. 24.↵
    1. Jiao P,
    2. Song H,
    3. Liu X,
    4. Song Y,
    5. Cui J,
    6. Wu S,
    7. Ye J,
    8. Qu N,
    9. Zhang T,
    10. Liao M
    . 2016. Pathogenicity, transmission and antigenic variation of H5N1 highly pathogenic avian influenza viruses. Front Microbiol 7:635. doi:10.3389/fmicb.2016.00635.
    OpenUrlCrossRef
  25. 25.↵
    1. Teo SM,
    2. Mok D,
    3. Pham K,
    4. Kusel M,
    5. Serralha M,
    6. Troy N,
    7. Holt BJ,
    8. Hales BJ,
    9. Walker ML,
    10. Hollams E,
    11. Bochkov YA,
    12. Grindle K,
    13. Johnston SL,
    14. Gern JE,
    15. Sly PD,
    16. Holt PG,
    17. Holt KE,
    18. Inouye M
    . 2015. The infant nasopharyngeal microbiome impacts severity of lower respiratory infection and risk of asthma development. Cell Host Microbe 17:704–715. doi:10.1016/j.chom.2015.03.008.
    OpenUrlCrossRefPubMed
  26. 26.↵
    1. Hooper LV,
    2. Littman DR,
    3. Macpherson AJ
    . 2012. Interactions between the microbiota and the immune system. Science 336:1268–1273. doi:10.1126/science.1223490.
    OpenUrlAbstract/FREE Full Text
  27. 27.↵
    1. Munster VJ,
    2. Baas C,
    3. Lexmond P,
    4. Waldenström J,
    5. Wallensten A,
    6. Fransson T,
    7. Rimmelzwaan GF,
    8. Beyer WEP,
    9. Schutten M,
    10. Olsen B,
    11. Albert DM,
    12. Fouchier RAM
    . 2007. Spatial, temporal, and species variation in prevalence of influenza A viruses in wild migratory birds. PLoS Pathog 3:e61. doi:10.1371/journal.ppat.0030061.
    OpenUrlCrossRefPubMed
  28. 28.↵
    1. Reeves AB,
    2. Pearce JM,
    3. Ramey AM,
    4. Meixell BW,
    5. Runstadler JA
    . 2011. Interspecies transmission and limited persistence of low pathogenic avian influenza genomes among Alaska dabbling ducks. Infect Genet Evol 11:2004–2010. doi:10.1016/j.meegid.2011.09.011.
    OpenUrlCrossRefPubMed
  29. 29.↵
    1. Franklin AB,
    2. VanDalen KK,
    3. Huyvaert KP
    . 2011. Avian influenza virus in aquatic environments: an ecological perspective, p 59–72. In Majumdar SK, Brenner FJ, Huffman JE, McLean RG, Panah AI, Pietrobon PJ, Keeler SP, Shive SE (ed), Pandemic influenza viruses: science, surveillance and public health. The Pennsylvania Academy of Science, Philadelphia, PA.
  30. 30.↵
    1. Hulse-Post DJ,
    2. Sturm-Ramirez KM,
    3. Humberd J,
    4. Seiler P,
    5. Govorkova EA,
    6. Krauss S,
    7. Scholtissek C,
    8. Puthavathana P,
    9. Buranathai C,
    10. Nguyen TD,
    11. Long HT,
    12. Naipospos TSP,
    13. Chen H,
    14. Ellis TM,
    15. Guan Y,
    16. Peiris JSM,
    17. Webster RG
    . 2005. Role of domestic ducks in the propagation and biological evolution of highly pathogenic H5N1 influenza viruses in Asia. Proc Natl Acad Sci U S A 102:10682–10687. doi:10.1073/pnas.0504662102.
    OpenUrlAbstract/FREE Full Text
  31. 31.↵
    1. Stallknecht DE,
    2. Kearney MT,
    3. Shane SM,
    4. Zwank PJ
    . 1990. Effects of pH, temperature, and salinity on persistence of avian influenza viruses in water. Avian Dis 34:412–418. doi:10.2307/1591429.
    OpenUrlCrossRefPubMedWeb of Science
  32. 32.↵
    1. Brown JD,
    2. Goekjian G,
    3. Poulson R,
    4. Valeika S,
    5. Stallknecht DE
    . 2009. Avian influenza virus in water: infectivity is dependent on pH, salinity and temperature. Vet Microbiol 136:20–26. doi:10.1016/j.vetmic.2008.10.027.
    OpenUrlCrossRefPubMedWeb of Science
  33. 33.↵
    1. Yamamoto Y,
    2. Nakamura K,
    3. Yamada M,
    4. Mase M
    . 2010. Persistence of avian influenza virus (H5N1) in feathers detached from bodies of infected domestic ducks. Appl Environ Microbiol 76:5496–5499. doi:10.1128/AEM.00563-10.
    OpenUrlAbstract/FREE Full Text
  34. 34.↵
    1. Dalziel AE,
    2. Delean S,
    3. Heinrich S,
    4. Cassey P
    . 2016. Persistence of low pathogenic influenza A virus in water: a systematic review and quantitative meta-analysis. PLoS One 11:e0161929. doi:10.1371/journal.pone.0161929.
    OpenUrlCrossRef
  35. 35.↵
    1. Vong S,
    2. Ly S,
    3. Mardy S,
    4. Holl D,
    5. Buchy P
    . 2008. Environmental contamination during influenza A virus (H5N1) outbreaks, Cambodia, 2006. Emerg Infect Dis 14:1303–1305. doi:10.3201/eid1408.070912.
    OpenUrlCrossRefPubMedWeb of Science
  36. 36.↵
    1. Krauss S,
    2. Walker D,
    3. Pryor SP,
    4. Niles L,
    5. Li CH,
    6. Hinshaw VS,
    7. Webster RG
    . 2004. Influenza A viruses of migrating wild aquatic birds in North America. Vector-Borne Zoonotic Dis 4:177–189. doi:10.1089/vbz.2004.4.177.
    OpenUrlCrossRefPubMedWeb of Science
  37. 37.↵
    1. Latorre-Margalef N,
    2. Grosbois V,
    3. Wahlgren J,
    4. Munster VJ,
    5. Tolf C,
    6. Fouchier RAM,
    7. Osterhaus ADME,
    8. Olsen B,
    9. Waldenström J
    . 2013. Heterosubtypic immunity to influenza A virus infections in mallards may explain existence of multiple virus subtypes. PLoS Pathog 9:e1003443. doi:10.1371/journal.ppat.1003443.
    OpenUrlCrossRefPubMed
  38. 38.↵
    1. Tolf C,
    2. Latorre-Margalef N,
    3. Wille M,
    4. Bengtsson D,
    5. Gunnarsson G,
    6. Grosbois V,
    7. Hasselquist D,
    8. Olsen B,
    9. Elmberg J,
    10. Waldenström J
    . 2013. Individual variation in influenza A virus infection histories and long-term immune responses in mallards. PLoS One 8:e61201. doi:10.1371/journal.pone.0061201.
    OpenUrlCrossRefPubMed
  39. 39.↵
    1. Costa TP,
    2. Brown JD,
    3. Howerth EW,
    4. Stallknecht DE
    . 2010. The effect of age on avian influenza viral shedding in mallards (Anas platyrhynchos). Avian Dis 54:581–585. doi:10.1637/8692-031309-ResNote.1.
    OpenUrlCrossRefPubMedWeb of Science
  40. 40.↵
    1. Costa TP,
    2. Brown JD,
    3. Howerth EW,
    4. Stallknecht DE
    . 2011. Variation in viral shedding patterns between different wild bird species infected experimentally with low-pathogenicity avian influenza viruses that originated from wild birds. Avian Pathol 40:119–124. doi:10.1080/03079457.2010.540002.
    OpenUrlCrossRefPubMedWeb of Science
  41. 41.↵
    1. Weiss B,
    2. Aksoy S
    . 2011. Microbiome influences on insect host vector competence. Trends Parasitol 27:514–522. doi:10.1016/j.pt.2011.05.001.
    OpenUrlCrossRefPubMedWeb of Science
  42. 42.↵
    1. Almand EA,
    2. Moore MD,
    3. Jaykus L-A
    . 2017. Virus-bacteria interactions: an emerging topic in human infection. Viruses 9:E58. doi:10.3390/v9030058.
    OpenUrlCrossRef
  43. 43.↵
    1. Grond K,
    2. Lanctot RB,
    3. Jumpponen A,
    4. Sandercock BK
    . 2017. Recruitment and establishment of the gut microbiome in arctic shorebirds. FEMS Microbiol Ecol 93. doi:10.1093/femsec/fix142.
    OpenUrlCrossRef
  44. 44.↵
    1. Dewar ML,
    2. Arnould JPY,
    3. Dann P,
    4. Trathan P,
    5. Groscolas R,
    6. Smith S
    . 2013. Interspecific variations in the gastrointestinal microbiota in penguins. Microbiologyopen 2:195–204. doi:10.1002/mbo3.66.
    OpenUrlCrossRefPubMedWeb of Science
  45. 45.↵
    1. Hird SM,
    2. Carstens BC,
    3. Cardiff SW,
    4. Dittmann DL,
    5. Brumfield RT
    . 2014. Sampling locality is more detectable than taxonomy or ecology in the gut microbiota of the brood-parasitic Brown-headed Cowbird (Molothrus ater). PeerJ 2:e321. doi:10.7717/peerj.321.
    OpenUrlCrossRef
  46. 46.↵
    1. Risely A,
    2. Waite DW,
    3. Ujvari B,
    4. Hoye BJ,
    5. Klaassen M
    . 2018. Active migration is associated with specific and consistent changes to gut microbiota in Calidris shorebirds. J Anim Ecol 87:428–437. doi:10.1111/1365-2656.12784.
    OpenUrlCrossRef
  47. 47.↵
    1. Ganz HH,
    2. Doroud L,
    3. Firl AJ,
    4. Hird SM,
    5. Eisen JA,
    6. Boyce WM
    . 2017. Community-level differences in the microbiome of healthy wild mallards and those infected by influenza A viruses. mSystems 2:e00188-16. doi:10.1128/mSystems.00188-16.
    OpenUrlAbstract/FREE Full Text
  48. 48.↵
    1. Collins MD,
    2. Lawson PA,
    3. Willems A,
    4. Cordoba JJ,
    5. Fernandez-Garayzabal J,
    6. Garcia P,
    7. Cai J,
    8. Hippe H,
    9. Farrow JA
    . 1994. The phylogeny of the genus Clostridium: proposal of five new genera and eleven new species combinations. Int J Syst Bacteriol 44:812–826. doi:10.1099/00207713-44-4-812.
    OpenUrlCrossRefPubMed
  49. 49.↵
    1. Lu J,
    2. Santo Domingo JW,
    3. Hill S,
    4. Edge TA
    . 2009. Microbial diversity and host-specific sequences of Canada goose feces. Appl Environ Microbiol 75:5919–5926. doi:10.1128/AEM.00462-09.
    OpenUrlAbstract/FREE Full Text
  50. 50.↵
    1. Hird SM,
    2. Sánchez C,
    3. Carstens BC,
    4. Brumfield RT
    . 2015. Comparative gut microbiota of 59 neotropical bird species. Front Microbiol 6:1403. doi:10.3389/fmicb.2015.01403.
    OpenUrlCrossRef
  51. 51.↵
    1. Escallón C,
    2. Becker MH,
    3. Walke JB,
    4. Jensen RV,
    5. Cormier G,
    6. Belden LK,
    7. Moore IT
    . 2017. Testosterone levels are positively correlated with cloacal bacterial diversity and the relative abundance of Chlamydiae in breeding male rufous‐collared sparrows. Funct Ecol 31:192–203. doi:10.1111/1365-2435.12696.
    OpenUrlCrossRef
  52. 52.↵
    1. Gunnarsson G,
    2. Latorre-Margalef N,
    3. Hobson KA,
    4. Van Wilgenburg SL,
    5. Elmberg J,
    6. Olsen B,
    7. Fouchier RAM,
    8. Waldenström J
    . 2012. Disease dynamics and bird migration–linking mallards Anas platyrhynchos and subtype diversity of the influenza A virus in time and space. PLoS One 7:e35679. doi:10.1371/journal.pone.0035679.
    OpenUrlCrossRefPubMed
  53. 53.↵
    1. Lindsay LL,
    2. Kelly TR,
    3. Plancarte M,
    4. Schobel S,
    5. Lin X,
    6. Dugan VG,
    7. Wentworth DE,
    8. Boyce WM
    . 2013. Avian influenza: mixed infections and missing viruses. Viruses 5:1964–1977. doi:10.3390/v5081964.
    OpenUrlCrossRef
  54. 54.↵
    1. Caporaso JG,
    2. Lauber CL,
    3. Walters WA,
    4. Berg-Lyons D,
    5. Lozupone CA,
    6. Turnbaugh PJ,
    7. Fierer N,
    8. Knight R
    . 2011. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc Natl Acad Sci U S A 108:4516–4522. doi:10.1073/pnas.1000080107.
    OpenUrlAbstract/FREE Full Text
  55. 55.↵
    1. Frank JA,
    2. Reich CI,
    3. Sharma S,
    4. Weisbaum JS,
    5. Wilson BA,
    6. Olsen GJ
    . 2008. Critical evaluation of two primers commonly used for amplification of bacterial 16S rRNA genes. Appl Environ Microbiol 74:2461–2470. doi:10.1128/AEM.02272-07.
    OpenUrlAbstract/FREE Full Text
  56. 56.↵
    1. Lang JM,
    2. Eisen JA,
    3. Zivkovic AM
    . 2014. The microbes we eat: abundance and taxonomy of microbes consumed in a day’s worth of meals for three diet types. PeerJ 2:e659. doi:10.7717/peerj.659.
    OpenUrlCrossRef
  57. 57.↵
    1. Magoč T,
    2. Salzberg SL
    . 2011. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27:2957–2963. doi:10.1093/bioinformatics/btr507.
    OpenUrlCrossRefPubMedWeb of Science
  58. 58.↵
    1. Caporaso JG,
    2. Kuczynski J,
    3. Stombaugh J,
    4. Bittinger K,
    5. Bushman FD,
    6. Costello EK,
    7. Fierer N,
    8. Peña AG,
    9. Goodrich JK,
    10. Gordon JI,
    11. Huttley GA,
    12. Kelley ST,
    13. Knights D,
    14. Koenig JE,
    15. Ley RE,
    16. Lozupone CA,
    17. McDonald D,
    18. Muegge BD,
    19. Pirrung M,
    20. Reeder J,
    21. Sevinsky JR,
    22. Turnbaugh PJ,
    23. Walters WA,
    24. Widmann J,
    25. Yatsunenko T,
    26. Zaneveld J,
    27. Knight R
    . 2010. QIIME allows analysis of high-throughput community sequencing data. Nat Methods 7:335–336. doi:10.1038/nmeth.f.303.
    OpenUrlCrossRefPubMedWeb of Science
  59. 59.↵
    1. McMurdie PJ,
    2. Holmes S
    . 2013. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One 8:e61217. doi:10.1371/journal.pone.0061217.
    OpenUrlCrossRefPubMed
  60. 60.↵
    1. Lozupone C,
    2. Knight R
    . 2005. UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol 71:8228–8235. doi:10.1128/AEM.71.12.8228-8235.2005.
    OpenUrlAbstract/FREE Full Text
  61. 61.↵
    1. Love MI,
    2. Huber W,
    3. Anders S
    . 2014. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol 15:550. doi:10.1186/s13059-014-0550-8.
    OpenUrlCrossRefPubMed
  62. 62.↵
    1. Liaw A,
    2. Wiener M
    . 2002. Classification and regression by Random Forest. R News doi:10.1057/9780230509993.
    OpenUrlCrossRef
PreviousNext
Back to top
Download PDF
Citation Tools
The Cloacal Microbiome of Five Wild Duck Species Varies by Species and Influenza A Virus Infection Status
Sarah M. Hird, Holly Ganz, Jonathan A. Eisen, Walter M. Boyce
mSphere Oct 2018, 3 (5) e00382-18; DOI: 10.1128/mSphere.00382-18

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print
Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this mSphere article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
The Cloacal Microbiome of Five Wild Duck Species Varies by Species and Influenza A Virus Infection Status
(Your Name) has forwarded a page to you from mSphere
(Your Name) thought you would be interested in this article in mSphere.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
The Cloacal Microbiome of Five Wild Duck Species Varies by Species and Influenza A Virus Infection Status
Sarah M. Hird, Holly Ganz, Jonathan A. Eisen, Walter M. Boyce
mSphere Oct 2018, 3 (5) e00382-18; DOI: 10.1128/mSphere.00382-18
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

avian microbiome
ducks
evolutionary biology
microbial ecology
microbiome

Related Articles

Cited By...

About

  • About mSphere
  • Board of Editors
  • Policies
  • For Reviewers
  • For the Media
  • Embargo Policy
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Author Warranty
  • Types of Articles
  • Getting Started
  • Ethics
  • Contact Us

Follow #mSphereJ

@ASMicrobiology

       

 

Website feedback

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

 

American Society for Microbiology
1752 N St. NW
Washington, DC 20036
Phone: (202) 737-3600

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

Online ISSN: 2379-5042