Genomic diversity and antimicrobial resistance of Prevotella spp. isolated from chronic lung disease airways
============================================================================================================

* Kasey A. Webb
* Olusola Olagoke
* Timothy Baird
* Jane Neill
* Amy Pham
* Timothy J. Wells
* Kay A. Ramsay
* Scott C. Bell
* Derek S. Sarovich
* Erin P. Price

## Abstract

Cystic fibrosis (CF) and chronic obstructive pulmonary disease (COPD) are characterised by increasingly frequent acute pulmonary exacerbations that reduce life quality and length. Human airways are home to a rich polymicrobial environment, including members of the obligately anaerobic genus, *Prevotella*. Despite their commonness, surprisingly little is known about the prevalence, role, genomic diversity, and antimicrobial resistance (AMR) potential of *Prevotella* species/strains in healthy and diseased airways. Here, we used comparative genomics to develop a real-time PCR assay to permit rapid *Prevotella* spp. quantification from cultures and clinical specimens. Assay specificity was validated across a panel of *Prevotella* and non-*Prevotella* species, followed by PCR screening of CF and COPD respiratory-derived cultures. Next, 35 PCR-positive isolates were subjected to whole-genome sequencing. Of eight identified species, *P. histicola, P. melaninogenica, P. nanceiensis, P. salivae* and *P. denticola* overlapped between participant cohorts. Phylogenomic analysis revealed considerable interhost but limited intrahost diversity, suggesting patient-specific lineages in the lower airways, probably from oral cavity aspirations. Correlation of phenotypic AMR profiles with AMR gene presence identified excellent correlation between *tetQ* presence and decreased doxycycline susceptibility, and *ermF* presence and decreased azithromycin susceptibility and clindamycin resistance. AMR rates were higher in the CF isolates, reflecting greater antibiotic use in this cohort. All tested *Prevotella* isolates were tobramycin-resistant, providing a potential selection method to improve *Prevotella* culture retrieval rates. Our addition of 35 airway-derived *Prevotella* genomes to public databases will enhance ongoing efforts to unravel the role of this diverse and enigmatic genus in both chronic respiratory diseases and healthy lungs.

**Data summary** Thirty-five *Prevotella spp*. genomes generated in this study are available in the Sequence Read Archive (SRA) and GenBank databases under BioProject accession PRJNA742126.

Key Words
*   *Prevotella*
*   genomics
*   cystic fibrosis
*   chronic obstructive pulmonary disease
*   respiratory infections
*   polymicrobial infections
*   real-time PCR
*   diagnostics
*   sputum

## Introduction

Chronic lung diseases are a leading cause of global healthcare burden and early mortality1,2. Although cystic fibrosis (CF) and chronic obstructive pulmonary disease (COPD) have fundamentally different aetiologies, both are characterised by impaired airway mucociliary clearance, which provides a favourable environment for microbial colonisation and persistence, and a challenging environment for pathogen eradication. Infections are the primary contributor in CF and COPD respiratory exacerbations (i.e. acute worsening of symptoms), during which inflammation and irreversible airway damage can occur3. In CF, 80-95% of deaths are associated with acute exacerbations and chronic bacterial infections, which eventually lead to respiratory failure4. In COPD, bacterial, viral, and fungal respiratory infections trigger ∼70-80% exacerbation events5. Despite their crucial role in reducing quality and length of life, there remain major gaps in understanding the collective microbial drivers associated with CF and COPD pathophysiology.

Routine culture-based diagnosis of CF and COPD exacerbations often overlooks more fastidious or uncultivable members within the airway microbiome, leaving a gap in fully understanding the complex composition and activity of CF and COPD microbiomes6,7. Culture-independent techniques like metataxonomics (e.g. 16S rRNA amplicon sequencing), metagenomics, and metatranscriptomics are gaining traction as adjuncts or potentially even alternatives to culture due to their ability to identify a much broader spectrum of microbes, including uncultured taxa, resulting in better diagnosis, more informed treatment strategies, and improved clinical outcomes8-10. In addition to identifying well-characterised microbes associated with respiratory exacerbations, molecular methods have unveiled previously unidentified microbial communities in chronic lung disease airways, including a diverse anaerobic microbiome11.

The *Prevotella* genus, which comprises ≥51 species12,13, is a common and often abundant obligate anaerobe in the human body, including in CF and COPD airways11,14-19. Despite its commonness, the role of *Prevotella* in disease progression is incompletely understood in CF11, and even less so in COPD. Several studies investigating *Prevotella* spp. in airway disease, particularly CF, have yielded conflicting evidence. Some studies suggest a likely pathogenic role18,20, whereas others have shown no role21, or even a benefit when *Prevotella* spp. are abundant22,23. A major contributor to these equivocal findings is that culture-dependent and metaxonomic methods largely fail to identify microbes to the species and strain level7,13,24,25. As such, only a handful of *Prevotella* studies have described potential species-or strain-specific mechanisms in airway disease pathogenesis.

The number of *Prevotella* genomes has been steadily increasing in recent years. However, in comparison to more intensely studied CF and COPD airway organisms, there is a relative paucity of respiratory *Prevotella* genomic data in public databases. Most *Prevotella* species contain only a single representative (and many were generated from a single metagenomic study)26, leading to issues with correct species assignments from metagenome-assembled genomes. Further, most genomic studies have concentrated on known pathogens in oral (e.g. *P. intermedia* and *P. nigrescens*), gut (e.g. *P. copri*), and urogenital (e.g. *P. bivia*) microbiomes27-29. There is a lack of knowledge about the *Prevotella* spp. commonly found in airways, particularly in COPD. Insufficient resolution to the species and strain level has limited our understanding of *Prevotella* spp. diversity, function, evolution, and antimicrobial resistance (AMR) potential within CF and COPD airways during exacerbations and stability13,24.

Here, we first used a comparative genomics approach to design and optimise a novel real-time PCR assay to permit the rapid, cost-effective, and specific detection of *Prevotella* spp. from both culture and clinical specimens. Next, 43 CF-COPD sputa and bronchial washings were cultured to isolate individual *Prevotella* strains for whole-genome sequencing (WGS) characterisation, followed by genomic analyses to determine species identity, strain relatedness, and within-host evolution. Cultured isolates were subject to AMR testing across 11 clinically-relevant antibiotics, and correlated with AMR determinants identified by WGS. Finally, *Prevotella* morphological characteristics were assessed for their diagnostic value.

## Methods

### Ethics statement

Collection and analysis of clinical samples was approved by The Prince Charles Hospital (TPCH) Human Research Ethics Committee (HREC), project IDs HREC/13/QPCH/127 (CF samples) and HREC/2019/QPCH/48013 (COPD samples). Site-specific approvals were obtained for CF recruitment at TPCH, Brisbane, Australia, and COPD recruitment at Sunshine Coast Hospital and Health Service, Sunshine Coast, Australia. All participants provided written consent.

### Sample collection and storage

A total of 36 CF sputa were obtained from 21 expectorating adults with CF presenting at TPCH Adult CF Centre between 2017 and 2019. Seven participants provided multiple sequential sputa (**Table S1**). Sputa were stored non-anaerobically at -80°C in Luria-Bertani broth containing 20% glycerol, and stored for 2-24 months prior to microbiological assessment. A total of 32 COPD sputa (*n*=21) or bronchial washings (BWs; *n*=11) were obtained from 21 current/former tobacco smokers during community or in-hospital consult at Sunshine Coast Hospital and Health Service in 2020 (**Table S1**). Five participants provided multiple sequential sputa (**Table S1**). BWs were taken from participants with suspected lung cancer, with most collected as paired specimens from both the tumour site (biopsied after BW collection) and approximately the same position in the non-tumour contralateral (contra) lung. BWs were pelleted at 4°C and ∼9,300 x g for 15 min, with supernatant removed and pellet resuspended in 800 µL phosphate-buffered saline. COPD samples were stored within 12-36 hours of collection using the same parameters as described for the CF specimens, and microbiologically processed within ∼2 weeks of collection. Clinical and antibiotic treatment regimens were recorded (**Table S1**) to permit correlation with *Prevotella* AMR profiles.

### Microbiological preparation and storage

