Identifying the plant-associated microbiome across aquatic and terrestrial environments: the effects of amplification method on taxa discovery

Structure of plant-associated microeukaryotes in roots and leaves of aquatic and terrestrial plants revealed by blocking peptide-nucleic acid (PNA) amplification

Studies of plant-microbe interactions, including mutualistic, antagonistic, parasitic, or commensal microbes, have greatly benefited our understanding of ecosystem functioning. New molecular identification tools have increasingly revealed the association patterns between microorganisms and plants. Here, we integrated long-read PacBio single-molecule sequencing technology with a blocking protein-nucleic acid (PNA) approach to minimise plant amplicons in a survey of plant-eukaryotic microbe relationships in roots and leaves of different aquatic and terrestrial plants to determine patterns of organ, host, and habitat preferences. The PNA approach reduced the samples' relative amounts of plant reads and did not distort the fungal and other microeukaryotic composition. Our analyses revealed that the eukaryotic microbiomes associated with leaves and roots of aquatic plants exhibit a much larger proportion of non-fungal microorganisms than terrestrial plants, and leaf and root microbiomes are similar. Terrestrial plants had much stronger differentiation of leaf and root microbiomes and stronger partner specificity than aquatic plants.

Fungi-Bacteria Associations in Wilt Diseased Rhizosphere and Endosphere by Interdomain Ecological Network Analysis

In the plant rhizosphere and endosphere, some fungal and bacterial species regularly co-exist, however, our knowledge about their co-existence patterns is quite limited, especially during invasion by bacterial wilt pathogens. In this study, the fungal communities from soil to endophytic compartments were surveyed during an outbreak of tobacco wilt disease caused by Ralstonia solanacearum. It was found that the stem endophytic fungal community was significantly altered by pathogen invasion in terms of community diversity, structure, and composition. The associations among fungal species in the rhizosphere and endosphere infected by R. solanacearum showed more complex network structures than those of healthy plants. By integrating the bacterial dataset, associations between fungi and bacteria were inferred by Inter-Domain Ecological Network (IDEN) approach. It also revealed that infected samples, including both the rhizosphere and endosphere, had more complex interdomain networks than the corresponding healthy samples. Additionally, the bacterial wilt pathogenic Ralstonia members were identified as the keystone genus within the IDENs of both root and stem endophytic compartments. Ralstonia members was negatively correlated with the fungal genera Phoma, Gibberella, and Alternaria in infected roots, as well as Phoma, Gibberella, and Diaporthe in infected stems. This suggested that those endophytic fungi may play an important role in resisting the invasion of R. solanacearum.

Elucidating biofilm diversity on water lily leaves through 16S rRNA amplicon analysis: Comparison of four DNA extraction kits

PREMISE

Within a broader study on leaf fossilization in freshwater environments, a long-term study on the development and microbiome composition of biofilms on the foliage of aquatic plants has been initiated to understand how microbes and biofilms contribute to leaf decay and preservation. Here, water lily leaves are employed as a study model to investigate the relationship between bacterial microbiomes, biodegradation, and fossilization. We compare four DNA extraction kits to reduce biases in interpretation and to identify the most suitable kit for the extraction of DNA from bacteria associated with biofilms on decaying water lily leaves for 16S rRNA amplicon analysis.

METHODS

We extracted surface-associated DNA from Nymphaea leaves in early stages of decay at two water depth levels using four commercially available kits to identify the most suitable protocol for bacterial extraction, applying a mock microbial community standard to enable a reliable comparison of the kits.

RESULTS

Kit 4, the FastDNA Spin Kit for Soil, resulted in high DNA concentrations with better quality and yielded the most accurate depiction of the mock community. Comparison of the leaves at two water depths showed no significant differences in community composition.

DISCUSSION

The success of Kit 4 may be attributed to its use of bead beating with a homogenizer, which was more efficient in the lysis of Gram-positive bacteria than the manual vortexing protocols used by the other kits. Our results show that microbial composition on leaves during early decay remains comparable and may change only in later stages of decomposition.

Host-associated microbe PCR (hamPCR) enables convenient measurement of both microbial load and community composition

