Plant species richness belowground: higher richness and new patterns revealed by next-generation sequencing

Environmental DNA: The next chapter

Molecular tools are an indispensable part of ecology and biodiversity sciences and implemented across all biomes. About a decade ago, the use and implementation of environmental DNA (eDNA) to detect biodiversity signals extracted from environmental samples opened new avenues of research. Initial eDNA research focused on understanding population dynamics of target species. Its scope thereafter broadened, uncovering previously unrecorded biodiversity via metabarcoding in both well-studied and understudied ecosystems across all taxonomic groups. The application of eDNA rapidly became an established part of biodiversity research, and a research field by its own. Here, we revisit key expectations made in a land-mark special issue on eDNA in Molecular Ecology in 2012 to frame the development in six key areas: (1) sample collection, (2) primer development, (3) biomonitoring, (4) quantification, (5) behaviour of DNA in the environment and (6) reference database development. We pinpoint the success of eDNA, yet also discuss shortfalls and expectations not met, highlighting areas of research priority and identify the unexpected developments. In parallel, our retrospective couples a screening of the peer-reviewed literature with a survey of eDNA users including academics, end-users and commercial providers, in which we address the priority areas to focus research efforts to advance the field of eDNA. With the rapid and ever-increasing pace of new technical advances, the future of eDNA looks bright, yet successful applications and best practices must become more interdisciplinary to reach its full potential. Our retrospect gives the tools and expectations towards concretely moving the field forward.

Niche types and community assembly

Studies of niche differentiation and biodiversity often focus on a few niche dimensions due to the methodological challenge of describing hyperdimensional niche space. However, this may limit our understanding of community assembly processes. We used the full spectrum of realized niche types to study arbuscular mycorrhizal fungal communities: distinguishing abiotic and biotic, and condition and resource, axes. Estimates of differentiation in relation to different niche types were only moderately correlated. However, coexisting taxon niches were consistently less differentiated than expected, based on a regional null model, indicating the importance of habitat filtering at that scale. Nonetheless, resource niches were relatively more differentiated than condition niches, which is consistent with the effect of a resource niche-based coexistence mechanism. Considering niche types, and in particular distinguishing resource and condition niches, provides a more complete understanding of community assembly, compared with studying individual niche axes or the full niche.

Soil community composition in dynamic stages of semi-natural calcareous grassland

European dry thin-soil calcareous grasslands (alvars) are species-rich semi-natural habitats. Cessation of traditional management, such as mowing and grazing, leads to shrub and tree encroachment and the local extinction of characteristic alvar species. While soil microbes are known to play a critical role in driving vegetation and ecosystem dynamics, more information is needed about their composition and function in grasslands of different dynamic stages. Here we assess the composition of soil fungal, prokaryotic, and plant communities using soil environmental DNA from restored alvar grasslands in Estonia. The study areas included grasslands that had experienced different degrees of woody encroachment prior to restoration (woody plant removal and grazing), as well as unmanaged open grasslands. We found that, in general, different taxonomic groups exhibited correlated patterns of between-community variation. Previous forest sites, which had prior to restoration experienced a high degree of woody encroachment by ectomycorrhizal Scots pine, were compositionally most distinct from managed open grasslands, which had little woody vegetation even prior to restoration. The functional structure of plant and fungal communities varied in ways that were consistent with the representation of mycorrhizal types in the ecosystems prior to restoration. Compositional differences between managed and unmanaged open grasslands reflecting the implementation of grazing without further management interventions were clearer among fungal, and to an extent prokaryotic, communities than among plant communities. While previous studies have shown that during woody encroachment of alvar grassland, plant communities change first and fungal communities follow, our DNA-based results suggest that microbial communities reacted faster than plant communities during the restoration of grazing management in alvar grassland. We conclude that while the plant community responds faster to cessation of management, the fungal community responds faster to restoration of management. This may indicate hysteresis, where the eventual pathway back to the original state (grazed ecosystem) differs from the pathway taken towards the alternative state (abandoned semi-natural grassland ecosystem).

Ecological drivers of fine-scale distribution of arbuscular mycorrhizal fungi in a semiarid Mediterranean scrubland

BACKGROUND AND AIMS

Arbuscular mycorrhizal (AM) fungi enhance the uptake of water and minerals by the plant hosts, alleviating plant stress. Therefore, AM fungal-plant interactions are particularly important in drylands and other stressful ecosystems. We aimed to determine the combined and independent effects of above- and below-ground plant community attributes (i.e. diversity and composition), soil heterogeneity and spatial covariates on the spatial structure of the AM fungal communities in a semiarid Mediterranean scrubland. Furthermore, we evaluated how the phylogenetic relatedness of both plants and AM fungi shapes these symbiotic relationships.

METHODS

We characterized the composition and diversity of AM fungal and plant communities in a dry Mediterranean scrubland taxonomically and phylogenetically, using DNA metabarcoding and a spatially explicit sampling design at the plant neighbourhood scale.

KEY RESULTS

The above- and below-ground plant community attributes, soil physicochemical properties and spatial variables explained unique fractions of AM fungal diversity and composition. Mainly, variations in plant composition affected the AM fungal composition and diversity. Our results also showed that particular AM fungal taxa tended to be associated with closely related plant species, suggesting the existence of a phylogenetic signal. Although soil texture, fertility and pH affected AM fungal community assembly, spatial factors had a greater influence on AM fungal community composition and diversity than soil physicochemical properties.

CONCLUSIONS

Our results highlight that the more easily accessible above-ground vegetation is a reliable indicator of the linkages between plant roots and AM fungi. We also emphasize the importance of soil physicochemical properties in addition to below-ground plant information, while accounting for the phylogenetic relationships of both plants and fungi, because these factors improve our ability to predict the relationships between AM fungal and plant communities.

Mapping the root systems of individual trees in a natural community using genotyping-by-sequencing

The architecture of root systems is an important driver of plant fitness, competition and ecosystem processes. However, the methodological difficulty of mapping roots hampers the study of these processes. Existing approaches to match individual plants to belowground samples are low throughput and species specific. Here, we developed a scalable sequencing-based method to map the root systems of individual trees across multiple species. We successfully applied it to a tropical dry forest community in the Brazilian Caatinga containing 14 species. We sequenced all 42 individual shrubs and trees in a 14 × 14 m plot using double-digest restriction site-associated sequencing (ddRADseq). We identified species-specific markers and individual-specific haplotypes from the data. We matched these markers to the ddRADseq data from 100 mixed root samples from across the centre (10 × 10 m) of the plot at four different depths using a newly developed R package. We identified individual root samples for all species and all but one individual. There was a strong significant correlation between belowground and aboveground size measurements, and we also detected significant species-level root-depth preference for two species. The method is more scalable and less labour intensive than the current techniques and is broadly applicable to ecology, forestry and agricultural biology.

