Citizen science monitoring reveals links between honeybee health, pesticide exposure and seasonal availability of floral resources.

We use a national citizen science monitoring scheme to quantify how agricultural intensification affects honeybee diet breadth (number of plant species). To do this we used DNA metabarcoding to identify the plants present in 527 honey samples collected in 2019 across Great Britain. The species richness of forage plants was negatively correlated with arable cropping area, although this was only found early in the year when the abundance of flowering plants was more limited. Within intensively farmed areas, honeybee diets were dominated by Brassica crops (including oilseed rape). We demonstrate how the structure and complexity of honeybee foraging relationships with plants is negatively affected by the area of arable crops surrounding hives. Using information collected from the beekeepers on the incidence of an economically damaging bee disease (Deformed Wing Virus) we found that the occurrence of this disease increased where bees foraged in agricultural land where there was a high use of foliar insecticides. Understanding impacts of land use on resource availability is fundamental to assessing long-term viability of pollinator populations. These findings highlight the importance of supporting temporally timed resources as mitigation strategies to support wider pollinator population viability.

Integration of DNA extraction, metabarcoding and an informatics pipeline to underpin a national citizen science honey monitoring scheme.

Worldwide honeybees (Apis mellifera L.) are one of the most widely kept domesticated animals, supporting domestic and commercial livelihoods through the production of honey and wax, as well as in the delivery of pollination services to crops. Quantifying which plant species are foraged upon by honeybees provides insights into their nutritional status as well as patterns of landscape scale habitat utilization. Here we outline a rapid and reproducible methodology for identifying environmental DNA (eDNA) originating principally from pollen grains suspended within honey. The process is based on a DNA extraction incorporating vacuum filtration prior to universal eukaryotic internal transcribed spacer 2 region (ITS2) amplicon generation, sequencing and identification. To provide a pre-cursor to sequence phylotyping, we outline systems for error-corrected processing amplicon sequence variant abundance tables that removes chimeras. This methodology underpins the new UK National Honey Monitoring Scheme.•We compare the efficacy and speed of centrifugation and filtration systems for removing pollen from honey samples as a precursor to plant DNA barcoding.•We introduce the 'HONEYPI' informatics pipeline, an open access resource implemented in python 2.7, to ensure long-term reproducibility during the process of amplicon sequence variant classification.

Identification of microbial signatures linked to oilseed rape yield decline at the landscape scale.

BACKGROUND: The plant microbiome plays a vital role in determining host health and productivity. However, we lack real-world comparative understanding of the factors which shape assembly of its diverse biota, and crucially relationships between microbiota composition and plant health. Here we investigated landscape scale rhizosphere microbial assembly processes in oilseed rape (OSR), the UK's third most cultivated crop by area and the world's third largest source of vegetable oil, which suffers from yield decline associated with the frequency it is grown in rotations. By including 37 conventional farmers' fields with varying OSR rotation frequencies, we present an innovative approach to identify microbial signatures characteristic of microbiomes which are beneficial and harmful to the host. RESULTS: We show that OSR yield decline is linked to rotation frequency in real-world agricultural systems. We demonstrate fundamental differences in the environmental and agronomic drivers of protist, bacterial and fungal communities between root, rhizosphere soil and bulk soil compartments. We further discovered that the assembly of fungi, but neither bacteria nor protists, was influenced by OSR rotation frequency. However, there were individual abundant bacterial OTUs that correlated with either yield or rotation frequency. A variety of fungal and protist pathogens were detected in roots and rhizosphere soil of OSR, and several increased relative abundance in root or rhizosphere compartments as OSR rotation frequency increased. Importantly, the relative abundance of the fungal pathogen Olpidium brassicae both increased with short rotations and was significantly associated with low yield. In contrast, the root endophyte Tetracladium spp. showed the reverse associations with both rotation frequency and yield to O. brassicae, suggesting that they are signatures of a microbiome which benefits the host. We also identified a variety of novel protist and fungal clades which are highly connected within the microbiome and could play a role in determining microbiome composition. CONCLUSIONS: We show that at the landscape scale, OSR crop yield is governed by interplay between complex communities of both pathogens and beneficial biota which is modulated by rotation frequency. Our comprehensive study has identified signatures of dysbiosis within the OSR microbiome, grown in real-world agricultural systems, which could be used in strategies to promote crop yield. Video abstract.