Clinical samples were directly subcultured onto anaerobic blood medium containing nalidixic acid and vancomycin (Thermo Fisher Scientific, Seventeen Mile Rocks, QLD, Australia; cat. no. PP2065). Plates were incubated anaerobically at 37°C in an AnaeroPack 2.5L jar with an AnaeroGen 2.5L gas pack (Thermo Fisher) for 2-5 days. Isolates were subcultured based on different morphotypes, with a second subculture for purity onto non-selective anaerobic blood medium (Thermo Fisher; cat. no. PP2039). *Prevotella* real-time PCR-positive isolates were again sub-cultured onto non-selective anaerobic blood medium for purity and stored non-anaerobically at -80°C in 1mL Luria-Bertani broth containing 20% glycerol. Additionally, WGS-confirmed, representative *Prevotella* isolates derived from CF samples were cultured from glycerol stocks onto both polyvitaminic-supplemented chocolate agar (PVCA; Edwards Group, Murrarie, QLD, Australia; cat. no. 1065) and non-selective anaerobic blood medium to assess phenotypic differences over a 10-day period.

### DNA extractions

DNA of presumptive *Prevotella* spp. isolates was initially extracted using 5% chelex-10030 (Bio-Rad, Gladesville, NSW, Australia), heat-treated at 95°C for 20 min, and diluted to 1:50 with molecular-grade H2O prior to PCR. DNA for WGS was extracted using the Gram-positive bacteria protocol of the DNeasy Blood and Tissue Kit (Qiagen, Chadstone, VIC, Australia). Qiagen extractions were quality-checked and quantified using the NanoDrop 2000 spectrophotometer (Thermo Fisher).

Clinical specimens were extracted using 3 parts TRI Reagent LS (Sigma-Aldrich, North Ryde, NSW, Australia) to 1 part specimen, 30µL 0.1M 2-mercaptoethanol, and ∼ 100µL equal mixture of 0.1 and 0.5mm zirconia beads (Daintree Scientific, St Helens, TAS, Australia). Homogenisation comprised eight rounds of 30 sec bead beating (each round followed by 30-60 sec on wet ice) using the ‘medium’ setting of the Minilys tissue homogeniser (Bertin Instruments, Montigny-le-Bretonneux, France). After removal of the RNA aqueous phase, back extraction was performed on the interphase and proteinaceous phase31 to isolate total DNA.

### *Prevotella* spp. real-time PCR assay design