What Hides in the Heights? The Case of the Iberian Endemism Bromus picoeuropeanus

Bromus picoeuropeanus is a recently described species belonging to a complex genus of grasses. It inhabits stony soils at heights ranging from 1600 to 2200 m in Picos de Europa (Cantabrian Mountains, northern Spain). This species is morphologically very similar to B. erectus, partially sharing its presumed distribution range. We aim to determine the relationship between these species and their altitudinal ranges in Picos de Europa and the Cantabrian Mountains by conducting phylogenetic analyses based on nuclear (ETS and ITS) and chloroplastic (trnL) markers. Phylogenetic trees were inferred by Maximum Likelihood and Bayesian Inference. Haplotype networks were estimated based on the plastid marker. Although the ITS topologies could not generate exclusive clades for these species, the ETS analyses generated highly supported B. picoeuropeanus exclusive clades, which included locations outside its altitudinal putative range. The ETS-ITS and ETS-ITS-trnL topologies generated B. picoeuropeanus exclusive clades, whereas the trnL-based trees and haplotype networks were unable to discriminate B. erectus and B. picoeuropeanus. This evidence suggests that B. picoeuropeanus is a separate species with a larger distribution than previously thought, opening new questions regarding the evolution of B. erectus and other similar species in European mountainous systems. However, more information is needed regarding B. picoeuropeanus susceptibility to temperature rises.

Metabarcoding of soil environmental DNA to estimate plant diversity globally

Introduction

Traditional approaches to collecting large-scale biodiversity data pose huge logistical and technical challenges. We aimed to assess how a comparatively simple method based on sequencing environmental DNA (eDNA) characterises global variation in plant diversity and community composition compared with data derived from traditional plant inventory methods.

Methods

We sequenced a short fragment (P6 loop) of the chloroplast trnL intron from from 325 globally distributed soil samples and compared estimates of diversity and composition with those derived from traditional sources based on empirical (GBIF) or extrapolated plant distribution and diversity data.

Results

Large-scale plant diversity and community composition patterns revealed by sequencing eDNA were broadly in accordance with those derived from traditional sources. The success of the eDNA taxonomy assignment, and the overlap of taxon lists between eDNA and GBIF, was greatest at moderate to high latitudes of the northern hemisphere. On average, around half (mean: 51.5% SD 17.6) of local GBIF records were represented in eDNA databases at the species level, depending on the geographic region.

Discussion

eDNA trnL gene sequencing data accurately represent global patterns in plant diversity and composition and thus can provide a basis for large-scale vegetation studies. Important experimental considerations for plant eDNA studies include using a sampling volume and design to maximise the number of taxa detected and optimising the sequencing depth. However, increasing the coverage of reference sequence databases would yield the most significant improvements in the accuracy of taxonomic assignments made using the P6 loop of the trnL region.

Root hemiparasitic plants are associated with more even communities across North America

Root hemiparasitic plants both compete with and extract resources from host plants. By reducing the abundance of dominant plants and releasing subordinates from competitive exclusion, they can have an outsized impact on plant communities. Most research on the ecological role of hemiparasites is manipulative and focuses on a small number of hemiparasitic taxa. Here, we ask whether patterns in natural plant communities match the expectation that hemiparasites affect the structure of plant communities. Our data were collected on 129 national park units spanning the continental United States. The most common hemiparasite genera were Pedicularis, Castilleja, Krameria, and Comandra. We used null models and linear mixed models to determine whether hemiparasites were associated with changes in community richness and evenness. Hemiparasite presence did not affect community metrics. Hemiparasite abundance was positively associated with increasing evenness of herbaceous species, but not with species richness. The associations that we observed on a continental scale are consistent with evidence that the impacts of root hemiparasitic plants on evenness can be substantial and abundance dependent but that effects on richness are less pronounced. Hemiparasites mediate competitive exclusion in communities to facilitate species coexistence and merit consideration of inclusion in ecological theories of coexistence.

A Critical Assessment of the Congruency between Environmental DNA and Palaeoecology for the Biodiversity Monitoring and Palaeoenvironmental Reconstruction

The present study suggests that standardized methodology, careful site selection, and stratigraphy are essential for investigating ancient ecosystems in order to evaluate biodiversity and DNA-based time series. Based on specific keywords, this investigation reviewed 146 publications using the SCOPUS, Web of Science (WoS), PUBMED, and Google Scholar databases. Results indicate that environmental deoxyribose nucleic acid (eDNA) can be pivotal for assessing and conserving ecosystems. Our review revealed that in the last 12 years (January 2008–July 2021), 63% of the studies based on eDNA have been reported from aquatic ecosystems, 25% from marine habitats, and 12% from terrestrial environments. Out of studies conducted in aquatic systems using the environmental DNA (eDNA) technique, 63% of the investigations have been reported from freshwater ecosystems, with an utmost focus on fish diversity (40%). Further analysis of the literature reveals that during the same period, 24% of the investigations using the environmental DNA technique were carried out on invertebrates, 8% on mammals, 7% on plants, 6% on reptiles, and 5% on birds. The results obtained clearly indicate that the environmental DNA technique has a clear-cut edge over other biodiversity monitoring methods. Furthermore, we also found that eDNA, in conjunction with different dating techniques, can provide better insight into deciphering eco-evolutionary feedback. Therefore, an attempt has been made to offer extensive information on the application of dating methods for different taxa present in diverse ecosystems. Last, we provide suggestions and elucidations on how to overcome the caveats and delineate some of the research avenues that will likely shape this field in the near future. This paper aims to identify the gaps in environmental DNA (eDNA) investigations to help researchers, ecologists, and decision-makers to develop a holistic understanding of environmental DNA (eDNA) and its utility as a palaeoenvironmental contrivance.

Plastid Genome of Equisetum xylochaetum from the Atacama Desert, Chile and the Relationships of Equisetum Based on Frequently Used Plastid Genes and Network Analysis

