Abstract
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.