The ratio of microbial population size relative to the amount of host tissue, or 'microbial load', is a fundamental metric of colonization and infection, but it cannot be directly deduced from microbial amplicon data such as 16S rRNA gene counts. Because existing methods to determine load, such as serial dilution plating, quantitative PCR, and whole metagenome sequencing add substantial cost and/or experimental burden, they are only rarely paired with amplicon sequencing. We introduce host-associated microbe PCR (hamPCR), a robust strategy to both quantify microbial load and describe interkingdom microbial community composition in a single amplicon library. We demonstrate its accuracy across multiple study systems, including nematodes and major crops, and further present a cost-saving technique to reduce host overrepresentation in the library prior to sequencing. Because hamPCR provides an accessible experimental solution to the well-known limitations and statistical challenges of compositional data, it has far-reaching potential in culture-independent microbiology.

Characterizing Microbiomes via Sequencing of Marker Loci: Techniques To Improve Throughput, Account for Cross-Contamination, and Reduce Cost

New approaches to characterizing microbiomes via high-throughput sequencing provide impressive gains in efficiency and cost reduction compared to approaches that were standard just a few years ago. However, the speed of method development has been such that staying abreast of the latest technological advances is challenging. Moreover, shifting laboratory protocols to include new methods can be expensive and time consuming. To facilitate adoption of new techniques, we provide a guide and review of recent advances that are relevant for single-locus sequence-based study of microbiomes—from extraction to library preparation—including a primer regarding the use of liquid-handling automation in small-scale academic settings. Additionally, we describe several amendments to published techniques to improve throughput, track contamination, and reduce cost. Notably, we suggest adding synthetic DNA molecules to each sample during nucleic acid extraction, thus providing a method of documenting incidences of cross-contamination. We also describe a dual-indexing scheme for Illumina sequencers that allows multiplexing of many thousands of samples with minimal PhiX input. Collectively, the techniques that we describe demonstrate that laboratory technology need not impose strict limitations on the scale of molecular microbial ecology studies.

IMPORTANCE New methods to characterize microbiomes reduce technology-imposed limitations to study design, but many new approaches have not been widely adopted. Here, we present techniques to increase throughput and reduce contamination alongside a thorough review of current best practices.

Evaluation of Established Methods for DNA Extraction and Primer Pairs Targeting 16S rRNA Gene for Bacterial Microbiota Profiling of Olive Xylem Sap

Next-generation sequencing has revolutionized our ability to investigate the microbiota composition of diverse and complex environments. However, a number of factors can affect the accuracy of microbial community assessment, such as the DNA extraction method, the hypervariable region of 16S rRNA gene targeted, or the PCR primers used for amplification. The aim of this study was to assess the influence of commercially available DNA extraction kits and different primer pairs to provide a non-biased vision of the composition of bacterial communities present in olive xylem sap. For that purpose, branches from "Picual" and "Arbequina" olive cultivars were used for xylem sap extraction using a Scholander chamber device. The DNA extraction protocol significantly affected xylem sap bacterial community assessment. That resulted in significant differences in alpha (Richness) and beta diversity (UniFrac distances) metrics among DNA extraction protocols, with the 12 DNA extraction kits evaluated being clustered in four groups behaving differently. Although the core number of taxa detected by all DNA extraction kits included four phyla, seven classes, 12 orders, 16 or 21 families, and 12 or 14 genera when using the Greengenes or Silva database for taxonomic assignment, respectively, some taxa, particularly those identified at low frequency, were detected by some DNA extraction kits only. The most accurate depiction of a bacterial mock community artificially inoculated on sap samples was generated when using the PowerPlant DNA extraction kit, the combination of 799F/1193R primers amplifying the hypervariable V5-V7 region, and the Silva 132 database for taxonomic assignment. The DESeq2 analysis displayed significant differences among genera abundance between the different PCR primer pairs tested. Thus, Enterobacter, Granulicatella, Prevotella, and Brevibacterium presented a significant higher abundance in all PCR protocols when compared with primer pair 799F/1193R, while the opposite was true for Pseudomonas and Pectobacterium. The methodological approach followed in this study can be useful to optimize plant-associated microbiome analysis, especially when exploring new plant niches. Some of the DNA extraction kits and PCR primers selected in this study will contribute to better characterize bacterial communities inhabiting the xylem sap of olives or other woody crop species.

DNA Metabarcoding for the Characterization of Terrestrial Microbiota-Pitfalls and Solutions