The modern pteridophyte genus Equisetum is the only survivor of Sphenopsida, an ancient clade known from the Devonian. This genus, of nearly worldwide distribution, comprises approximately 15 extant species. However, genomic information is limited. In this study, we assembled the complete chloroplast genome of the giant species Equisetum xylochaetum from a metagenomic sequence and compared the plastid genome structure and protein-coding regions with information available for two other Equisetum species using network analysis. Equisetum chloroplast genomes showed conserved traits of quadripartite structure, gene content, and gene order. Phylogenetic analysis based on plastome protein-coding regions corroborated previous reports that Equisetum is monophyletic, and that E. xylochaetum is more closely related to E. hyemale than to E. arvense. Single-gene phylogenetic estimation and haplotype analysis showed that E. xylochaetum belonged to the subgenus Hippochaete. Single-gene haplotype analysis revealed that E. arvense, E. hyemale, E. myriochaetum, and E. variegatum resolved more than one haplotype per species, suggesting the presence of a high diversity or a high mutation rate of the corresponding nucleotide sequence. Sequences from E. bogotense appeared as a distinct group of haplotypes representing the subgenus Paramochaete that diverged from Hippochaete and Equisetum. In addition, the taxa that were frequently located at the joint region of the map were E. scirpoides and E. pratense, suggesting the presence of some plastome characters among the Equiseum subgenera.

A starting guide to root ecology: strengthening ecological concepts and standardising root classification, sampling, processing and trait measurements

In the context of a recent massive increase in research on plant root functions and their impact on the environment, root ecologists currently face many important challenges to keep on generating cutting-edge, meaningful and integrated knowledge. Consideration of the below-ground components in plant and ecosystem studies has been consistently called for in recent decades, but methodology is disparate and sometimes inappropriate. This handbook, based on the collective effort of a large team of experts, will improve trait comparisons across studies and integration of information across databases by providing standardised methods and controlled vocabularies. It is meant to be used not only as starting point by students and scientists who desire working on below-ground ecosystems, but also by experts for consolidating and broadening their views on multiple aspects of root ecology. Beyond the classical compilation of measurement protocols, we have synthesised recommendations from the literature to provide key background knowledge useful for: (1) defining below-ground plant entities and giving keys for their meaningful dissection, classification and naming beyond the classical fine-root vs coarse-root approach; (2) considering the specificity of root research to produce sound laboratory and field data; (3) describing typical, but overlooked steps for studying roots (e.g. root handling, cleaning and storage); and (4) gathering metadata necessary for the interpretation of results and their reuse. Most importantly, all root traits have been introduced with some degree of ecological context that will be a foundation for understanding their ecological meaning, their typical use and uncertainties, and some methodological and conceptual perspectives for future research. Considering all of this, we urge readers not to solely extract protocol recommendations for trait measurements from this work, but to take a moment to read and reflect on the extensive information contained in this broader guide to root ecology, including sections I-VII and the many introductions to each section and root trait description. Finally, it is critical to understand that a major aim of this guide is to help break down barriers between the many subdisciplines of root ecology and ecophysiology, broaden researchers' views on the multiple aspects of root study and create favourable conditions for the inception of comprehensive experiments on the role of roots in plant and ecosystem functioning.

Exploring plant diversity through soil DNA in Thai national parks for influencing land reform and agriculture planning

Background

The severe deforestation, as indicated in national forest data, is a recurring problem in many areas of Northern Thailand, including Doi Suthep-Pui National Park. Agricultural expansion in these areas, is one of the major drivers of deforestation, having adverse consequences on local plant biodiversity. Conserving biodiversity is mainly dependent on the biological monitoring of species distribution and population sizes. However, the existing conventional approaches for monitoring biodiversity are rather limited.

Methods

Here, we explored soil DNA at four forest types in Doi Suthep-Pui National Park in Northern Thailand. Three soil samples, composed of different soil cores mixed together, per sampling location were collected. Soil biodiversity was investigated through eDNA metabarcoding analysis using primers targeting the P6 loop of the plastid DNA trnL (UAA) intron.

Results

The distribution of taxa for each sample was found to be similar between replicates. A strong congruence between the conventional morphology- and eDNA-based data of plant diversity in the studied areas was observed. All species recorded by conventional survey with DNA data deposited in the GenBank were detected through the eDNA analysis. Moreover, traces of crops, such as lettuce, maize, wheat and soybean, which were not expected and were not visually detected in the forest area, were identified. It is noteworthy that neighboring land and areas in the studied National Park were once used for crop cultivation, and even to date there is still agricultural land within a 5-10 km radius from the forest sites where the soil samples were collected. The presence of cultivated area near the forest may suggest that we are now facing agricultural intensification leading to deforestation. Land reform for agriculture usage necessitates coordinated planning in order to preserve the forest area. In that context, the eDNA-based data would be useful for influencing policies and management towards this goal.

Looking for Hidden Enemies of Metabarcoding: Species Composition, Habitat and Management Can Strongly Influence DNA Extraction while Examining Grassland Communities

Despite the raising preoccupation, the critical question of how the plant community is composed belowground still remains unresolved, particularly for the conservation priority types of vegetation. The usefulness of metabarcoding analysis of the belowground parts of the plant community is subjected to a considerable bias, that often impedes detection of all species in a sample due to insufficient DNA quality or quantity. In the presented study we have attempted to find environmental factors that determine the amount and quality of DNA extracted from total plant tissue from above- and belowground samples (1000 and 10,000 cm2). We analyzed the influence of land use intensity, soil properties, species composition, and season on DNA extraction. The most important factors for DNA quality were vegetation type, soil conductometry (EC), and soil pH for the belowground samples. The species that significantly decreased the DNA quality were Calamagrostis epigejos, Coronilla varia, and Holcus lanatus. For the aboveground part of the vegetation, the season, management intensity, and certain species-with the most prominent being Centaurea rhenana and Cirsium canum-have the highest influence. Additionally, we found that sample size, soil granulation, MgO, organic C, K2O, and total soil N content are important for DNA extraction effectiveness. Both low EC and pH reduce significantly the yield and quality of DNA. Identifying the potential inhibitors of DNA isolation and predicting difficulties of sampling the vegetation plots for metabarcoding analysis will help to optimize the universal, low-cost multi-stage DNA extraction procedure in molecular ecology studies.

User-friendly bioinformatics pipeline gDAT (graphical downstream analysis tool) for analysing rDNA sequences

High-throughput sequencing (HTS) of multiple organisms in parallel (metabarcoding) has become a routine and cost-effective method for the analysis of microbial communities in environmental samples. However, careful data treatment is required to identify potential errors in HTS data, and the large volume of data generated by HTS requires in-house experience with command line tools for downstream analysis. This paper introduces a pipeline that incorporates the most common command line tools into an easy-to-use graphical interface-gDAT. By using the Python scripting language, the pipeline is compatible with the latest Windows, macOS and Linux operating systems. The pipeline supports analysis of Sanger, 454, IonTorrent, Illumina and PacBio sequences, allows custom modification of quality filtering steps, and implements both open and closed-reference operational taxonomic unit-picking for sequence identification. Predefined parameters are optimized for analysis of small subunit (SSU) rRNA gene amplicons from arbuscular mycorrhizal fungi, but the pipeline is widely applicable to metabarcoding studies targeting a broad range of organisms. The pipeline was additionally tested with data using general eukaryotic primers from the SSU gene region and fungal primers from the internal transcribed spacer (ITS) marker region. We describe the pipeline design and evaluate its performance and speed by conducting analysis of example data sets using different marker regions sequenced on Illumina platforms. The graphical interface, with the option to use the command line if needed, provides an accessible tool for rapid data analysis with repeatability and logging capabilities. Keeping the software open-source maximizes code accessibility, allowing scrutiny and bug fixes by the community.

