An updated LSU database and pipeline for environmental DNA identification of arbuscular mycorrhizal fungi.

Recent work established a backbone reference tree and phylogenetic placement pipeline for identification of arbuscular mycorrhizal fungal (AMF) large subunit (LSU) rDNA environmental sequences. Our previously published pipeline allowed any environmental sequence to be identified as putative AMF or within one of the major families. Despite this contribution, difficulties in implementation of the pipeline remain. Here, we present an updated database and pipeline with (1) an expanded backbone tree to include four newly described genera and (2) several changes to improve ease and consistency of implementation. In particular, packages required for the pipeline are now installed as a single folder (conda environment) and the pipeline has been tested across three university computing clusters. This updated backbone tree and pipeline will enable broadened adoption by the community, advancing our understanding of these ubiquitous and ecologically important fungi.

Environmental identification of arbuscular mycorrhizal fungi using the LSU rDNA gene region: an expanded database and improved pipeline.

Arbuscular mycorrhizal fungi (AMF; Glomeromycota) are difficult to culture; therefore, establishing a robust amplicon-based approach to taxa identification is imperative to describe AMF diversity. Further, due to low and biased sampling of AMF taxa, molecular databases do not represent the breadth of AMF diversity, making database matching approaches suboptimal. Therefore, a full description of AMF diversity requires a tool to determine sequence-based placement in the Glomeromycota clade. Nonetheless, commonly used gene regions, including the SSU and ITS, do not enable reliable phylogenetic placement. Here, we present an improved database and pipeline for the phylogenetic determination of AMF using amplicons from the large subunit (LSU) rRNA gene. We improve our database and backbone tree by including additional outgroup sequences. We also improve an existing bioinformatics pipeline by aligning forward and reverse reads separately, using a universal alignment for all tree building, and implementing a BLAST screening prior to tree building to remove non-homologous sequences. Finally, we present a script to extract AMF belonging to 11 major families as well as an amplicon sequencing variant (ASV) version of our pipeline. We test the utility of the pipeline by testing the placement of known AMF, known non-AMF, and Acaulospora sp. spore sequences. This work represents the most comprehensive database and pipeline for phylogenetic placement of AMF LSU amplicon sequences within the Glomeromycota clade.