Soil-borne microbes are major ecological players in terrestrial environments since they cycle organic matter, channel nutrients across trophic levels and influence plant growth and health. Therefore, the identification, taxonomic characterization and determination of the ecological role of members of soil microbial communities have become major topics of interest. The development and continuous improvement of high-throughput sequencing platforms have further stimulated the study of complex microbiota in soils and plants. The most frequently used approach to study microbiota composition, diversity and dynamics is polymerase chain reaction (PCR), amplifying specific taxonomically informative gene markers with the subsequent sequencing of the amplicons. This methodological approach is called DNA metabarcoding. Over the last decade, DNA metabarcoding has rapidly emerged as a powerful and cost-effective method for the description of microbiota in environmental samples. However, this approach involves several processing steps, each of which might introduce significant biases that can considerably compromise the reliability of the metabarcoding output. The aim of this review is to provide state-of-the-art background knowledge needed to make appropriate decisions at each step of a DNA metabarcoding workflow, highlighting crucial steps that, if considered, ensures an accurate and standardized characterization of microbiota in environmental studies.

Jumping the green wall: The use of PNA-DNA clamps to enhance microbiome sampling depth in wildlife microbiome research

As microbiome research moves away from model organisms to wildlife, new challenges for microbiome high-throughput sequencing arise caused by the variety of wildlife diets. High levels of contamination are commonly observed emanating from the host (mitochondria) or diet (chloroplast). Such high contamination levels affect the overall sequencing depth of wildlife samples thus decreasing statistical power and leading to poor performance in downstream analysis. We developed an amplification protocol utilizing PNA-DNA clamps to maximize the use of resources and to increase the sampling depth of true microbiome sequences in samples with high levels of plastid contamination. We chose two study organisms, a bat (Leptonyteris yerbabuenae) and a bird (Mimus parvulus), both relying on heavy plant-based diets that sometimes lead to traces of plant-based fecal material producing high contamination signals from chloroplasts and mitochondria. On average, our protocol yielded a 13-fold increase in bacterial sequence amplification compared with the standard protocol (Earth Microbiome Protocol) used in wildlife research. For both focal species, we were able to increase significantly the percentage of sequences available for downstream analyses after the filtering of plastids and mitochondria. Our study presents the first results obtained by using PNA-DNA clamps to block the PCR amplification of chloroplast and mitochondrial DNA from the diet in the gut microbiome of wildlife. The method involves a cost-effective molecular technique instead of the filtering out of unwanted sequencing reads. As 33% and 26% of birds and bats, respectively, have a plant-based diet, the tool that we present here will optimize the sequencing and analysis of wild microbiomes.

Heterogeneity of the rice microbial community of the Chinese centuries-old Honghe Hani rice terraces system

The Honghe Hani rice terraces system (HHRTS) is a traditional rice cultivation system where Hani people cultivate remarkably diverse rice varieties. Recent introductions of modern rice varieties to the HHRTS have significantly increased the severity of rice diseases within the terraces. Here, we determine the impacts of these recent introductions on the composition of the rice-associated microbial communities. We confirm that the HHRTS contains a range of both traditional HHRTS landraces and introduced modern rice varieties and find differences between the microbial communities of these two groups. However, this introduction of modern rice varieties has not strongly impacted the overall diversity of the HHRTS rice microbial community. Furthermore, we find that the rice varieties (i.e. groups of closely related genotypes) have significantly structured the rice microbial community composition (accounting for 15%-22% of the variance) and that the core microbial community of HHRTS rice plants represents less than 3.3% of all the microbial taxa identified. Collectively, our study suggests a highly diverse HHRTS rice holobiont (host with its associated microbes) where the diversity of rice hosts mirrors the diversity of their microbial communities. Further studies will be needed to better determine how such changes might impact the sustainability of the HHRTS.

Engineering CRISPR/Cas9 to mitigate abundant host contamination for 16S rRNA gene-based amplicon sequencing

BACKGROUND

High-throughput sequencing of bacterial 16S rRNA gene (16S-seq) is a useful and common method for studying bacterial community structures. However, contamination of the 16S rRNA genes from the mitochondrion and plastid hinders the sensitive bacterial 16S-seq in plant microbiota profiling, especially for some plant species such as rice. To date, efficiently mitigating such host contamination without a bias is challenging in 16S rRNA gene-based amplicon sequencing.

RESULTS

We developed Cas-16S-seq method to reduce abundant host contamination for plant microbiota profiling. This method utilizes the Cas9 nuclease and specific guide RNA (gRNA) to cut 16S rRNA targets during library construction, thereby removing host contamination in 16S-seq. We used rice as an example to validate the feasibility and effectiveness of Cas-16S-seq. We established a bioinformatics pipeline to design gRNAs that specifically target rice 16S rRNA genes without bacterial 16S rRNA off-targets. We compared the effectiveness of Cas-16S-seq with that of the commonly used 16S-seq method for artificially mixed 16S rRNA gene communities, paddy soil, rice root, and phyllosphere samples. The results showed that Cas-16S-seq substantially reduces the fraction of rice 16S rRNA gene sequences from 63.2 to 2.9% in root samples and from 99.4 to 11.6% in phyllosphere samples on average. Consequently, Cas-16S-seq detected more bacterial species than the 16S-seq in plant samples. Importantly, when analyzing soil samples, Cas-16S-seq and 16S-seq showed almost identical bacterial communities, suggesting that Cas-16S-seq with host-specific gRNAs that we designed has no off-target in rice microbiota profiling.