Contrasting community assembly processes structure lotic bacteria metacommunities along the river continuum.

The heterogeneous nature of lotic habitats plays an important role in the complex ecological and evolutionary processes that structure the microbial communities within them. Due to such complexity, our understanding of lotic microbial ecology still lacks conceptual frameworks for the ecological processes that shape these communities. We explored how bacterial community composition and underlying ecological assembly processes differ between lotic habitats by examining community composition and inferring community assembly processes across four major habitat types (free-living, particle-associated, biofilm on benthic stones and rocks, and sediment). This was conducted at 12 river sites from headwater streams to the main river in the River Thames, UK. Our results indicate that there are distinct differences in the bacterial communities between four major habitat types, with contrasting ecological processes shaping their community assembly processes. While the mobile free-living and particle-associated communities were consistently less diverse than the fixed sediment and biofilm communities, the latter two communities displayed higher homogeneity across the sampling sites. This indicates that the relative influence of deterministic environmental filtering is elevated in sediment and biofilm communities compared with free-living and particle-associated communities, where stochastic processes play a larger role.

Lung function and microbiota diversity in cystic fibrosis.

BACKGROUND: Chronic infection and concomitant airway inflammation is the leading cause of morbidity and mortality for people living with cystic fibrosis (CF). Although chronic infection in CF is undeniably polymicrobial, involving a lung microbiota, infection surveillance and control approaches remain underpinned by classical aerobic culture-based microbiology. How to use microbiomics to direct clinical management of CF airway infections remains a crucial challenge. A pivotal step towards leveraging microbiome approaches in CF clinical care is to understand the ecology of the CF lung microbiome and identify ecological patterns of CF microbiota across a wide spectrum of lung disease. Assessing sputum samples from 299 patients attending 13 CF centres in Europe and the USA, we determined whether the emerging relationship of decreasing microbiota diversity with worsening lung function could be considered a generalised pattern of CF lung microbiota and explored its potential as an informative indicator of lung disease state in CF. RESULTS: We tested and found decreasing microbiota diversity with a reduction in lung function to be a significant ecological pattern. Moreover, the loss of diversity was accompanied by an increase in microbiota dominance. Subsequently, we stratified patients into lung disease categories of increasing disease severity to further investigate relationships between microbiota characteristics and lung function, and the factors contributing to microbiota variance. Core taxa group composition became highly conserved within the severe disease category, while the rarer satellite taxa underpinned the high variability observed in the microbiota diversity. Further, the lung microbiota of individual patient were increasingly dominated by recognised CF pathogens as lung function decreased. Conversely, other bacteria, especially obligate anaerobes, increasingly dominated in those with better lung function. Ordination analyses revealed lung function and antibiotics to be main explanators of compositional variance in the microbiota and the core and satellite taxa. Biogeography was found to influence acquisition of the rarer satellite taxa. CONCLUSIONS: Our findings demonstrate that microbiota diversity and dominance, as well as the identity of the dominant bacterial species, in combination with measures of lung function, can be used as informative indicators of disease state in CF. Video Abstract.

Auricular transcutaneous vagus nerve stimulation improves memory persistence in naïve mice and in an intellectual disability mouse model.

BACKGROUND: Vagus nerve stimulation (VNS) using non-invasive approaches have attracted great attention due to their anti-epileptic, anti-depressive and pro-cognitive effects. It has been proposed that auricular transcutaneous VNS (atVNS) could benefit intellectual disability disorders, but preclinical data supporting this idea is limited. OBJECTIVE: To develop an atVNS device for mice and to test its efficacy on memory performance in naïve mice and in a mouse model for intellectual disability. METHODS: Naïve outbreed CD-1 mice and a model for fragile X syndrome, the Fmr1 knockout (Fmr1KO), were used to assess the effect of atVNS in the novel object-recognition memory performance. RESULTS: We found that atVNS significantly improves memory persistence in naïve mice. Notably, atVNS was efficacious in normalizing the object-recognition memory deficit in the Fmr1KO model. CONCLUSION: Our data show that atVNS improves memory persistence in naïve mice and in a model of intellectual disability and support further studies taking advantage of preclinical mouse models of cognitive disorders.

Cannabinoid type-1 receptor blockade restores neurological phenotypes in two models for Down syndrome.