Due to high levels of genetic variability among *Prevotella* spp., the highly conserved 16S-23S rRNA genes were targeted. A custom BLAST32 2.10.2+ database was created using 60 *Prevotella* spp. genomes encompassing 26 assigned and nine unassigned species (**Table S2**). Discontiguous MegaBLAST was carried out against the 4.8kb *P. intermedia* ATCC 25611 23S-16S rRNA locus (GenBank accession [CP019300.1](http://medrxiv.org/lookup/external-ref?link_type=GEN&access_num=CP019300.1&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom), coordinates 63001-67800) using the custom BLAST database. BLAST .xml outputs were aligned and visualised using Geneious Prime vR10. A highly conserved area within the *Prevotella* spp. 16S rRNA gene was assessed further to identify potential PCR oligos, with conserved regions possessing several single nucleotide polymorphisms (SNPs) in other microbes exploited during probe and primer design. Primers and probes were initially assessed using Primer Express v3.0 (Applied Biosystems) to test melting temperature (*T*m) suitability in Black Hole Quencher format (primer *T*m=60°C; probe *T*m=70°C). All candidate primers and probes were assessed for self- and hetero-dimer formation using Beacon Designer ([http://www.premierbiosoft.com/qpcr/](http://www.premierbiosoft.com/qpcr/)), with oligo sequences requiring a ΔG value of >-4 kcal/mol. Candidate oligos were subsequently subjected to NCBI Microbial Nucleotide BLAST and Nucleotide BLAST searches to confirm specificity for *Prevotella* spp.. As the 16S region is highly conserved in all bacteria, the forward primer incorporated a penultimate mismatch to prevent non-specific binding to closely related non-*Prevotella* taxa. Additionally, a BHQplus probe (CAL Fluor Gold 540-pdCApdUApdUpdCGApdCGGGGpgUGG-BHQ1; Biosearch Technologies, Petaluma, California, USA) was included to enhance specificity and to permit assay use on clinical specimens. Forward (5’-TTCTCCAGtCCAGGATGTGtC-3’) and reverse (5’-CGGGGATAACAGGCTGATCC-3’) primers were manufactured by Macrogen Inc. (Geumcheon-gu, Seoul, South Korea), with lowercase nucleotides representing intentional mismatches.

Real-time PCRs were carried out on a Bio-Rad CFX96 thermocycler using 0.35µM primers and probe, 1X Sso Advanced Universal mastermix (Bio-Rad), molecular-grade H2O, and 1µL template, to a 5µL final volume. Thermocycling parameters were: 95°C for 2 min, followed by 45 cycles of 95°C for 5 sec and 60°C for 20 sec. DNA integrity for all clinical samples was assessed via amplification towards β-globin33 (human) and universal 16S rRNA34 (bacteria); the universal 16S rRNA assay was also used to verify bacterial culture DNA integrity, and a panfungal assay (O. Olagoke and E. Price, unpublished) was used to verify *Candida* spp. DNA integrity.

### Real-time PCR assay optimisation and specificity testing

The final optimised real-time PCR assay was assessed for accuracy by running against a diversity panel of six *Prevotella* spp. (one each of *P. copri, P. denticola, P. histicola, P. melaninogenica, P. nigrescens*, and *P. salivae*), 1 human DNA sample, and 40 non-*Prevotella* spp., with a focus on common respiratory tract microbes (*n*=1 unless otherwise stated): *Achromobacter xylosoxidans, Bordetella bronchiseptica, Burkholderia cepacia, Burkholderia territorii, Burkholderia thailandensis, Candida albicans, Candida dubliniensis, Cupriavidus metallidurans, Enterobacter aerogenes, Enterobacter cloacae* (*n*=3), *Enterococcus raffinosus, Klebsiella oxytoca, Klebsiella pneumoniae* (*n*=2), *Lacticaseibacillus paracasei* (*n*=2), *Pseudomonas aeruginosa* (*n*=12), *Staphylococcus aureus* (*n*=3), *Staphylococcus haemolyticus, Stenotrophomonas maltophilia, Veillonella atypica* (*n=*3), *Veillonella parvula*, and *Veilonella dispar*. All DNA was normalised spectrophotometrically to 2ng/µL prior to PCR.

Limits of detection (LoD) and quantification (LoQ) for the optimised assay were determined using *P. melaninogenica* LMG 28911 DNA, using serial dilutions of 4ng/µL to 4 × 10−6 ng/µL across 8 replicates per dilution. The upper and lower LoD/LoQ limits were determined as described elsewhere35. Genome equivalence was determined based on an average *Prevotella* genome size of 3.29Mbp and a single copy of the 16S rRNA gene.

### 16S rRNA gene diversity profiling (metataxonomics)

Metataxonomics was performed by the Australian Genome Research Facility (AGRF; Melbourne, VIC, Australia) for total CF DNA and ID Genomics (Seattle, WA, USA) for total COPD DNA using standard protocols for V3-V4 (i.e. 341F and 806R) primers.

### Antimicrobial resistance (AMR) profiling

Twenty-eight glycerol-stocked *Prevotella* cultures were subcultured onto non-selective anaerobic blood medium, or PVCA where no growth was obtained on the anaerobic medium (*n*=2), and incubated at 37°C under anaerobic conditions for 2-3 days. Cultures were then resuspended into phosphate-buffered saline to OD600=0.75 (Biochrom, Cambridge, UK). Six supplemented chocolate plates per isolate were lawn-inoculated with the OD-normalised cells as previously described36, and antibiotic discs (Edwards Group) towards 11 antibiotics (amoxicillin-clavulanate [AMC] 30µg; azithromycin [AZM] 15ug; ceftazidime [CAZ] 30µg; chloramphenicol [CHL] 30µg; clindamycin [CLI] 2µg; doxycycline [DOX] 30µg; imipenem [IPM] 30µg; meropenem [MEM] 10µg; metronidazole [MZ] 5µg; piperacillin-tazobactam [TZP] P:100µg, TZ:10µg; tobramycin [TOB] 10µg) were placed onto inoculated plates. Plates were incubated at 37°C under anaerobic conditions for 48h prior to measurement. As no disc diffusion reference ranges exist for *Prevotella* spp., growth up to the disc (i.e. 6mm diameter) was denoted as AMR, which underestimates the true rate of AMR but provides a clear cut-off for unambiguous interpretation. Diameters ≤30mm for any antibiotic were denoted as reduced susceptibility.

### Whole-Genome Sequencing (WGS)

PCR-confirmed *Prevotella* isolates were sequenced using the Illumina HiSeq 2500 platform with 150bp seqWell plexWell paired-end read chemistry at AGRF.

### Genome assemblies

*Prevotella* genomes were *de novo* assembled using default Microbial Genome Assembly Pipeline v1.1 parameters ([https://github.com/dsarov/MGAP\---|Microbial-Genome-Assembler-Pipeline](https://github.com/dsarov/MGAP\---|Microbial-Genome-Assembler-Pipeline)) or SPAdes37.

### Speciation and mixture analysis

*Prevotella* spp. assemblies were initially assessed for species identification and mixtures using the Ribosomal MLST38 tool. Isolates containing mixtures of *Prevotella* and non-*Prevotella* species (*n*=8), or multiple *Prevotella* species (*n*=1) (**Table S3**), were excluded from further investigation.

### Phylogenetic analysis

A custom BLAST database of *Prevotella* genomes (*n*=234) comprising taxa of known (*n*=56) and unassigned (*n*=40) species, and including non-mixed isolate genomes sequenced in this study (*n*=35 genomes), was created using BLAST 2.10.2+. BLAST analysis was performed against the full-length 1,534bp *P. melaninogenica* 16S rRNA gene (RefSeq ID: NC_014371.1; bases 962,784-964,317). Hits were filtered to remove redundancies (final *n*=199), and reads aligned using ClustalX2 v2.139. A dendogram was constructed using the neighbour-joining method and the Tamura-Nei model and viewed in FigTree v1.4.0 ([http://tree.bio.ed.ac.uk/software/figtree/](http://tree.bio.ed.ac.uk/software/figtree/)) using midpoint rooting.

### Phylogenomics

Variant identification of our CF-COPD genomes (‘SCHI’ prefix) was performed using SPANDx v4.040. Reference genomes for *P. denticola, P. histicola, P. melaninogenica, P. nanceiensis*, and *P. salivae* were downloaded from GenBank (IDs GCA\_000193395.1, GCA\_000234055.1, GCA\_000144405.1, GCA\_000163055.2, and GCA_000477535.1, respectively) and used for read mapping. Synthetic reads were generated from GenBank genome assemblies for *P. denticola* (GCA\_000191765.2, GCA\_000759205.1, GCA\_000421205.1, GCA\_003515045.1, and GCA_900454835.1), *P. histicola* (GCA_000613925.1), *P. melaninogenica* (GCA\_000163035.1, GCA\_000759305.1, GCA\_002208725.2, GCA\_003609775.1, and GCF_000517825.1), *P. nanceiensis* (GCA\_000517825.1 and GCA_000613985.1), and *P. salivae* (GCA_000185845.1 and GCA_902374235.1) using ART v2.3.741.

Illumina reads for *P. histicola* T05-04T (SRR1518629) and *P. nanceiensis* ADL-403 (SRR755359) were downloaded from the Sequence Read Archive database and quality filtered using Trimmomatic v0.3942 parameters described elsewhere43. The resulting biallelic, orthologous SNP matrix generated by SPANDx was phylogenetically reconstructed using the maximum parsimony heuristic function of PAUP* v4.0a (build 168). Phylogenies were visualised in FigTree; all trees were midpoint-rooted.

### AMR determinant analysis

The Resistance Gene Identifier v5.1.1, which employs the Comprehensive Antibiotic Resistance Database (CARD) v3.1.044, was used to identify putative AMR determinants.

## Results and Discussion

Molecular typing approaches have been instrumental in unveiling previously untapped microbial diversity in both healthy and diseased lungs. Using these techniques, *Prevotella* spp. have consistently been identified as one of the most prevalent taxa in CF and COPD airways10,14,19. However, most studies to date have employed 16S metataxonomics (usually targeting the V3-V4 gene region), which only characterises *Prevotella* to the genus level, despite species and strain resolution being fundamental to unravelling the role of these anaerobes in chronic lung disease. To address this knowledge gap, our study used a combination of culture-dependent and -independent methods to examine the genomic and phenotypic diversity, and prevalence, of *Prevotella* spp. present in CF and COPD airways.

We first used a comparative genomics-guided approach to develop an accurate real-time PCR assay to enable the rapid detection of *Prevotella* spp., to the exclusion of all other organisms, both from cultures and clinical specimens. Our PCR assay was reasonably sensitive, with LoD and LoQ both being 400fg/µL, or 111 genome equivalents (**Figure 1**). When tested against a diverse panel of 47 microbial and human DNA, the PCR demonstrated 100% specificity, amplifying all six tested *Prevotella* species, but none of the non-*Prevotella* spp. DNA and 10 no-template controls. On clinical specimens, the PCR identified *Prevotella* in 81% CF and 71% COPD participants, and in 81% CF and 66% COPD specimens (**Table 1**). This difference in detection rates between diseases was expected given the generally lower total bacterial load in COPD airways1. Metataxonomics had higher *Prevotella* sensitivity than PCR, identifying *Prevotella* in six PCR-negative samples, although there was one false-negative (SCHI0049.C.3, which was PCR-and culture-positive); there was otherwise near-perfect correlation between PCR and metataxonomics (**Table S3**). In the six PCR-negative samples, *Prevotella* relative abundance was very low (0.004-0.3%; median=0.05%), indicating that DNA levels were below the PCR LoD threshold. Thus, culture-independent methods are highly suitable for accurate *Prevotella* spp. detection from polymicrobial specimens, particularly real-time PCR, which is often performed upon clinical sample receipt, thereby permitting same-day detection of *Prevotella* in respiratory samples.

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/31/2021.08.30.21262864/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2021/08/31/2021.08.30.21262864/F1)

Figure 1. 
Limits of detection (LoD) and quantification (LoQ) for the *Prevotella* real-time PCR assay. The LoQ and LoD were both 400fg/µL according to previously defined parameters35, or 111 genome equivalents.

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/31/2021.08.30.21262864/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2021/08/31/2021.08.30.21262864/F2)

Figure 2. 
Full-length 16S ribosomal RNA gene phylogenetic analysis of the 35 *Prevotella* strains retrieved in this study (identified by ‘SCHI’ prefix) against 199 non-redundant public *Prevotella* genomes comprising 56 taxa of known species and 40 taxa of unassigned species. Species clades containing isolates identified from our cystic fibrosis (CF) and chronic obstructive pulmonary disease (COPD) specimens are coloured.

View this table:
[Table 1.](http://medrxiv.org/content/early/2021/08/31/2021.08.30.21262864/T1)

Table 1. 
*Prevotella* spp. prevalence in cystic fibrosis (CF) and chronic obstructive pulmonary disease (COPD) according to microbiological culture, real-time PCR, and 16S rRNA gene metataxonomics.

Culture of sputa and BWs onto three different agars (anaerobic blood medium containing nalidixic acid and vancomycin, non-selective anaerobic blood medium, and PVCA) yielded a much lower *Prevotella* positivity rate compared with PCR and metataxonomics. Only 44% CF and 28% COPD clinical specimens, respectively, were *Prevotella*-positive via culture, compared with 81% and 66% using PCR, and 84% and 91% using metataxonomics (**Table 1**). The reduced retrieval rate in COPD vs. CF can be explained by lower relative *Prevotella* abundance according to metataxonomics (0.3-74.3% [median=9.6%) for CF vs. 0.002%-36.4% [median=0.6%] for COPD). The lower prevalence of *Prevotella* using culture is likely due to multiple factors, such as an inability of culture-independent methods to distinguish viable from non-viable cells, nutritional fastidiousness, cells being in a viable but non-culturable state, and outcompetition from faster-growing taxa9,19,22. A prior study found that four different culture media were required to recover 42 *Prevotella* isolates from 26 samples, and only 20 isolates were recoverable on more than one medium19. Moreover, all four media were required to recover *P. melaninogenica* from all sputum samples19, suggesting that growth requirements may be species- and strain-specific. In contrast to this prior work, we used aerobic collection and storage conditions, which likely precluded recovery of *Prevotella* taxa with poor/no aerotolerance45. Consistent with their fastidiousness, one *P. histicola* and five *P. melaninogenica* initially isolated in our study were unable to be resuscitated from glycerol stocks upon repeated attempts, despite previously being culturable. Taken together, it is probable that several *Prevotella* species or strains were not recovered by culture in our study due to multiple factors, thereby underestimating species and strain diversity within our cohorts.

To gain greater insight into *Prevotella* residing in CF and COPD airways, we performed WGS and phylogenetic analysis on cultured isolates to determine species identity and strain relatedness. From 36 CF sputa, 27 PCR-positive *Prevotella* colonies were retrieved, of which 18 were *Prevotella* spp.; the remaining nine contained a mixture of *Prevotella* and other taxa, or a mixture of multiple *Prevotella* spp., due to insufficient subculturing for purity prior to WGS (**Table S3**). To avoid known bioinformatic analysis issues associated with variant calling from mixtures46, these mixed data were excluded from phylogenetic analysis. Among the 18 isolates, WGS identified seven *Prevotella* spp.: *P. histicola* (*n*=5), *P. melaninogenica* (*n*=5), *P. denticola* (*n*=3), *P. nanceiensis* (*n*=2), *P. jejuni* (*n*=1), *P. nigrescens* (*n*=1), and *P. salivae* (*n*=1). Of the 17 COPD-derived isolates, none were mixed due to a greater number of passages carried out to obtain purity. WGS identified six *Prevotella* spp.: *P. melaninogenica* (*n*=8), *P. histicola* (*n*=3), *P. salivae* (*n*=3), *P. buccae* (*n*=1), *P. denticola* (*n*=1), and *P. nanceiensis* (*n*=1) (**Table S3**). Despite being fundamentally different diseases, there was little difference in *Prevotella* species found in CF and COPD, with only *P. jejuni* and *P. nigrescens* being unique to CF, and *P. buccae* being unique to COPD. One explanation for this overlap is that the oral cavity and oropharynx are common habitats for all retrieved species, except *P. jejuni*47; it is thus possible that some strains we isolated were upper airway residents. Whilst oral cavity contamination is unavoidable during respiratory secretion collection, especially for sputum, microbes can also enter lower airways via small-volume oral cavity aspirations; indeed, the risk of microbial colonisation and persistence is greatly increased in mucociliary-impaired CF and COPD airways24,48-52. In support of lower airway origin, we identified *Prevotella* spp. in both sputa and BWs. Unlike sputa, BWs are collected via bronchoscopy and thus far less prone to oral/oropharynx contamination, providing greater confidence that many of our retrieved *Prevotella* isolates originated from the lower airways.

The origin of *P. jejuni* in one CF participant is somewhat perplexing. This species was first documented in the small intestine of a child with coeliac disease in 201253, and little is currently known about the distribution of this organism in healthy and diseased niches. One possibility is that *P. jejuni* entered the CF airways via inhalation of gastroesophageal aspirations. Alternatively, microbial inter-compartment cross-talk between the lungs and gut (the ‘Gut-Lung axis’) via the mesenteric lymphatic system may explain its presence in the lungs54. As a relatively new taxon, little is known about possible mechanisms of *P. jejuni* pathogenesis, with more work needed to determine the prevalence and virulence potential of this bacterium in both healthy and diseased airways.

The most frequently isolated species in our study – *P. melaninogenica* (37%), *P. histicola* (23%), *P. salivae* (11%), *P. denticola* (11%), and *P. nanceiensis* (9%) – are common in CF14,15,19 and COPD25 airways. However, their precise role in airway disease is poorly understood due to conflicting reports about their beneficial vs. pathogenic potential. For instance, *P. melaninogenica, P. salivae*, and *P. nanceiensis* may inhibit non-typeable *Haemophilus influenzae*-induced pro-inflammatory cytokine production by lung dendritic cells52, increased *P. melaninogenica* prevalence in CF airways has been associated with improved lung function15, and *P. histicola* can mute pro-inflammatory pathways triggered by *P. aeruginosa* in CF airways55. In contrast, *P. melaninogenica* produced the greatest variety of short-chain fatty acids (associated with inflammation) amongst five tested anaerobes18, and *P. denticola* shares similarities in its lipid A moieties to *P. intermedia* (a known pathogenic species), pointing towards pro-inflammatory traits56. Finally, certain strains of the periodontal disease-associated species, *P. nigrescens*, a close relative of *P. intermedia*, are thought to encode an enzymatic cell hydrolysis mechanism that destroys tissues. Given the reported conflicting roles of *Prevotella* spp. in chronic lung disease, these collective results suggest that virulence mechanisms vary among *Prevotella* strains, which requires further exploration57,58. Our genomic data will assist these future efforts.

In our experience, metagenomic analysis tools can struggle with *Prevotella* species-level identification from respiratory specimens due to a lack of accurate representative species in public databases. Exemplifying this problem, few publicly available genomic data are available for *P. nanceiensis* – only three were available at the time this study was undertaken, one of which is misclassified as *P. melaninogenica*, and another as *Prevotella* sp. (**Figure 3**). Thus, our addition of three *P. nanceiensis* genomes doubles the available *P. nanceiensis* data. Similarly, there are currently only a handful of nonredundant *P. jejuni* genomes in public databases, and none are from lower airway specimens; our addition *P. jejuni* SCHI0011.S.13 genome is therefore valuable for future respiratory and gastroenterological studies. Taken together, our *P. nanceiensis* and *P. jejuni* genomes will help public database curation efforts, and in turn, improve species identification using metagenomics.

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/31/2021.08.30.21262864/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2021/08/31/2021.08.30.21262864/F3)

Figure 3. Maximum parsimony within-species phylogenomic reconstruction of *P. denticola, P. histicola, P. melaninogenica, P. nanceiensis* and *P. salivae*.
Strains identified from COPD and CF specimens from this study are labelled with the ‘SCHI’ prefix; the remaining strains are public genomes. The number of SNPs separating each strain is labelled on branches; where two or more strains were retrieved from a single participant, coloured branch labels represent the number of SNPs separating these strains. CI, consistency index. Incorrectly speciated public genomes are denoted by an asterisk.

To determine strain relatedness within and between participant airways, we undertook within-species analyses of our *P. denticola, P. histicola, P. melaninogenica, P. nanceiensis* and *P. salivae* genomes (**Figure 3**). Whole-genome SNP analysis revealed high levels of strain and interhost diversity in all examined *Prevotella* species. For example, genetic variation was as high as ∼72k SNPs among *P. histicola* strains, an impressive feat for an organism encoding a 3Mbp genome, and *P. nanceiensis* strains SCHI0010.S.3 and SCHI0011.S.12 were more closely related to strains isolated from France and the US, respectively, than to other strains retrieved from other Australian participants. In contrast, intrahost diversity was relatively limited (0-297 SNPs), with two instances of genetically indistinguishable strains: *P. histicola* SCHI0034.S.5 and SCHI0034.S.6, and *P. melaninogenica* SCHI0034.S.7 and SCHI0034.S.9 (**Figure 3**); the latter pair were retrieved from a BW and sputum collected on the same day (**Table S1; Table S3**). The only exception was *P. melaninogenica* from SCHI0027, where two genetically distinct lineages separated by >20k SNPs were retrieved from this participant. There was no evidence of genetically-related strains between any of the Australian participants for any *Prevotella* spp. (**Figure 3**). These findings indicate that interhost diversity is very high, whereas within-host *Prevotella* diversity is limited, possibly due to frequent microaspiration of *Prevotella* from upper to lower airways (or via the gut-lung axis) rather than chronic infection; however, more in-depth sampling is required to investigate *Prevotella* persistence.

Antibiotics are routinely used in the maintenance and treatment of acute respiratory infections in CF and COPD15,59. Given its uncertain pathogenic role, *Prevotella* spp. are not targeted by these treatment regimens60. Despite this, AMR rates in *Prevotella* are trending upwards59, possibly due to bystander selection. To examine this phenomenon, we determined AMR profiles for all culturable *Prevotella* strains (*n*=27) using disc diffusion against 11 antibiotics commonly used to treat lung infections15,59, or which have known activity against anaerobic taxa47. AMR (AZM, CAZ, CLI, MZ, TOB) or reduced susceptibility (AMC, AZM, CAZ, DOX, MEM) was observed in at least one strain, whereas all were susceptible towards CHL, TZP, and IPM, an observation supported by previous studies47,59,61. Reflecting much greater lifetime antimicrobial use in the CF cohort62, the 18 tested CF strains had much higher AMR incidence (57% vs. 8%) and AMR gene presence (83% vs. 35%) compared with COPD strains, and this rate increased to 100% (69% for COPD) when including strains with reduced susceptibility towards one or more antibiotics (**Table 2**). In addition, almost all CF participants were undergoing antimicrobial treatment with at least two different antibiotic classes at the time of sample collection, vs. almost no antibiotic use in the COPD participants (**Table S1**). One strain, CF isolate *P. histicola* SCHI0009.S.6, could be classed as multidrug-resistant due to exhibiting AMR towards three antibiotic classes63 – β-lactams, lincosamides, and imidazoles – along with reduced susceptibility towards the tetracycline-class antibiotic, DOX (**Table 2**). Our findings confirm bystander AMR selection in off-target *Prevotella* as a direct result of heavy antimicrobial use in CF (and to a much lesser extent, COPD), which may have ramifications for *Prevotella* persistence as a potential driver of disease pathogenesis.

View this table:
[Table 2.](http://medrxiv.org/content/early/2021/08/31/2021.08.30.21262864/T2)

Table 2. 
Disc diffusion and antimicrobial resistance gene results for genome-sequenced *Prevotella* isolates from cystic fibrosis (CF) or chronic obstructive pulmonary disease (COPD) respiratory secretions, and antibiotic treatment at time of isolate collection.

TOB was the only antibiotic tested that showed AMR in all tested isolates, including the type strain, *P. melaninogenica* LMG 28911 (**Table 2**). This finding that concurs with prior work reporting 100% TOB resistance among 157 *Prevotella* isolates59. A CF microbiome study also found that *Prevotella* abundance was largely unchanged upon TOB treatment64. These collective findings indicate that all respiratory *Prevotella* species are intrinsically resistant to TOB. Notably, we observed that the median relative *Prevotella* abundance according to metataxonomics was quite different between CF (9.6%) vs. COPD (0.6%) specimens (**Table S3**). Although the basis of this higher relative *Prevotella* abundance in CF is unknown, one possibility is that most of our CF participants were receiving TOB at the time of sample collection (**Table S1**; **Table 2**), whereas no COPD participant was receiving this antibiotic. The predominance of *Prevotella* spp. in the CF airways might thus be attributed to much more frequent TOB use in this cohort65-67. Although not assessed in this study, adding TOB to selective media may improve *Prevotella* recovery rates from respiratory specimens by suppressing growth of other taxa. After TOB, CLI resistance was most common (6/14 CF isolates), followed by AZM (4/14 CF isolates), MZ (1/13 COPD and 3/14 CF isolates), and CAZ (2/14 CF isolates). Decreased susceptibility towards AMC, CAZ, and DOX was common in both CF and COPD isolates; however, only CF isolates exhibited decreased AZM susceptibility (6/14 isolates). These results provide further evidence of bystander selection in *Prevotella* spp., with this effect especially noticeable in CF-derived isolates.

AMR genes represent a major concern in infection control as these genes can rapidly transfer AMR among different genera and species. To determine AMR gene potential, CARD analysis was performed on the 35 CF-COPD *Prevotella* genomes and used to correlate horizontally acquired AMR-conferring genes with AMR/decreased susceptibility phenotypes. Consistent with their higher AMR rates, CF isolates demonstrated greater prevalence of AMR-conferring genes compared with COPD (**Table 2**). Strikingly, 67% CF strains encoded *ermF*, a horizontally-acquired gene that confers resistance towards macrolide, lincosamide, and streptogramin-class antibiotics such as AZM and CLI47,61,68. In contrast, *ermF* was absent in all COPD strains, a somewhat unexpected finding given that COPD participants are frequently given long-term macrolides (e.g. azithromycin, erythromycin) as prophylaxis to reduce exacerbation frequency69. None of the CF or COPD participants were receiving these antibiotics at the time of sample collection, suggesting that cross-resistance or prior CF antibiotic treatment may be driving *ermF* acquisition and maintenance in *Prevotella. ermF* presence correlated well with AZM and CLI AMR/decreased susceptibility (87% and 75%, respectively). *ermF* rates in our Australian isolates exceed those previously reported (23-68% in UK, USA, and Italian isolates from both healthy and diseased participants)68,70. This discrepancy is likely due to only respiratory disease-associated isolates being tested in our study, but may also be due to differences in AMR cut-offs between studies, *in vivo* antibiotic exposure levels, and potentially species or strain composition, noting that one study only identified most isolates to the genus level68.

*tetQ*, which confers resistance to tetracycline-class antibiotics (e.g. DOX)68,70, was present in 54% of our strains, and more common in CF-derived strains (67% vs. 41%). This rate is much higher than *tetQ* prevalence reported in UK CF (21%)68 and human-derived Dutch (30%)61 strains.. Our *tetQ*-harbouring strains consistently showed much smaller DOX disc diameters (12-30mm; median 22mm) compared with *tetQ-*negative strains (40-64mm; median 54mm) with 100% correlation between *tetQ* presence and decreased DOX susceptibility (**Table 2**). A prior study reported that 83% of *tetQ*-positive isolates demonstrated intermediate or full resistance towards tetracycline68, suggesting that *tetQ* is a reliable marker of reduced susceptibility towards tetracycline-class antibiotics. Taken together, these results confirm that *tetQ* presence robustly confers decreased DOX susceptibility in *Prevotella* spp. Although our dataset is small, we propose *Prevotella* spp. DOX disc diameters of ≥40mm to denote strains lacking *tetQ*, and ≤30mm to denote strains that may be harbouring *tetQ*, with the caveat that other enigmatic determinants conferring reduced DOX susceptibility (e.g. chromosomal mutations) may be present in lieu of *tetQ*. These criteria also differentiate *Prevotella* spp. exhibiting decreased DOX susceptibility from DOX-sensitive strains.

The *cfxA*-type β-lactamase gene is the most clinically important AMR mechanism in *Prevotella* spp.71 as it confers AMR towards cephamycins61 and some other β-lactam antibiotics68, which are commonly used in respiratory infection treatment59. Concerningly, *cfxA* is often encoded on mobile elements, meaning that AMR can potentially be transmitted from *Prevotella* spp. to known pathogenic taxa in the respiratory microbiome72. We found *cfxA* alleles 2, 3, 4, and/or 5 in 54% strains (72% CF and 35% COPD). *cfxA*3 and *cfxA*2 were most common, being present in 37% and 14% strains, respectively. In comparison, a study of UK and US isolates reported much higher *cfxA* carriage rates (86%), with 29% and 52% strains encoding *cfxA*3 and *cfxA*2, respectively. Their study also found that *cfxA*3-encoding isolates exhibited significantly higher CAZ resistance levels68, which we also observed in our *cfxA*3-harbouring strains (median disc diameter=12mm vs. 32mm in non-*cfxA*3 strains) (**Table 2**). In *cfxA*3-encoding strains displaying CAZ sensitivity, reduced *cfxA*3 expression73 may explain this phenotype, although this was not explored in our study. In addition, we did not test *Prevotella* susceptibility towards cephamycins so were unable to determine whether *cfxA* variants correlated with these β-lactam antibiotics.

MZ resistance was evident in 4/28 tested strains (3 CF; 1 COPD), despite none of the participants receiving nitroimidazoles at the time of sample collection. Other than TOB, MZ was the only antibiotic that exhibited AMR in any of the COPD strains. The *nim* genes, which encode for nitroimidazole reductase enzymes and have previously been associated with MZ resistance in *P. baroniae*74 and *P. bivia*75, were not detected in any isolate. Therefore, the basis of MZ resistance in our isolates remains unknown.

Of the eight species identified in our study, four are conventionally considered dark-pigmented – *P. denticola, P. histicola, P. melaninogenica*, and *P. nigrescens*. Pigmented *Prevotella* reportedly harbour a greater cache of virulence factors and are more commonly associated with infection47. However, we found limited phenotypic differences within and between species on non-selective anaerobic blood medium at Day 3 (**Figure 4**), with most presenting as convex, shiny, white colonies; only *P. nigrescens* demonstrated tan pigmentation in the primary streak at this time point (**Figure 4F**). In contrast, a dark brown to black pigment was produced by *P. nigrescens* on PVCA by Day 3 (**Figure 5D**), but not on selective or non-selective anaerobic media, even by Day 10, with only medium-brown colonies present (**Figure 5B and 5C**). *P. histicola* demonstrated tan-coloured pigmentation on both non-selective anaerobic medium and PVCA, but only after five days (**Figure 6**), and unlike *P. nigrescens*, did not develop black pigment on PVCA at Days 5 or 10. *P. melaninogenica* reportedly produces a characteristic dark pigment that is commonly associated with the species76,77; yet, no pigmentation was observed in the six *P. melaninogenica* isolates on any medium by Day 10. A likely explanation for this phenomenon is that different strains require different media for pigment production78. Collectively, the limited pigmentation observed among *Prevotella* spp. suggests that phenotypic characteristics alone should not be used for *Prevotella* genus or species identification.

![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/31/2021.08.30.21262864/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2021/08/31/2021.08.30.21262864/F4)

Figure 4. 
Three-day growth of *Prevotella* spp. on non-selective anaerobic blood medium. (A) *P. melaninogenica*; (B) *P. histicola*; (C) *P. denticola*; (D) *P. nanceiensis*; (E) *P. jejuni*; (F) *P. nigrescens*; (G) *P. salivae*; and (H) *P. buccae*.

![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/31/2021.08.30.21262864/F5.medium.gif)

[Figure 5.](http://medrxiv.org/content/early/2021/08/31/2021.08.30.21262864/F5)

Figure 5. 
Phenotype comparison of *Prevotella nigrescens* on non-selective anaerobic blood medium (top) and polyvitaminic-supplemented chocolate agar (bottom) at (A and D) 3 days; (B and E) 5 days; and (C and F) 10 days.

![Figure 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/31/2021.08.30.21262864/F6.medium.gif)

[Figure 6.](http://medrxiv.org/content/early/2021/08/31/2021.08.30.21262864/F6)

Figure 6. 
Phenotype comparison of *Prevotella histicola* on non-selective anaerobic blood medium (top) and polyvitaminic-supplemented chocolate agar (bottom) at (A and C) 3 days and (B and D) 5 days.

We acknowledge several study limitations. First, small sample size (43 enrolled participants; 69 samples) and inclusion of only COPD and CF samples makes it difficult to generalise our results across a larger cohorts. Nevertheless, the broad agreement of our findings with previous literature10,14,19 demonstrates a reasonable estimation of species diversity in these cohorts. Second, due to non-anaerobic sample handling and storage, non-aerotolerant *Prevotella* species may have been rendered non-viable and thus not recoverable45. Future establishment of disc diffusion guidelines in *Prevotella* will enable more accurate assessment of AMR phenotypes in this genus. Third, due to our use of a strict AMR cut-off (i.e. 6mm), AMR rates in our *Prevotella* are almost certainly underestimated. Fourth, we did not assess healthy airway microbiota due to difficulties associated with obtaining respiratory secretions from healthy participants. Future work should continue to seek strain-level resolution of *Prevotella* strains in CF and COPD airways, including the design of species-specific assays towards commonly isolated and known pathogenic species, assessment of *Prevotella* from healthy airways, and incorporation of higher-resolution molecular technologies such as metagenomics and metatranscriptomics to gain a clearer understanding of *Prevotella* function within complex microbiomes. Co-culture experiments between *Prevotella* species and known respiratory pathogens, such as *P. aeruginosa*, would be beneficial for elucidating the role of *Prevotella* in potentiating or abrogating pathogen virulence. Finally, phylogeographic analysis of *Prevotella* genomes may permit a deeper understanding of *Prevotella* spp. distribution and strain diversity in chronic lung disease globally, and population genomic approaches will yield a better understanding of *Prevotella* origin and capacity for airway-to-airway transmission.

## Supporting information

Merged supplemental files [[supplements/262864_file09.docx]](pending:yes)

## Data Availability

Thirty-five Prevotella spp. genomes generated in this study are available in the Sequence Read Archive (SRA) and GenBank databases under BioProject accession PRJNA742126.

[https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA742126](https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA742126) 

## Acknowledgements

We thank Tamieka Fraser (University of the Sunshine Coast) for laboratory assistance, Peter Timms (University of the Sunshine Coast) and Laura Sherrard (Queen’s University Belfast) for insightful feedback on this manuscript, Tania Duarte (University of Queensland) for helpful discussions about the challenges of *Prevotella* identification using current metagenomic pipelines, and the Belgian Co-ordinated Collection of Micro-organisms (BCCM) for providing type strains. This work was funded by Advance Queensland (awards AQRF13016-17RD2 and AQIRF0362018), the Wishlist Sunshine Coast Hospital Foundation (award 2019-24), and the University of the Sunshine Coast.

*   Received August 30, 2021.
*   Revision received August 30, 2021.
*   Accepted August 31, 2021.


*   © 2021, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## References

1.  1.Sze MA, Hogg JC, Sin DD. Bacterial microbiome of lungs in COPD. Int J Chron Obstruct Pulmon Dis 2014; 9: 229–38.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24591822&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 

2.  2.Fenker DE, McDaniel CT, Panmanee W, et al. A comparison between two pathophysiologically different yet microbiologically similar lung diseases: cystic fibrosis and chronic obstructive pulmonary disease. Int J Respir Pulm Med 2018; 5(2).
    
    

3.  3.Wedzicha JA, Seemungal TA. COPD exacerbations: defining their cause and prevention. Lancet 2007; 370(9589): 786–96.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(07)61382-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17765528&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000249209300033&link_type=ISI) 

4.  4.Lyczak JB, Cannon CL, Pier GB. Lung infections associated with cystic fibrosis. Clin Microbiol Rev 2002; 15(2): 194–222.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiY21yIjtzOjU6InJlc2lkIjtzOjg6IjE1LzIvMTk0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDgvMzEvMjAyMS4wOC4zMC4yMTI2Mjg2NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

5.  5.Jacobs DM, Pandit U, Sethi S. Acute exacerbations in chronic obstructive pulmonary disease: should we use antibiotics and if so, which ones? Curr Opin Infect Dis 2019; 32(2): 143–51.
    
    

6.  6.Larsen JM. The immune response to Prevotella bacteria in chronic inflammatory disease. Immunology 2017; 151(4): 363–74.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/imm.12760&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28542929&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 

7.  7.Gilpin DF, Nixon KA, Bull M, et al. Evidence of persistence of Prevotella spp. in the cystic fibrosis lung. J Med Microbiol 2017; 66(6): 825–32.
    
    

8.  8.Ditz B, Christenson S, Rossen J, et al. Sputum microbiome profiling in COPD: beyond singular pathogen detection. Thorax 2020; 75(4): 338–44.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToidGhvcmF4am5sIjtzOjU6InJlc2lkIjtzOjg6Ijc1LzQvMzM4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDgvMzEvMjAyMS4wOC4zMC4yMTI2Mjg2NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

9.  9.Pattison SH, Rogers GB, Crockard M, Elborn JS, Tunney MM. Molecular detection of CF lung pathogens: current status and future potential. J Cyst Fibros 2013; 12(3): 194–205.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jcf.2013.01.007&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23402821&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 

10. 10.Lee SW, Kuan CS, Wu LS, Weng JT. Metagenome and metatranscriptome profiling of moderate and severe COPD sputum in Taiwanese Han males. PLoS One 2016; 11(7): e0159066.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0159066&link_type=DOI) 

11. 11.Muhlebach MS, Hatch JE, Einarsson GG, et al. Anaerobic bacteria cultured from cystic fibrosis airways correlate to milder disease: a multisite study. Eur Respir J 2018; 52(1).
    
    

12. 12.Amat S, Lantz H, Munyaka PM, Willing BP. Prevotella in pigs: the positive and negative associations with production and health. Microorganisms 2020; 8(10).
    
    

13. 13.Wybo I, Soetens O, De Bel A, et al. Species identification of clinical Prevotella isolates by matrix-assisted laser desorption ionization-time of flight mass spectrometry. J Clin Microbiol 2012; 50(4): 1415–8.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNtIjtzOjU6InJlc2lkIjtzOjk6IjUwLzQvMTQxNSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzA4LzMxLzIwMjEuMDguMzAuMjEyNjI4NjQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

14. 14.Tunney MM, Field TR, Moriarty TF, et al. Detection of anaerobic bacteria in high numbers in sputum from patients with cystic fibrosis. Am J Respir Crit Care Med 2008; 177(9): 995–1001.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1164/rccm.200708-1151OC&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18263800&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000255486100013&link_type=ISI) 

15. 15.Cuthbertson L, Walker AW, Oliver AE, et al. Lung function and microbiota diversity in cystic fibrosis. Microbiome 2020; 8(1): 45.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s40168-020-00810-3&link_type=DOI) 

16. 16.Han MK, Huang YJ, Lipuma JJ, et al. Significance of the microbiome in obstructive lung disease. Thorax 2012; 67(5): 456–63.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToidGhvcmF4am5sIjtzOjU6InJlc2lkIjtzOjg6IjY3LzUvNDU2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDgvMzEvMjAyMS4wOC4zMC4yMTI2Mjg2NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

17. 17.Wang Z, Maschera B, Lea S, et al. Airway host-microbiome interactions in chronic obstructive pulmonary disease. Respir Res 2019; 20(1): 113.
    
    

18. 18.Mirkovic B, Murray MA, Lavelle GM, et al. The role of short-chain fatty acids, produced by anaerobic bacteria, in the cystic fibrosis airway. Am J Respir Crit Care Med 2015; 192(11): 1314–24.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1164/rccm.201505-0943OC&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26266556&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 

19. 19.Field TR, Sibley CD, Parkins MD, Rabin HR, Surette MG. The genus Prevotella in cystic fibrosis airways. Anaerobe 2010; 16(4): 337–44.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.anaerobe.2010.04.002&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20412865&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000282076300005&link_type=ISI) 

20. 20.Ulrich M, Beer I, Braitmaier P, et al. Relative contribution of Prevotella intermedia and Pseudomonas aeruginosa to lung pathology in airways of patients with cystic fibrosis. Thorax 2010; 65(11): 978–84.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToidGhvcmF4am5sIjtzOjU6InJlc2lkIjtzOjk6IjY1LzExLzk3OCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzA4LzMxLzIwMjEuMDguMzAuMjEyNjI4NjQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

21. 21.Jorth P, Ehsan Z, Rezayat A, et al. Direct lung sampling indicates that established pathogens dominate early infections in children with cystic fibrosis. Cell Rep 2019; 27(4): 1190–204.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.celrep.2019.03.086&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 

22. 22.Zemanick ET, Wagner BD, Sagel SD, Stevens MJ, Accurso FJ, Harris JK. Reliability of quantitative real-time PCR for bacterial detection in cystic fibrosis airway specimens. PLoS One 2010; 5(11): e15101.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0015101&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21152087&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 

23. 23.O’Neill K, Bradley JM, Johnston E, et al. Reduced bacterial colony count of anaerobic bacteria is associated with a worsening in lung clearance index and inflammation in cystic fibrosis. PLoS One 2015; 10(5): e0126980.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0126980&link_type=DOI) 

24. 24.Einarsson GG, Zhao J, LiPuma JJ, Downey DG, Tunney MM, Elborn JS. Community analysis and co-occurrence patterns in airway microbial communities during health and disease. ERJ Open Res 2019; 5(3).
    
    

25. 25.Wang Z, Liu H, Wang F, et al. A refined view of airway microbiome in chronic obstructive pulmonary disease at species and strain-levels. Front Microbiol 2020; 11: 1758.
    
    

26. 26.Human Microbiome Jumpstart Reference Strains Consortium. A catalog of reference genomes from the human microbiome. Science 2010; 328(5981): 994–9.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzMjgvNTk4MS85OTQiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wOC8zMS8yMDIxLjA4LjMwLjIxMjYyODY0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

27. 27.Gupta VK, Chaudhari NM, Iskepalli S, Dutta C. Divergences in gene repertoire among the reference Prevotella genomes derived from distinct body sites of human. BMC Genomics 2015; 16: 153.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12864-015-1350-6&link_type=DOI) 