CONCLUSION

Our Cas-16S-seq can efficiently remove abundant host contamination without a bias for 16S rRNA gene-based amplicon sequencing, thereby enabling deeper bacterial community profiling with a low cost and high flexibility. Thus, we anticipate that this method would be a useful tool for plant microbiomics. Video Abstract.

Effects of Earthworms on the Microbiomes and Antibiotic Resistomes of Detritus Fauna and Phyllospheres

Our understanding of the influences of animals on microbial communities and antibiotic resistance genes (ARGs) throughout the food web is limited, particularly in the soil ecosystem. We performed a microcosm experiment using the combination of 16S rRNA gene high-throughput sequencing and high-throughput qPCR to study the effects of earthworms on the microbiomes and resistomes of the phyllospheres and soil detritus fauna. A full factorial design was employed in this study, including two earthworm treatments (control and earthworm addition) and two soil treatments (control and manure amendment). Earthworms changed the bacterial composition of the soil detritus fauna and the phyllosphere by modifying the shared habitat. Earthworms also reduced the number and abundance of ARGs in the soil fauna and phyllosphere microbiomes by changing their microbial communities, suggesting that earthworms may be able to provide a sustainable and natural solution to address the ARG crisis in the soil ecosystem. This study provides a new understanding on the effects of soil animals on microbial communities, driving shifts in the resistome of other animals and the phyllosphere, by altering the shared soil habitat.

Evidence for Co-evolutionary History of Early Diverging Lycopodiaceae Plants With Fungi

Lycopods are tracheophytes in the Kingdom Plantae and represent one of the oldest lineages of living vascular plants. Symbiotic interactions between these plants with fungi and bacteria, including fine root endophytes in Endogonales, have been hypothesized to have helped early diverging plant lineages colonize land. However, attempts to study the lycopod rhizobiome in its natural environment are still limited. In this study, we used Illumina amplicon sequencing to characterize fungal and bacterial diversity in nine Lycopodiaceae (club moss) species collected in New Zealand. This was done with generic fungal ITS rDNA primers, as well as Endogonales- and arbuscular mycorrhizal fungi (AMF)-selective primer sets targeting the 18S rDNA, and generic bacterial primers targeting the V4 region of the 16S rDNA. We found that the Lycopodiaceae rhizobiome was comprised of an unexpected high frequency of Basidiomycota and Ascomycota coincident with a low abundance of Endogonales and Glomerales. The distribution and abundance of Endogonales varied with host lycopod, and included a novel taxon as well as a single operational taxonomic unit (OTU) that was detected across all plant species. The Lycopodiaceae species with the greatest number and also most unique OTUs was Phlegmariurus varius, while the plant species that shared the most fungal OTUs were Lycopodiella fastigiatum and Lycopodium scariosum. The bacterial OTU distribution was generally not consistent with fungal OTU distribution. For example, community dissimilarity analysis revealed strong concordance between the evolutionary histories of host plants with the fungal community but not with the bacterial community, indicating that Lycopodiaceae have evolved specific relationships with their fungal symbionts. Notably, nearly 16% of the ITS rDNA fungal diversity detected in the Lycopodiaceae rhizobiome remained poorly classified, indicating there is much plant-associated fungal diversity left to describe in New Zealand.

Kelp beds and their local effects on seawater chemistry, productivity, and microbial communities