Plant Taxonomy: A Historical Perspective, Current Challenges, and Perspectives

Taxonomy is the science that explores, describes, names, and classifies all organisms. In this introductory chapter, we highlight the major historical steps in the elaboration of this science, which provides baseline data for all fields of biology and plays a vital role for society but is also an independent, complex, and sound hypothesis-driven scientific discipline.In a first part, we underline that plant taxonomy is one of the earliest scientific disciplines that emerged thousands of years ago, even before the important contributions of the Greeks and Romans (e.g., Theophrastus, Pliny the Elder, and Dioscorides). In the fifteenth-sixteenth centuries, plant taxonomy benefited from the Great Navigations, the invention of the printing press, the creation of botanic gardens, and the use of the drying technique to preserve plant specimens. In parallel with the growing body of morpho-anatomical data, subsequent major steps in the history of plant taxonomy include the emergence of the concept of natural classification , the adoption of the binomial naming system (with the major role of Linnaeus) and other universal rules for the naming of plants, the formulation of the principle of subordination of characters, and the advent of the evolutionary thought. More recently, the cladistic theory (initiated by Hennig) and the rapid advances in DNA technologies allowed to infer phylogenies and to propose true natural, genealogy-based classifications.In a second part, we put the emphasis on the challenges that plant taxonomy faces nowadays. The still very incomplete taxonomic knowledge of the worldwide flora (the so-called taxonomic impediment) is seriously hampering conservation efforts that are especially crucial as biodiversity has entered its sixth extinction crisis. It appears mainly due to insufficient funding, lack of taxonomic expertise, and lack of communication and coordination. We then review recent initiatives to overcome these limitations and to anticipate how taxonomy should and could evolve. In particular, the use of molecular data has been era-splitting for taxonomy and may allow an accelerated pace of species discovery. We examine both strengths and limitations of such techniques in comparison to morphology-based investigations, we give broad recommendations on the use of molecular tools for plant taxonomy, and we highlight the need for an integrative taxonomy based on evidence from multiple sources.

msGBS: A new high-throughput approach to quantify the relative species abundance in root samples of multispecies plant communities

Plant interactions are as important belowground as aboveground. Belowground plant interactions are however inherently difficult to quantify, as roots of different species are difficult to disentangle. Although for a couple of decades molecular techniques have been successfully applied to quantify root abundance, root identification and quantification in multispecies plant communities remains particularly challenging. Here we present a novel methodology, multispecies genotyping by sequencing (msGBS), as a next step to tackle this challenge. First, a multispecies meta‐reference database containing thousands of gDNA clusters per species is created from GBS derived High Throughput Sequencing (HTS) reads. Second, GBS derived HTS reads from multispecies root samples are mapped to this meta‐reference which, after a filter procedure to increase the taxonomic resolution, allows the parallel quantification of multiple species. The msGBS signal of 111 mock‐mixture root samples, with up to 8 plant species per sample, was used to calculate the within‐species abundance. Optional subsequent calibration yielded the across‐species abundance. The within‐ and across‐species abundances highly correlated (R² range 0.72–0.94 and 0.85–0.98, respectively) to the biomass‐based species abundance. Compared to a qPCR based method which was previously used to analyse the same set of samples, msGBS provided similar results. Additional data on 11 congener species groups within 105 natural field root samples showed high taxonomic resolution of the method. msGBS is highly scalable in terms of sensitivity and species numbers within samples, which is a major advantage compared to the qPCR method and advances our tools to reveal hidden belowground interactions.

see also the Perspective by Josep Piñol

The role of root community attributes in predicting soil fungal and bacterial community patterns

Roots are assumed to play a major role in structuring soil microbial communities, but most studies exploring the relationships between microbes and plants at the community level have only used aboveground plant distribution as a proxy. However, a decoupling between belowground and aboveground plant components may occur due to differential spreading of plant canopies and root systems. Thus, soil microbe-plant links are not completely understood. Using a combination of DNA metabarcoding and spatially explicit sampling at the plant neighbourhood scale, we assessed the influence of the plant root community on soil bacterial and fungal diversity (species richness, composition and β-diversity) in a dry Mediterranean scrubland. We found that root composition and biomass, but not richness, predict unique fractions of variation in microbial richness and composition. Moreover, bacterial β-diversity was related to root β-diversity, while fungal β-diversity was related to aboveground plant β-diversity, suggesting that plants differently influence both microbial groups. Our study highlights the role of plant distribution both belowground and aboveground, soil properties and other spatially structured factors in explaining the heterogeneity in soil microbial diversity. These results also show that incorporating data on both plant community compartments will further our understanding of the relationships between soil microbial and plant communities.

Plant evenness modulates the effect of plant richness on soil bacterial diversity

Understanding the relationships between aboveground and belowground biodiversity will help to expand our knowledge on how ecological communities and processes are interactively determined, and thus provide new perspectives for the conservation of biodiversity. Despite the theoretical analyses generally predicting a positive relationship between plant richness and soil microbial diversity, the results from empirical studies have been mixed, probably due to the effect of plant evenness. To investigate this relationship, we conducted field experiments in two geographically distinct sites (Linhai and Shenmu, >1400km apart), by simultaneously manipulating plant richness (2, 4, and 8 species) and evenness (homogeneous versus non-homogeneous). After one year, the bacterial response to plant richness with different plant evenness levels was evaluated using terminal restriction fragment length polymorphism (T-RFLP) analysis. Our results showed that plant evenness modulated plant richness effects on bacterial community, as reflected by the more pronounced positive correlations between bacterial richness and plant richness in homogeneous plant community than in the non-homogeneous treatment. Additionally, plant community structure significantly affected bacterial communities only in the homogeneous treatment in Shenmu, but not in the non-homogeneous treatments. Our results demonstrate that plant evenness could regulate plant richness effects on bacterial alpha- and beta-diversity and thus provide valuable insights into the association between aboveground and belowground biodiversity.

Structural equation modeling of a winnowed soil microbiome identifies how invasive plants re-structure microbial networks

