Microbial Communities in the Upper Respiratory Tract of Patients with Asthma and Chronic Obstructive Pulmonary Disease

Respiratory infections are well-known triggers of chronic respiratory diseases. Recently, culture-independent tools have indicated that lower airway microbiota may contribute to pathophysiologic processes associated with asthma and chronic obstructive pulmonary disease (COPD). However, the relationship between upper airway microbiota and chronic respiratory diseases remains unclear. This study was undertaken to define differences of microbiota in the oropharynx of asthma and COPD patients relative to those in healthy individuals. To account for the qualitative and quantitative diversity of the 16S rRNA gene in the oropharynx, the microbiomes of 18 asthma patients, 17 COPD patients, and 12 normal individuals were assessed using a high-throughput next-generation sequencing analysis. In the 259,572 total sequence reads, α and β diversity measurements and a generalized linear model revealed that the oropharynx microbiota are diverse, but no significant differences were observed between asthma and COPD patients. Pseudomonas spp. of Proteobacteria and Lactobacillus spp. of Firmicutes were highly abundant in asthma and COPD. By contrast, Streptococcus, Veillonella, Prevotella, and Neisseria of Bacteroidetes dominated in the healthy oropharynx. These findings are consistent with previous studies conducted in the lower airways and suggest that oropharyngeal airway microbiota are important for understanding the relationships between the various parts of the respiratory tract with regard to bacterial colonization and comprehensive assessment of asthma and COPD.


Introduction
Asthma and chronic obstructive pulmonary disease (COPD) are common chronic inflammatory diseases of the airways, and their prevalence and mortality rates continue to increase, establishing them as leading causes of disability worldwide [1]. The World Health Organization (WHO) estimates that 235 million people are annually affected by asthma and approximately 250,000 people die per year worldwide [2]. COPD also affects 329 million people worldwide, resulting in approximately 3 million deaths annually [3,4].
Clinical studies of chronic respiratory diseases have been focused on their association with common respiratory viruses. Respiratory viral infections are common and usually self-limiting illnesses in healthy adults but are a major cause of exacerbations in patients with asthma and COPD [5]. One of the triggers for development of asthma and COPD is impaired lung function, caused due to bacterial infections. In asthma, Streptococcus pneumoniae, Haemophilus influenzae, and Moraxella catarrhalis are clinically relevant contributors to asthma exacerbations, particularly when sinusitis is present [1,2]. Atypical bacteria such as Chlamydophila pneumoniae and Mycoplasma pneumoniae may also be associated with asthma [6][7][8]. The role that bacterial pathogens play in acute exacerbations of COPD is controversial. Pseudomonas aeruginosa, Haemophilus influenzae, Moraxella catarrhalis, and Streptococcus pneumoniae are known to be associated with a significantly increased risk of exacerbation of COPD. However, the prevalence of bacterial pathogens in the lower airways was found to be the same during acute exacerbations and stable disease; thus, changes in bacterial load are unlikely to be an important mechanism underlying exacerbation of COPD [9][10][11].
Human microbiota is diverse across individuals and different sites in the body. Oral microbiota are known to be involved in upper and lower respiratory infections, which can also develop into atopic airway diseases such as allergic rhinitis and bronchial asthma [12]. The oropharynx is constantly exposed to both inhaled and ingested microbes, including those cleared by mucociliary mechanisms from the upper and lower respiratory tracts and those contained in saliva [13]. The clinical relevance of the microbiota of the oropharyngeal wall is unknown. However, it is known that the diversity of the oral microbiota can be affected by antibiotics, probiotics, diet, and gut microbiota [14].
Oral microbiota can trigger and strike as a critical component in conditions like Bronchial asthma and COPD [15], particularly considering that the normal healthy lung is not a sterile organ as previously assumed [16]. To date, the relationship between airway microbiota and chronic respiratory diseases has been investigated only in lower airway samples such as sputum and bronchoalveolar lavage (BAL) samples [17,18,19]. However, the microbial inhabitants of the upper respiratory tract in asthma and COPD patients has been insufficiently characterized. The induction or suppression of systemic immune tolerance to challenge antigens can be attributed to the link between oral tolerance and airway tolerance. Although strong evidence exists to implicate bacterial infections in the course and pathogenesis of airway diseases, systematic studies of organisms in the airways are lacking [12]. So it becomes crucial to define the characteristics of the oropharyngeal microbiota, unfolding the pathogenic roles of these microorganisms. This in turn can lead to identification of potential targets and the discovery of novel therapies for the prevention and treatment of oral complications [20].
Traditionally identification of bacteria in the oral cavity has been done by culturing methods. Analysis of microbiota through culture-dependent methodologies is both expensive and timeconsuming in comparison with many of the available cultureindependent methods. It is evident to assess minor or nondominant populations, culture-independent strategies are important for analysis of microbial ecology. Previous studies demonstrate culture-independent approaches have revealed more diverse microbiota compared to culture-dependent microbiota, of which nearly 40-60% bacterial taxa in the mouth were uncultured or validly described [21].
Recent high throughput methods like next-generation sequencing have revolutionized the way diversity of microbial communities in varied environments, including the human body is analyzed [22]. In this study, we aimed to comprehensively characterize the bacterial communities of the oropharyngeal airway in asthma and COPD patients relative to healthy individuals using a highthroughput 16S rRNA gene-pyrosequencing analysis.