Kelp forests are known as key habitats for species diversity and macroalgal productivity; however, we know little about how these biogenic habitats interact with seawater chemistry and phototroph productivity in the water column. We examined kelp forest functions at three locales along the Olympic Peninsula of Washington state by quantifying carbonate chemistry, nutrient concentrations, phytoplankton productivity, and seawater microbial communities inside and outside of kelp beds dominated by the canopy kelp species Nereocystis luetkeana and Macrocystis pyrifera. Kelp beds locally increased the pH, oxygen, and aragonite saturation state of the seawater, but lowered seawater inorganic carbon content and total alkalinity. Although kelp beds depleted nitrate and phosphorus concentrations, ammonium and dissolved organic carbon (DOC) concentrations were enhanced. Kelp beds also decreased chlorophyll concentrations and carbon fixed by phytoplankton, although kelp carbon fixation more than compensated for any difference in phytoplankton production. Kelp beds entrained distinct microbial communities, with higher taxonomic and phylogenetic diversity compared to seawater outside of the kelp bed. Kelp forests thus had significant effects on seawater chemistry, productivity and the microbial assemblages in their proximity. Thereby, the diversity of pathways for carbon and nitrogen cycling was also enhanced. Overall, these observations suggest that the contribution of kelp forests to nearshore carbon and nitrogen cycling is greater than previously documented.

Interdomain ecological networks between plants and microbes

While macroscopic interkingdom relationships have been intensively investigated in various ecosystems, the above-belowground ecology in natural ecosystems has been poorly understood, especially for the plant-microbe associations at a regional scale. In this study, we proposed a workflow to construct interdomain ecological networks (IDEN) between multiple plants and various microbes (bacteria and archaea in this study). Across 30 latitudinal forests in China, the regional IDEN showed particular topological features, including high connectance, nested structure, asymmetric specialization and modularity. Also, plant species exhibited strong preference to specific microbial groups, and the observed network was significantly different from randomly rewired networks. Network module analysis indicated that a majority of microbes associated with plants within modules rather than across modules, suggesting specialized associations between plants and microorganisms. Consistent plant-microbe associations were captured via IDENs constructed within individual forest locations, which reinforced the validity of IDEN analysis. In addition, the plant-forest link distribution showed the geographical distribution of plants had higher endemicity than that of microorganisms. With cautious experimental design and data processing, this study shows interdomain species associations between plants and microbes in natural forest ecosystems and provides new insights into our understanding of meta-communities across different domain species.

Aedes spp. and Their Microbiota: A Review

Aedes spp. are a major public health concern due to their ability to be efficient vectors of dengue, Chikungunya, Zika, and other arboviruses. With limited vaccines available and no effective therapeutic treatments against arboviruses, the control of Aedes spp. populations is currently the only strategy to prevent disease transmission. Host-associated microbes (i.e., microbiota) recently emerged as a promising field to be explored for novel environmentally friendly vector control strategies. In particular, gut microbiota is revealing its impact on multiple aspects of Aedes spp. biology, including vector competence, thus being a promising target for manipulation. Here we describe the technological advances, which are currently expanding our understanding of microbiota composition, abundance, variability, and function in the two main arboviral vectors, the mosquitoes Aedes aegypti and Aedes albopictus. Aedes spp. microbiota is described in light of its tight connections with the environment, with which mosquitoes interact during their various developmental stages. Unraveling the dynamic interactions among the ecology of the habitat, the mosquito and the microbiota have the potential to uncover novel physiological interdependencies and provide a novel perspective for mosquito control.

The Origin, Succession, and Predicted Metabolism of Bacterial Communities Associated with Leaf Decomposition

Community ecologists have traditionally treated individuals within a species as uniform, with individual-level biodiversity rarely considered as a regulator of community and ecosystem function. In our study system, we have documented clear evidence of within-species variation causing local ecosystem adaptation to fluxes across ecosystem boundaries. In this striking pattern of a “home-field advantage,” leaves from individual trees tend to decompose most rapidly when immediately adjacent to their parent tree. Here, we merge community ecology experiments with microbiome approaches to describe how bacterial communities adjust to within-species variation in leaves over spatial scales of less than a kilometer. The results show that bacterial community compositional changes facilitate rapid ecosystem responses to environmental change, effectively maintaining high rates of carbon and nutrient cycling through ecosystems.

Intraspecific variation in plant nutrient and defensive traits can regulate ecosystem-level processes, such as decomposition and transformation of plant carbon and nutrients. Understanding the regulatory mechanisms of ecosystem functions at local scales may facilitate predictions of the resistance and resilience of these functions to change. We evaluated how riverine bacterial community assembly and predicted gene content corresponded to decomposition rates of green leaf inputs from red alder trees into rivers of Washington State, USA. Previously, we documented accelerated decomposition rates for leaves originating from trees growing adjacent to the site of decomposition versus more distant locales, suggesting that microbes have a “home-field advantage” in decomposing local leaves. Here, we identified repeatable stages of bacterial succession, each defined by dominant taxa with predicted gene content associated with metabolic pathways relevant to the leaf characteristics and course of decomposition. “Home” leaves contained bacterial communities with distinct functional capacities to degrade aromatic compounds. Given known spatial variation of alder aromatics, this finding helps explain locally accelerated decomposition. Bacterial decomposer communities adjust to intraspecific variation in leaves at spatial scales of less than a kilometer, providing a mechanism for rapid response to changes in resources such as range shifts among plant genotypes. Such rapid responses among bacterial communities in turn may maintain high rates of carbon and nutrient cycling through aquatic ecosystems.