28. 28.Zhang Y, Zhen M, Zhan Y, Song Y, Zhang Q, Wang J. Population-genomic insights into variation in Prevotella intermedia and Prevotella nigrescens isolates and its association with periodontal disease. Front Cell Infect Microbiol 2017; 7: 409.
    
    

29. 29.Tett A, Huang KD, Asnicar F, et al. The Prevotella copri complex comprises four distinct clades underrepresented in Westernized populations. Cell Host Microbe 2019; 26(5): 666–79.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.chom.2019.08.018&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31607556&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 

30. 30.de Lamballerie X, Zandotti C, Vignoli C, Bollet C, de Micco P. A one-step microbial DNA extraction method using “Chelex 100” suitable for gene amplification. Res Microbiol 1992; 143(8): 785–90.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0923-2508(92)90107-Y&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1298031&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1992JX00400006&link_type=ISI) 

31. 31.Triant DA, Whitehead A. Simultaneous extraction of high-quality RNA and DNA from small tissue samples. J Hered 2009; 100(2): 246–50.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jhered/esn083&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18840898&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000263420000014&link_type=ISI) 

32. 32.Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol 1990; 215(3): 403–10.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1006/jmbi.1990.9999&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=2231712&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1990ED16700008&link_type=ISI) 

33. 33.Klaassen CH, Jeunink MA, Prinsen CF, et al. Quantification of human DNA in feces as a diagnostic test for the presence of colorectal cancer. Clin Chem 2003; 49(7): 1185–7.
    
    [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiY2xpbmNoZW0iO3M6NToicmVzaWQiO3M6OToiNDkvNy8xMTg1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDgvMzEvMjAyMS4wOC4zMC4yMTI2Mjg2NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

34. 34.Nadkarni MA, Martin FE, Jacques NA, Hunter N. Determination of bacterial load by real-time PCR using a broad-range (universal) probe and primers set. Microbiology (Reading) 2002; 148(Pt 1): 257–66.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1099/00221287-148-1-257&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11782518&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000173355100027&link_type=ISI) 