The development of microbial networks is central to ecosystem functioning and is the hallmark of complex natural systems. Characterizing network development over time and across environmental gradients is hindered by the millions of potential interactions among community members, limiting interpretations of network evolution. We developed a feature selection approach using data winnowing that identifies the most ecologically influential microorganisms within a network undergoing change. Using a combination of graph theory, leave-one-out analysis, and statistical inference, complex microbial communities are winnowed to identify the core organisms responding to external gradients or functionality, and then network development is evaluated against these externalities. In a plant invasion case study, the winnowed microbial network became more influential as the plant invasion progressed as a result of direct plant-microbe links rather than the expected indirect plant-soil-microbe links. This represents the first use of structural equation modeling to predict microbial network evolution, which requires identification of keystone taxa and quantification of the ecological processes underpinning community structure and function patterns.

Belowground plant parts are crucial for comprehensively estimating total plant richness in herbaceous and woody habitats

Most studies consider aboveground plant species richness as a representative biodiversity measure. This approach inevitably assumes that the partitioning of total plant species richness into above- and belowground components is constant or at least consistent within and across vegetation types. However, with studies considering belowground plant richness still scarce and completely absent along vegetation gradients, this assumption lacks experimental support. Novel DNA sequencing techniques allow economical, high-throughput species identification of belowground environmental samples, enabling the measurement of the contributions of both above- and belowground plant components to total plant richness. We investigated above- and belowground plant species richness in four vegetation types (birch forest, heath, low alpine tundra, high alpine tundra) at the scale of herbaceous plant neighborhoods (dm) using 454 sequencing of the chloroplast trnL (UAA) intron to determine the plant species richness of environmental root samples and combined it with aboveground data from vegetation surveys to obtain total plant species richness. We correlated the measured plant species richness components with each other and with their respective plant biomass components within and across vegetation types. Total plant species richness exceeded aboveground richness twice on average and by as much as three times in low alpine tundra, indicating that a significant fraction of belowground plant richness cannot be recorded aboveground. More importantly, no consistent relationship among richness components (above- and belowground) was found within or across vegetation types, indicating that aboveground richness alone cannot predict total plant richness in contrasting vegetation types. Finally, no consistent relationship between plant richness and the corresponding biomass component was found. Our results clearly show that aboveground plant richness alone is a poor estimator of total plant species richness within and across different vegetation types. Consequently, it is crucial to account for belowground plant richness in future plant ecological studies in order to validate currently accepted plant richness patterns, as well as to measure potential changes in plant community composition in a changing environment.

Multiscale patterns and drivers of arbuscular mycorrhizal fungal communities in the roots and root-associated soil of a wild perennial herb

Arbuscular mycorrhizal (AM) fungi form diverse communities and are known to influence above-ground community dynamics and biodiversity. However, the multiscale patterns and drivers of AM fungal composition and diversity are still poorly understood. We sequenced DNA markers from roots and root-associated soil from Plantago lanceolata plants collected across multiple spatial scales to allow comparison of AM fungal communities among neighbouring plants, plant subpopulations, nearby plant populations, and regions. We also measured soil nutrients, temperature, humidity, and community composition of neighbouring plants and nonAM root-associated fungi. AM fungal communities were already highly dissimilar among neighbouring plants (c. 30 cm apart), albeit with a high variation in the degree of similarity at this small spatial scale. AM fungal communities were increasingly, and more consistently, dissimilar at larger spatial scales. Spatial structure and environmental drivers explained a similar percentage of the variation, from 7% to 25%. A large fraction of the variation remained unexplained, which may be a result of unmeasured environmental variables, species interactions and stochastic processes. We conclude that AM fungal communities are highly variable among nearby plants. AM fungi may therefore play a major role in maintaining small-scale variation in community dynamics and biodiversity.

Predicting the structure of soil communities from plant community taxonomy, phylogeny, and traits

There are numerous ways in which plants can influence the composition of soil communities. However, it remains unclear whether information on plant community attributes, including taxonomic, phylogenetic, or trait-based composition, can be used to predict the structure of soil communities. We tested, in both monocultures and field-grown mixed temperate grassland communities, whether plant attributes predict soil communities including taxonomic groups from across the tree of life (fungi, bacteria, protists, and metazoa). The composition of all soil community groups was affected by plant species identity, both in monocultures and in mixed communities. Moreover, plant community composition predicted additional variation in soil community composition beyond what could be predicted from soil abiotic characteristics. In addition, analysis of the field aboveground plant community composition and the composition of plant roots suggests that plant community attributes are better predictors of soil communities than root distributions. However, neither plant phylogeny nor plant traits were strong predictors of soil communities in either experiment. Our results demonstrate that grassland plant species form specific associations with soil community members and that information on plant species distributions can improve predictions of soil community composition. These results indicate that specific associations between plant species and complex soil communities are key determinants of biodiversity patterns in grassland soils.

Plant DNA metabarcoding of lake sediments: How does it represent the contemporary vegetation

Metabarcoding of lake sediments have been shown to reveal current and past biodiversity, but little is known about the degree to which taxa growing in the vegetation are represented in environmental DNA (eDNA) records. We analysed composition of lake and catchment vegetation and vascular plant eDNA at 11 lakes in northern Norway. Out of 489 records of taxa growing within 2 m from the lake shore, 17–49% (mean 31%) of the identifiable taxa recorded were detected with eDNA. Of the 217 eDNA records of 47 plant taxa in the 11 lakes, 73% and 12% matched taxa recorded in vegetation surveys within 2 m and up to about 50 m away from the lakeshore, respectively, whereas 16% were not recorded in the vegetation surveys of the same lake. The latter include taxa likely overlooked in the vegetation surveys or growing outside the survey area. The percentages detected were 61, 47, 25, and 15 for dominant, common, scattered, and rare taxa, respectively. Similar numbers for aquatic plants were 88, 88, 33 and 62%, respectively. Detection rate and taxonomic resolution varied among plant families and functional groups with good detection of e.g. Ericaceae, Roseaceae, deciduous trees, ferns, club mosses and aquatics. The representation of terrestrial taxa in eDNA depends on both their distance from the sampling site and their abundance and is sufficient for recording vegetation types. For aquatic vegetation, eDNA may be comparable with, or even superior to, in-lake vegetation surveys and may therefore be used as an tool for biomonitoring. For reconstruction of terrestrial vegetation, technical improvements and more intensive sampling is needed to detect a higher proportion of rare taxa although DNA of some taxa may never reach the lake sediments due to taphonomical constrains. Nevertheless, eDNA performs similar to conventional methods of pollen and macrofossil analyses and may therefore be an important tool for reconstruction of past vegetation.

Archaea and bacteria mediate the effects of native species root loss on fungi during plant invasion