Intellectual disability is the most limiting hallmark of Down syndrome, for which there is no gold-standard clinical treatment yet. The endocannabinoid system is a widespread neuromodulatory system involved in multiple functions including learning and memory processes. Alterations of this system contribute to the pathogenesis of several neurological and neurodevelopmental disorders. However, the involvement of the endocannabinoid system in the pathogenesis of Down syndrome has not been explored before. We used the best-characterized preclinical model of Down syndrome, the segmentally trisomic Ts65Dn model. In male Ts65Dn mice, cannabinoid type-1 receptor (CB1R) expression was enhanced and its function increased in hippocampal excitatory terminals. Knockdown of CB1R in the hippocampus of male Ts65Dn mice restored hippocampal-dependent memory. Concomitant with this result, pharmacological inhibition of CB1R restored memory deficits, hippocampal synaptic plasticity and adult neurogenesis in the subgranular zone of the dentate gyrus. Notably, the blockade of CB1R also normalized hippocampal-dependent memory in female Ts65Dn mice. To further investigate the mechanisms involved, we used a second transgenic mouse model overexpressing a single gene candidate for Down syndrome cognitive phenotypes, the dual specificity tyrosine-phosphorylation-regulated kinase 1A (DYRK1A). CB1R pharmacological blockade similarly improved cognitive performance, synaptic plasticity and neurogenesis in transgenic male Dyrk1A mice. Our results identify CB1R as a novel druggable target potentially relevant for the improvement of cognitive deficits associated with Down syndrome.

Soil bacterial networks are less stable under drought than fungal networks.

Soil microbial communities play a crucial role in ecosystem functioning, but it is unknown how co-occurrence networks within these communities respond to disturbances such as climate extremes. This represents an important knowledge gap because changes in microbial networks could have implications for their functioning and vulnerability to future disturbances. Here, we show in grassland mesocosms that drought promotes destabilising properties in soil bacterial, but not fungal, co-occurrence networks, and that changes in bacterial communities link more strongly to soil functioning during recovery than do changes in fungal communities. Moreover, we reveal that drought has a prolonged effect on bacterial communities and their co-occurrence networks via changes in vegetation composition and resultant reductions in soil moisture. Our results provide new insight in the mechanisms through which drought alters soil microbial communities with potential long-term consequences, including future plant community composition and the ability of aboveground and belowground communities to withstand future disturbances.

Links between soil microbial communities and plant traits in a species-rich grassland under long-term climate change.

Climate change can influence soil microorganisms directly by altering their growth and activity but also indirectly via effects on the vegetation, which modifies the availability of resources. Direct impacts of climate change on soil microorganisms can occur rapidly, whereas indirect effects mediated by shifts in plant community composition are not immediately apparent and likely to increase over time. We used molecular fingerprinting of bacterial and fungal communities in the soil to investigate the effects of 17 years of temperature and rainfall manipulations in a species-rich grassland near Buxton, UK. We compared shifts in microbial community structure to changes in plant species composition and key plant traits across 78 microsites within plots subjected to winter heating, rainfall supplementation, or summer drought. We observed marked shifts in soil fungal and bacterial community structure in response to chronic summer drought. Importantly, although dominant microbial taxa were largely unaffected by drought, there were substantial changes in the abundances of subordinate fungal and bacterial taxa. In contrast to short-term studies that report high resistance of soil fungi to drought, we observed substantial losses of fungal taxa in the summer drought treatments. There was moderate concordance between soil microbial communities and plant species composition within microsites. Vector fitting of community-weighted mean plant traits to ordinations of soil bacterial and fungal communities showed that shifts in soil microbial community structure were related to plant traits representing the quality of resources available to soil microorganisms: the construction cost of leaf material, foliar carbon-to-nitrogen ratios, and leaf dry matter content. Thus, our study provides evidence that climate change could affect soil microbial communities indirectly via changes in plant inputs and highlights the importance of considering long-term climate change effects, especially in nutrient-poor systems with slow-growing vegetation.

Helminth burden and ecological factors associated with alterations in wild host gastrointestinal microbiota.