35. 35.Price EP, Dale JL, Cook JM, et al. Development and validation of Burkholderia pseudomallei-specific real-time PCR assays for clinical, environmental or forensic detection applications. PLoS One 2012; 7(5): e37723.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0037723&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22624061&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 

36. 36.Hudzicki J. Kirby-Bauer Disk Diffusion Susceptibility Test Protocol. 2016. [https://asm.org/Protocols/Kirby-Bauer-Disk-Diffusion-Susceptibility-Test-Pro](https://asm.org/Protocols/Kirby-Bauer-Disk-Diffusion-Susceptibility-Test-Pro) (accessed 11Dec20).
    
    

37. 37.Bankevich A, Nurk S, Antipov D, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol 2012; 19(5): 455–77.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1089/cmb.2012.0021&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22506599&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 

38. 38.Jolley KA, Bliss CM, Bennett JS, et al. Ribosomal multilocus sequence typing: universal characterization of bacteria from domain to strain. Microbiology (Reading) 2012; 158(Pt 4): 1005–15.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1099/mic.0.055459-0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22282518&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000303633200015&link_type=ISI) 

39. 39.Larkin MA, Blackshields G, Brown NP, et al. Clustal W and Clustal X version 2.0. Bioinformatics 2007; 23(21): 2947–8.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btm404&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17846036&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000251197700021&link_type=ISI) 