Although invasive plants can drive ecosystem change, little is known about the directional nature of belowground interactions between invasive plants, native roots, bacteria, archaea and fungi. We used detailed bioinformatics and a recently developed root assay on soils collected in fescue grassland along a gradient of smooth brome (Bromus inermis Leyss) invasion to examine the links between smooth brome shoot litter and root, archaea, bacteria and fungal communities. We examined (1) aboveground versus belowground influences of smooth brome on soil microbial communities, (2) the importance of direct versus microbe-mediated impacts of plants on soil fungal communities, and (3) the web of roots, shoots, archaea, bacteria and fungi interactions across the A and B soil horizons in invaded and non-invaded sites. Archaea and bacteria influenced fungal composition, but not vice versa, as indicated by redundancy analyses. Co-inertia analyses suggested that bacterial-fungal variance was driven primarily by 12 bacterial operational taxonomic units (OTUs). Brome increased bacterial diversity via smooth brome litter in the A horizon and roots in the B horizon, which then reduced fungal diversity. Archaea increased abundance of several bacterial OTUs, and the key bacterial OTUs mediated changes in the fungi's response to invasion. Overall, native root diversity loss and bacterial mediation were more important drivers of fungal composition than were the direct effects of increases in smooth brome. Critically, native plant species displacement and root loss appeared to be the most important driver of fungal composition during invasion. This causal web likely gives rise to the plant-fungi feedbacks, which are an essential factor determining plant diversity in invaded grassland ecosystems.

The best of both worlds: A combined approach for analyzing microalgal diversity via metabarcoding and morphology-based methods

An increasing number of studies use next generation sequencing (NGS) to analyze complex communities, but is the method sensitive enough when it comes to identification and quantification of species? We compared NGS with morphology-based identification methods in an analysis of microalgal (periphyton) communities. We conducted a mesocosm experiment in which we allowed two benthic grazer species to feed upon benthic biofilms, which resulted in altered periphyton communities. Morphology-based identification and 454 (Roche) pyrosequencing of the V4 region in the small ribosomal unit (18S) rDNA gene were used to investigate the community change caused by grazing. Both the NGS-based data and the morphology-based method detected a marked shift in the biofilm composition, though the two methods varied strongly in their abilities to detect and quantify specific taxa, and neither method was able to detect all species in the biofilms. For quantitative analysis, we therefore recommend using both metabarcoding and microscopic identification when assessing the community composition of eukaryotic microorganisms.

Large-Scale Monitoring of Plants through Environmental DNA Metabarcoding of Soil: Recovery, Resolution, and Annotation of Four DNA Markers

In a rapidly changing world we need methods to efficiently assess biodiversity in order to monitor ecosystem trends. Ecological monitoring often uses plant community composition to infer quality of sites but conventional aboveground surveys only capture a snapshot of the actively growing plant diversity. Environmental DNA (eDNA) extracted from soil samples, however, can include taxa represented by both active and dormant tissues, seeds, pollen, and detritus. Analysis of this eDNA through DNA metabarcoding provides a more comprehensive view of plant diversity at a site from a single assessment but it is not clear which DNA markers are best used to capture this diversity. Sequence recovery, annotation, and sequence resolution among taxa were evaluated for four established DNA markers (matK, rbcL, ITS2, and the trnL P6 loop) in silico using database sequences and in situ using high throughput sequencing of 35 soil samples from a remote boreal wetland. Overall, ITS2 and rbcL are recommended for DNA metabarcoding of vascular plants from eDNA when not using customized or geographically restricted reference databases. We describe a new framework for evaluating DNA metabarcodes and, contrary to existing assumptions, we found that full length DNA barcode regions could outperform shorter markers for surveying plant diversity from soil samples. By using current DNA barcoding markers rbcL and ITS2 for plant metabarcoding, we can take advantage of existing resources such as the growing DNA barcode database. Our work establishes the value of standard DNA barcodes for soil plant eDNA analysis in ecological investigations and biomonitoring programs and supports the collaborative development of DNA barcoding and metabarcoding.

Evaluating a multigene environmental DNA approach for biodiversity assessment

Background

There is an increasing demand for rapid biodiversity assessment tools that have a broad taxonomic coverage. Here we evaluate a suite of environmental DNA (eDNA) markers coupled with next generation sequencing (NGS) that span the tree of life, comparing them with traditional biodiversity monitoring tools within ten 20×20 meter plots along a 700 meter elevational gradient.

Results

From six eDNA datasets (one from each of 16S, 18S, ITS, trnL and two from COI) we identified sequences from 109 NCBI taxonomy-defined phyla or equivalent, ranging from 31 to 60 for a given eDNA marker. Estimates of alpha and gamma diversity were sensitive to the number of sequence reads, whereas beta diversity estimates were less sensitive. The average within-plot beta diversity was lower than between plots for all markers. The soil beta diversity of COI and 18S markers showed the strongest response to the elevational variation of the eDNA markers (COI: r=0.49, p<0.001; 18S: r=0.48, p<0.001). Furthermore pairwise beta diversities for these two markers were strongly correlated with those calculated from traditional vegetation and invertebrate biodiversity measures.

Conclusions

Using a soil-based eDNA approach, we demonstrate that standard phylogenetic markers are capable of recovering sequences from a broad diversity of eukaryotes, in addition to prokaryotes by 16S. The COI and 18S eDNA markers are the best proxies for aboveground biodiversity based on the high correlation between the pairwise beta diversities of these markers and those obtained using traditional methods.

Electronic supplementary material

The online version of this article (doi:10.1186/s13742-015-0086-1) contains supplementary material, which is available to authorized users.

Fine-scale belowground species associations in temperate grassland

Evaluating how belowground processes contribute to plant community dynamics is hampered by limited information on the spatial structure of root communities at the scale that plants interact belowground. In this study, roots were mapped to the nearest one mm and molecularly identified by species on vertical (0-15 cm deep) surfaces of soil blocks excavated from dry and mesic grasslands in Yellowstone National Park (YNP) to examine the spatial relationships among species at the scale that roots interact. Our results indicated that average interspecific root - root distances for the majority of species were within a distance (3 mm) that roots have been shown to compete for resources. Most species placed their roots at random, although low root numbers for many species probably led to overestimating the occurrence of random patterns. According to theory, we expected that most of the remaining species would segregate their root systems to avoid competition. Instead we found that more species aggregated than segregated from others. Based on previous investigations, we hypothesize that species aggregate to increase uptake of water, nitrogen and/or phosphorus made available by neighbouring roots, or as a consequence of a reduction in the pathogenicity of soil biota growing in multispecies mixtures. Our results indicate that YNP grassland root communities are organized as closely interdigitating networks of species that potentially can support strong interactions among many species combinations. Future root research should address the prevalence and functional consequences of species aggregation across plant communities.

