Yousri Abdelhafiz1
Jorge M. O. Fernandes1
Simone Larger2
Davide Albanese2
Claudio Donati2
Omid Jafari1,3
Artem V. Nedoluzhko1
Viswanath Kiron1*- 1Faculty of Biosciences and Aquaculture, Nord University, Bodø, Norway
- 2Unit of Computational Biology, Fondazione Edmund Mach, San Michele all’Adige, Italy
- 3International Sturgeon Research Institute, Iranian Fisheries Science Research Institute, Agricultural Research, Education and Extension Organization, Rasht, Iran
In industrial animal production, breeding strategies are essential to produce offspring of better quality and vitality. It is also known that host microbiome has a bearing on its health. Here, we report for the first time the influence of crossbreeding strategy, inbreeding or outbreeding, on the buccal and intestinal bacterial communities in female Nile tilapia (Oreochromis niloticus). Crossbreeding was performed within a family and between different fish families to obtain the inbred and outbred study groups, respectively. The genetic relationship and structure analysis revealed significant genetic differentiation between the inbred and outbred groups. We also employed a 16S rRNA gene sequencing technique to understand the significant differences between the diversities of the bacterial communities of the inbred and outbred groups. The core microbiota composition in the mouth and the intestine was not affected by the crossbreeding strategy but their abundance varied between the two groups. Furthermore, opportunistic bacteria were abundant in the buccal cavity and intestine of the outbred group, whereas beneficial bacteria were abundant in the intestine of the inbred group. The present study indicates that crossbreeding can influence the abundance of beneficial bacteria, core microbiome and the inter-individual variation in the microbiome.
Introduction
Animals are bred for food, fibers, transport, protection, company as well as for other purposes such as scientific research (Flint and Woolliams, 2008). Domestication of different animals, mainly livestock species started several years ago and presently crossbreeding programs are essential tools to improve the productivity, efficiency, and sustainability of domesticated animals (Hill, 2014, 2016). Initially, livestock were selected based on desired phenotypic traits. Over the past 50years, there has been a remarkable increase in livestock production due to the improvement in breeding practices and better understanding of genetics. Genetics plays an important role in modern breeding programs, which combine basic breeding concepts and emerging technologies (Schultz et al., 2020).
Crossbreeding of farmed animals and agricultural plants is well-established compared to those of farmed aquatic animals (D’ambrosio et al., 2019; Gratacap et al., 2019). However, the production of fish based on crossbreeding programs is expected to increase as the farming of fish such as Nile tilapia (Oreochromis niloticus) and Atlantic salmon (Salmo salar) is expanding rapidly (Gjedrem et al., 2012; Gjedrem and Rye, 2018; Mehar et al., 2019). Several strategies such as selective breeding have been implemented to increase the production of fast-growing fish species and their disease resistance (Lind et al., 2012; Ina-Salwany et al., 2019). Nevertheless, outbreak of many diseases such as Tenacibaculosis (yellow mouth), Streptococcosis and Vibriosis has led to high mortality in fish farms and the industry has suffered huge economic losses (Jantrakajorn et al., 2014; Ina-Salwany et al., 2019; Wynne et al., 2020). The industry has hardly taken steps to selectively breed fishes in order to shape the microbiota as an indicator of health. It has been reported that selective breeding can produce fishes with microbiota that can be manipulated to improve disease resistance (Piazzon et al., 2020).
Currently, there are many genetically improved tilapia and GIFT (Genetically Improved Farmed Tilapia) is the most known breed. Although many studies have employed genetically improved tilapia (Bolivar and Newkirk, 2002; Romana-Eguia et al., 2005; Santos et al., 2013; Mehar et al., 2019), to our knowledge there are only a couple of reports about the microbiome composition in selectively bred fish (Kokou et al., 2018; Brown et al., 2019). In mouse, selective breeding is known to increase the inter-individual gut microbiota similarity (Pang et al., 2012); variation is less in the case of inbred animals compared to their outbred counterparts (Hufeldt et al., 2010). Researchers have also succeeded in producing outbred mice with stable gut microbiota (Hart et al., 2018). Furthermore, the association between the gut microbiome and breeding was studied in mouse models by analysing the effect of the gut microbiome on different breeds (Pang et al., 2012; Kreisinger et al., 2014; Ericsson et al., 2015; Oriá et al., 2018). This link was also explored in plants by examining the impact of the microorganisms on host phenotype (Wagner et al., 2020). Moreover, the microbial taxa that is widespread among the host population is vertically transmitted, and host factors provide them with the optimum ecosystem for colonization (Risely, 2020).
Selective breeding affects host genetic selection, which in turn shapes the gut microbiome (Kokou et al., 2018) that has an important role in, among others, maintaining the host health. The paucity of information regarding the mating strategy-caused changes in fish microbiome that can signal disease propensity led us to examine the differences in the bacteria associated with inbred and outbred Nile tilapia using next-generation sequencing technology.
Materials and Methods
Fish Husbandry and Sample Collection
Fertilized eggs (n=180) of Nile tilapia, were obtained from wild fish captured from the Nile river, Luxor, Egypt (location GPS: 25°39'56'' N, 32°37'07'' E). These eggs were disinfected with hydrogen peroxide for 10min and placed in egg rockers (Cobalt Aquatics, Rock Hill, South Carolina, United States) installed in a 60L tank with UV treated water, containing 5% NaCl. Around 85% of the eggs were hatched at 28°C within 4days. The hatched larvae were placed in fish transport bags filled with UV treated and 100% oxygen saturated water. These larvae were shipped, within approximately 18h, to the Research Station of Nord University, Bodø, Norway via air and their survival rate exceeded 95%. The transported larvae were reared at a maximum density of 27 fish/m3 for 5months in a freshwater recirculating system. The rearing conditions were: dissolved oxygen – 100%, water temperature –28°C, photoperiod – LD 13:11. The fish were fed Amber Neptun pellets (0.15–0.8mm, Skretting, Stavanger, Norway) during the rearing period. These fish were designated as the F0 generation and were used for the breeding study.
We randomly chose males and females and produced the inbred and outbred groups. When the fish reached 3,570 degree·days, we anesthetized and PIT-tagged them for tracing the individual families.
Prior to sampling, fish were not fed for 48h. They were sacrificed by immersion in an emulsion containing 12ml of clove oil (Sigma Aldrich, St. Louis, Missouri, United States), 96% ethanol (1:10v/v) and 10L of water (Simões et al., 2011; Konstantinidis et al., 2020). Female fish were used for the study as they are maternal mouthbrooders. Twenty fish each from the inbred and outbred groups were used in this study, and three body sites (mouth, anterior and posterior intestine) of female Nile tilapia were targeted for examining the bacterial communities. Mucus samples from the buccal cavity were taken using swabs (Copan Italia, Brescia, Italy), which were transferred to cryotubes and immediately frozen in liquid nitrogen. Then, the same fish were aseptically dissected to collect the anterior and posterior intestine. The intestine samples were also transferred to cryotubes and snap-frozen in liquid nitrogen. The collected samples were stored at −80°C until further use.
DNA Extraction for Whole-Genome Sequencing
DNA was extracted from fast muscle using DNeasy Blood and Tissue Kit based on the guidelines provided by the manufacturer (Qiagen, Hilden, Germany). The Invitrogen Qubit 3.0 fluorometer (ThermoFisher Scientific, Waltham, Massachusetts, United States) was used to quantify the concentration of DNA in the samples. Quality (based on 260/280 and 260/230 absorbance ratios) and integrity (based on DIN values) of the extracted DNA samples were checked using Nanodrop 1000 Spectrophotometer (ThermoFisher Scientific) and TapeStation 2200 DNA screen (Agilent Technologies, Santa Clara, California, United States), respectively.
DNA Extraction for 16S Amplicon Analysis
All the procedures mentioned here were performed under sterile conditions. Before extracting the DNA, intestine samples were transferred to a sterile Petri dish and placed on a cool-pack on dry ice. The intestine was opened and transferred to a 5ml tube containing 1.4mm Zirconium oxide beads (Cayman Chemical, Ann Arbor, Michigan, United States) and 2ml of InhibitEX buffer (Qiagen). Thereafter, DNA was extracted immediately using QIAamp DNA stool Mini Kit (Qiagen) according to the manufacturer’s protocol. The final elution volume was 75μl (ATE buffer). The same extraction method was employed for the mouth samples. The quality and quantity of the extracted DNA were checked with NanoDrop spectrophotometer ND-8000 (ThermoFisher Scientific).
Libraries Preparation and Sequencing
Whole-Genome Sequences
The Nextera DNA Flex library preparation kit with dual indices was used to prepare whole genome libraries based on the manufacturer’s protocol (Illumina, San Diego, California, United States). After tagmentation of the extracted gDNA samples using Bead-linked transposomes at 55°C for 15min, the sheared and tagmented gDNA was washed at 30°C for 15min. Amplification of the tagmented gDNA was performed using a 5-cycle PCR programme wherein the index 1 (i7) and index 2 (i5) adapters were added for sequencing cluster formation. The PCR program was started with an incubation at 68°C for 3min and a subsequent pre-denaturation at 98°C for 3min. In the following step, 5cycles of denaturation at 98°C for 45s, annealing at 62°C for 30s and extension at 68°C for 2min were first performed, followed by a final extension at 68°C for 1min. In the final step of the library preparation, the amplified libraries were purified through a double-sided bead (Bead-linked transposome; Illumina) purification procedure. The quality and normality of the libraries were assessed with the Agilent Tapestation instrument using High Sensitivity D1000 screen tape. After normalization based on the minimum observed molarity, the barcoded samples were pooled before the sequencing run. The 75bp paired-end sequencing was done on a NextSeq 500 sequencer (Illumina) at the sequencing platform of Nord University.
Bacterial 16S Sequences
Under sterile conditions, 16S rRNA gene libraries were constructed from DNA extracts using the specific bacterial primers 341F (5'CCTACGGGNGGCWGCAG 3') and 805R (5'GACTACNVGGGTWTCTAATCC 3'; Klindworth et al., 2013) flanked by overhang Illumina adapters targeting the hypervariable V3–V4 region (~460bp). PCR reactions were performed for each sample in 25μl, using Q5® High-Fidelity 2X Master Mix (New England Biolabs, Ipswich, Massachusetts, United States) and 2.5μl of DNA template (5ng/μl). PCR conditions consisted of an initial denaturation step at 95°C for 10min (1cycle), 30cycles at 95°C for 30s, 57°C for 30s, 72°C for 1min, and a final extension step at 72°C for 7min (1cycle).
An agarose gel (1.5%) was employed to check the amplified products. The PCR products were purified using the CleanNGS system (CleanNA, Waddinxveen, Netherlands) following the manufacturer’s instructions. The purified product was subjected to a second PCR (8cycles, 16S Metagenomic Sequencing Library Preparation, Illumina); this step was done to add dual indices and Illumina sequencing adapters Nextera XT Index Primer (Illumina). CleanNGS (CleanNA) was used to purify the obtained amplicon libraries. The quality of the libraries was checked on a Tapestation 2200 platform (Agilent Technologies). Thereafter, the libraries were quantified using the Quant-IT PicoGreen dsDNA assay kit (ThermoFisher Scientific) by the Synergy2 microplate reader (Biotek, Winooski, Vermont, United States). Next, the pooled libraries were quantified using the KAPA Library quantification kit (Roche, Basel, Switzerland). The libraries were checked by realtime qPCR LightCycler 480 (Roche) and then sequenced on an Illumina® MiSeq (PE300) platform (MiSeq Control Software 2.5.0.5 and Real-Time Analysis software 1.18.54.0).
Sequence Analysis
Whole-Genome Sequences
In order to perform demultiplexing and obtain the fastq files, the Illumina Experiment Manager v1.18.1 along with bcl2fastq v2.20.0.422 was used. Thereafter, dual adapter indexes and Ns from the 3' end of the raw reads were trimmed and the quality of the cleaned fastq files was assessed employing Trime_galore v0.4.4 (Babraham Bioinformatics; http://www.bioinformatics.babraham.ac.uk/projects/trim_galore/). The clean reads were then aligned to the reference genome O_niloticus_UMD_NMBU, GCA_001858045.3 (Conte et al., 2017) using Bowtie2 v0.12.8 with the --very-sensitive option (Langmead and Salzberg, 2012). The bcftools pipeline was applied for variant calling (Li, 2011), and the generated SAM files were converted to the binary format and sorted based on coordinates using samtools v1.9. Also, the samtools markedup command was used to mark duplicate reads. Then variants were called using bcftools mpileup command (bcftools 1.9) with the minimum base and mapping quality of 20 (−q 20 −Q 20). Using bcftools filter command accompanied by the options --SnpGap 5 -i ‘MQ>20 and QUAL>20 and DP>100 and DP<450 and TYPE=“snp,” only SNP variants were kept in the Variant Call Format (VCF). The missing genotypes were imputed using imp-states=1,600 option in Beagle v5.0 (Browning and Browning, 2016). Thereafter, using vcftools, the non-biallelic SNP variants were omitted so that the generated VCF file had only the biallelic SNPs (Danecek et al., 2011). This VCF file was read by vcfR package (Knaus and Grünwald, 2017).
Bacterial 16S Amplicon Sequences
The generated reads were truncated at 270bp using VSEARCH (Rognes et al., 2016), and then processed using MICCA pipeline (v1.7.2; Albanese et al., 2015). Sequences with a minimum overlap length of 60bp and a maximum mismatch of 20bp were merged. Next, the forward and reverse primers were trimmed off the merged reads and reads which did not contain the primers were discarded. Thereafter, the sequences with an expected error rate (Edgar and Flyvbjerg, 2015) >0.75 were filtered out and shorter than 400bp sequences were discarded. Filtered reads were denoised using the “de novo unoise” method implemented in MICCA, which utilise UNOISE3 algorithm (Edgar, 2016). The denoising method generates amplicon sequence variants (ASVs) which is based on correcting sequencing errors and determining true biological sequences at single-nucleotide resolution. The taxonomic assignment of the representative bacterial ASVs was performed using RDP classifier. The sequences were aligned using the NAST (Desantis et al., 2006) multiple sequence aligner, and a phylogenetic tree was prepared using the FastTree software available in the MICCA pipeline.
Statistical Analysis of Host Genetic Data
To quantify the genetic diversity of the inbred and outbred groups, we first determined the genetic diversity within members of the crossbred groups, and then the between groups genetic diversity. For this, we quantified the level of heterozygosity, using the population package of the Stacks 2.3b. Next, to assess the level of genetic differentiation based on allele frequencies between different groups, the Fst index was calculated using the StAMPP package (Pembleton et al., 2013). In order to quantify the genetic relationship between the inbred and outbred groups, Nei-based genetic distance between individuals was estimated using poppr (Kamvar et al., 2015) and adegenet (Jombart, 2008) packages and visualized using pheatmap package (Kolde and Kolde, 2015). Then the genetic relationship between the crossbred groups was assessed by PCoA (employing the abovementioned Nei-based genetic distance), also using the ape package (Paradis and Schliep, 2019). PERMANOVA (Permutational Multivariate Analysis of Variance) was performed to decipher the significance of genetic differences between the inbred and outbred groups. To further analyze the population structure of the inbred and outbred groups, admixture analysis was performed in adegenet for values of ancestries (K) from 1 to 10 with 10 repeats for each value of K, decided based on Bayesian Information Criteria. Four samples were removed due to the low quality of sequences.
Statistical Analysis of 16S Amplicon Data
Statistical analysis was conducted using R (version 3.6.3) software. The packages phyloseq (McMurdie and Holmes, 2013) and vegan (Oksanen et al., 2013) were employed to analyse the data. All plots were made using ggplot2 package (Wickham, 2011).
To understand the differences between the proportions of different bacteria in the inbred and outbred groups, we performed chi-square test and the associated post hoc analyses. A subset of the most dominant phyla was employed for this analysis. The similarities/differences in α-diversity were checked by Wilcoxon rank-sum test. Bacterial β-diversity was determined using unweighted and weighted UniFrac distances (Lozupone and Knight, 2005). Differences between the bacterial communities of the two groups were visualized by PCoA. After checking the dispersions within the data set of each group, statistically significant differences between the groups were assessed using PERMANOVA (Anderson, 2001; with 9,999 permutations), implemented in adonis function of the vegan R-package (Oksanen et al., 2013). DESeq2 (Love et al., 2014) package was employed to detect the differentially abundant ASVs in the non-rarefied data (McMurdie and Holmes, 2014). The core microbiota was analysed using the packages microbiome and microbiome utilities; at a detection level of 0.2% and prevalence level of 90%. The differences in the core bacterial community in the two crossbred groups were analysed by performing PERMANOVA on weighted and unweighted UniFrac distances. The abundances in the different ASVs which made up the core microbiome were analyzed using Spearman test (Zar, 2014) and correlation plot package (Wei and Simko, 2017).
Results
Genetic Background-Associated Changes in the Microbiome of Nile Tilapia
A total of 11,578,530 SNPs were obtained after the initial SNP calling. Bcftools was employed to first calculate genotype likelihoods for each position and then filter out every position with actual sequence variant. Thus, 4,693,720 SNPs were filtered out and finally after biallelic filtration, 6,825,083 SNP variants with an average coverage of 1.74 per sample were used in the final VCF file.
The genetic diversity analysis based on nucleotide sequences revealed that the observed heterozygosity (Ho) values were slightly higher compared to the expected heterozygosity values (He; Table 1).
Table 1. Observed (Ho) and expected (He) heterozygosity of the crossbred female Nile tilapia.
The fixation index (Fst) value within groups was 0.04 for both Inbred-S1 vs. Inbred-C6 and Outbred-S3 vs. Outbred-C9 comparisons. On the other hand, the Fst values between crossbred groups were in the range 0.06–0.08 (Table 2).
Table 2. Genetic differentiation, based on Fst index, of the crossbred female Nile tilapia.
The Nei-based genetic distances between the inbred and outbred groups were employed to generate a heatmap to understand their genetic relationships; differences between the groups are seen in Supplementary Figure 1. Principal Coordinates Analysis (PCoA) based on the Nei-based genetic distance indicated that the first two components captured 17.7 and 7.8% of the variation in the data set (Figure 1A). Furthermore, a PERMANOVA test based on the same genetic distance showed that the inbred and outbred groups were significantly different (p=0.001). The genetic sub-population clustering based on admixture analysis revealed that K=2 was the optimal number to explain the genetic structure of the inbred and outbred groups (Supplementary Figure 2). The results also indicated that 4 inbred individuals are genetically similar to the outbred population.
Figure 1. Genetic differentiation and microbiome in the inbred and outbred groups of Nile tilapia. (A) Principal coordinates analysis (PCoA) plot based on 6,825,083SNPs of the inbred and outbred groups. The ellipses were generated assuming that the data are from a multivariate normal distribution. (B) Phylum-level relative abundance of the microbial composition in the inbred and outbred groups. (C) Relative abundance of top 12 genera in the inbred and outbred groups.
To delineate the effect of genetic selection on gut microbiota composition, the inbred and outbred Nile tilapia were reared in a common garden and the environmental and nutritional factors that affect the microbiota were kept constant throughout the experimental period. The amplicon library of 16S rRNA gene, generated 12,034,190 high-quality reads with an average coverage of 54,016 reads per sample. Due to the variation in sample size, the reads were rarefied to 18,000 reads per sample (without replacement). Out of the 120 samples, six libraries with a number of reads below the cut off were discarded. After normalization we obtained 14,228 ASVs, distributed among 30 phyla and 695 genera.
First, we investigated the dominant communities in the two groups. In their order of dominance, the most dominant bacterial phyla in both the inbred and outbred fish groups were Proteobacteria, Fusobacteria, Firmicutes, Bacteroidetes, and Actinobacteria (Figure 1B). This order of dominance was reflected in the microbial composition at the genus level also. Most of the dominant genera in the two crossbred groups belonged to the phylum Proteobacteria (Acinetobacter, Curvibacter, Enhydrobacter, Escherichia/Shigella, Plesiomonas, Pseudomonas, Psychrobacter, and Undibacterium). The most abundant genus was Cetobacterium which belongs to the phylum Fusobacteria (Figure 1C).
To understand the differences in proportions of the dominant communities in each study group, we performed chi-square test. The analyses revealed that the abundances of the most dominant phyla in both the inbred and outbred groups were significantly different (Supplementary Table 1).
To characterize the microbial diversity within the samples, we calculated three ecological indexes, namely the Chao1 estimator of the number of species, which is a measure of richness, the Shannon diversity which measures the evenness of the microbial populations and the Simpson diversity, which measures the importance of dominant species (Marcon and Hérault, 2015; Hsieh et al., 2016). Shannon diversity analysis showed that the microbial diversity in the mouth of the inbred group was lower compared to the outbred group (Figure 2, p=0.01). The Simpson diversity analysis indicated that there were fewer dominant ASVs in the posterior intestine of the inbred group (Figure 3, p=0.04). Although there were no significant differences in species richness of the communities associated with the two groups, in each body site (Supplementary Table 2), there was an increasing trend (p=0.08; inbred higher richness) in the case of the anterior intestine (Figure 4). Furthermore, the diversity analysis of dominant bacteria (Simpson diversity) in the mouth and anterior intestine revealed a trend in differences (p=0.08 and 0.06, respectively; Figures 2, 4).
Figure 2. Alpha diversity of the bacteria in the mouth of the inbred and outbred groups of Nile tilapia. Species richness of the groups is not significantly different. Shannon diversity is higher in the outbred group (p=0.007, indicated with an asterisk). Simpson diversity indicated an increasing trend in the dominant ASVs of the outbred group (p=0.08). The boxplots show minimum, lower quartile, median, upper quartile and maximum values.
Figure 3. Alpha diversity of the bacteria in the posterior intestine of the inbred and outbred groups of Nile tilapia. Simpson diversity analysis showed that the dominant ASVs are higher in the outbred groups (p=0.04, indicated with an asterisk). The boxplots show minimum, lower quartile, median, upper quartile and maximum values.
Figure 4. Alpha diversity of the bacteria in the anterior intestine of the inbred and outbred groups of Nile tilapia. There is an increasing trend in the species richness of the inbred group (p=0.07). Simpson diversity shows an increasing trend in the dominant ASVs of the outbred group (p=0.06). The boxplots show minimum, lower quartile, median, upper quartile and maximum values.
Beta diversity analysis was performed to evaluate the overall dissimilarity between the two crossbred groups (Figure 5). The results of PERMANOVA on the unweighted UniFrac distances showed a significant difference between the bacterial composition in the posterior intestine of the inbred and outbred groups (p=0.003). There was no significant difference between the communities in the mouth or the anterior intestine of the two groups (p=0.082 and 0.311, respectively). In the mouth, there may exist a difference in composition between the two groups, based on the observed trend (Table 3).
Figure 5. Principal coordinates analyses (PCoA) using unweighted and weighted UniFrac distance matrices of the bacteria in the different body sites of the inbred and outbred groups of Nile tilapia. The ellipses were generated assuming that the data are from a multivariate normal distribution.
Table 3. Results of the analysis of homogeneity of group dispersions and PERMANOVA using distance (unweighted and weighted UniFrac) matrices.
Considering the weighted UniFrac distance, there was a significant difference in the community composition of the anterior intestine (p=0.001). In addition, there was a significant difference in the community of posterior intestine (p=0.003), but not in the case of mouth (p=0.37; Table 3).
Differential Abundance of ASVs: Outbred Group vs. Inbred Group
The package DESeq2 was used to identify the ASVs with a significantly different abundance in the outbred group compared to the inbred group. In the mouth, the bacteria belonging to Actinobacteria, Armatimonadetes, Firmicutes, and Proteobacteria were differentially abundant. There were six genera that belonged to the phylum Proteobacteria. Bacteria belonging to two genera (Psychrobacter and Polaromonas) were 5-fold higher in the outbred group, while those of Pseudomonas and Acinetobacter were 20-fold higher in the same group. Furthermore, an ASV of the genus Limnohabitans was about 9-fold lower and Comamonas was 20-fold lower in the outbred group. Lachnospiracea_incertae_sedis were about 5-fold higher in the outbred group, whereas the genus Bacillus was 20-fold lower. These two genera belong to the phylum Firmicutes. Also, Armatimonadetes_gp5 was 20-fold lower in the outbred group (Supplementary Figure 3). In the anterior intestine, the majority of the ASVs that were differentially abundant in the outbred group had fold changes between −5 and−15 (Supplementary Figures 4, 5) compared to the inbred group. However, the fold changes of the differentially abundant ASVs of Bacteroidetes, Fusobacteria, and Proteobacteria in the anterior intestine were between −5 and−10 (Supplementary Figure 4 shows selected differentially abundant ASVs; Supplementary Figure 5 shows all the differentially abundant ASVs). Similarly, in the posterior intestine, out of 31 ASVs that were differentially abundant, 30 ASVs had fold changes between −5 and−28 in the outbred group, while only one ASV that belongs to Acinetobacter was 20-fold higher in the outbred group (Supplementary Figure 6). Moreover, ASVs of Pediococcus and Bifidobacterium which belong to Firmicutes and Actinobacteria, respectively, were lower (log fold change; −5 and−8, respectively) in the posterior intestine of the outbred group (Supplementary Figure 6).
Core Microbiome and Variability in Taxa
In the mouth, 9 ASVs of the core microbiota belonged to the genera Staphylococcus, Curvibacter, Undibacterium, Escherichia/Shigella, Enhydrobacter, Propionibacterium, and Cetobacterium. However, two bacteria were classified only up to the order level – Actinomycetales, Sphingobacteriales (Figure 6). Taking all the 9 ASVs together, we observed a significant difference in the core microbiome in the inbred and outbred groups; only for unweighted UniFrac distance (R2=0.073, p=0.043; weighted UniFrac distance showed no significant difference; R2=0.024, p=0.445; Table 4). In the anterior and posterior intestine, the core ASVs were Staphylococcus, Plesiomonas, Undibacterium, Enhydrobacter, Propionibacterium, and Cetobacterium (Figures 7A,B). One extra genus was a member of the core microbiota in the anterior intestine (Escherichia/Shigella). One ASV in the anterior and posterior intestine was not classified up to the genus level, but was annotated as Actinomycetales (Figures 7A,B). The core microbiota in the anterior intestine of the inbred group was different from that of the outbred group; the weighted UniFrac distances-based assessment indicated the significant difference (PERMANOVA test; R2=0.155, p=0.001) between the two crossbred groups. As for the posterior intestine, we cannot specify that there is a significant difference between the crossbred groups (Table 4). The inter-individual variation in the abundance of the core microbiota in the intestine samples of the inbred group was less pronounced compared to the outbred groups (Supplementary Figure 7). On the other hand, the inter-individual variation in the abundances was more pronounced in the mouth of the inbred compared to the outbred group (Supplementary Figure 8).
Figure 6. Core microbiota in the mouth of the inbred and outbred groups of Nile tilapia. NAs: Not classified at the genus level, but at the order level, they are classified as Actinomycetales, Sphingobacteriales, in both groups.
Table 4. Results of the analysis of homogeneity of group dispersions and PERMANOVA using distance (unweighted and weighted UniFrac) matrices of the core microbiota.
Figure 7. Core microbiota in the anterior and posterior intestine of the inbred and outbred groups of Nile tilapia. (A) anterior intestine and (B) posterior intestine. NA: at the order level is classified as Actinomycetales.
Discussion
The genetic structure of wild/domestic/experimental animals can be altered through breeding to retain desired phenotypic and genotypic traits across generations. It is known that selective breeding can preserve desired traits, which can affect the bacterial profile that is highly correlated to host health.
Gut microbiota in fish has been studied extensively in recent years considering mainly its importance in host health. In the present study, we used genetically distinct (based on SNP analysis) inbred and outbred Nile tilapia to investigate the impact of crossbreeding on the composition of the mouth and intestine bacteria.
Mouth and Intestine Bacterial Community Composition and Diversity in the Inbred and Outbred Groups
Although male Nile tilapia are widely farmed because of their higher growth rate, in the present study, we analyzed the microbial community in females, which are mouthbrooders. Hence, we believe that studying the microbial communities in its mouth will yield interesting results. In humans, microbiota is transferred from different body sites of mothers to infants (Ferretti et al., 2018). Moreover, microbial symbionts from discus (Symphysodon aequifasciata, another fish of the family Cichlidae) parents are vertically transferred to fry through feeding of a cutaneous mucus secretion (Sylvain and Derome, 2017). The most dominant phyla found in our samples were Proteobacteria, Fusobacteria, Firmicutes, Bacteroidetes, and Actinobacteria (Supplementary Figure 9). These are known to be the most represented phyla in model fishes such as zebrafish and threespine stickleback (Legrand et al., 2020). They are also dominant in farmed fishes like Nile tilapia even though many factors including diet (Ray et al., 2017; Souza et al., 2020), rearing systems (Giatsis et al., 2015; Yukgehnaish et al., 2020), and salinity (Zhang et al., 2016; Yukgehnaish et al., 2020) affect the abundance of these phyla in the gut. However, the role of crossbreeding in shaping microbial communities has not yet been reported in fish although it is studied in mice (Pang et al., 2012; Kreisinger et al., 2014), mammals (Alessandri et al., 2019), and plants (Wagner et al., 2020).
The dominant phyla were the same in both the inbred and outbred groups of Nile tilapia. Proteobacteria are facultative anaerobes, and they are the most abundant bacterial phylum in fish gut (Egerton et al., 2018). Furthermore, bacteria such as Escherichia and Enhydrobacter belonging to this phylum have the ability to make the gut environment conducive to strict anaerobes which colonize healthy gut (Shin et al., 2015). Although the aforementioned genera were present in the mouth and intestine of both the outbred and inbred fish, their abundances in the two groups were different. In addition, the genus Curvibacter which was present in both groups is known to have a critical role in colonization in freshwater invertebrates (Wein et al., 2018).
Alpha diversity analysis revealed that our crossbreeding strategy increased the microbial evenness in the mouth of the outbred group, in which we observed apparently higher species richness. The increasing trend in the dominant bacteria in the mouth and the anterior intestine of the outbred group along with the significant increase in the posterior intestine suggests that the dominant bacteria in the outbred groups are more diverse compared to the inbred group. On the other hand, the increasing trend in the species richness in the anterior intestine of the inbred group suggests that the bacterial community is more diverse in this intestinal segment of the inbred group compared to the outbred group. The abovementioned findings are similar to the results of the PCoA analysis that used UniFrac distances. Microbial diversity is believed to have a positive correlation with host health (Deng et al., 2019). However, Reese and Dunn (2018) have stated that “understanding diversity in host-associated microbial communities will not be as simple as ‘more diversity is better’.” Hence, it is not ideal to correlate host health with the diversity in the outbred group. Studies in Nile tilapia have not reported a significant difference in the diversity of gut microbiota as a pathogenic effect (Suphoronski et al., 2019; Silva et al., 2020). On the other hand, while diet was shown to increase the species richness of bacteria in the gut, another environmental factor, salinity, was found to decrease the richness of bacteria in Nile tilapia (Zhang et al., 2016). The implication of the increasing trend in diversity in the anterior intestine of the inbred group should be clarified by conducting studies on the bacteria in this segment and their effect on nutritional physiology (Hallali et al., 2018). Thus, in addition to the aforementioned factors, we suggest that crossbreeding is a determinant of both the mouth and intestine bacterial diversity in female Nile tilapia.
Significant Differences Between the ASV Abundance of the Inbred and Outbred Groups
Fish gut harbors complex and diverse microbial communities, and the site is a reservoir of many opportunistic pathogens belonging to the genera Acinetobacter, Aeromonas, Psychrobacter, Flavobacterium, Pseudomonas, and Pleisomonas. Many commensal bacteria including Cetobacterium, Methylobacterium, Sphingomonas, and Propionibacterium (Suphoronski et al., 2019; Legrand et al., 2020; Silva et al., 2020) that colonise the fish gut are essential for the production of vitamin B12 and antimicrobial metabolites (Suphoronski et al., 2019; Legrand et al., 2020), protection against pathogens such as Flavobacterium (Boutin et al., 2014), and improving host health (Boutin et al., 2013). The differential ASV analysis revealed that the abundances of some of these opportunistic pathogens (Psychrobacter, Pseudomonas, and Acinetobacter) were more than 5-fold in the mouth of the outbred group compared to the inbred group. In the anterior and posterior intestine of the outbred group, although the opportunistic pathogens belonging to the genera Acinetobacter, Aeromonas, Pleisomonas, Psychrobacter, Pseudomonas, and Flavobacterium were differentially abundant, their fold changes were less than 5-fold. The bacterial community in the mouth is extensively exposed to the external environment, and we found that the opportunistic pathogens in the mouth are more abundant in the outbred group. On the other hand, the abundance of potential pathogens was lower in the intestine of the inbred group. Pseudomonas sp. are opportunistic pathogens and they cause high mortality in farmed fishes (Oh et al., 2019). Moreover, bacteria belonging to Flavobacterium were reported to cause acute bacteremia primarily in small fishes or more chronic disease in larger fishes (Semple et al., 2020). Although the outbred fish had a more diverse microbiome, they appear to harbor potential opportunistic bacteria also.
Interestingly, the abundance of potential beneficial bacteria (Cetobacterium, Methylobacterium, Sphingomonas, and Propionibacterium; Boutin et al., 2013, 2014; Suphoronski et al., 2019; Legrand et al., 2020) was higher in the inbred group. Many studies report that commensal microbiota in the gut plays an important role in regulating the growth of other microbes by competing for space and nutrition. The mouth of the inbred fish had higher abundance of Aeromonas sp. which was found to compete for nutrients and play a negative role during infection (Wiles et al., 2016; Legrand et al., 2020). On the other hand, the bacteria that had higher abundance in the posterior intestine of the inbred tilapia, namely Enhydrobacter sp., is a commensal microbe in rainbow trout (Oncorhynchus mykiss), which is known to produce entericidin, and this antitoxin peptide inhibits the growth of certain pathogens such as those belonging to Flavobacterium (Legrand et al., 2020). Furthermore, Pediococcus and Bifidobacterium which were found to be more abundant in the anterior and posterior intestine of the inbred groups compared to the outbred group are known to outcompete some invasive pathogens, associated with tilapia intestinal mucosa (Ferguson et al., 2010; Standen et al., 2013) and promote fish growth (Ayyat et al., 2014). Thus, the inbred group had a higher abundance of potential beneficial commensal bacteria.
Changes in Core Microbiome
The transient allochthonous microbiome of fish is associated with digesta and is usually expelled after some period as they are predominantly influenced by diet. On the other hand, the resident microbes that belong to the autochthonous microbiome colonise the mucus surface in the gut and make up the core microbiome (Egerton et al., 2018). These microbial communities, which are known to be vertically transmitted (Risely, 2020), associate with the host’s cells (Egerton et al., 2018; Legrand et al., 2020). In the present study, the core microbiome in each body site was determined based on the ASVs present in all samples in each group. However, the inter-individual variation in abundance that we observed is similar to the learning from studies on zebrafish (Burns et al., 2016) and mice (Pang et al., 2012). In mice, inbreeding was found to reduce the inter-individual variation (Pang et al., 2012). The inter-individual variation in the core microbiome in the intestine of the inbred group is much lesser compared to the outbred group. In contrast, such similarity was not observed in the mouth of the inbred fish; this was attributed to the effect of external environment in other studies (Lokesh and Kiron, 2016; Krotman et al., 2020). However, in the present study, environmental factors were kept constant throughout the study period. In humans, the initial oral colonizers from the vagina and mother’s milk and mouth can be perturbed by environmental factors (Kilian, 2018).
The most dominant bacterial phylum in the two study groups was Proteobacteria. Nevertheless, Cetobacterium (phylum Fusobacteria) was found to be dominant in the anterior and posterior intestine of the inbred group, while its proportion was reduced in the outbred group. Previous studies conducted on Nile tilapia showed that the composition of Cetobacterium spp., the most prevalent genera in tilapia gut, was not affected by diets (Ray et al., 2017) or presence of pathogens (Suphoronski et al., 2019; Silva et al., 2020). Other reports that studied the influence of factors including rearing environment (Giatsis et al., 2015), and salinity (Zhang et al., 2016) on the gut microbial composition substantiates our finding that Cetobacterium is a core member of the bacterial community. Based on the present study, it appears that the crossbreeding strategy does not impact the presence of this core member in the mouth and intestine of Nile tilapia.
Some of the commonly reported bacteria in the intestine of Nile tilapia (Staphylococcus, Cetobacterium, Plesiomonas, Enhydrobacter, Undibacterium, and Propionibacterium) were present in both groups. However, some core microbiome members such as Pseudomonas and Curvibacter were present only in the mouth of both groups. A study employing turbot (Scophthalmus maximus) showed that a similar microbiome community was present in the intestine of different breeds fed with different diets and reared in different water environments. In addition, it was reported that core microbiome could colonize fish gut for a long term and it could have a vital physiological significance to the host (Zhang et al., 2020). This suggests that fishes preserve their core microbiome community despite differences in environmental factors.
Host Genetics and Intestine Microbiome
Growing evidence shows that host genetics plays a key role in shaping the gut microbiome of mammals (Hufeldt et al., 2010; Miller et al., 2018; Alessandri et al., 2019), but not to the same degree as that of environmental factors (Davenport, 2016). While there are many reports on diet-based microbiota differences in fish, evidences of fish genetics-associated microbiota are sparse (Li et al., 2014; Kokou et al., 2018).
Our genetic diversity analysis indicated a small but significant difference between the inbred and outbred fish. Unexpectedly, the observed heterozygosity was slightly higher than the expected heterozygosity, probably arising from the low genetic diversity values in both the inbred and outbred groups. The Ho, He, and Fst results that we obtained are likely due to small number of founders with a similar genetic background since the F0 generation of the fish were caught from the same area. The F0 itself may have lost considerable genetic diversity, as noted for birds; a small number of founders in a population increased the probability of inbreeding and associated gene diversity loss (Jamieson, 2011).
Wild Nile tilapia populations in West Africa are reported to have low diversity, especially, the species within a particular region; for example in Gambia River and the far western region of the Niger River (Lind et al., 2019). Nile tilapia is seen as a range-limited species in these areas, and founder effect was reported to be the reason for their genetic diversity reduction (Lind et al., 2019). In addition, Fst results also indicated the low genetic differentiation within the inbred groups as well as the outbred groups.
Anthropogenic needs not only alter species behavior, feeding habits, rearing environment, and traits within the host genotype but also reshape the gut microbiota of domesticated/captivated animals (Li et al., 2014; Alessandri et al., 2019). A study on blue tilapia, which was selectively bred to retain a host genotype, has reported that gut microbiome was linked to host genotype as well as specific bacteria such as Cetobacterium somerae (Kokou et al., 2018). This bacterium is a cobalamin producer (Tsuchiya et al., 2008; Degnan et al., 2014) and fishes with high abundance of C. somerae do not require dietary vitamin B12 (Sugita et al., 1991; Tsuchiya et al., 2008).
In order to analyse the genetic effect (by controlling the mating strategy) on the mouth and gut microbiota, the fish were kept in the same environmental conditions and fed the same diet, since both these factors are determinants of host microbial communities. Thus, crossing strategy influenced the microbial alpha diversity and composition in Nile tilapia. A similar effect on the midgut microbiota composition was observed in selectively bred trout (Brown et al., 2019). In addition, a study conducted on mice suggested that the alpha diversity of the gastrointestinal tract microbiota is slightly decreased in the inbred individuals (Kreisinger et al., 2014). Thus, the differences in the diversities of the microbial communities of the two groups could be attributed to crossbreeding strategy.
The differences in abundance of the microbial composition of the core microbiome in the individual samples from the mouth of the inbred fish were more pronounced compared to the outbred groups. In the mouth, influence of an external environmental factor (water) appears to surpass that of the host genetics. On the other hand, there was more similarity in the abundance of the bacterial communities in the individual intestine samples of the inbred group compared to the outbred group of Nile tilapia. Host genetics is known to have a long-lasting effect on the gut microbial communities and this is due to maternal transfer during early development (Kreisinger et al., 2014). A core microbiota is heritable in several species (Hauffe and Barelli, 2019), including cichlids (Baldo et al., 2017). The similarities in the abundances of the taxa in the inbred group of Nile tilapia, which is also a cichlid fish, suggest that the microbial composition in the gut is more established without being affected by the external environment. A study conducted in mice showed that the inter-individual variation in the gut microbiome of the inbred group is lower compared to the outbred animals (Hufeldt et al., 2010). Furthermore, in humans the similarity of the gut microbiome is higher among closer relatives in families (Zoetendal et al., 2001). Therefore, this finding suggests that the genetic factor is more prominent in the intestine of the inbred groups and the effect is likely the inheritance of the microbial profile to the offspring of the fish, especially the core microbiome.
We report for the first time the effect of inbreeding and outbreeding on the mouth and intestine microbiome in Nile tilapia. The genetic relationship and structure analysis indicated the genetic differentiation between the inbred and outbred groups. Differential ASV analysis revealed the abundance of the potential opportunistic pathogens such as Flavobacterium in the outbred group and beneficial bacteria like Bifidobacterium and Pediococcus in the inbred group. We also found that Cetobacterium is the core member in both groups, but its abundance was higher in the intestine of the inbred group. The inbred fish which has less inter-individual microbiome variability, could be a better choice for controlled studies that examine the maternal transfer of intestine microbiome to offspring. We highlight that crossbreeding can influence Nile tilapia bacterial communities.
Data Availability Statement
The data used in this study is available at European Nucleotide Archive (ENA) with accession no PRJEB40093. Whole-genome data of the host are available at Sequence Read Archive (SRA) with accession no PRJNA719847.
Ethics Statement
The animal study was reviewed and approved by Norwegian Animal Research Authority, FOTS ID 1042.
Author Contributions
VK, JF, and YA designed the study. YA carried out the sampling and lab work, analysed the data, and wrote the manuscript. CD also analysed the data. SL helped in sequencing and data generation. DA was involved in initial data analyses. OJ and AN prepared the whole genome shotgun sequencing libraries. OJ performed the SNPs analysis of the host. YA, JF, CD, and VK interpreted the data. VK, JF, and CD reviewed and edited the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This study was funded by the European Research Council (ERC) under the European Union’s Horizon 2020 Research and Innovation Programme (grant agreement no. 683210), the Research Council of Norway under the Toppforsk Programme (grant agreement no. 250548/F20), and the Faculty of Biosciences and Aquaculture (Nord University, Norway).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We are grateful to Hilde Ribe, Steinar Johnsen, Øivind Torslett, particularly to Kaspar Klaudiussen at Nord University’s research station (Bodø, Norway) for their assistance with fish husbandry and welfare of the fish. The authors are thankful to Massimo Pindo (Edmund Mach Foundation, San Michele, Italy) who performed formal analysis and data curation. We also thank Erika Stefani and Stefano Piazza at the sequencing platform of Edmund Mach Foundation for their significant assistance in library preparation and sequencing. Lastly, we express our gratitude to Bisa Saraswathy for her assistance in data analysis, scientific input and preparation of the manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.709611/full#supplementary-material
References
Albanese, D., Fontana, P., De Filippo, C., Cavalieri, D., and Donati, C. (2015). MICCA: a complete and accurate software for taxonomic profiling of metagenomic data. Sci. Rep. 5:9743. doi: 10.1038/srep09743
Alessandri, G., Milani, C., Mancabelli, L., Mangifesta, M., Lugli, G. A., Viappiani, A., et al. (2019). The impact of human-facilitated selection on the gut microbiota of domesticated mammals. FEMS Microbiol. Ecol. 95:fiz121. doi: 10.1093/femsec/fiz121
Anderson, M. J. (2001). A new method for non-parametric multivariate analysis of variance. Austral. Ecol. 26, 32–46. doi: 10.1111/j.1442-9993.2001.01070.pp.x
Ayyat, M., Labib, H. M., and Mahmoud, H. K. (2014). A probiotic cocktail as a growth promoter in Nile tilapia (Oreochromis niloticus). J. Appl. Aquac. 26, 208–215. doi: 10.1080/10454438.2014.934164
Baldo, L., Pretus, J. L., Riera, J. L., Musilova, Z., Nyom, A. R. B., and Salzburger, W. (2017). Convergence of gut microbiotas in the adaptive radiations of African cichlid fishes. ISME J. 11, 1975–1987. doi: 10.1038/ismej.2017.62
Bolivar, R. B., and Newkirk, G. F. (2002). Response to within family selection for body weight in Nile tilapia (Oreochromis niloticus) using a single-trait animal model. Aquaculture 204, 371–381. doi: 10.1016/S0044-8486(01)00824-9
Boutin, S., Bernatchez, L., Audet, C., and Derôme, N. (2013). Network analysis highlights complex interactions between pathogen, host and commensal microbiota. PLoS One 8:e84772. doi: 10.1371/journal.pone.0084772
Boutin, S., Sauvage, C., Bernatchez, L., Audet, C., and Derome, N. (2014). Inter individual variations of the fish skin microbiota: host genetics basis of mutualism? PLoS One 9:e102649. doi: 10.1371/journal.pone.0102649
Brown, R. M., Wiens, G. D., and Salinas, I. (2019). Analysis of the gut and gill microbiome of resistant and susceptible lines of rainbow trout (Oncorhynchus mykiss). Fish Shellfish Immunol. 86, 497–506. doi: 10.1016/j.fsi.2018.11.079
Browning, B. L., and Browning, S. R. (2016). Genotype imputation with millions of reference samples. Am. J. Hum. Genet. 98, 116–126. doi: 10.1016/j.ajhg.2015.11.020
Burns, A. R., Stephens, W. Z., Stagaman, K., Wong, S., Rawls, J. F., Guillemin, K., et al. (2016). Contribution of neutral processes to the assembly of gut microbial communities in the zebrafish over host development. ISME J. 10, 655–664. doi: 10.1038/ismej.2015.142
Conte, M. A., Gammerdinger, W. J., Bartie, K. L., Penman, D. J., and Kocher, T. D. (2017). A high quality assembly of the Nile tilapia (Oreochromis niloticus) genome reveals the structure of two sex determination regions. BMC Genomics 18:341. doi: 10.1186/s12864-017-3723-5
D’ambrosio, J., Phocas, F., Haffray, P., Bestin, A., Brard-Fudulea, S., Poncet, C., et al. (2019). Genome-wide estimates of genetic diversity, inbreeding and effective size of experimental and commercial rainbow trout lines undergoing selective breeding. Genet. Sel. Evol. 51:26. doi: 10.1186/s12711-019-0468-4
Danecek, P., Auton, A., Abecasis, G., Albers, C. A., Banks, E., Depristo, M. A., et al. (2011). The variant call format and VCFtools. Bioinformatics 27, 2156–2158. doi: 10.1093/bioinformatics/btr330
Davenport, E. R. (2016). Elucidating the role of the host genome in shaping microbiome composition. Gut Microbes 7, 178–184. doi: 10.1080/19490976.2016.1155022
Degnan, P. H., Taga, M. E., and Goodman, A. L. (2014). Vitamin B12 as a modulator of gut microbial ecology. Cell Metab. 20, 769–778. doi: 10.1016/j.cmet.2014.10.002
Deng, F., Li, Y., and Zhao, J. (2019). The gut microbiome of healthy long-living people. Aging (Albany NY) 11, 289–290. doi: 10.18632/aging.101771
Desantis, T., Hugenholtz, P., Keller, K., Brodie, E., Larsen, N., Piceno, Y., et al. (2006). NAST: a multiple sequence alignment server for comparative analysis of 16S rRNA genes. Nucleic Acids Res. 34, W394–W399. doi: 10.1093/nar/gkl244
Edgar, R. C. (2016). UNOISE2: improved error-correction for Illumina 16S and ITS amplicon sequencing. bioRxiv [Preprint]. doi: 10.1101/081257
Edgar, R. C., and Flyvbjerg, H. (2015). Error filtering, pair assembly and error correction for next-generation sequencing reads. Bioinformatics 31, 3476–3482. doi: 10.1093/bioinformatics/btv401
Egerton, S., Culloty, S., Whooley, J., Stanton, C., and Ross, R. P. (2018). The gut microbiota of marine fish. Front. Microbiol. 9:873. doi: 10.3389/fmicb.2018.00873
Ericsson, A. C., Davis, J. W., Spollen, W., Bivens, N., Givan, S., Hagan, C. E., et al. (2015). Effects of vendor and genetic background on the composition of the fecal microbiota of inbred mice. PLoS One 10:e0116704. doi: 10.1371/journal.pone.0116704
Ferguson, R., Merrifield, D. L., Harper, G. M., Rawling, M. D., Mustafa, S., Picchietti, S., et al. (2010). The effect of Pediococcus acidilactici on the gut microbiota and immune status of on-growing red tilapia (Oreochromis niloticus). J. Appl. Microbiol. 109, 851–862. doi: 10.1111/j.1365-2672.2010.04713.x
Ferretti, P., Pasolli, E., Tett, A., Asnicar, F., Gorfer, V., Fedi, S., et al. (2018). Mother-to-infant microbial transmission from different body sites shapes the developing infant gut microbiome. Cell Host Microbe 24, 133.e135–145.e135. doi: 10.1016/j.chom.2018.06.005
Flint, A., and Woolliams, J. (2008). Precision animal breeding. Philos. Trans. R. Soc. B Biol. Sci. 363, 573–590. doi: 10.1098/rstb.2007.2171
Giatsis, C., Sipkema, D., Smidt, H., Heilig, H., Benvenuti, G., Verreth, J., et al. (2015). The impact of rearing environment on the development of gut microbiota in tilapia larvae. Sci. Rep. 5:18206. doi: 10.1038/srep18206
Gjedrem, T., Robinson, N., and Rye, M. (2012). The importance of selective breeding in aquaculture to meet future demands for animal protein: a review. Aquaculture 350, 117–129. doi: 10.1016/j.aquaculture.2012.04.008
Gjedrem, T., and Rye, M. (2018). Selection response in fish and shellfish: a review. Rev. Aquac. 10, 168–179. doi: 10.1111/raq.12154
Gratacap, R. L., Wargelius, A., Edvardsen, R. B., and Houston, R. D. (2019). Potential of genome editing to improve aquaculture breeding and production. Trends Genet. 35, 672–684. doi: 10.1016/j.tig.2019.06.006
Hallali, E., Kokou, F., Chourasia, T. K., Nitzan, T., Con, P., Harpaz, S., et al. (2018). Dietary salt levels affect digestibility, intestinal gene expression, and the microbiome, in Nile tilapia (Oreochromis niloticus). PLoS One 13:e0202351. doi: 10.1371/journal.pone.0202351
Hart, M. L., Ericsson, A. C., Lloyd, K. K., Grimsrud, K. N., Rogala, A. R., Godfrey, V. L., et al. (2018). Development of outbred CD1 mouse colonies with distinct standardized gut microbiota profiles for use in complex microbiota targeted studies. Sci. Rep. 8:10107. doi: 10.1038/s41598-018-28448-0
Hauffe, H. C., and Barelli, C. (2019). Conserve the germs: the gut microbiota and adaptive potential. Conserv. Genet. 20, 19–27. doi: 10.1007/s10592-019-01150-y
Hill, W. G. (2014). Applications of population genetics to animal breeding, from wright, fisher and lush to genomic prediction. Genetics 196, 1–16. doi: 10.1534/genetics.112.147850
Hill, W. G. (2016). Is continued genetic improvement of livestock sustainable? Genetics 202, 877–881. doi: 10.1534/genetics.115.186650
Hsieh, T., Ma, K., and Chao, A. (2016). iNEXT: an R package for rarefaction and extrapolation of species diversity (Hill numbers). Methods Ecol. Evol. 7, 1451–1456. doi: 10.1111/2041-210X.12613
Hufeldt, M. R., Nielsen, D. S., Vogensen, F. K., Midtvedt, T., and Hansen, A. K. (2010). Variation in the gut microbiota of laboratory mice is related to both genetic and environmental factors. Comp. Med. 60, 336–347.
Ina-Salwany, M. Y., Al-Saari, N., Mohamad, A., Mursidi, F.-A., Mohd-Aris, A., Amal, M. N. A., et al. (2019). Vibriosis in fish: a review on disease development and prevention. J. Aquat. Anim. Health 31, 3–22. doi: 10.1002/aah.10045
Jamieson, I. G. (2011). Founder effects, inbreeding, and loss of genetic diversity in four avian reintroduction programs. Conserv. Biol. 25, 115–123. doi: 10.1111/j.1523-1739.2010.01574.x
Jantrakajorn, S., Maisak, H., and Wongtavatchai, J. (2014). Comprehensive investigation of Streptococcosis outbreaks in cultured Nile tilapia, Oreochromis niloticus, and red tilapia, Oreochromis sp., of Thailand. J. World Aquacult. Soc. 45, 392–402. doi: 10.1111/jwas.12131
Jombart, T. (2008). Adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics 24, 1403–1405. doi: 10.1093/bioinformatics/btn129
Kamvar, Z. N., Brooks, J. C., and Grünwald, N. J. (2015). Novel R tools for analysis of genome-wide population genetic data with emphasis on clonality. Front. Gent. 6:208. doi: 10.3389/fgene.2015.00208
Kilian, M. (2018). The oral microbiome–friend or foe? Eur. J. Oral Sci. 126, 5–12. doi: 10.1111/eos.12527
Klindworth, A., Pruesse, E., Schweer, T., Peplies, J., Quast, C., Horn, M., et al. (2013). Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 41:e1. doi: 10.1093/nar/gks808
Knaus, B. J., and Grünwald, N. J. (2017). vcfr: a package to manipulate and visualize variant call format data in R. Mol. Ecol. 17, 44–53. doi: 10.1111/1755-0998.12549
Kokou, F., Sasson, G., Nitzan, T., Doron-Faigenboim, A., Harpaz, S., Cnaani, A., et al. (2018). Host genetic selection for cold tolerance shapes microbiome composition and modulates its response to temperature. eLife 7:e36398. doi: 10.7554/eLife.36398
Konstantinidis, I., Sætrom, P., Mjelle, R., Nedoluzhko, A. V., Robledo, D., and Fernandes, J. M. (2020). Major gene expression changes and epigenetic remodelling in Nile tilapia muscle after just one generation of domestication. Epigenetics 15, 1052–1067. doi: 10.1080/15592294.2020.1748914
Kreisinger, J., Čížková, D., Vohánka, J., and Piálek, J. (2014). Gastrointestinal microbiota of wild and inbred individuals of two house mouse subspecies assessed using high-throughput parallel pyrosequencing. Mol. Ecol. 23, 5048–5060. doi: 10.1111/mec.12909
Krotman, Y., Yergaliyev, T. M., Shani, R. A., Avrahami, Y., and Szitenberg, A. (2020). Dissecting the factors shaping fish skin microbiomes in a heterogeneous inland water system. Microbiome 8:9. doi: 10.1186/s40168-020-0784-5
Langmead, B., and Salzberg, S. L. (2012). Fast gapped-read alignment with Bowtie 2. Nat. Methods 9:357. doi: 10.1038/nmeth.1923
Legrand, T. P., Wynne, J. W., Weyrich, L. S., and Oxley, A. P. (2020). A microbial sea of possibilities: current knowledge and prospects for an improved understanding of the fish microbiome. Rev. Aquac. 12, 1101–1134. doi: 10.1111/raq.12375
Li, H. (2011). A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics 27, 2987–2993. doi: 10.1093/bioinformatics/btr509
Li, J., Ni, J., Li, J., Wang, C., Li, X., Wu, S., et al. (2014). Comparative study on gastrointestinal microbiota of eight fish species with different feeding habits. J. Appl. Microbiol. 117, 1750–1760. doi: 10.1111/jam.12663
Lind, C. E., Agyakwah, S. K., Attipoe, F. Y., Nugent, C., Crooijmans, R. P., and Toguyeni, A. (2019). Genetic diversity of Nile tilapia (Oreochromis niloticus) throughout West Africa. Sci. Rep. 9:16767. doi: 10.1038/s41598-019-53295-y
Lind, C., Ponzoni, R., Nguyen, N., and Khaw, H. (2012). Selective breeding in fish and conservation of genetic resources for aquaculture. Reprod. Domest. Anim. 47, 255–263. doi: 10.1111/j.1439-0531.2012.02084.x
Lokesh, J., and Kiron, V. (2016). Transition from freshwater to seawater reshapes the skin-associated microbiota of Atlantic salmon. Sci. Rep. 6:19707. doi: 10.1038/srep19707
Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15:550. doi: 10.1186/s13059-014-0550-8
Lozupone, C., and Knight, R. (2005). UniFrac: a new phylogenetic method for comparing microbial communities. Appl. Environ. Microbiol. 71, 8228–8235. doi: 10.1128/AEM.71.12.8228-8235.2005
Marcon, E., and Hérault, B. (2015). entropart: an R package to measure and partition diversity. J. Stat. Softw. 67, 1–26. doi: 10.18637/jss.v067.i08
McMurdie, P. J., and Holmes, S. (2013). phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One 8:e61217. doi: 10.1371/journal.pone.0061217
McMurdie, P. J., and Holmes, S. (2014). Waste not, want not: why rarefying microbiome data is inadmissible. PLoS Comput. Biol. 10:e1003531. doi: 10.1371/journal.pcbi.1003531
Mehar, M., Mekkawy, W., Mcdougall, C., and Benzie, J. A. (2019). Fish trait preferences: a review of existing knowledge and implications for breeding programmes. Rev. Aquac. 12, 1273–1296. doi: 10.1111/raq.12382
Miller, E. T., Svanbäck, R., and Bohannan, B. J. (2018). Microbiomes as metacommunities: understanding host-associated microbes through metacommunity ecology. Trends Ecol. Evol. 33, 926–935. doi: 10.1016/j.tree.2018.09.002
Oh, W. T., Kim, J. H., Jun, J. W., Giri, S. S., Yun, S., Kim, H. J., et al. (2019). Genetic characterization and pathological analysis of a novel bacterial pathogen, Pseudomonas tructae, in rainbow trout (Oncorhynchus mykiss). Microorganisms 7:432. doi: 10.3390/microorganisms7100432
Oksanen, J., Blanchet, F. G., Kindt, R., Legendre, P., Minchin, P., O’hara, R., et al. (2013). Community ecology package. R package version, 2.0–2.
Oriá, R. B., Malva, J. O., Foley, P. L., Freitas, R. S., Bolick, D. T., and Guerrant, R. L. (2018). Revisiting inbred mouse models to study the developing brain: the potential role of intestinal microbiota. Front. Hum. Neurosci. 12:358. doi: 10.3389/fnhum.2018.00358
Pang, W., Stradiotto, D., Krych, L., Karlskov-Mortensen, P., Vogensen, F. K., Nielsen, D. S., et al. (2012). Selective inbreeding does not increase gut microbiota similarity in BALB/c mice. Lab. Anim. 46, 335–337. doi: 10.1258/la.2012.012040
Paradis, E., and Schliep, K. (2019). ape 5.0: an environment for modern phylogenetics and evolutionary analyses in R. Bioinformatics 35, 526–528. doi: 10.1093/bioinformatics/bty633
Pembleton, L. W., Cogan, N. O., and Forster, J. W. (2013). St AMPP: an R package for calculation of genetic differentiation and structure of mixed-ploidy level populations. Mol. Ecol. 13, 946–952. doi: 10.1111/1755-0998.12129
Piazzon, M. C., Naya-Català, F., Perera, E., Palenzuela, O., Sitjà-Bobadilla, A., and Pérez-Sánchez, J. (2020). Genetic selection for growth drives differences in intestinal microbiota composition and parasite disease resistance in gilthead sea bream. Microbiome 8:168. doi: 10.1186/s40168-020-00922-w
Ray, C., Bujan, N., Tarnecki, A., Davis, A. D., Browdy, C., and Arias, C. (2017). Analysis of the gut microbiome of Nile tilapia Oreochromis niloticus L. fed diets supplemented with Previda® and saponin. J. Fishscicom 11:36. doi: 10.21767/1307-234X.1000116
Reese, A. T., and Dunn, R. R. (2018). Drivers of microbiome biodiversity: a review of general rules, feces, and ignorance. MBio 9:e01294–18. doi: 10.1128/mBio.01294-18
Risely, A. (2020). Applying the core microbiome to understand host–microbe systems. J. Anim. Ecol. 89, 1549–1558. doi: 10.1111/1365-2656.13229
Rognes, T., Flouri, T., Nichols, B., Quince, C., and Mahé, F. (2016). VSEARCH: a versatile open source tool for metagenomics. PeerJ 4:e2584. doi: 10.7717/peerj.2584
Romana-Eguia, M. R. R., Ikeda, M., Basiao, Z. U., and Taniguchi, N. (2005). Genetic changes during mass selection for growth in Nile tilapia, Oreochromis niloticus (L.), assessed by microsatellites. Aquac. Res. 36, 69–78. doi: 10.1111/j.1365-2109.2004.01185.x
Santos, V. B. D., Mareco, E. A., and Dal Pai Silva, M. (2013). Growth curves of Nile tilapia (Oreochromis niloticus) strains cultivated at different temperatures. Acta Sci. Anim. Sci. 35, 235–242. doi: 10.4025/actascianimsci.v35i3.19443
Schultz, B., Serão, N., and Ross, J. W. (2020). “Genetic improvement of livestock, from conventional breeding to biotechnological approaches,” in Animal Agriculture. eds. F. W. Bazer, G. C. Lamb, and G. Wu (London, UK: Elsevier), 393–405.
Semple, S. L., Bols, N. C., Lumsden, J. S., and Dixon, B. (2020). Understanding the pathogenesis of Flavobacterium psychrophilum using the rainbow trout monocyte/macrophage-like cell line, RTS11, as an infection model. Microb. Pathog. 139:103910. doi: 10.1016/j.micpath.2019.103910
Shin, N.-R., Whon, T. W., and Bae, J.-W. (2015). Proteobacteria: microbial signature of dysbiosis in gut microbiota. Trends Biotechnol. 33, 496–503. doi: 10.1016/j.tibtech.2015.06.011
Silva, B. R. D. S., Derami, M. S., Paixão, D. A., Persinoti, G. F., Dias Da Silveira, W., and Maluta, R. P. (2020). Comparison between the intestinal microbiome of healthy fish and fish experimentally infected with Streptococcus agalactiae. Aquac. Res. 51, 3412–3420. doi: 10.1111/are.14676
Simões, L. N., Lombardi, D. C., Gomide, A. T. M., and Gomes, L. C. (2011). Efficacy of clove oil as anesthetic in handling and transportation of Nile tilapia, Oreochromis niloticus (Actinopterygii: Cichlidae) juveniles. Zoologia 28, 285–290. doi: 10.1590/S1984-46702011000300001
Souza, F. P. D., Lima, E. C. S. D., Urrea-Rojas, A. M., Suphoronski, S. A., Facimoto, C. T., Bezerra Júnior, J. D. S., et al. (2020). Effects of dietary supplementation with a microalga (Schizochytrium sp.) on the hemato-immunological, and intestinal histological parameters and gut microbiota of Nile tilapia in net cages. PLoS One 15:e0226977. doi: 10.1371/journal.pone.0226977
Standen, B., Rawling, M., Davies, S., Castex, M., Foey, A., Gioacchini, G., et al. (2013). Probiotic Pediococcus acidilactici modulates both localised intestinal-and peripheral-immunity in tilapia (Oreochromis niloticus). Fish Shellfish Immunol. 35, 1097–1104. doi: 10.1016/j.fsi.2013.07.018
Sugita, H., Miyajima, C., and Deguchi, Y. (1991). The vitamin B12-producing ability of the intestinal microflora of freshwater fish. Aquaculture 92, 267–276. doi: 10.1016/0044-8486(91)90028-6
Suphoronski, S., Chideroli, R., Facimoto, C., Mainardi, R., Souza, F., Lopera-Barrero, N., et al. (2019). Effects of a phytogenic, alone and associated with potassium diformate, on tilapia growth, immunity, gut microbiome and resistance against Francisellosis. Sci. Rep. 9:6045. doi: 10.1038/s41598-019-42480-8
Sylvain, F. -É., and Derome, N. (2017). Vertically and horizontally transmitted microbial symbionts shape the gut microbiota ontogenesis of a skin-mucus feeding discus fish progeny. Sci. Rep. 7:5263. doi: 10.1038/s41598-017-05662-w
Tsuchiya, C., Sakata, T., and Sugita, H. (2008). Novel ecological niche of Cetobacterium somerae, an anaerobic bacterium in the intestinal tracts of freshwater fish. Lett. Appl. Microbiol. 46, 43–48. doi: 10.1111/j.1472-765X.2007.02258.x
Wagner, M. R., Roberts, J. H., Balint-Kurti, P. J., and Holland, J. B. (2020). Microbiome composition differs in hybrid and inbred maize. bioRxiv [Preprint]. doi: 10.1101/2020.01.13.904979
Wei, T., and Simko, V. (2017). R package “corrplot”: Visualization of a correlation matrix (Version 0.84). Available at: https://github.com/taiyun/corrplot (Accessed April 28, 2021).
Wein, T., Dagan, T., Fraune, S., Bosch, T. C., Reusch, T. B., and Hülter, N. F. (2018). Carrying capacity and colonization dynamics of Curvibacter in the hydra host habitat. Front. Microbiol. 9:443. doi: 10.3389/fmicb.2018.00443
Wickham, H. (2011). ggplot2. Wiley Interdiscip. Rev. Comput. Stat. 3, 180–185. doi: 10.1002/wics.147
Wiles, T. J., Jemielita, M., Baker, R. P., Schlomann, B. H., Logan, S. L., Ganz, J., et al. (2016). Host gut motility promotes competitive exclusion within a model intestinal microbiota. PLoS Biol. 14:e1002517. doi: 10.1371/journal.pbio.1002517
Wynne, J. W., Thakur, K. K., Slinger, J., Samsing, F., Milligan, B., Powell, J. F. F., et al. (2020). Microbiome profiling reveals a microbial dysbiosis during a natural outbreak of Tenacibaculosis (yellow mouth) in Atlantic Salmon. Front. Microbiol. 11:586387. doi: 10.3389/fmicb.2020.586387
Yukgehnaish, K., Kumar, P., Sivachandran, P., Marimuthu, K., Arshad, A., Paray, B. A., et al. (2020). Gut microbiota metagenomics in aquaculture: factors influencing gut microbiome and its physiological role in fish. Rev. Aquac. 12, 1903–1927. doi: 10.1111/raq.12416
Zar, J. H. (2014). Spearman rank correlation: overview. Wiley StatsRef: Statistics Reference Online.
Zhang, M., Sun, Y., Liu, Y., Qiao, F., Chen, L., Liu, W.-T., et al. (2016). Response of gut microbiota to salinity change in two euryhaline aquatic animals with reverse salinity preference. Aquaculture 454, 72–80. doi: 10.1016/j.aquaculture.2015.12.014
Zhang, Z., Yu, Y., Jiang, Y., Wang, Y., Liao, M., Rong, X., et al. (2020). The intestine of artificially bred larval turbot (Scophthalmus maximus) contains a stable core group of microbiota. Arch. Microbiol. 202, 2619–2628. doi: 10.1007/s00203-020-01984-y
Keywords: breeding, Nile tilapia, microbiome, 16S amplicon, whole-genome sequencing, core microbiome
Citation: Abdelhafiz Y, Fernandes JMO, Larger S, Albanese D, Donati C, Jafari O, Nedoluzhko AV and Kiron V (2021) Breeding Strategy Shapes the Composition of Bacterial Communities in Female Nile Tilapia Reared in a Recirculating Aquaculture System. Front. Microbiol. 12:709611. doi: 10.3389/fmicb.2021.709611
Edited by:
Haidong Yao, Karolinska Institutet (KI), SwedenReviewed by:
Hetron Mweemba Munang’andu, Norwegian University of Life Sciences, NorwayHien Van Doan, Chiang Mai University, Thailand
Copyright © 2021 Abdelhafiz, Fernandes, Larger, Albanese, Donati, Jafari, Nedoluzhko and Kiron. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Viswanath Kiron, kiron.viswanath@nord.no