Characterization of the nasopharyngeal microbiota in health and during rhinovirus challenge
© Allen et al.; licensee BioMed Central Ltd. 2014
Received: 1 November 2013
Accepted: 28 May 2014
Published: 25 June 2014
The bacterial communities of the nasopharynx play an important role in upper respiratory tract infections (URTIs). Our study represents the first survey of the nasopharynx during a known, controlled viral challenge. We aimed to gain a better understanding of the composition and dynamics of the nasopharyngeal microbiome during viral infection.
Rhinovirus illnesses were induced by self-inoculation using the finger to nose or eye natural transmission route in ten otherwise healthy young adults. Nasal lavage fluid samples (NLF) samples were collected at specific time points before, during, and following experimental rhinovirus inoculation. Bacterial DNA from each sample (N = 97 from 10 subjects) was subjected to 16S rRNA sequencing by amplifying the V1-V2 hypervariable region followed by sequencing using the 454-FLX platform.
This survey of the nasopharyngeal microbiota revealed a highly complex microbial ecosystem. Taxonomic composition varied widely between subjects and between time points of the same subject. We also observed significantly higher diversity in not infected individuals compared to infected individuals. Two genera – Neisseria and Propionibacterium – differed significantly between infected and not infected individuals. Certain phyla, including Firmicutes, Actinobacteria, and Proteobacteria, were detected in all samples.
Our results reveal the complex and diverse nature of the nasopharyngeal microbiota in both healthy and viral-challenged adults. Although some phyla were common to all samples, differences in levels of diversity and selected phyla were detected between infected and uninfected participants. Deeper, species-level metagenomic sequencing in a larger sample is warranted.
KeywordsMicrobiota nasopharynx longitudinal sequencing rhinovirus illness
The nasopharynx is frequently colonized by both commensal and pathogenic bacteria. Many pathogenic species, including Streptococcus pneumonia, Haemophilus influenza, Moraxella catarrahlis, Staphylococcus aureus, and Neisseria meningitidis, exist in the nasopharynx of apparently healthy individuals[1–4]. Ling et al. found the nasopharyngeal microbiota to be distinct from other body sites surveyed (saliva, dominant hands, and feces) in healthy Chinese adults.
Understanding the relationships between microbes of the upper respiratory tract (URT) during perturbations is anticipated to provide insights into the pathogenesis of URT infections. Studies of the nasopharyngeal microbiota in children have observed changes due to season (winter/fall versus spring) and treatment with antimicrobials or the heptavalent conjugated pneumococcal polysaccharide vaccine[6, 7]. Specific commensal taxa have been negatively associated with colonization of known pathogenic bacteria and with acute otitis media in children, and these relationships changed depending on antibiotic usage.
Like pathogenic bacteria, viruses, including rhinoviruses, enteroviruses, coronaviruses, and adenoviruses, have been found in asymptomatic, healthy individuals[9, 10]. Infections of the URT display complex interactions between bacterial pathogens and viruses, both synergistic and competitive[11–14]. Viruses in the URT can alter bacterial adherence, bacterial colonization, and host immune response[17–19]. A study of the nasopharyngeal microbiota of children with severe bronchitis showed microbiotal shifts depending on the viral infection (human rhinovirus (HRV) only, respiratory syncytial virus (RSV) only, or co-infection). The lung microbiota in chronic obstructive pulmonary disease (COPD) patients and healthy individuals during a rhinovirus challenge showed significant changes after infection.
We previously used microarray technology to investigate bacteria in the adult nasopharynx prior to, during, and after experimental rhinovirus infection. The microarray approach allowed us to detect bacteria at the species level, but we were only able to test the presence or absence of about 60 clinically informative bacteria. The goal of the present study was to characterize the relative abundance of microbial communities of the nasopharynx at the phylum and genus levels, by sequencing the 16S rRNA gene of bacteria in nasal lavage fluid (NLF) samples from the same adults during rhinovirus challenge.
Assessment of infection
This protocol has previously been described. Briefly, a subject was considered infected if the challenge rhinovirus was detected in the NLF at least once during the five days following inoculation. Additionally, the serum antibody response to the challenge rhinovirus type 39 was examined in sera obtained prior to inoculation and three weeks following inoculation by standard methods. Subjects with a 4-fold increase in antibody titer to the challenge virus were considered to be infected. If either of these tests were positive, the subject was considered infected. Assessment of illness of the participants is described in Additional file1.
DNA was isolated from NLF using bead beating and phenol-chloroform extraction methods. This protocol was optimized from Ren T, et al. to increase DNA output from a more difficult sample type (NLF). Briefly, 1.8 mL of nasal wash was transferred to a 2-mL BeadBeater tube (BioSpec, Bartlesville, OK) and centrifuged at 2,348 g (5,000 rpm) for 10 minutes. Supernatant (1.3 mL) was aspirated, leaving all pelleted bacterial cells plus 500 μL supernatant for DNA extraction and subsequent sequencing. Then, 20 μL Proteinase K (20 mg/mL) (Roche, Basel, Switzerland) was added and incubated at 60˚C for two hours then 95˚C for 10 minutes. A total of 200 μL 10% SDS (Fischer Scientific, Waltham, MA), 400 μL of 0.1 mm zirconia/silica beads (BioSpec, Bartlesville, OK), 100 μL of buffer (10 mM Tris-HCl pH 7.5, 5 mM EDTA, 100 mM NaCl), and 500 μL of phenol:chloroform:isoamyl alcohol (25:24:1) (Fischer Scientific, Waltham, MA) were added to the tube. The sample was homogenized on a Biospec Mini-BeadBeater-8 for 3 minutes and then centrifuged for 10 minutes at 10,000 g. The top aqueous layer was transferred to a new Phase-Lock-Gel tube (VWR, Radnor, PA) with an equal volume of chloroform:isoamyl alcohol (24:1) (Acros Organics, Fischer Scientific, Waltham, MA) and centrifuged at 20,000 g for 10 minutes. This step was repeated twice to ensure removal of phenol. The sample was then ethanol precipitated, and the DNA pellet was resuspended in 50 μL 10 mM Tris-HCl pH 8.0. Each set of extractions was accompanied by a negative control to ensure no contaminants in the extracted DNA.
Tag-PCR amplification of the V1-V2 regions of the bacterial 16S rRNA gene and pyrosequencing
This protocol has been previously described. Briefly, the V1-V2 hypervariable regions of the 16S rRNA gene were amplified from extracted DNA samples using two primers containing the universal sequences 27 F (5′-AGRGTTTGATCMTGGCTCAG-3′) and 534R (5′-TTACCGCGGCTGCTGGCAC-3′). Each sample was tagged with a unique barcode (10 bp) added to the 5′ end of the forward primer. Conditions for amplification were 94˚C for 3 minutes, then 30 cycles of 94˚C for 30 seconds, 57˚C for 45 seconds, and 72˚C for 60 seconds with a final extension step of 72˚C for 5 minutes. Each set of samples was amplified with a negative control to ensure no contaminants from reagents. PCR products were then run on a gel to verify amplification of the correct size. Amplicons from all samples were quantified using the Qubit® 2.0 Fluorometer and then pooled in equal molar ratios with 49 samples in pool 1 and 48 samples in pool 2. To achieve a final pool of 20 μL of sample with a 10 ng/μL concentration, 40.8 ng of each sample were added to pool 1 and 41.67 ng were added from each sample to pool 2. The pooled samples were then gel-purified and the concentration of each sample was measured using the Qubit 2.0 Fluorometer to ensure the minimum of 20 μL of 10 ng/μL sample for sequencing. Each pooled sample was sequenced using titanium chemistry on a 454 Life Science Genome Sequencer FLX platform at the University of Virginia Department of Biology Genome Core Facility using the standard protocol for sequencing. The samples were run on an Agilent 2100 Bioanalyzer High-Sensitivity chip to determine the size distribution of the library samples. Then, a KAPA Biosystems qPCR assay was run to determine the effective concentration of the samples.
Quality assessment and filtering of sequences
This quality of the read sequences was assessed using FastQC v0.10.1. This assessment revealed that a significant fraction of our reads were identical short (<100 bp) reads (homopolymers). Based on this assessment, and the distribution of read lengths, we filtered the sequences by length, removing any reads shorter than 100 bp in length, or greater than 600 bp. This filter alone reduced our initial dataset from 4.32 Mbp, to 3.32 Mbp (approximately 23.1% reduction). In addition, reads with an average quality score lower than 25 (phred scale) were removed, as were any reads that contained 50-base windows with average quality scores below that threshold. Qiime’s default filters also removed any reads with mismatches in the primer or barcode sequences, and any reads with more than six ambiguous bases.
Clustering and filtering of operational taxonomic units (OTUs)
Sequences were clustered into operational taxonomic units (OTUs) (97% identity cutoff) using the uclust algorithm v1.2.22 (de novo clustering) in QIIME v1.8, and assigned taxonomy using RDP classifier v2.2, trained with the Greengenes (gg_13_8) dataset. OTUs were removed from the dataset if they could not be classified at the domain level, or were classified as archaea, or were classified as chloroplast sequences. OTUs estimated to be putatively chimeric by the uchime algorithm v4.2.40 were also filtered out, as were singleton OTUs. After filtering, we were left with a mean of 1021.6 sequences per sample (SD 630.84, median 913, minimum 134, maximum 4,501). Samples containing fewer sequences than three SD below the mean count (146 sequences/sample) (n = 1) were eliminated from downstream analyses.
Analyses of taxonomic composition and diversity were performed using QIIME, version 1.8. Default parameters were used for the core analyses. OTU representatives were aligned with pyNAST v1.2.2, and phylogenies were generated using FastTree v2.1.3. Alpha diversity (within-population diversity) of samples was measured using counts of observed species (OTUs), the chao1 estimator for species richness, and the Shannon diversity index, which estimates total diversity taking into account both species richness and evenness. Beta diversity (diversity between populations) was calculated using UniFrac distances (unweighted) between samples (based on the relative abundance of OTUs), and visualized using principal coordinates analysis (PCoA). In order to compare diversity of OTUs in the nasopharynx to that of other body sites, we also performed a combined beta diversity analysis of our nasopharyngeal samples with samples from several different body sites (gut, oral cavity, external auditory canal, nostril, and hair) from Costello et al., as well as adenoid samples from Ren et al..
Analyses of phylum abundance between inoculation states (before, during, and after) and between subjects, was performed using analysis of variance (ANOVA). For significant results, the Tukey Kramer test was subsequently conducted. ANOVA was also used for the analyses of class, order, and family abundances between inoculation states (before, during, and after). Comparisons of alpha diversity between the infected and non-infected sample groups were done using Mann-Whitney Wilcoxon tests, as were the abundances of selected genera between the infected and non-infected samples, and the comparison of inter-individual and intra-individual Unifrac distances.
Route of inoculation
Criteria for infection status
Analysis of nasopharyngeal bacterial richness and composition by sequencing
The most abundant genus in the nasopharyngeal microbiota across samples was Alloiococcus (21.2% of total abundance), followed by Corynebacterium (14.7%), Staphylococcus (8.9%), Haemophilus (5.3%), Propionibacterium (4.6%), and Streptococcus (3.2%). All other classified OTUs belonged to genera comprising less than 2% of the total abundance. It should be noted that these abundances reflect total abundance across all samples. Taxonomic composition varied widely between samples (Additional file 2).
Bacterial community structure of infected and not infected subjects
Comparisons of genus abundance between infected and uninfected samples, using the Mann-Whintey Wilcoxon test
Mean Relative Abundance
This study is the first to survey and describe the microbiota of the nasopharynx in healthy adults during a rhinovirus challenge. We found a highly complex microbiome that varied between individuals, and significant differences in two genera (Neisseria and Propionibacterium) between infected and not infected individuals. Strengths of our study include a controlled rhinovirus study, longitudinal analyses, and use of high-throughput sequencing technology. Limitations include the modest sample size and sequencing depth.
The results from this survey have been compared to the results of the microarray study  using the same NLF samples. Due to the design of the microarray study, it can only detect the presence or absence of bacterial species on the panel included, while the 16S method is able to determine relative abundance of all bacteria in the sample. We can only compare the number of mismatches in those genera that the microarray was able to detect (n = 8). Out of the total 96 samples, the percent of samples which were detected in one study but not the other ranged from 4% to 53% in each genus. These results can be due to two aliquots from the same NLF sample extracted using two different methods. It could also be due to the microarray detecting species in very low abundance because of enrichment due to PCR.
Due to the level of complexity of the nasopharyngeal microbiota, it is possible that significant changes in the bacterial community before, during, and after infection could likely be detected with deeper sequencing and/or recruiting more subjects. Some large changes in relative abundances of phyla within individuals at close time points were observed. It will be important to investigate whether these reflect true biological fluctuations in this dynamic ecological niche or low microbial density via greater sequencing depth. Additionally, using true shotgun metagenomics and taxonomic classification  to distinguish bacteria at the species level may allow significant changes in pathogenic and opportunistic bacteria during infection to be discovered. Information on the species level is especially important because, in numerous cases, there are pathogenic and commensal bacteria within the same genus, for example, Haemophilius influenzae, a known pathogen, and its closest phylogenetic relative Haemophilus haemolyticus, a commensal bacterium [37, 38].
PCoA of the nasopharyngeal microbiota samples along with samples from various other body habitats shows logical clustering with similar locations throughout the URT. The nasopharynx, adenoid, ear, hair, and nostril samples all formed a single large cluster. While this cluster was distinct from the separate clusters formed by the gut and oral samples, there was not a clear differentiation between the other body sites. Bacterial niches in close proximity have some sharing of bacteria due to physiological events like sneezing and nose blowing [39–41], which could help to explain the overlap between the microbiota of the nasopharynx and other proximal body habitats. Differences in age and health status of subjects along with season of sample collection may impact results. Participants in the Costello, et al. study included healthy adults, with two samples collected during June and two during the following September . Adenoids in the Ren, et al. study were from children undergoing adenoidectomy, collected year round . Demographic data was not accessible because these were discarded surgical specimens. NLF samples in our study were collected during the fall season in healthy adults. Future studies will need to sample complementary niches in the URT to understand coordinated changes of bacterial communities in response to viral infection.
In our study, we surveyed the microbiota of the healthy nasopharynx during viral infection and cataloged a highly diverse population of microbes. The nine dominant phyla of the adenoid microbiota  contained the dominant phyla found in our study including Firmicutes, Proteobacteria, and Actinobacteria. We found the same dominant phyla as a study of nasopharyngeal microbiota of healthy Chinese adults . We did not observe an effect of viral infection on microbial composition, which is likely due to a small population, insufficient sequencing depth to detect the full complexity of the microbiota, and the inability to classify OTUs to the species level. In addition to deeper sequencing to achieve species level information, shotgun metatranscriptomics could be used as a complementary approach, to characterize the changes in gene expression that the invasion of a virus can elicit in a microbial community . This information would help us to address the open question of the functional role that the microbiota plays in viral infections. It is uncertain, for instance, whether the microbiota can help to protect the host from viral infection by providing a barrier or by initiating an immune response. Alternatively, it is also possible that the microbiota can aid a virus in invading the host mucosa. One metagenomic study of the airway DNA virome of cystic fibrosis (CF) and non-CF individuals showed that the metabolic profiles of CF and non-CF individuals were distinctly different . This study used the approach of finding target pathways that have altered expression levels in the diseased state, rather than taxonomic composition differences of the microbiota to discover novel therapeutic strategies. Our study has highlighted the highly diverse microbial composition of the adult nasopharynx in both healthy individuals and during viral infection.
Our finding that each individual had a unique nasopharyngeal microbiota broadened our views on how to study and treat URTIs. Most prior studies have looked at a population of individuals with a specific URTI or complication (that is, sinusitis) and sampled each individual for specific pathogens. From the findings of such studies, conclusions may be drawn as to which pathogens are likely involved and targeted treatments may be created. However, with a new understanding of the complexity of the NP microbiota at the individual level, we can apply this understanding by recognizing that future studies need to sample patients with URTIs or complications at baseline comprehensively, and then during acute disease to determine which bacteria respond, expand, disappear. This approach is anticipated to lead to improved understanding of disease and more tailored approaches to treatments of these diseases.
In conclusion, we have conducted a longitudinal survey of the nasopharyngeal microbiota of healthy adults during rhinovirus challenge. This microbiota is highly diverse and varies greatly between individuals. The microbiota of each individual differs between time points throughout the viral challenge, though variability within individuals was less than between individuals. No significant changes in bacterial presence or relative abundance between time points were found due to viral infection. Further studies are needed to further characterize this microbiota and the response of the microbiota to viral infection down to the species level using deeper sequencing. Additionally, metagenomic sequencing will give us insight into transcriptional changes of the host and microbiota during viral infection. With this knowledge, we will better understand the pathogenesis of URT infections.
Availability of supporting data
The additional files and data set supporting the results of this article are available in Figshare, http://dx.doi.org/10.6084/m9.figshare.832471.
analysis of variance
chronic obstructive pulmonary disease
nasal lavage fluid - the sample type used to comprehensively survey the microbiota of the nasopharynx
operational taxonomic unit - approximations of bacterial species based on a percent similarity threshold for DNA sequence identity at a specified locus (usually 16S rRNA). PCoA: principal coordinates analysis - method used to visualize beta diversity between samples using UniFrac distances
upper respiratory tract - URT refers to the parts of the respiratory tract including the nasal cavity the pharynx, and the larynx
upper respiratory tract infection - general term used to refer to acute infections of the URT which are usually caused by viral infection.
We thank Kathleen Ashe for her technical assistance throughout this study. The rhinovirus inoculum pool was kindly provided by RB Turner. This study was supported in part by the Pendleton Pediatric Infectious Disease Laboratory at the University of Virginia (extraction, PCR, and sequencing) and by an unrestricted grant from NexBio, Inc (rhinovirus challenge study and analysis of sequencing data). Ms Allen received support from the National Institute of General Medical Sciences of the National Institutes of Health under Award Number T32GM008715.
- Chi DH, Hendley JO, French P, Arango P, Hayden FG, Winther B: Nasopharyngeal reservoir of bacterial otitis media and sinusitis pathogens in adults during wellness and viral respiratory illness. Am J Rhinol. 2003, 17: 209-214.PubMedGoogle Scholar
- Hendley JO, Hayden FG, Winther B: Weekly point prevalence of Streptococcus pneumoniae, Hemophilus influenzae and Moraxella catarrhalis in the upper airways of normal young children: effect of respiratory illness and season. APMIS. 2005, 113: 213-220.View ArticlePubMedGoogle Scholar
- Puhakka T, Mäkelä MJ, Malmström K, Uhari M, Savolainen J, Terho EO, Pulkkinen M, Ruuskanen O: The common cold: effects of intranasal fluticasone propionate treatment. J Allergy Clin Immunol. 1998, 101: 726-731.View ArticlePubMedGoogle Scholar
- Bogaert D, Keijser B, Huse S, Rossen J, Veenhoven R, van Gils E, Bruin J, Montijn R, Bonten M, Sanders E: Variability and diversity of nasopharyngeal microbiota in children: a metagenomic analysis. PLoS One. 2011, 6: e17035-View ArticlePubMedPubMed CentralGoogle Scholar
- Ling Z, Liu X, Luo Y, Yuan L, Nelson KE, Wang Y, Xiang C, Li L: Pyrosequencing analysis of the human microbiota of healthy Chinese undergraduates. BMC Genomics. 2013, 14: 390-View ArticlePubMedPubMed CentralGoogle Scholar
- Hilty M, Qi W, Brugger SD, Frei L, Agyeman P, Frey PM, Aebi S, Muhlemann K: Nasopharyngeal microbiota in infants with acute otitis media. J Infect Dis. 2012, 205: 1048-1055.View ArticlePubMedGoogle Scholar
- Biesbroek G, Wang X, Keijser BJ, Eijkemans RM, Trzcinski K, Rots NY, Veenhoven RH, Sanders EA, Bogaert D: Seven-valent pneumococcal conjugate vaccine and nasopharyngeal microbiota in healthy children. Emerg Infect Dis. 2014, 20: 201-210.View ArticlePubMedPubMed CentralGoogle Scholar
- Pettigrew MM, Laufer AS, Gent JF, Kong Y, Fennie KP, Metlay JP: Upper respiratory tract microbial communities, acute otitis media pathogens, and antibiotic use in healthy and sick children. Appl Environ Microbiol. 2012, 78: 6262-6270.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang Q, Garrity GM, Tiedje JM, Cole JR: Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007, 73: 5261-5267.View ArticlePubMedPubMed CentralGoogle Scholar
- van der Zalm MM, van Ewijk BE, Wilbrink B, Uiterwaal CS, Wolfs TF, van der Ent CK: Respiratory pathogens in children with and without respiratory symptoms. J Pediatr. 2009, 154: 396-400. 400 e391View ArticlePubMedGoogle Scholar
- Murphy TF, Bakaletz LO, Smeesters PR: Microbial interactions in the respiratory tract. Pediatr Infect Dis J. 2009, 28: S121-S126.View ArticlePubMedGoogle Scholar
- Peltola VT, McCullers JA: Respiratory viruses predisposing to bacterial infections: role of neuraminidase. Pediatr Infect Dis J. 2004, 23: S87-S97.View ArticlePubMedGoogle Scholar
- van den Bergh MR, Biesbroek G, Rossen JW, de Steenhuijsen Piters WA, Bosch AA, van Gils EJ, Wang X, Boonacker CW, Veenhoven RH, Bruin JP: Associations between pathogens in the upper respiratory tract of young children: interplay between viruses and bacteria. PLoS One. 2012, 7: e47711-View ArticlePubMedPubMed CentralGoogle Scholar
- Brogden KA, Guthmiller JM, Taylor CE: Human polymicrobial infections. Lancet. 2005, 365: 253-255.View ArticlePubMedGoogle Scholar
- Avadhanula V, Rodriguez CA, Devincenzo JP, Wang Y, Webby RJ, Ulett GC, Adderson EE: Respiratory viruses augment the adhesion of bacterial pathogens to respiratory epithelium in a viral species- and cell type-dependent manner. J Virol. 2006, 80: 1629-1636.View ArticlePubMedPubMed CentralGoogle Scholar
- Tong HH, Weiser JN, James MA, DeMaria TF: Effect of influenza A virus infection on nasopharyngeal colonization and otitis media induced by transparent or opaque phenotype variants of Streptococcus pneumoniae in the chinchilla model. Infect Immun. 2001, 69: 602-606.View ArticlePubMedPubMed CentralGoogle Scholar
- McNamee LA, Harmsen AG: Both influenza-induced neutrophil dysfunction and neutrophil-independent mechanisms contribute to increased susceptibility to a secondary Streptococcus pneumoniae infection. Infect Immun. 2006, 74: 6707-6721.View ArticlePubMedPubMed CentralGoogle Scholar
- Colamussi ML, White MR, Crouch E, Hartshorn KL: Influenza A virus accelerates neutrophil apoptosis and markedly potentiates apoptotic effects of bacteria. Blood. 1999, 93: 2395-2403.PubMedGoogle Scholar
- Abramson JS, Hudnor HR: Role of the sialophorin (CD43) receptor in mediating influenza A virus-induced polymorphonuclear leukocyte dysfunction. Blood. 1995, 85: 1615-1619.PubMedGoogle Scholar
- Hyde ER, Petrosino JF, Piedra PA, Camargo CA, Espinola JA, Mansbach JM: Nasopharyngeal Proteobacteria are associated with viral etiology and acute wheezing in children with severe bronchiolitis. J Allergy Clin Immunol. 2014, 133: 1220-1222.View ArticlePubMedPubMed CentralGoogle Scholar
- Molyneaux PL, Mallia P, Cox MJ, Footitt J, Willis-Owen SA, Homola D, Trujillo-Torralbo MB, Elkin S, Kon OM, Cookson WO, Moffatt MF, Johnston S: Outgrowth of the bacterial airway microbiome after rhinovirus exacerbation of chronic obstructive pulmonary disease. Am J Respir Crit Care Med. 2013, 188: 1224-1231.View ArticlePubMedPubMed CentralGoogle Scholar
- Allen EK, Pitkaranta A, Maki M, Hendley JO, Laakso S, Sale MM, Winther B: Bacteria in the nose of young adults during wellness and rhinovirus colds: detection by culture and microarray methods in 100 nasal lavage specimens. Int Forum Allergy Rhinol. 2013, 3: 731-739.View ArticlePubMedGoogle Scholar
- Hamparian VV: Rhinoviruses. Diagnostic procedures for viral, rickettsial and chlamydial infections. Edited by: Lennette EH, Schmidt NJ. 1979, Washington, D.C: American Public Health Association, 562-Google Scholar
- Ren T, Glatt DU, Nguyen TN, Allen EK, Early SV, Sale M, Winther B, Wu M: 16S rRNA survey revealed complex bacterial communities and evidence of bacterial interference on human adenoids. Environ Microbiol. 2013, 15: 535-547.View ArticlePubMedGoogle Scholar
- FastQC A quality control tool for high throughput sequence data.http://www.bioinformatics.babraham.ac.uk/projects/fastqc/,
- Edgar RC: Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010, 26: 2460-2461.View ArticlePubMedGoogle Scholar
- Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Pena AG, Goodrich JK, Gordon JI, Huttley GA, Kelley ST, Knights D, Koenig JE, Ley RE, Lozupone CA, McDonald D, Muegge BD, Pirrung M, Reeder J, Sevinsky JR, Turnbaugh PJ, Walters WA, Widmann J, Yatsunenko T, Zanveld J, Knight R:QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010, 7: 335-336.View ArticlePubMedPubMed CentralGoogle Scholar
- DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, Huber T, Dalevi D, Hu P, Andersen GL: Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006, 72: 5069-5072.View ArticlePubMedPubMed CentralGoogle Scholar
- Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R: UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011, 27: 2194-2200.View ArticlePubMedPubMed CentralGoogle Scholar
- Caporaso JG, Bittinger K, Bushman FD, DeSantis TZ, Andersen GL, Knight R: PyNAST: a flexible tool for aligning sequences to a template alignment. Bioinformatics. 2010, 26: 266-267.View ArticlePubMedPubMed CentralGoogle Scholar
- Price MN, Dehal PS, Arkin AP: FastTree 2–approximately maximum-likelihood trees for large alignments. PLoS One. 2010, 5: e9490-View ArticlePubMedPubMed CentralGoogle Scholar
- Lozupone C, Knight R: UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol. 2005, 71: 8228-8235.View ArticlePubMedPubMed CentralGoogle Scholar
- Ramette A: Multivariate analyses in microbial ecology. FEMS Microbiol Ecol. 2007, 62: 142-160.View ArticlePubMedPubMed CentralGoogle Scholar
- Costello EK, Lauber CL, Hamady M, Fierer N, Gordon JI, Knight R: Bacterial community variation in human body habitats across space and time. Science. 2009, 326: 1694-1697.View ArticlePubMedPubMed CentralGoogle Scholar
- Biesbroek G, Sanders EA, Roeselers G, Wang X, Caspers MP, Trzcinski K, Bogaert D, Keijser BJ: Deep sequencing analyses of low density microbial communities: working at the boundary of accurate microbiota detection. PLoS One. 2012, 7: e32942-View ArticlePubMedPubMed CentralGoogle Scholar
- Segata N, Waldron L, Ballarini A, Narasimhan V, Jousson O, Huttenhower C: Metagenomic microbial community profiling using unique clade-specific marker genes. Nat Methods. 2012, 9: 811-814.View ArticlePubMedPubMed CentralGoogle Scholar
- Murphy TF, Brauer AL, Sethi S, Kilian M, Cai X, Lesse AJ: Haemophilus haemolyticus: a human respiratory tract commensal to be distinguished from Haemophilus influenzae. J Infect Dis. 2007, 195: 81-89.View ArticlePubMedGoogle Scholar
- McCrea KW, Xie J, LaCross N, Patel M, Mukundan D, Murphy TF, Marrs CF, Gilsdorf JR: Relationships of nontypeable Haemophilus influenzae strains to hemolytic and nonhemolytic Haemophilus haemolyticus strains. J Clin Microbiol. 2008, 46: 406-416.View ArticlePubMedPubMed CentralGoogle Scholar
- Gwaltney JM, Hendley JO, Phillips CD, Bass CR, Mygind N, Winther B: Nose blowing propels nasal fluid into the Paranasal Sinuses. Clin Infect Dis. 2000, 30: 387-391.View ArticlePubMedGoogle Scholar
- Schwartz RH, Pitkaranta A: Computed tomography imaging of the Maxillary and Ethmoid Sinuses in children with short-duration purulent Rhinorrhea. Otolaryngol Head Neck Surg. 2001, 124: 160-163.View ArticlePubMedGoogle Scholar
- Winther B, Gwaltney JM, Phillips CD, Hendley JO: Radiopaque contrast dye in nasopharynx reaches the middle ear during swallowing and/or yawning. Acta Otolaryngol. 2005, 125: 625-628.View ArticlePubMedGoogle Scholar
- Segata N, Boernigen D, Tickle TL, Morgan XC, Garrett WS, Huttenhower C: Computational meta’omics for microbial community studies. Mol Syst Biol. 2013, 9: 666-View ArticlePubMedPubMed CentralGoogle Scholar
- Willner D, Furlan M, Haynes M, Schmieder R, Angly FE, Silva J, Tammadoni S, Nosrat B, Conrad D, Rohwer F: Metagenomic analysis of respiratory tract DNA viral communities in cystic fibrosis and non-cystic fibrosis individuals. PLoS One. 2009, 4: e7370-View ArticlePubMedPubMed CentralGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.