40. 40.Sarovich DS, Price EP. SPANDx: a genomics pipeline for comparative analysis of large haploid whole genome re-sequencing datasets. BMC Res Notes 2014; 7: 618.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1756-0500-7-618&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25201145&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 

41. 41.Huang W, Li L, Myers JR, Marth GT. ART: a next-generation sequencing read simulator. Bioinformatics 2012; 28(4): 593–4.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btr708&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22199392&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000300490500023&link_type=ISI) 

42. 42.Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 2014; 30(15): 2114–20.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btu170&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24695404&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000340049100004&link_type=ISI) 

43. 43.Price EP, Viberg LT, Kidd TJ, Bell SC, Currie BJ, Sarovich DS. Transcriptomic analysis of longitudinal Burkholderia pseudomallei infecting the cystic fibrosis lung. Microb Genom 2018; 4(8).
    
    

44. 44.Alcock BP, Raphenya AR, Lau TTY, et al. CARD 2020: antibiotic resistome surveillance with the comprehensive antibiotic resistance database. Nucleic Acids Res 2020; 48(D1): D517–D25.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkz935&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31665441&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 

45. 45.Silva VL, Carvalho MA, Nicoli JR, Farias LM. Aerotolerance of human clinical isolates of Prevotella spp. J Appl Microbiol 2003; 94(4): 701–7.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12631206&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 