Infection by gastrointestinal helminths of humans, livestock and wild animals is common, but the impact of such endoparasites on wild hosts and their gut microbiota represents an important overlooked component of population dynamics. Wild host gut microbiota and endoparasites occupy the same physical niche spaces with both affecting host nutrition and health. However, associations between the two are poorly understood. Here we used the commonly parasitized European shag (Phalacrocorax aristotelis) as a model wild host. Forty live adults from the same colony were sampled. Endoscopy was employed to quantify helminth infection in situ. Microbiota from the significantly distinct proventriculus (site of infection), cloacal and faecal gastrointestinal tract microbiomes were characterised using 16S rRNA gene-targeted high-throughput sequencing. We found increasingly strong associations between helminth infection and microbiota composition progressing away from the site of infection, observing a pronounced dysbiosis in microbiota when samples were partitioned into high- and low-burden groups. We posit this dysbiosis is predominately explained by helminths inducing an anti-inflammatory environment in the proventriculus, diverting host immune responses away from themselves. This study, within live wild animals, provides a vital foundation to better understand the mechanisms that underpin the three-way relationship between helminths, microbiota and hosts.

Soil Fungal:Bacterial Ratios Are Linked to Altered Carbon Cycling.

Despite several lines of observational evidence, there is a lack of consensus on whether higher fungal:bacterial (F:B) ratios directly cause higher soil carbon (C) storage. We employed RNA sequencing, protein profiling and isotope tracer techniques to evaluate whether differing F:B ratios are associated with differences in C storage. A mesocosm (13)C labeled foliar litter decomposition experiment was performed in two soils that were similar in their physico-chemical properties but differed in microbial community structure, specifically their F:B ratio (determined by PLFA analyses, RNA sequencing and protein profiling; all three corroborating each other). Following litter addition, we observed a consistent increase in abundance of fungal phyla; and greater increases in the fungal dominated soil; implicating the role of fungi in litter decomposition. Litter derived (13)C in respired CO2 was consistently lower, and residual (13)C in bulk SOM was higher in high F:B soil demonstrating greater C storage potential in the F:B dominated soil. We conclude that in this soil system, the increased abundance of fungi in both soils and the altered C cycling patterns in the F:B dominated soils highlight the significant role of fungi in litter decomposition and indicate that F:B ratios are linked to higher C storage potential.

Environmental DNA metabarcoding of lake fish communities reflects long-term data from established survey methods.

Organisms continuously release DNA into their environments via shed cells, excreta, gametes and decaying material. Analysis of this 'environmental DNA' (eDNA) is revolutionizing biodiversity monitoring. eDNA outperforms many established survey methods for targeted detection of single species, but few studies have investigated how well eDNA reflects whole communities of organisms in natural environments. We investigated whether eDNA can recover accurate qualitative and quantitative information about fish communities in large lakes, by comparison to the most comprehensive long-term gill-net data set available in the UK. Seventy-eight 2L water samples were collected along depth profile transects, gill-net sites and from the shoreline in three large, deep lakes (Windermere, Bassenthwaite Lake and Derwent Water) in the English Lake District. Water samples were assayed by eDNA metabarcoding of the mitochondrial 12S and cytochrome b regions. Fourteen of the 16 species historically recorded in Windermere were detected using eDNA, compared to four species in the most recent gill-net survey, demonstrating eDNA is extremely sensitive for detecting species. A key question for biodiversity monitoring is whether eDNA can accurately estimate abundance. To test this, we used the number of sequence reads per species and the proportion of sampling sites in which a species was detected with eDNA (i.e. site occupancy) as proxies for abundance. eDNA abundance data consistently correlated with rank abundance estimates from established surveys. These results demonstrate that eDNA metabarcoding can describe fish communities in large lakes, both qualitatively and quantitatively, and has great potential as a complementary tool to established monitoring methods.

Respiratory microbiota resistance and resilience to pulmonary exacerbation and subsequent antimicrobial intervention.

Pulmonary symptoms in cystic fibrosis (CF) begin in early life with chronic lung infections and concomitant airway inflammation leading to progressive loss of lung function. Gradual pulmonary function decline is interspersed with periods of acute worsening of respiratory symptoms known as CF pulmonary exacerbations (CFPEs). Cumulatively, CFPEs are associated with more rapid disease progression. In this study multiple sputum samples were collected from adult CF patients over the course of CFPEs to better understand how changes in microbiota are associated with CFPE onset and management. Data were divided into five clinical periods: pre-CFPE baseline, CFPE, antibiotic treatment, recovery, and post-CFPE baseline. Samples were treated with propidium monoazide prior to DNA extraction, to remove the impact of bacterial cell death artefacts following antibiotic treatment, and then characterised by 16S rRNA gene-targeted high-throughput sequencing. Partitioning CF microbiota into core and rare groups revealed compositional resistance to CFPE and resilience to antibiotics interventions. Mixed effects modelling of core microbiota members revealed no significant negative impact on the relative abundance of Pseudomonas aeruginosa across the exacerbation cycle. Our findings have implications for current CFPE management strategies, supporting reassessment of existing antimicrobial treatment regimens, as antimicrobial resistance by pathogens and other members of the microbiota may be significant contributing factors.

