Assessing anaerobic gut fungal (Neocalliamstigomycota) diversity using PacBio D1/D2 LSU rRNA amplicon sequencing and multi-year isolation

Abstract

The anaerobic gut fungi (AGF, Neocallimastigomycota) reside in the alimentary tracts of herbivores where they play a central role in the breakdown of ingested plant material. Accurate assessment of AGF diversity has been hampered by inherent deficiencies of the internal transcribed spacer1 (ITS1) region as a phylogenetic marker. Here, we report on the development and implementation of the D1/D2 region of the large ribosomal subunit (D1/D2 LSU) as a novel marker for assessing AGF diversity in culture-independent surveys. Sequencing a 1.4-1.5 Kbp amplicon encompassing the ITS1-5.8S rRNA-ITS2-D1/D2 LSU region in the ribosomal RNA locus from fungal strains and environmental samples generated a reference D1/D2 LSU database for all cultured AGF genera, as well as the majority of candidate genera encountered in prior ITS1-based diversity surveys. Subsequently, a D1/D2 LSU-based diversity survey using long read PacBio SMRT sequencing technology was conducted on fecal samples from 21 wild and domesticated herbivores. Twenty-eight genera and candidate genera were identified in the 17.7 K sequences obtained, including multiple novel lineages that were predominantly, but not exclusively, identified in wild herbivores. Association between certain AGF genera and animal lifestyles, or animal host family was observed. Finally, to address the current paucity of AGF isolates, concurrent isolation efforts utilizing multiple approaches to maximize recovery yielded 216 isolates belonging to twelve different genera, several of which have no prior cultured-representatives. Our results establish the utility of D1/D2 LSU and PacBio sequencing for AGF diversity surveys, and the culturability of a wide range of AGF taxa, and demonstrate that wild herbivores represent a yet-untapped reservoir of AGF diversity.