Sequencing summary by 16S rRNA high-throughput sequencing
A total of 1,146,486 tags were obtained. After quality control, a total of 1,119,154 effective tags were retrieved (Table S1), and these samples were assigned to a total of 3310 OTUs. The rarefaction analyses indicated that the diversity of all samples was well represented by the number of sequences analyzed (Figure S1). The bacterial sequences in the female and male microbial communities were clustered into 1335 and 1666 OTUs at a 97% sequence similarity, respectively. The bacterial sequences in the microbial community in the surrounding seawater were clustered into 1453 OTUs with 97% sequence similarity (Figure S2).
Sexual differences in α-diversity and β-diversity
The results of α-diversity and β-diversity showed that there was no significant gender difference (p > 0.05) in the diversity of Porphyra haitanensis (Figs. 1 and 2), but the epiphytic bacterial communities were clustered separately according to gender (Fig. 3).
Community structure of the epiphytic bacterial community
A total of 3310 OTUs were identified in 9 samples (Figure S2), which were distributed in 34 phyla and 84 classes and 410 genera.
In this study, 34 phyla were identified (Fig. 4A). Proteobacteria was the most abundant phylum in most samples. Other dominant taxa included Bacteroides, Actinomycetes, and Cyanobacteria. The abundance of Proteobacteria was higher in seawater (61.59%) than in male Porphyra haitanensis (41.25%) and female Porphyra haitanensis (32.43%).
At the class level, a total of 84 classes were identified. The bacterial taxa with high relative abundance at the class level (Fig. 4B) were Bacteroidia, Alphaproteobacteria, and Gammaproteobacteria. The average relative abundance of Bacteroidia were male Porphyra haitanensis (41.61%), female Porphyra haitanensis (49.82%), and seawater (8.48%); the average relative abundance of Alphaproteobacteria were male Porphyra haitanensis (26.30%), female Porphyra haitanensis (25.96%), and seawater (25.07%); the average relative abundance of Gammaproteobacteria were male Porphyra haitanensis (14.42%), female Porphyra haitanensis (5.87%), and seawater (34.64%).
At the genus level, a total of 410 genera were identified. The bacterial taxa with high relative abundance at the genus level (Fig. 4C) were Tenacibaculum, Maribacter, and Granulosicoccus. The average relative abundance of Tenacibaculum were male Porphyra haitanensis (1.44%), female Porphyra haitanensis (25.78%) and seawater (0.06%); the average relative abundance of Maribacter were male Porphyra haitanensis (14.87%), female Porphyra haitanensis (1.66%), and seawater (0.22%); the average relative abundance of Granulosicoccus were male Porphyra haitanensis (12.20%), female Porphyra haitanensis (3.04%), and seawater (0.15%).
Core microbiota in male and female Porphyra haitanensis
A total of 752 OTUs were common in male and female Porphyra haitanensis, being mapped to 17 phyla and 31 classes and 147 genera. Proteobacteria and Bacteroidetes were the core microbiota in male and female Porphyra haitanensis. Proteobacteria was represented by Granulosicoccus, Hellea, Litorimonas, Marinicella, and Sphingorhabdus, which were present in all male and female Porphyra haitanensis samples. The phylum Bacteroidetes was represented by Maribacter, Lewinella, Euzebyella, Tunicatimonas, and Wenyingzhuangia, which were present in all male and female Porphyra haitanensis samples.
Indicative species of the epiphytic bacterial communities on male and female Porphyra haitanensis
There were 33 and 31 indicative epiphytic bacteria taxa for the male and female Porphyra haitanensis, respectively (Fig. 5) and among them, bacteria from the Maribacter (male 14.87%, female 1.66%) and the Tenacibaculum (male 1.44%, female 25.78%) were the most indicative for male and female Porphyra haitanensis, respectively, according the negative linear discriminant analysis (LDA) score (Figure S3).