Genomic and Environmental Factors Shape the Active Gill Bacterial Community of an Amazonian Teleost Holobiont

Fish bacterial communities provide functions critical for their host’s survival in contrasting environments. These communities are sensitive to environmental-specific factors (i.e., physicochemical parameters, bacterioplankton), and host-specific factors (i.e., host genetic background). The relative contribution of these factors shaping Amazonian fish bacterial communities is largely unknown. Here, we investigated this topic by analyzing the gill bacterial communities of 240 wild flag cichlids (Mesonauta festivus) from 4 different populations (genetic clusters) distributed across 12 sites in 2 contrasting water types (ion-poor/ acidic black water and ion-rich/circumneutral white water). Transcriptionally active gill bacterial communities were characterized by a 16S rRNA metabarcoding approach carried on RNA extractions. They were analyzed using comprehensive data sets from the hosts genetic background (Genotyping-By-Sequencing), the bacterioplankton (16S rRNA) and a set of 34 environmental parameters. Results show that the taxonomic structure of 16S rRNA gene transcripts libraries were significantly different between the 4 genetic clusters and also between the 2 water types. However, results suggest that the contribution of the host’s genetic background was relatively weak in comparison to the environment-related factors in structuring the relative abundance of different active gill bacteria species. This finding was also confirmed by a mixed-effects modeling analysis, which indicated that the dissimilarity between the taxonomic structure of bacterioplanktonic communities possessed the best explicative power regarding the dissimilarity between gill bacterial communities’ structure, while pairwise fixation indexes (FST) from the hosts’ genetic data only had a weak explicative power. We discuss these results in terms of bacterial community assembly processes and flag cichlid fish ecology.