Study subjects
All of the samples were collected from individuals with newly diagnosed asthma and COPD in accordance with the ethical guidelines mandated by the commission a protocol (number #12-0015) that was approved by the University's Human Subjects Institutional Review Board (IRB) of Chung-Ang University College of Medicine, Seoul, Korea. Human subjects remain anonymous and the informed consent was written. All participants gave written informed consent for participation in the study. The current diagnosis of asthma was based on a history of wheezing, shortness of breath, and cough, which were variable in severity and frequency [23,24]. COPD patients were further classified into four distinct stages according to Global Initiative on Obstructive Lung Disease (GOLD) guidelines [25]. Patients with bronchial asthma had mild intermittent disease and were treated as being naive. For the oropharyngeal swabs, 3M quick swabs (3M Microbiology Products; St. Paul, MN, USA) were inserted into the oral cavity until resistance was met at the oropharynx. The swabs were then rotated 180 degrees and withdrawn while taking care not to contaminate the swabs with oral bacteria. After swabbing, the swabs were transported immediately to the lab for identification and subsequent microbiomic analysis. Each swab sample was aseptically placed in a microtube and centrifuged at 14,000 rpm for 10 min, and the supernatant was removed. To detect possible contamination, negative controls were prepared and subjected to the same procedures.

DNA extraction from oropharyngeal swab samples
The first step in the analysis was extraction of DNA from the bacterial pellet that was obtained from oropharyngeal swab samples, followed by polymerase chain reaction (PCR) analysis. Samples were extracted individually by using the cetyltrimethylammonium bromide method [26]. Purified DNA was dissolved in sterile water with 40 mg/mL RNase A and quantified using an Infinite 200 NanoQuant microplate reader (Tecan; Mä nnedorf, Switzerland) at a wavelength of 260 nm.  Table 3. Abundance table of normal, asthma and COPD. GS-FLX 454 pyrosequencing by using 16S rRNA gene amplification PCR amplification of the 16S rRNA gene was performed following established procedures [27]. The V1-V3 region of the 16S rRNA gene was amplified with primers 8F (59-CTGCTGCCTYCCGTA- 39) and 530R (59-GTAT-TACCGCGGCTGCTG-39) with the MID sequence. PCR was performed in final reaction mixtures of 25 mL containing 5-25 ng of genomic DNA, each primer at 0.4 mM, 0.2 mM dNTPs (Takara Bio; Shiga, Japan), 1.5 mM MgCl 2 (Invitrogen; Carlsbad, CA, USA), 2.0 U of Hot Start Taq polymerase (Takara), and 1.0 mL of reaction buffer (Invitrogen). PCR amplification was performed on a GeneAmp PCR system 9700 thermocycler (Applied Biosystems; Foster City, CA, USA) under the following conditions: initial denaturation for 5 min at 94uC, followed by 30 cycles of 30 s at 94uC, 30 s at 55uC, and 1 min at 72uC, with a final extension for 10 min at 72uC and cooling to 4uC. The PCR products were resolved on a 1.2% Seakem LE agarose gel (FMC Bioproducts; Rockland, ME, USA) and visualized after ethidium bromide staining. 454 sequencing was carried out on the pooled PCR products of equal amounts. Amplicon pyrosequencing was performed using a 454/Roche GS-FLX Titanium instrument (Roche Molecular Systems; Branchburg, NJ, USA). After sequencing, the sequence file (.fna) and quality file (.qual) were recovered by MID groups. Quality control of the post-processing sequences was performed to eliminate sequences with ambiguous base calls, homopolymers .6, quality scores ,Q25, read lengths ,250 bp, and low-quality sequence ends and those without a 100% match to barcodes. Random sampling was used to analyze the diversity of the low reads of the sequence and the high reads of the sequence set. Tests of the sequence analysis methods are described in PRINSEQ [28]. The raw data were submitted to the National Center for Biotechnology Information-Sequence Read Archive (NCBI-SRA; http://www.ncbi.nlm.nih.gov/sra) with accession number SRX316595.

16S rRNA gene sequence analysis
The nucleotide sequence homology of the amplified 16S rRNA gene V1-V3 regions was determined using the nucleotide BLAST (BLASTn) program, which is available at the SILVA website (http://www.arb-silva.de/, Release 115) and the BLAST webpage of the NCBI. Accurate identification of organisms through comparative analysis of 16S rRNA gene sequences is strongly dependent on the quality of the database that is used. The curated Ribosomal Database Project (RDP; http://rdp.cme.msu.edu, Release10, update 32) provides bacterial and archaeal small subunit rRNA gene sequences in an aligned and annotated format, and has achieved major improvements in the detection of sequence anomalies [29]. Among all of the online tools provided by the RDP website, the RDP pyrosequencing pipeline includes information of Shannon-Wiener, Chao1, and ACE and operational taxonomic units (OTUs) through the rarefaction analysis tool. A phylogenetic tree was built using the neighbor-joining method in Clustal_X [30] with 1,000 bootstrap replicates. Fast UniFrac was used to evaluate the relatedness of samples [31,32] (http://unifrac.colorado.edu/static/welcome.html). Weighted and normalized PCoA and Jackknife Environment Clusters were used to evaluate similarity among samples, where each sample represented an environment. Hierarchical cluster analysis was performed using complete linkage and UniFrac distances as a measure of similarity [33]. Weighted and normalized Principal Coordinate Analysis (PCoA) was performed to evaluate similarity among the samples, where each sample represented an environment. R software (version 2.15.2) with gplots and the RColor-Brewer package was used for all other analytical processes.

Statistical analysis
Factor analysis was applied to track the important variables from an observed set of variables which are influencing the final response. Principal components were identified with over 70% of the cumulative variance. The component matrix of the chosen principal components was calculated, and variables whose factor loadings lower than 0.6 were deleted. A generalized linear model (GLM) was constructed using the final set of chosen variables. The following equation (a) represents the GLM: where t is the total number of occurrences,

x~1
, Normal group 0, Patients group , and g(m)~log m=t ð Þ in the pooled group of samples from asthmatic patients, COPD patients, and healthy controls.
The parameter beta was tested at a 5% significance level by using (a).

Study subjects
The 47 characteristics of the subjects are summarized in Table 1. For the current microbiota study, subjects were segregated into three groups: 18 patients with asthma (10 males, mean age 53.4617.1 years; 8 females, mean age 55.0613.0 years). All patients had a history of wheezing, cough, and dyspnea, with predicted post-bronchodilator FEV1/FVC 75.466.0, FEV1% 108.1625.6, and FVC% 105.4623.7. COPD subjects (15 males, mean age 68.967 years; 22 females, mean age 74.560.7 years) received a diagnosis of stage I (n = 6), stage II (n = 10), stage III (n = 1), or stage IV (n = 0) based on Global Initiative on Obstructive Lung Disease (GOLD) guidelines [25], with predicted post-bronchodilator FEV1/FVC 60.3613.3, FEV1% 78.2626.0, and FVC% 105.4623.7. Healthy control subjects (six males, mean age 50.2622.6 years; six females, mean age 59.3617.7 years) exhibited no evidence of underlying respiratory disease.

Sequence reads of the 16S rRNA gene
The diversity of microbiota in the oropharyngeal swab sample was characterized by sequencing the V1-V3 region of the 16S rRNA gene, which was amplified from genomic DNA samples of 18 asthmatic patients, 17 COPD patients, and 12 healthy controls. In total, 259,572 bacterial sequence reads were identified for three group samples with average read lengths of 273-476 bp (mean 6 standard deviation, 398.9641.9) ( Table 2). The 191,616 sequence reads obtained after trimming of low-quality sequences were then analyzed (mean 6 standard deviation, control: 321461136; asthma: 414063723; and COPD: 461965710).
The Chao1 estimator of total species richness sequences in healthy controls and asthmatic and COPD patients is described in Table 2. Ecological diversity was also estimated by using the Shannon-Wiener index, which takes into account relative abundance. This index is widely applied in ecological studies as a measurement of biodiversity, as it accounts for both the number of different taxa and their relative abundance. For this purpose, the 16S rRNA fragments were identified from the normalized pyrosequencing datasets and classified using the RDP classifier. Only 16S rRNA sequences of genus level with at least 80% confidence was used for Shannon-Wiener index. The resulting Shannon-Wiener index value demonstrated that the normal group was more diverse than the COPD and asthmatic groups because the sequence data obtained from the 454/Roche GS FLX Titanium sequencer provided a high estimate of biodiversity.
Species richness, diversity, and rarefaction calculations were used to analyze whether the diversity of the oropharyngeal microbial community was sufficiently covered by the sequence data. Species coverage, richness, evenness, and diversity were estimated for the combined data set for each of the bacterial sequences from the oropharyngeal samples. In addition, OTUs (5% difference) were calculated for the healthy control group and patients with asthma and COPD. Rarefaction analyses with RDP analysis curves at the OTUs reached saturation for the subjects with lower microbial diversity when the sequences of each sample type and the total number of OTUs were estimated. The diversity ranged from 29 to 410 OTUs with 95% similarity in each of the 47 samples (mean 6 standard deviation, control: 207699; asthma: 128673; and COPD: 1576100) ( Table 2). The calculation of rarefaction curves was used for all the samples, which showed that additional sampling would not have provided any additional OTUs ( Figure S1) with very few exceptions. The samples from the asthma, COPD, and healthy control groups contained a total of 190 different bacterial genera representing 12 different bacterial phyla (Table S1).

Phylum identification of bacterial communities
Seven bacterial phyla were detected, but most sequences were assigned to three phyla: Firmicutes, Proteobacteria, and Bacteroidetes. In 17 of 18 asthmatic patients, the majority of bacteria belonged to Firmicutes (59.8%), at a proportion higher than that in the healthy controls (48.4%). In addition, in 17 of 17 COPD patients, 61.6% of the bacterial phyla belonged to Firmicutes. Table 4. Statistical analysis of the healthy control, asthma, and chronic obstructive pulmonary disease (COPD) groups. Therefore, an increase in the prevalence of Firmicutes from 59.8% in the asthma sample to 61.6% in the COPD sample was observed. In the 12 healthy control samples, Proteobacteria and Bacteroidetes composed 29.02% and 16.20% of the bacteria community, respectively. In asthmatic and COPD patients, Proteobacteria accounted for 35.61% and 28.04% of the bacteria community, respectively. However, the percentage of Bacteroidetes was drastically decreased in the asthma (2.62%) and COPD (7.17%) groups ( Table 3).

Phylum
The relative abundance of the main phyla identified in the normal control group and in the asthmatic and COPD patients was analyzed using an R script with gplots and the RColorBrewer package ( Figure 1). The Firmicutes and Proteobacteria phyla were highly abundant in the asthma and COPD groups. The Bacteroidetes and Actinobacteria phyla had low abundance in the COPD group but were highly abundant in the healthy control group. These results indicate that the increase in microbial communities in the asthma and COPD groups was associated with an abundance of Proteobacteria and Firmicutes along with a reduction in the number of Bacteroidetes.

Genus identification of bacterial communities
To diversity of the microbial communities present in the oropharyngeal swab samples, was determined from BLASTn analysis, where the assembled sequences showed closest relationship with their neighbouring genera. Although the three data sets showed broad taxonomic similarity, there was notable variation in each sample at the genus level. In total, 190 different bacterial genera were identified in this study: 105, 131, and 133 different genera were present in the healthy control, asthma, and COPD groups, respectively (Table 3). Most sequences in the microbiome library from the healthy control group were Streptococcus spp. In the healthy control group, Neisseria spp. (Proteobacteria) and Prevotella spp. (Bacteroidetes) constituted 21.69% and 13.42% of the sequences, respectively, and thus comprised more than 40% of the bacterial community in the healthy control group compared to just over 10% in the asthmatic and COPD patients. When compared to the healthy controls, Lactobacillus spp. (29.67% and 34.50%) and Pseudomonas spp. (15.80% and 10.13%) were much more prevalent in patients with asthma or COPD, respectively. Lactobacillus spp. was present in 14 of 18 asthmatic patients and 15 of 17 patients with COPD, and Pseudomonas spp. was present in 10 of 18 asthmatic patients and eight of 17 patients with COPD. This was in contrast to the low frequencies of Lactobacillus spp. (0.07%) and Pseudomonas spp. (0.06%) in the healthy control population, i.e., six and four, respectively, of a total of 12 individuals. The proportion of Pseudomonas spp. decreased from 15.80% in the asthma sample to 10.13% in the COPD sample, and that of Lactobacillus increased from 29.67% in the asthma sample to 34.50% in the COPD sample. Among the phylum Bacteroidetes, populations of Prevotella spp. that were observed in the healthy control group were replaced with Pseudomonas spp. of Proteobacteria in the asthma and COPD groups (Table 3). It should be noted that this analysis was coverage-limited.

Identification by statistical analysis
We found that the three phyla and six genera in the normal controls and patients with asthma and COPD that were identified by sequencing and BLASTn analysis were significantly different in terms of frequency among samples, as indicated by the GLM statistical analysis (Table 4). In the statistical analysis that estimated the number of bacterial genera in each sample, if the value was less than 0, the genus was considered to occur at a greater frequency in the healthy sample. If the value was greater than 0, the genus was considered to occur at a greater frequency in the patient samples. Pseudomonas spp. (5.471) and Lactobacillus spp. (6.0006) were found with greater frequencies in the asthmatic and COPD patients. Streptococcus spp., Neisseria spp., Veillonella spp., and Prevotella spp. were found with greater frequencies in the healthy controls (p,0.0001). Jackknife clustering and PCoA were performed using Fast UniFrac [34].
The identification of microbiota in the healthy control, asthma, and COPD groups (Figure 2A), proved to be fairly robust with the Jackknife clustering method (75% bootstrap on all nodes but one). This is also supported by Fast UniFrac PCoA plot results ( Figure 2B). The dendrogram and PCoA plots depict well, the differences in overall microbial community structure in each group, highlighting the differences between three groups. It also apparent from these figures that the microbiotas in asthma samples were closer to outlying COPD samples. On the contrary in PCoA plots, asthma samples clustered with both sides of normal and COPD samples. There is a distant overlap between asthma and COPD subjects with many outliers including the controls, emphasizing absence of tight clustering between the disease states. A derangement of some of the samples was observed with asthma no 18, no 25, no 30, and COPD no 85 falling in the normal control cluster, while COPD no 88, no 93, no 94 clustered into the asthma group. Normal control sample no 3, no 5, no 15 and asthma sample no 17, no 21, no 22, no 27, no 83, no 84 were found to be a part of the COPD group.

Discussion
Respiratory diseases include the common cold, influenza, pharyngitis, sinusitis, bronchitis, pneumonia, asthma, and COPD, which account for more than 10% of visits to primary care physicians [35]. Among these, conditions like asthma and COPD are heterogeneous syndromes of chronic airway inflammation with undetermined pathogenesis and high mortality rates [36]. There are many guidelines for the diagnosis and management of asthma and COPD, but to date both diseases have often been misdiagnosed and undertreated [36,37]. The key to confront this problem is to understand and characterize the microbiome of these disease conditions in comparison to the healthy state. The human microbiome studies can identify the link between healthy and diseased states and have the capacity to reveal disease mechanisms.
So far the studies of the microbial community in asthma and COPD of human airways, have targeted the lower respiratory tracts using samples such as BAL and the sputum, showing controversial results about the role of the microbiota. There are still gaps when compared to previous studies in this aspect, for instance trying to understand the role of the microbial populations in the lower respiratory tract of oropharyngeal airway in conditions like asthma and COPD which is still unclear. So it becomes imperative to get a true representation of the microbiome in these disease states.
In present decade it has become possible to capture the true microbial diversity with modern high throughput culture-independent techniques which aid in comparing the microbiome between the healthy and diseased conditions. These studies can provide new information on the oropharyngeal microbiota in asthma and COPD patients. Previously, several studies investigated the links between airway bacterial concentrations and COPD, but with targeted molecular typing methods aimed at detecting specific bacterial species in the sputum [10,11]. These studies demonstrated that changes in the bacterial load in the sputum were unlikely to be an important mechanism in the exacerbation of COPD. These results were also supported by a study of the microbiome of healthy smokers and individuals with COPD using 454 pyrosequencing analysis, which showed that relatively low numbers of diverse types of bacteria existed in the lung [19].
In the current study, we gathered comprehensive information concentrating on the whole microbial diversity of the lower respiratory tract as opposed to previous studies with a nontargeted culture independent approach. We have three main conclusions from this study. First in terms of overall diversity, the oropharyngeal airway microbiota are quite diverse but reveal only few differences between asthma and COPD patients. Secondly, Pseudomonas spp. of the phylum Proteobacteria and Lactobacillus spp. of the phylum Firmicutes prove to be the most dominant populations in asthma and COPD patients but however were not detected in the healthy oropharynx. The final conclusion from this study is the dominance of groups like, Streptococcus, Veillonella, Prevotella, and Neisseria, belonging to the phylum Bacteroidetes, in the healthy oropharynx.
Briefly highlighting the methodology, sequences were generated using a high-throughput GS-FLX 454 pyrosequencer and aligned using tools like PRINSEQ, RDP, and Clustal_X. A GLM analysis of the 259,572 total sequence reads reveals a rich microbial environment in the oropharyngeal airway in asthma and COPD individuals. Distance-based statistical UniFrac tests were used for detecting the changes in abundance and rare lineages [33]. Minimal differences were observed between asthma and COPD individuals (Proteobacteria, Firmicutes; p,0.0001). As observed Pseudomonas (estimated asthma: 5.471; COPD: 4.890) and Lactobacillus spp. (estimated asthma: 6.006; COPD: 6.203) were found to be in abundance in asthma and COPD subjects when compared to healthy individuals. This is consistent with the findings of previous diversity studies in the bronchi of asthma and COPD patients, by 16S rRNA gene-based Sanger sequencing and by the 16S rRNA PhyloChip methods in endotracheal aspirates [38,39].
Hilty et al. [39] analyzed 5,054 16S rRNA bacterial clones in the bronchi and found that the members of Proteobacteria were in dominance in asthmatic children, whereas Bacteroidetes, particularly Prevotella, were more frequent in controls than in adult or pediatric asthmatic or COPD patients. Sze et al. [40] evaluated the lung tissue microbiomes of patients with severe COPD at the time of lung transplantation and noted increased bacterial diversity of the phylum Firmicutes, attributable to an increase in Lactobacillus species. Using 16S rRNA-based culture-independent profiling methods, Huang et al. [41] reported increased microbiota at exacerbation of COPD primarily belonging to phylum Proteobacteria, including some non-typical COPD pathogens. In another study a combination of high-density 16S ribosomal RNA microarray and parallel clone library-sequencing analysis in asthmatic patients revealed a rich microbial environment in the oropharyngeal airway [42].
In the oropharyngeal airway of asthma and COPD patients, Pseudomonas, Stenotrophomonas, and Lactobacillus were found to be dominant. Whereas the main phyla reported in the bronchoalvelolar lavage fluid of COPD patients are Actinobacteria, Firmicutes, and Proteobacteria, with dominance of Pseudomonas, Haemophilus, or Stenotrophomonas. in the respective phyla [19]. However, Haemophilus was not found in our study of the oropharynx. Pseudomonas and Stenotrophomonas maltophilia are the most common pathogenic bacteria associated with cystic fibrosis [43][44][45]. In COPD, Pseudomonas was also reported to a risk factor for severe acute exacerbation [19] and Lactobacillus is known to be involved in severe COPD and fatal lung abscess [19,40]. Therefore, further study is necessary to understand the pathogenesis of diseases at microanatomic sites of the respiratory tract.
In our study the healthy controls showed higher frequencies of organisms from Streptococcus, Veillonella, Prevotella, and Neisseria of phylum Firmicutes when compared to asthmatic and COPD patients. Normally, Streptococcus species are considered to be commensal bacteria of the oral cavity; however, S. pneumoniae, S. mitis, and S. oralis are members of the viridians streptococci and are associated with pneumonia and endocarditis and are considered as opportunistic pathogens [46][47][48][49]. Veillonella is a type of subgingival microflora found at periodontal disease sites and is known to be a periodontopathogenic bacteria [50]. Recently, Prevotella copri was reported to be involved in inflammation [51].
Previously, several studies reported less abundance of Neisseria in smokers [16,52,53]. This study also showed that Neisseria was less dominant in asthmatic and COPD patients, implying that Neisseria is an important indicator of the health of the oral and upper respiratory tract. Therefore, the potential relationship between these oropharynx microbiota and chronic respiratory diseases should be evaluated.
Previously exquisite study of the bronchial microbiome obtained from lower respiratory samples avoiding upper airway contami-nation in the healthy adults revealed overlapping of certain groups of bacteria similar to upper airways in six subjects [16]. Therefore, our results suggest that the distribution and profile of microbiota are consistent between the oropharynx and bronchi in asthma and COPD patients and that the pathogenic role of oropharyngeal commensal bacteria changes with n the progression of chronic respiratory diseases.