Rearing and foraging affects bumblebee (Bombus terrestris) gut microbiota.

Bumblebees are ecologically and economically important as pollinators of crop and wild plants, especially in temperate systems. Species, such as the buff-tailed bumblebee (Bombus terrestris), are reared commercially to pollinate high-value crops. Their highly specific gut microbiota, characterized by low diversity, may affect nutrition and immunity and are likely to be important for fitness and colony health. However, little is known about how environmental factors affect bacterial community structure. We analysed the gut microbiota from three groups of worker bumblebees (B. terrestris) from distinct colonies that varied in rearing and foraging characteristics: commercially reared with restricted foraging (RR); commercially reared with outside foraging (RF); and wild-caught workers (W). Contrary to previous studies, which indicate that bacterial communities are highly conserved across workers, we found that RF individuals had an intermediate community structure compared with RR and W types. Further, this was shaped by differences in the abundances of common operational taxonomic units (OTUs) and the diversity of rare OTUs present, which we propose results from an increase in the variety of carbohydrates obtained through foraging.

The spatial organization and microbial community structure of an epilithic biofilm.

Microbial biofilms are common on lithic surfaces, including stone buildings. However, the ecology of these communities is poorly understood. Few studies have focused on the spatial characteristics of lithobiontic biofilms, despite the fact that spatial structure has been demonstrated to influence ecosystem function (and hence biodegradation) and community diversity. Furthermore, relatively few studies have utilized molecular techniques to characterize these communities, even though molecular methods have revealed unexpected microbial diversity in other habitats. This study investigated (1) the spatial structure and (2) the taxonomic composition of an epilithic biofilm using molecular techniques, namely amplicon pyrosequencing and terminal restriction fragment length polymorphism. Dispersion indices and Mantel correlograms were used to test for the presence of spatial structure in the biofilm. Diversity metrics and rank-abundance distributions (RADs) were also generated. The study revealed spatial structure on a centimetre scale in eukaryotic microbes (fungi and algae), but not the bacteria. Fungal and bacterial communities were highly diverse; algal communities much less so. The RADs were characterized by a distinctive 'hollow' (concave up) profile and long tails of rare taxa. These findings have implications for understanding the ecology of epilithic biofilms and the spatial heterogeneity of stone biodeterioration.

PIPITS: an automated pipeline for analyses of fungal internal transcribed spacer sequences from the Illumina sequencing platform.

Studying fungal biodiversity using data generated from Illumina MiSeq sequencing platforms poses a number of bioinformatic challenges with the analysis typically involving a large number of tools for each analytical step from quality filtering to generating identified operational taxonomic unit (OTU) abundance tables.Here, we introduce PIPITS, an open-source stand-alone suite of software for automated processing of Illumina MiSeq sequences for fungal community analysis. PIPITS exploits a number of state of the art applications to process paired-end reads from quality filtering to producing OTU abundance tables.We provide detailed descriptions of the pipeline and show its utility in the analysis of 9 396 092 sequences generated on the MiSeq platform from Illumina MiSeq. PIPITS is the first automated bioinformatics pipeline dedicated for fungal ITS sequences which incorporates ITSx to extract subregions of ITS and exploits the latest RDP Classifier to classify sequences against the curated UNITE fungal data set.

Implications of multiple freeze-thawing on respiratory samples for culture-independent analyses.

BACKGROUND: Best practice when performing culture-independent microbiological analysis of sputum samples involves their rapid freezing and storage at -80°C. However, accessing biobanked collections can mean that material has been passed through repeated freeze-thaw cycles. The aim of this study was to determine the impact of these cycles on microbial community profiles. METHODS: Sputum was collected from eight adults with cystic fibrosis, and each sample was subjected to six freeze-thaw cycles. Following each cycle, an aliquot was removed and treated with propidium monoazide (PMA) prior to DNA extraction and 16S rRNA gene pyrosequencing. RESULTS: The impact of freeze-thaw cycles was greatest on rare members of the microbiota, with variation beyond that detected with within-sample repeat analysis observed after three cycles. CONCLUSION: Four or more freeze thaw cycles result in a significant distortion of microbiota profiles from CF sputum.