Spatial heterogeneity of plant-soil feedback affects root interactions and interspecific competition

Plant-soil feedback is receiving increasing interest as a factor influencing plant competition and species coexistence in grasslands. However, we do not know how spatial distribution of plant-soil feedback affects plant below-ground interactions. We investigated the way in which spatial heterogeneity of soil biota affects competitive interactions in grassland plant species. We performed a pairwise competition experiment combined with heterogeneous distribution of soil biota using four grassland plant species and their soil biota. Patches were applied as quadrants of 'own' and 'foreign' soils from all plant species in all pairwise combinations. To evaluate interspecific root responses, species-specific root biomass was quantified using real-time PCR. All plant species suffered negative soil feedback, but strength was species-specific, reflected by a decrease in root growth in own compared with foreign soil. Reduction in root growth in own patches by the superior plant competitor provided opportunities for inferior competitors to increase root biomass in these patches. These patterns did not cascade into above-ground effects during our experiment. We show that root distributions can be determined by spatial heterogeneity of soil biota, affecting plant below-ground competitive interactions. Thus, spatial heterogeneity of soil biota may contribute to plant species coexistence in species-rich grasslands.

Biases underlying species detection using fluorescent amplified-fragment length polymorphisms yielded from roots

BACKGROUND

Roots of different plant species are typically morphologically indistinguishable. Of the DNA-based techniques, fluorescent amplified-fragment length polymorphisms (FAFLPs) are considered reliable, high throughput, inexpensive methods to identify roots from mixed species samples. False-negatives, however, are not uncommon and their underlying causes are poorly understood. We investigated several sources of potential biases originating in DNA extraction and amplification. Specifically, we examined the effects of sample storage, tissue, and species on DNA yield and purity, and the effects of DNA concentration and fragment size on amplification of three non-coding chloroplast regions (trnT-trnL intergenic spacer, trnL intron, and trnL-trnF intergenic spacer).

RESULTS

We found that sample condition, tissue and species all affected DNA yield. A single freeze-thaw reduces DNA yield, DNA yield is less for roots than shoots, and species vary in the amount of DNA yielded from extractions. The effects of template DNA concentration, species identity, and their interaction on amplicon yield differed across the three chloroplast regions tested. We found that the effect of species identity on amplicon production was generally more pronounced than that of DNA concentration. Though these factors influenced DNA yield, they likely do not have a pronounced effect on detection success of fragments and only underscore the restriction on the use of FAFLPs for measuring species presence rather than their abundance. However, for two of the regions tested-the trnT-trnL intergenic spacer and the trnL intron-size-based fragment competition occurred and the likelihood of detection was higher for smaller than larger fragments. This result reveals a methodological bias when using FAFLPs.

CONCLUSIONS

To avoid potential bias with the use of FAFLPs, we recommend users check for the disproportionate absence of species detected belowground versus aboveground as a function of fragment size, and explore other regions, aside from the trnT-trnL intergenic spacer and trnL intron, for amplification.

Aggregated and complementary: symmetric proliferation, overyielding, and mass effects explain fine-root biomass in soil patches in a diverse temperate deciduous forest landscape

Few studies describe root distributions at the species level in diverse forests, although belowground species interactions and traits are often assumed to affect fine-root biomass (FRB). We used molecular barcoding to study how FRB of trees relates to soil characteristics, species identity, root diversity, and root traits, and how these relationships are affected by proximity to ecotones in a temperate forest landscape. We found that soil patch root biomass increased in response to soil resources across all species, and there was little belowground vertical or horizontal spatial segregation among species. Root traits and species relative abundance did not explain significant variation in FRB after correcting for soil fertility. A positive relationship between phylogenetic diversity and FRB indicated significant belowground overyielding attributable to local root diversity. Finally, variation in FRB explained by soil fertility and diversity was reduced near ecotones, but only because of a reduction in biomass in periodically anoxic areas. These results suggest that symmetric responses to soil properties are coupled with complementary species traits and interactions to explain variation in FRB among soil patches. In addition, landscape-level dispersal among habitats and across ecotones helps explain variation in the strength of these relationships in complex landscapes.

A molecular method to identify species of fine roots and to predict the proportion of a species in mixed samples in subtropical forests

Understanding of belowground interactions among tree species and the fine root (≤2 mm in diameter) contribution of a species to forest ecosystem production are mostly restricted by experimental difficulties in the quantification of the species composition. The available approaches have various defects. By contrast, DNA-based methods can avoid these drawbacks. Quantitative real-time polymerase chain reaction (PCR) is an advanced molecular technology, but it is difficult to develop specific primer sets. The method of next-generation sequencing has several limitations, such as inaccurate sequencing of homopolymer regions, as well as being time-consuming, and requiring special knowledge for data analysis. This study evaluated the potential of the DNA-sequence-based method to identify tree species and to quantify the relative proportion of each species in mixed fine root samples. We discriminated the species by isolating DNA from individual fine roots and amplifying the plastid trnL(UAA; i.e., tRNA-Leu-UAA) intron using the PCR. To estimate relative proportions, we extracted DNA from fine root mixtures. After the plastid trnL(UAA) intron amplification and TA-cloning, we sequenced the positive clones from each mixture. Our results indicated that the plastid trnL(UAA) intron spacer successfully distinguished tree species of fine roots in subtropical forests. In addition, the DNA-sequence-based approach could reliably estimate the relative proportion of each species in mixed fine root samples. To our knowledge, this is the first time that the DNA-sequence-based method has been used to quantify tree species proportions in mixed fine root samples in Chinese subtropical forests. As the cost of DNA-sequencing declines and DNA-sequence-based methods improve, the molecular method will be more widely used to determine fine root species and abundance.

Universal and blocking primer mismatches limit the use of high-throughput DNA sequencing for the quantitative metabarcoding of arthropods

The quantification of the biological diversity in environmental samples using high-throughput DNA sequencing is hindered by the PCR bias caused by variable primer-template mismatches of the individual species. In some dietary studies, there is the added problem that samples are enriched with predator DNA, so often a predator-specific blocking oligonucleotide is used to alleviate the problem. However, specific blocking oligonucleotides could coblock nontarget species to some degree. Here, we accurately estimate the extent of the PCR biases induced by universal and blocking primers on a mock community prepared with DNA of twelve species of terrestrial arthropods. We also compare universal and blocking primer biases with those induced by variable annealing temperature and number of PCR cycles. The results show that reads of all species were recovered after PCR enrichment at our control conditions (no blocking oligonucleotide, 45 °C annealing temperature and 40 cycles) and high-throughput sequencing. They also show that the four factors considered biased the final proportions of the species to some degree. Among these factors, the number of primer-template mismatches of each species had a disproportionate effect (up to five orders of magnitude) on the amplification efficiency. In particular, the number of primer-template mismatches explained most of the variation (~3/4) in the amplification efficiency of the species. The effect of blocking oligonucleotide concentration on nontarget species relative abundance was also significant, but less important (below one order of magnitude). Considering the results reported here, the quantitative potential of the technique is limited, and only qualitative results (the species list) are reliable, at least when targeting the barcoding COI region.