46. 46.Aziz A, Currie BJ, Mayo M, Sarovich DS, Price EP. Comparative genomics confirms a rare melioidosis human-to-human transmission event and reveals incorrect phylogenomic reconstruction due to polyclonality. Microb Genom 2020; 6(2).
    
    

47. 47.Xie Y, Chen J, He J, et al. Antimicrobial resistance and prevalence of resistance genes of obligate anaerobes isolated from periodontal abscesses. J Periodontol 2014; 85(2): 327–34.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1902/jop.2013.130081&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23659425&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000331140100017&link_type=ISI) 

48. 48.Goddard AF, Staudinger BJ, Dowd SE, et al. Direct sampling of cystic fibrosis lungs indicates that DNA-based analyses of upper-airway specimens can misrepresent lung microbiota. Proc Natl Acad Sci U S A 2012; 109(34): 13769–74.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTA5LzM0LzEzNzY5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDgvMzEvMjAyMS4wOC4zMC4yMTI2Mjg2NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

49. 49.Pragman AA, Lyu T, Baller JA, et al. The lung tissue microbiota of mild and moderate chronic obstructive pulmonary disease. Microbiome 2018; 6(1): 7.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s40168-017-0381-4&link_type=DOI) 

50. 50.Dickson RP, Erb-Downward JR, Freeman CM, et al. Bacterial topography of the healthy human lower respiratory tract. mBio 2017; 8(1).
    
    

