Comparative analysis of the endophytic bacteria inhabiting the phyllosphere of aquatic fern Azolla species by high-throughput sequencing

Background Azolla is a small floating fern living in symbiosis with nitrogen-fixing cyanobacteria and provides a variety of important ecosystem benefits. Previous studies have presented that Azolla harbors diverse bacteria that may play a key role in host fitness and productivity. However, the characteristics of endophytic bacteria inhabiting the phyllosphere of different species of Azolla have not yet been fully understood. Results In this study, the 16S ribosomal DNA (rDNA) V5-V7 region of bacteria was determined by Illumina high-throughput sequencing platform to study the diversity and richness of endophytic bacterial communities in the phyllosphere of five Azolla species collected from different countries. A total of 1150 operational taxonomic units (OTUs) were detected for the endophytic bacteria community. According to the alpha diversity indices, the diversity of bacteria was ordered as Azolla imbricata > A. pinnata > A. filiculoides > A. mexicana > A. caroliniana. The PCoA results displayed that the bacterial communities of A. mexicana and A. caroliniana shared the highest similarity, followed by the similarity between A. pinnata and A. imbricata, and they were significantly distinct from the community of A. filiculoides. The dominant bacteria of Azolla mainly belonged to the phylum of Proteobacteria, followed by Actinobacteria, Chlorobillobacteria, and Firmicutes. In detail, the relative abundance of Proteobacteria in A. imbricata was 52.23%, whereas it was more than 80.00% in the other four species of Azolla. Notably, Herbaspirillum (45.91%, 44.08%) and Methylophilus (29.97%, 37.96%) were the main genera inhabiting A. mexicana and A. caroliniana respectively. Ferrovibrio (18.54%) and Rhizobium (16.68%) were the dominant genera inhabiting A. filiculoides. The group of unidentified genera (41.63%, 44.92%) consisted most of the bacteria in A. imbricata and A. pinnata respectively. Further analysis suggested that the significant different bacteria identified in LDA Effect Size analysis existed Azolla species-specific patterns. Conclusions In summary, all results suggested that the diversity and composition of the endophytic bacterial communities were different in Azolla species.