Time between collection and storage significantly influences bacterial sequence composition in sputum samples from cystic fibrosis respiratory infections.

Spontaneously expectorated sputum is traditionally used as the sampling method for the investigation of lower airway infections. While guidelines exist for the handling of these samples for culture-based diagnostic microbiology, there is no comparable consensus on their handling prior to culture-independent analysis. The increasing incorporation of culture-independent approaches in diagnostic microbiology means that it is of critical importance to assess potential biases. The aim of this study was to assess the impact of delayed freezing on culture-independent microbiological analyses and to identify acceptable parameters for sample handling. Sputum samples from eight adult cystic fibrosis (CF) patients were collected and aliquoted into sterile Bijou bottles. Aliquots were stored at room temperature before being frozen at -80 °C for increasing intervals, up to a 72-h period. Samples were treated with propidium monoazide to distinguish live from dead cells prior to DNA extraction, and 16S rRNA gene pyrosequencing was used to characterize their bacterial compositions. Substantial variation was observed in samples with high-diversity bacterial communities over time, whereas little variation was observed in low-diversity communities dominated by recognized CF pathogens, regardless of time to freezing. Partitioning into common and rare species demonstrated that the rare species drove changes in similarity. The percentage abundance of anaerobes over the study significantly decreased after 12 h at room temperature (P = 0.008). Failure to stabilize samples at -80 °C within 12 h of collection results in significant changes in the detected community composition.

Marine bacterial communities are resistant to elevated carbon dioxide levels.

It is well established that the release of anthropogenic-derived CO2 into the atmosphere will be mainly absorbed by the oceans, with a concomitant drop in pH, a process termed ocean acidification. As such, there is considerable interest in how changes in increased CO2 and lower pH will affect marine biota, such as bacteria, which play central roles in oceanic biogeochemical processes. Set within an ecological framework, we investigated the direct effects of elevated CO2, contrasted with ambient conditions on the resistance and resilience of marine bacterial communities in a replicated temporal seawater mesocosm experiment. The results of the study strongly indicate that marine bacterial communities are highly resistant to the elevated CO2 and lower pH conditions imposed, as demonstrated from measures of turnover using taxa­time relationships and distance­decay relationships. In addition, no significant differences in community abundance, structure or composition were observed. Our results suggest that there are no direct effects on marine bacterial communities and that the bacterial fraction of microbial plankton holds enough flexibility and evolutionary capacity to withstand predicted future changes from elevated CO2 and subsequent ocean acidification.

Reducing bias in bacterial community analysis of lower respiratory infections.

High-throughput pyrosequencing and quantitative PCR (Q-PCR) analysis offer greatly improved accuracy and depth of characterisation of lower respiratory infections. However, such approaches suffer from an inability to distinguish between DNA derived from viable and non-viable bacteria. This discrimination represents an important step in characterising microbial communities, particularly in contexts with poor clearance of material or high antimicrobial stress, as non-viable bacteria and extracellular DNA can contribute significantly to analyses. Pre-treatment of samples with propidium monoazide (PMA) is an effective approach to non-viable cell exclusion (NVCE). However, the impact of NVCE on microbial community characteristics (abundance, diversity, composition and structure) is not known. Here, adult cystic fibrosis (CF) sputum samples were used as a paradigm. The effects of PMA treatment on CF sputum bacterial community characteristics, as analysed by pyrosequencing and enumeration by species-specific (Pseudomonas aeruginosa) and total bacterial Q-PCR, were assessed. At the local community level, abundances of both total bacteria and of P. aeruginosa were significantly lower in PMA-treated sample portions. Meta-analysis indicated no overall significant differences in diversity; however, PMA treatment resulted in a significant alteration in local community membership in all cases. In contrast, at the metacommunity level, PMA treatment resulted in an increase in community evenness, driven by an increase in diversity, predominately representing rare community members. Importantly, PMA treatment facilitated the detection of both recognised and emerging CF pathogens, significantly influencing 'core' and 'satellite' taxa group membership. Our findings suggest failure to implement NVCE may result in skewed bacterial community analyses.