Environmental DNA metabarcoding reveals the influence of human activities on microeukaryotic plankton along the Chinese coastline

Microeukaryotic plankton, with its extremely diverse taxa, is a key component in both the marine food web and biogeochemical cycling. Coastal seas, which are home to the numerous microeukaryotic plankton that underpin the functions of these aquatic ecosystems, are often impacted by human activities. However, understanding the biogeographical patterns of diversity and community structure of microeukaryotic plankton and the role that major shaping factors play at the continent scale is still a challenge in coastal ecology. Here, the biogeographic patterns of biodiversity, community structure, and co-occurrence patterns were investigated by environmental DNA (eDNA) based approaches. Unlike most eDNA studies, we combined several methods (in silico PCR, mock and environmental communities) to systematically evaluate the specificity and coverage of primers to overcome the limitation of marker selection on biodiversity recovery. The 1380F/1510R primer set showed the best performance for the amplification of coastal plankton with the highest coverage, sensitivity, and resolution. We showed a unimodal pattern for planktonic alpha diversity with latitude (P < 0.001), and nutrient-related factors (NO3N, NO2N, and NH4N) were the leading predictors for spatial patterning. Significant regional biogeographic patterns and potential drivers for planktonic communities were found across coastal regions. All communities generally fitted the regional distance-decay relationship (DDR) model with the strongest spatial turnover rate was found in the Yalujiang (YLJ) estuary (P < 0.001). The environmental factors, especially inorganic nitrogen and heavy metals (HMs), had the greatest impact on planktonic community similarity in the Beibu Bay (BB) and East China Sea (ECS). Furthermore, we observed spatial plankton co-occurrence patterns, and the networked topology and structure were strongly driven by potential anthropogenic activity factors (nutrients and HMs). Overall, our study provided a systematic approach for metabarcode primer selection in eDNA-based biodiversity monitoring and revealed that the spatial pattern of the microeukaryotic plankton community was mainly controlled by regional human activity-related factors.