51. 51.Bassis CM, Erb-Downward JR, Dickson RP, et al. Analysis of the upper respiratory tract microbiotas as the source of the lung and gastric microbiotas in healthy individuals. mBio 2015; 6(2): e00037.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1128/mBio.00037-15&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25736890&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 

52. 52.Larsen JM, Steen-Jensen DB, Laursen JM, et al. Divergent pro-inflammatory profile of human dendritic cells in response to commensal and pathogenic bacteria associated with the airway microbiota. PLoS One 2012; 7(2): e31976.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0031976&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22363778&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 

53. 53.Hedberg ME, Israelsson A, Moore ERB, et al. Prevotella jejuni sp. nov., isolated from the small intestine of a child with coeliac disease. Int J Syst Evol Microbiol 2013; 63(Pt 11): 4218–23.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1099/ijs.0.052647-0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23793857&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 

54. 54.Enaud R, Prevel R, Ciarlo E, et al. The gut-lung axis in health and respiratory diseases: a place for inter-organ and inter-kingdom crosstalks. Front Cell Infect Microbiol 2020; 10: 9.
    
    

55. 55.Bertelsen A, Elborn SJ, Schock BC. Toll like Receptor signalling by Prevotella histicola activates alternative NF-kappaB signalling in cystic fibrosis bronchial epithelial cells compared to P. aeruginosa. PLoS One 2020; 15(10): e0235803.
    
    

56. 56.Di Lorenzo F, Silipo A, Matier T, et al. Prevotella denticola lipopolysaccharide from a cystic fibrosis isolate possesses a unique chemical structure. European J Org Chem 2016; 9: 1732–8.
    
    

57. 57.Yanagisawa M, Kuriyama T, Williams DW, Nakagawa K, Karasawa T. Proteinase activity of Prevotella species associated with oral purulent infection. Curr Microbiol 2006; 52(5): 375–8.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16586023&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 

58. 58.Alauzet C, Marchandin H, Lozniewski A. New insights into Prevotella diversity and medical microbiology. Future Microbiol 2010; 5(11): 1695–718.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2217/fmb.10.126&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21133690&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000286290800013&link_type=ISI) 

59. 59.Sherrard LJ, Graham KA, McGrath SJ, et al. Antibiotic resistance in Prevotella species isolated from patients with cystic fibrosis. J Antimicrob Chemother 2013; 68(10): 2369–74.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jac/dkt191&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23696621&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000326973500033&link_type=ISI) 

60. 60.Brill SE, Law M, El-Emir E, et al. Effects of different antibiotic classes on airway bacteria in stable COPD using culture and molecular techniques: a randomised controlled trial. Thorax 2015; 70(10): 930–8.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToidGhvcmF4am5sIjtzOjU6InJlc2lkIjtzOjk6IjcwLzEwLzkzMCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzA4LzMxLzIwMjEuMDguMzAuMjEyNjI4NjQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

61. 61.Veloo ACM, Baas WH, Haan FJ, Coco J, Rossen JW. Prevalence of antimicrobial resistance genes in Bacteroides spp. and Prevotella spp. Dutch clinical isolates. Clin Microbiol Infect 2019; 25(9): 1156 e9–e13.
    
    

62. 62.Sherrard LJ, Tunney MM, Elborn JS. Antimicrobial resistance in the respiratory microbiota of people with cystic fibrosis. Lancet 2014; 384(9944): 703–13.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(14)61137-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25152272&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 

63. 63.Hahn A, Burrell A, Fanous H, et al. Antibiotic multidrug resistance in the cystic fibrosis airway microbiome is associated with decreased diversity. Heliyon 2018; 4(9): e00795.
    
    

64. 64.Nelson MT, Wolter DJ, Eng A, et al. Maintenance tobramycin primarily affects untargeted bacteria in the CF sputum microbiome. Thorax 2020; 75(9): 780–90.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToidGhvcmF4am5sIjtzOjU6InJlc2lkIjtzOjg6Ijc1LzkvNzgwIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDgvMzEvMjAyMS4wOC4zMC4yMTI2Mjg2NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

65. 65.Bittar F, Richet H, Dubus JC, et al. Molecular detection of multiple emerging pathogens in sputa from cystic fibrosis patients. PLoS One 2008; 3(8): e2908.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0002908&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18682840&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 

66. 66.Shakil S, Khan R, Zarrilli R, Khan AU. Aminoglycosides versus bacteria--a description of the action, resistance mechanism, and nosocomial battleground. J Biomed Sci 2008; 15(1): 5–14.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s11373-007-9194-y&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17657587&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000251371100001&link_type=ISI) 

67. 67.Thornton CS, Surette MG. Potential contributions of anaerobes in cystic fibrosis airways. J Clin Microbiol 2021; 59(3).
    
    

68. 68.Sherrard LJ, Schaible B, Graham KA, et al. Mechanisms of reduced susceptibility and genotypic prediction of antibiotic resistance in Prevotella isolated from cystic fibrosis (CF) and non-CF patients. J Antimicrob Chemother 2014; 69(10): 2690–8.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jac/dku192&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24917582&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 

69. 69.Parameswaran GI, Sethi S. Long-term macrolide therapy in chronic obstructive pulmonary disease. CMAJ 2014; 186(15): 1148–52.
    
    [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiY21haiI7czo1OiJyZXNpZCI7czoxMToiMTg2LzE1LzExNDgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wOC8zMS8yMDIxLjA4LjMwLjIxMjYyODY0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

70. 70.Arzese AR, Tomasetig L, Botta GA. Detection of tetQ and ermF antibiotic resistance genes in Prevotella and Porphyromonas isolates from clinical specimens and resident microbiota of humans. J Antimicrob Chemother 2000; 45(5): 577–82.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jac/45.5.577&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10797077&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000087238900004&link_type=ISI) 

71. 71.Fernández-Canigia L, Cejas D, Gutkind G, Radice M. Detection and genetic characterization of beta-lactamases in Prevotella intermedia and Prevotella nigrescens isolated from oral cavity infections and peritonsillar abscesses. Anaerobe 2015; 33: 8–13.
    
    

72. 72.Sherrard LJ, McGrath SJ, McIlreavey L, et al. Production of extended-spectrum β-lactamases and the potential indirect pathogenic role of Prevotella isolates from the cystic fibrosis respiratory microbiota. Int J Antimicrob Agents 2016; 47(2): 140–5.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijantimicag.2015.12.004&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26774156&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F31%2F2021.08.30.21262864.atom) 

73. 73.Jolivet-Gougeon A, Tamanai-Shacoori Z, Desbordes L, Burggraeve N, Cormier M, Bonnaure-Mallet M. Genetic analysis of an ambler class A extended-spectrum beta-lactamase from Capnocytophaga ochracea. J Clin Microbiol 2004; 42(2): 888–90.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNtIjtzOjU6InJlc2lkIjtzOjg6IjQyLzIvODg4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDgvMzEvMjAyMS4wOC4zMC4yMTI2Mjg2NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

74. 74.Alauzet C, Mory F, Teyssier C, et al. Metronidazole resistance in Prevotella spp. and description of a new nim gene in Prevotella baroniae. Antimicrob Agents Chemother 2010; 54(1): 60–4.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWFjIjtzOjU6InJlc2lkIjtzOjc6IjU0LzEvNjAiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wOC8zMS8yMDIxLjA4LjMwLjIxMjYyODY0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

75. 75.Veloo ACM, Chlebowicz M, Winter HLJ, Bathoorn D, Rossen JWA. Three metronidazole-resistant Prevotella bivia strains harbour a mobile element, encoding a novel nim gene, nimK, and an efflux small MDR transporter. J Antimicrob Chemother 2018; 73(10): 2687–90.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jac/dky236&link_type=DOI) 

76. 76.Allison HE, Hillman JD. Cloning and characterization of a Prevotella melaninogenica hemolysin. Infect Immun 1997; 65(7): 2765–71.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiaWFpIjtzOjU6InJlc2lkIjtzOjk6IjY1LzcvMjc2NSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzA4LzMxLzIwMjEuMDguMzAuMjEyNjI4NjQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

77. 77.Zhou X, Li Y. Subgingival microbes: from healthy microflora to disease. Atlas of Oral Microbiology: Elsevier; 2015.
    
    

78. 78.1.  Baron S
    
    Finegold SM. Anaerobic Gram-negative bacilli (Chapter 20). In: Baron S, ed. Medical Microbiology (4th edition). Galveston, TX: University of Texas Medical Branch at Galveston; 1996.