Estimating belowground plant abundance with DNA metabarcoding

Most work on plant community ecology has been performed above ground, neglecting the processes that occur in the soil. DNA metabarcoding, in which multiple species are computationally identified in bulk samples, can help to overcome the logistical limitations involved in sampling plant communities belowground. However, a major limitation of this methodology is the quantification of species' abundances based on the percentage of sequences assigned to each taxon. Using root tissues of five dominant species in a semi-arid Mediterranean shrubland (Bupleurum fruticescens, Helianthemum cinereum, Linum suffruticosum, Stipa pennata and Thymus vulgaris), we built pairwise mixtures of relative abundance (20%, 50% and 80% biomass), and implemented two methods (linear model fits and correction indices) to improve estimates of root biomass. We validated both methods with multispecies mixtures that simulate field-collected samples. For all species, we found a positive and highly significant relationship between the percentage of sequences and biomass in the mixtures (R-2 = .44-.66), but the equations for each species (slope and intercept) differed among them, and two species were consistently over- and under-estimated. The correction indices greatly improved the estimates of biomass percentage for all five species in the multispecies mixtures, and reduced the overall error from 17% to 6%. Our results show that, through the use of post-sequencing quantification methods on mock communities, DNA metabarcoding can be effectively used to determine not only species' presence but also their relative abundance in field samples of root mixtures. Importantly, knowledge of these aspects will allow us to study key, yet poorly understood, belowground processes.