Going underground: root traits as drivers of ecosystem processes

Ecologists are increasingly adopting trait-based approaches to understand how community change influences ecosystem processes. However, most of this research has focussed on aboveground plant traits, whereas it is becoming clear that root traits are important drivers of many ecosystem processes, such as carbon (C) and nutrient cycling, and the formation and structural stability of soil. Here, we synthesise emerging evidence that illustrates how root traits impact ecosystem processes, and propose a pathway to unravel the complex roles of root traits in driving ecosystem processes and their response to global change. Finally, we identify research challenges and novel technologies to address them.

A molecular identification protocol for roots of boreal forest tree species

PREMISE OF THE STUDY

Roots play a key role in many ecological processes, yet our ability to identify species from bulk root samples is limited. Molecular tools may be used to identify species from root samples, but they have not yet been developed for most systems. Here we present a PCR-based method previously used to identify roots of grassland species, modified for use in boreal forests. •

METHODS

We used repeatable interspecific size differences in fluorescent amplified fragment length polymorphisms of three noncoding chloroplast DNA regions to identify seven woody species common to boreal forests in Alberta, Canada. •

RESULTS

Abies balsamea, Alnus crispa, Betula papyrifera, Pinus contorta, and Populus tremuloides were identifiable to species, while Picea glauca and Picea mariana were identifiable to genus. In mixtures of known composition of foliar DNA, species were identified with 98% accuracy using one region. Mixed root samples of unknown composition were identified with 100% accuracy; four species were identified using one region, while three species were identified using two regions. •

DISCUSSION

This methodology is accurate, efficient, and inexpensive, and thus a valuable approach for ecological studies of roots. Furthermore, this method has now been validated for both grassland and boreal forest systems, and thus may also have applications in any plant community.

More closely related plants have more distinct mycorrhizal communities

Neighbouring plants are known to vary from having similar to dissimilar arbuscular mycorrhizal fungal (AMF) communities. One possibility is that closely related plants have more similar AMF communities than more distantly related plants, an indication of phylogenetic host specificity. Here, we investigated the structure of AMF communities among dominant grassland plants at three sites in the Northern Great Plains to test whether the pairwise phylogenetic distance among plant species was correlated with pairwise AMF community dissimilarity. For eight dominant and co-occurring grassland plant species, we reconstructed a phylogeny based on DNA data and characterized the AMF communities of their roots at each site. Community analyses revealed that AMF communities varied among sites and among plant species. Contrary to expectations for phylogenetic host specificity, we found that within a site more closely related plants had more distinct AMF communities despite their having similar phenologies. Associations with unique AMF communities may enhance the functional complementarity of related species and promote their coexistence.

Next generation restoration genetics: applications and opportunities

Restoration ecology is a young scientific discipline underpinning improvements in the rapid global expansion of ecological restoration. The application of molecular tools over the past 20 years has made an important contribution to understanding genetic factors influencing ecological restoration success. Here we illustrate how recent advances in next generation sequencing (NGS) methods are revolutionising the practical contribution of genetics to restoration. Novel applications include a dramatically enhanced capacity to measure adaptive variation for optimal seed sourcing, high-throughput assessment and monitoring of natural and restored biological communities aboveground and belowground, and gene expression analysis as a measure of genetic resilience of restored populations. Challenges remain in data generation, handling and analysis, and how best to apply NGS for practical outcomes in restoration.

Species richness of arbuscular mycorrhizal fungi: associations with grassland plant richness and biomass

Although experiments show a positive association between vascular plant and arbuscular mycorrhizal fungal (AMF) species richness, evidence from natural ecosystems is scarce. Furthermore, there is little knowledge about how AMF richness varies with belowground plant richness and biomass. We examined relationships among AMF richness, above- and belowground plant richness, and plant root and shoot biomass in a native North American grassland. Root-colonizing AMF richness and belowground plant richness were detected from the same bulk root samples by 454-sequencing of the AMF SSU rRNA and plant trnL genes. In total we detected 63 AMF taxa. Plant richness was 1.5 times greater belowground than aboveground. AMF richness was significantly positively correlated with plant species richness, and more strongly with below- than aboveground plant richness. Belowground plant richness was positively correlated with belowground plant biomass and total plant biomass, whereas aboveground plant richness was positively correlated only with belowground plant biomass. By contrast, AMF richness was negatively correlated with belowground and total plant biomass. Our results indicate that AMF richness and plant belowground richness are more strongly related with each other and with plant community biomass than with the plant aboveground richness measures that have been almost exclusively considered to date.

Environmental metabarcodes for insects: in silico PCR reveals potential for taxonomic bias

Studies of insect assemblages are suited to the simultaneous DNA-based identification of multiple taxa known as metabarcoding. To obtain accurate estimates of diversity, metabarcoding markers ideally possess appropriate taxonomic coverage to avoid PCR-amplification bias, as well as sufficient sequence divergence to resolve species. We used in silico PCR to compare the taxonomic coverage and resolution of newly designed insect metabarcodes (targeting 16S) with that of existing markers [16S and cytochrome oxidase c subunit I (COI)] and then compared their efficiency in vitro. Existing metabarcoding primers amplified in silico <75% of insect species with complete mitochondrial genomes available, whereas new primers targeting 16S provided >90% coverage. Furthermore, metabarcodes targeting COI appeared to introduce taxonomic PCR-amplification bias, typically amplifying a greater percentage of Lepidoptera and Diptera species, while failing to amplify certain orders in silico. To test whether bias predicted in silico was observed in vitro, we created an artificial DNA blend containing equal amounts of DNA from 14 species, representing 11 insect orders and one arachnid. We PCR-amplified the blend using five primer sets, targeting either COI or 16S, with high-throughput amplicon sequencing yielding more than 6 million reads. In vitro results typically corresponded to in silico PCR predictions, with newly designed 16S primers detecting 11 insect taxa present, thus providing equivalent or better taxonomic coverage than COI metabarcodes. Our results demonstrate that in silico PCR is a useful tool for predicting taxonomic bias in mixed template PCR and that researchers should be wary of potential bias when selecting metabarcoding markers.