Successional Dynamics and Seascape-Level Patterns of Microbial Communities on the Canopy-Forming Kelps Nereocystis luetkeana and Macrocystis pyrifera

Canopy-forming kelps create underwater forests that are among the most productive marine ecosystems. On the Pacific coast of North America, two canopy-forming kelps with contrasting life histories co-occur; Macrocystis pyrifera, a perennial species, and Nereocystis luetkeana, an annual species. Kelp blade-associated microbes were sampled from 12 locations across a spatial gradient in Washington, United States, from the outer Pacific Coast to Puget Sound. Microbial communities were characterized using next-generation Illumina sequencing of 16S rRNA genes. At higher taxonomic levels (bacterial phylum and class), canopy-forming kelps hosted remarkably similar microbial communities, but at the amplicon sequence variant level, microbial communities on M. pyrifera and N. luetkeana were host-specific and distinct from free-living bacteria in the surrounding seawater. Microbial communities associated with blades of each kelp species displayed significant geographic variation. The microbiome of N. luetkeana changed along the spatial gradient and was significantly correlated to salinity, with outer Pacific coast sites enriched in Bacteroidetes (family Saprospiraceae) and Gammaproteobacteria (Granulosicoccus sp.), and southern Puget Sound sites enriched in Alphaproteobacteria (family Hyphomonadaceae). We also examined microbial community development and succession on meristematic and apical N. luetkeana blade tissues throughout the summer growing season on Tatoosh Island, WA. Across all dates, microbial communities were less diverse on younger, meristematic blade tissue compared to the older, apical tissues. In addition, phylogenetic relatedness among microbial taxa increased from meristematic to apical blade tissues, suggesting that the addition of microbial taxa to the community was a non-random process that selected for certain phylogenetic groups of microbes. Microbial communities on older, apical tissues displayed significant temporal variation throughout the summer and microbial taxa that were differentially abundant over time displayed clear patterns of community succession. Overall, we report that host species identity, geographic location, and blade tissue age shape the microbial communities on canopy-forming kelps.

Chloroplast sequence variation and the efficacy of peptide nucleic acids for blocking host amplification in plant microbiome studies

BACKGROUND

The ability to efficiently characterize microbial communities from host individuals can be limited by co-amplification of host organellar sequences (mitochondrial and/or plastid), which share a common ancestor and thus sequence similarity with extant bacterial lineages. One promising approach is the use of sequence-specific peptide nucleic acid (PNA) clamps, which bind to, and block amplification of, host-derived DNA. Universal PNA clamps have been proposed to block host plant-derived mitochondrial (mPNA) and plastid (pPNA) sequences at the V4 16S rRNA locus, but their efficacy across a wide range of host plant species has not been experimentally tested.

RESULTS

Using the universal PNA clamps, we amplified and sequenced root microbial communities from replicate individuals of 32 plant species with a most recent common ancestor inferred at 140 MYA. We found the average rate of host plastid contamination across plant species was 23%, however, particular lineages exhibited much higher rates (62-94%), with the highest levels of contamination occurring in the Asteraceae. We investigated chloroplast sequence variation at the V4 locus across 500 land plant species (Embryophyta) and found six lineages with mismatches between plastid and the universal pPNA sequence, including all species within the Asteraceae. Using a modified pPNA for the Asteraceae sequence, we found (1) host contamination in Asteraceae species was reduced from 65 to 23%; and (2) host contamination in non-Asteraceae species was increased from 12 to 69%. These results demonstrate that even single nucleotide mismatches can lead to drastic reductions in pPNA efficacy in blocking host amplification. Importantly, we found that pPNA type (universal or modified) had no effect on the detection of individual bacterial taxa, or estimates of within and between sample bacterial diversity, suggesting that our modification did not introduce bias against particular bacterial lineages.

CONCLUSIONS

When high similarity exists between host organellar DNA and PCR target sequences, PNA clamps are an important molecular tool to reduce host contamination during amplification. Here, we provide a validated framework to modify universal PNA clamps to accommodate host variation in organellar sequences.