Journal of Oceanology and Limnology   2021, Vol. 39 issue(5): 1961-1971     PDF
Institute of Oceanology, Chinese Academy of Sciences

Article Information

ZHOU Lei, HAN Yaoquan, WANG Dapeng, LI Yusen, HUANG Xiande, HE Anyou
Comparison of fungal community composition within different intestinal segments of tilapia and bighead carp
Journal of Oceanology and Limnology, 39(5): 1961-1971

Article History

Received May. 30, 2020
accepted in principle Oct. 12, 2020
accepted for publication Nov. 6, 2020
Comparison of fungal community composition within different intestinal segments of tilapia and bighead carp
Lei ZHOU1,2, Yaoquan HAN1, Dapeng WANG1, Yusen LI1,2, Xiande HUANG2, Anyou HE1     
1 Guangxi Key Laboratory of Aquatic Genetic Breeding and Healthy Aquaculture, Guangxi Academy of Fishery Sciences, Nanning 530021, China;
2 Joint Laboratory of Guangdong Province and Hong Kong Region on Marine Bioresource Conservation and Exploitation, College of Marine Sciences, South China Agricultural University, Guangzhou 510642, China
Abstract: Although intestinal fungi play important roles in host health and disease, the composition and diversity of fungal communities remain poorly reported in fish. In this study, fungi in the fore-, mid-, and hindintestine of tilapia (Oreochromis mossambicus) and bighead carp (Aristichthys nobilis) from Hongchaojiang Reservoir in Guangxi, China were investigated by ITS sequencing. Based on this, we obtained 1 763 478 high-quality tags, which clustered into 1 089 operational taxonomic units (OTUs). In total, 404 OTUs were annotated, of which 310, 68, and 26 belonged to Ascomycota, Basidiomycota, and other, respectively. Results show significant differences in the community composition of intestinal fungi between tilapia and bighead carp but not within their different intestinal segments. Furthermore, 154 of the 404 annotated OTUs were considered reliable and were classified into three trophic modes and nine guilds. The three trophic modes consisted of 108 OTUs of saprotrophic fungi, 41 OTUs of pathotrophic fungi, and five OTUs of symbiotrophic fungi. The top three most abundant OTUs overall (i.e., Otu000002, Scopulariopsis acremonium; Otu000018, Alternaria palandui; Otu000034, Aureobasidium pullulans) showed lower abundance in the hind-intestinal segments of bighead carp, suggesting uneven distribution of these fungi in this species. In addition, saprotrophic and pathotrophic fungi were markedly decreased in the hindintestine. It is indicated that the fungal community was not only related to host species specificity but also to the respective physiological functions of different intestinal segments. These findings provide valuable information on the composition, structure, and potential function of the intestinal fungi community associated with different intestinal segments in tilapia and bighead carp under natural conditions, thus highlighting the importance of fungi as an integral part of the intestinal microbiota in maintaining host health.
Keywords: fungi    intestinal microbiome    tilapia    bighead carp    Internal Transcribed Spacer (ITS) sequencing    

Host-associated microorganisms are vital for host species growth and survival (Eichmiller et al., 2016). In particular, the vertebrate intestinal tract is a complex ecosystem that provides optimal habitat for dynamic microorganism communities. These communities, in turn, play critical roles in host digestion, energy, and health (Lynch and Pedersen, 2016; Rooks and Garrett, 2016) and display complex microbial-microbial and host-microbial relationships (Wang et al., 2018).

Representing more than half of all vertebrate species, fish are an important class of animal for testing mutualistic symbiosis of microorganisms with their hosts (Talwar et al., 2018). Intestinal microbe often participates in the physiological and biological functions of their fish hosts (Rooks and Garrett, 2016; Blander et al., 2017). Although such microorganisms include many categories, including bacteria, viruses, and fungi, most of our knowledge on the intestinal microbiota derives from studies on bacterial studies, with few reports on viruses and fungi. As primary decomposers of organic material, fungi including molds and yeasts are a diverse and ubiquitous group of eukaryotic organisms that contain cell walls but no chlorophyll or plastids. They inhabit both terrestrial and aquatic environments, including within the bodies of vertebrates (e.g. intestinal segments) and substantially contribute to global biogeochemical cycling (Gadd, 2007; Nazir et al., 2017). Just like intestinal bacteria, fungi are an important part of the fish microbial community and play important roles in host health and disease. In aquaculture, various yeasts have been identified and studied as a part of the normal microbiota in fish intestines (Gatesoupe, 2007; Nayak, 2010; Sarlin and Philip, 2011; Øverland et al., 2013; Caruffo et al., 2015). Thus, understanding the composition and diversity of intestinal fungal communities is essential to promote fish health and productivity. However, previous studies on intestinal fungi in fish have primarily focused on culture-based methods, with the composition and diversity of fungal communities rarely reported.

Tilapia (Oreochromis mossambicus) is a tropical, salt-tolerant, omnivorous fish native to Africa, and it shows a good resistance to hypoxia and strong adaptability. Bighead carp (Aristichthys nobilis) is a unique filter-feeding freshwater fish native to China. Both species are considered economically important and are commonly found in rivers, lakes, and reservoirs of China. Although, intestinal bacteria in tilapia and bighead carp have been researched intensively (Borsodi et al., 2017; Wang et al., 2017; Hassaan et al., 2018; Li et al., 2018a, b, c; Zheng et al., 2018; Yu et al., 2019), studies on the composition, diversity, and role of intestinal fungi in these two fish remain limited. Such studies would improve our understanding of the composition and structure of microbes in fish intestines, as well as the ecological function of fish in different water bodies and their healthy ecological farming.

In this study, the intestinal fungi of tilapia and bighead carp from the same habitat (Hongchaojiang Reservoir, Guangxi, China) were investigated by Internal Transcribed Spacer (ITS) high-throughput sequencing. We hypothesized that the intestinal fungi would exhibit significant differences in composition and structure between these two fish, i.e., species specificity, and among different intestinal segments related to their respective intestinal physiological functions. Thus, we aim to examine and compare the composition and diversity of fungal communities in different intestinal segments of tilapia and bighead carp; and clarify the functions and roles of fungi in the fish intestine. This research shall provide a theoretical basis for further studies on intestinal microbiota.

2 MATERIAL AND METHOD 2.1 Animal and sample preparation

Using net cages (mesh size 0.7 cm), three individual fish of each species (average weight of ~0.5 kg) were captured randomly. After the anesthetization with tricaine methane sulfonate (MS-222) (60 mg/L, immersion for 2–3 min), the intestines of each fish were removed aseptically from the abdominal cavity using a sterile scalpel. Microbial samples were collected from the fore-, mid- and hind-intestinal segments by gently abrading the intestinal wall with a sterile spatula, with the obtained contents placed into separate 2.0-mL cryo tubes. Based on the six individual fish and three intestinal samples, we obtained 18 microbial samples in total. Upon collection, all samples were rapidly frozen in liquid nitrogen followed by storage in a-80-℃ refrigerator awaiting DNA extraction. All experiments were conducted under the Guide for the Care and Use of Laboratory Animals in China. The procedures were approved by the Animal Ethics Committee at the South China Agricultural University, China (approval ID: SYXK-2019-0136).

2.2 DNA extraction, polymerase chain reaction (PCR), and sequencing

In accordance with manufacturer provided procedures, DNA was first extracted from the microbial samples using an EZNA stool DNA Kit (Omega Biotek, USA). From the resulting ITS rDNA, the ITS2 region was PCR-amplified with specific ITS primers (KYO2F: GATGAAGAACGYAGYRAA; ITS4R: TCCTCCGCTTATTGATATGC). The PCR procedures (i.e., 95 ℃ for 2 min; 27 cycles at 98 ℃ for 10 s, 62 ℃ for 30 s, and 68 ℃ for 30 s; final extension at 68 ℃ for 10 min) were conducted in 5 μL of 10× KOD buffer, 1 μL of KOD polymerase, 5 μL of 2.5 mmol/L dNTPs, 100 ng of template DNA, 1.5 μL of each primer (5 μmol/L), and ddH2O to a volume of 50 μL (Toyobo, Japan). The purified PCR amplicons were combined in equimolar amounts and sent to Gene Denovo Biological Technology Co., Ltd. (Guangzhou, China) for paired-end sequencing (2×250, Illumina HiSeq 2500, USA).

2.3 Sequence assembly and taxonomic classification

Under specific criteria i.e., 10-bp minimum overlap and 2% mismatch error rate, the raw sequences were cleaned and assembled into raw tags using FLASH (v1.2.11) (Magoč and Salzberg, 2011), followed by QIIME (v1.9.1)-based filtering to acquire clean tags (Caporaso et al., 2010; Bokulich et al., 2013). These clean tags were searched using the reference database, with chimeric tags identified and then removed using QIIME (v1.9.1). The resulting effective tags were then binned into operational taxonomic unit (OTU) clusters (≥97% similarity) with UPARSE (Edgar, 2013). The absolute and relative abundances of the obtained OTUs were determined for the following analyses. All raw data were submitted to the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) under accession number SRP239684.

2.4 ITS fungal community analysis and functional prediction

For each microbial sample, Chao1, Sobs, Shannon, and Simpson indices were calculated to determine the overall (mean) alpha diversity for each sample (with 95% confidence interval), which was then evaluated against all other microbial samples collected from the same site. Beta diversity of samples was also determined and used for principal component analysis (PCA). Nonmetric multidimensional scaling (NMDS) was evaluated and plotted using R software. The Bray-Curtis dissimilarity index was applied for analysis of similarities (ANOSIM) to establish statistical differences in the microbial community and community overlap (i.e., R > 0.75, well-separated community; 0.50 < R≤0.75, separated but overlapping community; 0.25 < R≤0.50, separated but strongly overlapping community; and 0.25≤R, non-separated community) (Buttigieg and Ramette, 2014). Significant differences between samples were determined at P < 0.05 and P < 0.01.

From the above OTU abundances, FUNGuild-based functional annotation of the identified fungi was performed ( In brief, the FUNGuild database was used for functional predictions of the OTUs of interest. From the database, annotations are classified as 'highly probable' or 'probable' and include trophic mode, guild, and growth morphology. The trophic fungal guilds include pathotrophic, symbiotrophic, and saprotrophic (i.e., procuring nutrients by damaging live host cells, exchanging resources with live host cells, or destroying dead host cells, respectively) (Nguyen et al., 2016). Species can be classified as single or multi-mode. Guild annotations include animal pathogens, fungi (e.g., arbuscular mycorrhizal, ectomycorrhizal, ericoid mycorrhizal, lichenicolous, lichenized, and foliar endophytic fungi), mycoparasites, and plant pathogens (e.g., root endophytes, saprotrophs, and wood saprotrophs). Species can be single or multi-guild. Growth morphology annotations include yeast, facultative yeast, and thallus.

3 RESULT 3.1 Fungal composition in tilapia and bighead carp intestines

A total of 1 763 478 high-quality sequences (tags) were obtained from the 18 samples. The number of sequences per sample ranged from 81 827 to 121 770, which clustered into 1 089 OTUs. As shown in the Venn diagrams, we obtained a total of 667 OTUs for tilapia and 482 OTUs for bighead carp (Fig. 1c & Supplementary Table S1), of which 283 are overlapped between the two species (Fig. 1c) and 44 overlapped in the fore-, mid-, and hind-intestine of tilapia and bighead carp (Fig. 1d). In detail, we identified 288, 439, and 496 OTUs for the fore-, mid-, and hind-intestine of tilapia (OsFF, OsMF, and OsHF), respectively, and 242, 245, and 374 OTUs for the fore-, mid-, and hind-intestine of bighead carp (AnFF, AnMF and AnHF), respectively (Fig. 1a & b). A summary of the OTUs and tags of all samples are shown in Supplementary Table S1. The Venn diagrams (Fig. 1c & d) indicated that many intestinal fungi were shared not only within the different intestinal segments of one fish but also among the different fish. Among the 44 overlapping OTUs (Fig. 1d), 14 were annotated, including 13 belonging to Ascomycota and one belonging to Rozellomycota.

Fig.1 Venn diagrams showing OTUs of intestinal fungi in tilapia and bighead carp a. OTUs in fore-, mid-, and hind-intestine in tilapia (OsFF, OsMF, and OsHF, respectively); b. OTUs in fore-, mid-, and hind-intestine in bighead carp (AnFF, AnMF, and AnHF, respectively); c. total intestinal fungal OTUs of tilapia (Os) and bighead carp (An); d. flower diagram displaying fore-, mid-, and hind-intestine OTU overlap in tilapia and bighead carp.

Taxonomy stack distributions were determined to analyze the composition of each sample at the phylum and genus levels. Among the annotated OTUs, the top 10 with abundances ≥2% in the sample were selected, with others classified into the "other" category and unannotated OTUs classified into the "unclassified" category (Fig. 2). At the phylum level (Fig. 2a), the taxonomic compositions of OsFF, OsMF, and OsHF showed relatively high average abundances of Ascomycota (48.93%±19.20%, 26.88%±2.02%, and 33.40%±9.91%, respectively), followed by Rozellomycota (0.26%±0.38%, 0.60%±0.59%, and 0.38%±0.33%, respectively) and Basidiomycota (0.82%±0.77%, 2.46%±1.79%, and 1.21%±0.61%, respectively), whereas the composition in AnFF, AnMF, and AnHF showed a very high dominance of Ascomycota (39.30%±15.32%, 31.39%±23.70%, and 8.15%±0.47%, respectively), followed by Rozellomycota (0.06%±0.06%, 3.26%±4.60%, and 4.26%±1.04%, respectively) and Basidiomycota (1.46%±1.70%, 0.27%±0.27 and, 0.06%±0.05%, respectively). At genus level (Fig. 2b), the taxonomic compositions of OsFF, OsMF, and OsHF and AnFF, AnMF, and AnHF showed relatively high average abundances of Scopulariopsis (12.68%±8.06%, 2.03%±1.59%, 5.62%±4.53%, 7.07%±4.71%, 8.57%±7.39%, and 0.39%±0.40%, respectively), Chaetomium (4.90%±3.51%, 1.84%±1.07%, 0.65%±0.71%, 1.08%±1.48%, 2.42%±4.19%, and 0.11%±0.17%, respectively), and Alternaria (2.76%±0.54%, 1.97%±1.57%, 0.49%±0.25%, 0.60%±0.82%, 1.61%±1.80%, and 0.37%±0.46%, respectively).

Fig.2 Relative fungal abundances and species composition at phylum (a) and genus (b) levels based on taxonomy stack distributions Only top 10 species (abundances ≥2% in at least one sample) are displayed, with remaining species categorized into "Other". Unannotated tags were categorized into "Unclassified".

It is worth noting that, Alternaria (Dothideomycetes), Scopulariopsis (Sordariomycetes), and Chaetomium (Sordariomycetes) all belong to Ascomycota. Among the 1 089 fungal OTU sequences, 404 were annotated, including 310 in Ascomycota, 68 in Basidiomycota, three in Chytridiomycota, one in Glomeromycota, six in Mortierellomycota, three in Mucoromycota, and 13 in Rozellomycota. The dominant fungi in the intestines of tilapia and bighead carp were Ascomycota at the phylum level, which dominated in most segments and accounted for approximately one third of relative fungal abundances, and Scopulariopsis (Ascomycota) at the genus level.

3.2 Diversity and potential functions of fungi in intestines of tilapia and bighead carp

For alpha diversity (Chao1, Sobs, Shannon and Simpson indices), no significant differences were observed between the two species (Fig. 3a) or within tilapia (OsFF, OsMF, and OsHF) (Fig. 3b). Within bighead carp (AnFF, AnMF, and AnHF), however, fungi in the hind-intestine differed markedly from those in the fore- and mid-intestine (Fig. 3b).

Fig.3 Differences in alpha diversity of intestinal fungi between (a) and within (b) tilapia and bighead carp Four alpha diversity indices (i.e., Chao1, Sobs, Shannon, and Simpson) were determined based on 97% identity thresholds. Groups were compared with t-tests.

For beta diversity, the NMDS (Fig. 4) results show that tilapia and bighead carp are closely clustered but exhibit considerable differences in intestinal fungi. Pairwise ANOSIM (Table 1) indicates no significant differences among the fore-, mid-, and hind-intestinal segments within tilapia and bighead carp (P > 0.05). The fungal communities of tilapia and bighead carp were significantly different but strongly overlapping (ANOSIM-R=0.415, P=0.001), which is consistent to the Venn diagram results, and shows that tilapia and bighead carp shared many intestinal fungi.

Fig.4 Beta diversity of intestinal fungi within and between tilapia and bighead carp Bray-Curtis dissimilarity index and nonmetric multidimensional scaling (NMDS) were used by dividing samples into tilapia and bighead carp groups. When stress=0.11, NMDS accurately reflected differences in fungal community between two groups.
Table 1 Pairwise ANOSIM of different groups based on Bray-Curtis dissimilarity index

Based on the FUNGuild database, only 154 of the 1 089 OTU sequences (404 annotated OTUs) were considered reliable (probable or highly probable) and were classified into three trophic modes and nine guilds (Supplementary Table S2). The three trophic modes consisted of 108 OTUs of saprotrophic fungi, 41 OTUs of pathotrophic fungi, and five OTUs of symbiotrophic fungi. Among the 108 OTUs of saprotrophic fungi, 80 belonged to 31 known genera of Ascomycota, 25 belonged to 22 known genera of Basidiomycota, and three belonged to two known genera of Mucoromycota. In addition, 83 pathotrophic fungi were undefined pathotrophic, 14 were wood saprotrophic, nine were dung saprotrophic, and two were soil saprotrophic. Among the 41 OTUs of pathotrophic fungi, 35 belonged to 24 known genera of Ascomycota and six belonged to five known genera of Basidiomycota. In addition, 23 OTUs of pathotrophic fungi were from plant pathogens and three were from animal pathogens, which may have entered the intestines via the food consumed by the fish. Among the five OTUs of symbiotrophic fungi, two belonged to two known genera of Ascomycota, two belonged to two known genera of Basidiomycota, and one belonged to one unclassified genus of Glomeromycota.

For the fore-, mid- and hind-intestinal segments of tilapia and bighead carp, the majority of taxa (relative abundance) were classified as dung saprotrophic-undefined saprotrophic-wood saprotrophic (17.57%±10.85%, 3.87%±2.64%, 6.27%±4.05%, 8.15%±4.19%, 10.98%±10.21%, and 0.50%±0.56%, respectively), undefined saprotrophic (7.08%±5.83%, 5.99%±5.42%, 5.27%±4.11%, 8.08%±0.67%, 6.84%±4.58%, 1.20%±0.53%, and respectively), and animal pathogen-endophyte-plant pathogen-wood saprotrophic guilds (2.76%±0.54%, 1.97%±1.57%, 0.49%±0.25%, 0.60%±0.82%, 1.61%±1.80%, and 0.37%±0.46%, respectively) (Supplementary Fig. S1). For distribution of the relative abundance mentioned above, a heatmap was constructed based on guild assignments for the top 15 fungal species determined using FUNGuild. From top to bottom, the abundance of guild assignments decreased gradually. As shown in Fig. 5a, there were relatively smaller differences among OsFF, OsMF, and OsHF, but relatively greater differences among AnFF, AnMF, and AnHF. This was especially obvious for the top four guild assignments (dung saprotrophic-undefined saprotrophic-wood saprotrophic, undefined saprotrophic, animal pathogen-endophyte-plant pathogen-wood saprotrophic guilds). The PCA results for the predicted function of fungi in tilapia and bighead carp showed that fungi in the fore-, mid- and hind-intestinal segments functioned differently (Fig. 5b & c). For specific annotated OTUs, the three most abundant (i.e., Otu000002, Ascomycota; Sordariomycetes; Microascales; Microascaceae; Scopulariopsis; Scopulariopsis acremonium; Otu000018, Ascomycota; Dothideomycetes; Pleosporales; Pleosporaceae; Alternaria; Alternaria palandui; Otu000034, Ascomycota; Dothideomycetes; Dothideales; Dothioraceae; Aureobasidium; Aureobasidium pullulans) were classified as dung saprotrophic (Otu000002) and pathotrophic-saprotrophic-symbiotrophic (Otu000018 and Otu000034), and all showed lower abundance in the hind intestine of bighead carp, in accordance with the taxonomy stack distributions. Notably, compared to the above three OTUs, many of the other annotated OTUs were detected in only one of the three segments in tilapia and bighead carp, suggesting that these corresponding fungi might exhibit specific colonization.

Fig.5 The heatmap and principal component analysis (PCA) of intestinal fungi in tilapia and bighead carp a. heatmap of relative abundances of top 15 functions of intestinal fungi in tilapia and bighead carp. Values are normalized by z-score method. Guild assignments were determined using FUNGuild. PCA plots of fore-, mid-, and hind-intestinal fungal functions of tilapia (b) and bighead carp (c).

In overall, our results show that, for reliable OTUs, more than two thirds were saprotrophic fungi, less than one third were pathotrophic fungi, and only 3.2% were symbiotrophic fungi; and the distribution of fungi in the intestine was uneven, especially in the hind-intestine of bighead carp, where saprotrophic and pathotrophic fungi decreased sharply, suggesting that they may be involved in digestion and absorption.


In this study, Ascomycota fungi (phylum level) were dominant in the intestines of tilapia (48.93%±19.20%) and bighead carp (39.30%±15.32%). Previous research on the wood-eating fish Panaque nigrolineatus also showed Ascomycota (i.e., Sordariomycetes and Dothideomycetes) to be dominant in the gut (Marden et al., 2017). Ascomycota fungi are heterotrophs that obtain nutrients from dead or living organisms, often forming symbiotic relationships with algae, plants, and even arthropods (Yang et al., 2012). The Ascomycota accounts for more than 75% of all fungi and contain a wide variety of organisms, including yeasts, penicillin (Penicillium chrysogenum), morels, truffles, and most animal and plant pathogens, and thus play a large role in the recycling of dead plant material (Schoch et al., 2009; Aranda, 2016; Nguyen et al., 2017; Teixeira et al., 2017; Wijayawardene et al., 2018). For the potential functions of fungi based on the FUNGuild database, most reliable OTUs (> 2/3) were saprotrophic fungi. Saprotrophic fungi primarily regulate carbon, nitrogen, and phosphorus in nutrient cycling, and represent key agents of nutrient redistribution in ecosystems (Crowther et al., 2012; Talbot et al., 2013). Notably, although S. acremonium is reported to cause invasive infection in humans (Ellison et al., 1998), it has been identified as a dominant species in fore-, mid- and hind-intestinal segments in tilapia and bighead carp, suggesting that it may play a role in nutrition and absorption in the intestine due to its saprophytic ability. Thus, given its dominance, further studies on the effects of S. acremonium on intestinal function are required.

Tilapia and bighead carp showed significant differences in the composition and structure of the intestinal fungal microbiota, which may be due to the different feeding habits or trophic level of these two fish. From a feeding habit/diet perspective, tilapia fish are omnivorous, inhabits the lower and middle water layers, and dominated by plant-based bait under natural conditions such as the lakes and reservoirs. In contrast, filter-feeding bighead carp reside in the upper and middle water layers, and mainly eat zooplankton and protozoa such as rotifers, branch horns, copepods. Therefore, it is not surprising that the tilapia fish had more OTU and higher fungal richness, which supports the view that intestinal microbial diversity usually increases sequentially corresponding to carnivores, herbivores and omnivores (Kashinskaya et al., 2018). From the perspective of function, different intestinal segments in fish exhibit different functions such as digestion and absorption (Denstadli et al., 2004; Le et al., 2019; Lin et al., 2019; Zhou et al., 2019). In this study, whether for alpha diversity or distribution of abundant specific OTUs, no significant differences in fungi was found among the three intestinal segments, but obvious differences were found between the hind-intestine and the fore- and mid-intestines in bighead carp. For example, compared to that in the fore- and mid-intestines, saprotroph abundance decreased drastically in the hind-intestine of bighead carp. These results suggest that fungal community was not only related to host species specificity but also to the respective physiological functions of different intestinal segments.

It is worth mentioning that 41 OTUs were annotated as probable or highly probable pathotrophic fungi in this study, implying that the intestinal mucosa may be a potential site of pathogenic invasion and colonization by fungi. Although the proportion of fungi/mycobiome in the microbiome is relatively small, i.e., less than 0.1% in the human microbiome (Qin et al., 2010), as microbial saprotrophs, pathogens and mutualists, the fungal community plays a vital role in the health and disease of animals, plants, and ecosystems (Cui et al., 2013; Seed, 2014; Peay et al., 2016; Nash et al., 2017). Thus, our study highlights the importance of investigating the diversity of fungi and their relationship with host health.


Community composition and structure of the intestinal fungi of tilapia and bighead carp were comprehensively analyzed based on ITS sequencing. Our results demonstrated significant differences in the community of intestinal fungi in tilapia and bighead carp but not among their different intestinal segments. For the potential functions of OTUs in the two fish species, saprotrophic fungi were dominant, followed by pathotrophic fungi, although distribution were uneven, especially in the hind-intestine of bighead carp. This study should improve our understanding of the composition and function of fungi as an important part of the intestinal microbiota of tilapia and bighead carp.


All raw data were submitted to the NCBI Sequence Read Archive (SRA) under accession number SRP239684. The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Electronic supplementary material

Supplementary material (Supplementary Tables S1–S2 & Fig.S1) is available in the online version of this article at

Aranda E. 2016. Promising approaches towards biotransformation of polycyclic aromatic hydrocarbons with Ascomycota fungi. Current Opinion in Biotechnology, 38: 1-8. DOI:10.1016/j.copbio.2015.12.002
Blander J M, Longman R S, Iliev I D, Sonnenberg G F, Artis D. 2017. Regulation of inflammation by microbiota interactions with the host. Nature Immunology, 18: 851-860. DOI:10.1038/ni.3780
Bokulich N A, Subramanian S, Faith J J, Gevers D, Gordon J I, Knight R, Mills D A, Caporaso J G. 2013. Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. Nature Methods, 10: 57-59. DOI:10.1038/nmeth.2276
Borsodi A K, Szabó A, Krett G, Felföldi T, Specziár A, Boros G. 2017. Gut content microbiota of introduced bigheaded carps (Hypophthalmichthys spp. ) inhabiting the largest shallow lake in Central Europe. Microbiological Research, 195: 40-50. DOI:10.1016/j.micres.2016.11.001
Buttigieg P L, Ramette A. 2014. A guide to statistical analysis in microbial ecology: a community-focused, living review of multivariate data analyses. FEMS Microbiology Ecology, 90(3): 543-550. DOI:10.1111/1574-6941.12437
Caporaso J G, Kuczynski J, Stombaugh J, Bittinger K, Bushman F D, Costello E K, Fierer N, Peña A G, Goodrich J K, Gordon J I, Huttley G A, Kelley S T, Knights D, Koenig J E, Ley R E, Lozupone C A, McDonald D, Muegge B D, Pirrung M, Reeder J, Sevinsky J R, Turnbaugh P J, Walters W A, Widmann J, Yatsunenko T, Zaneveld J, Knight R. 2010. QIIME allows analysis of high-throughput community sequencing data. Nature Methods, 7: 335-336. DOI:10.1038/nmeth.f.303
Caruffo M, Navarrete N, Salgado O, Díaz A, López P, García K, Feijóo C G, Navarrete P. 2015. Potential probiotic yeasts isolated from the fish gut protect zebrafish (Danio rerio) from a Vibrio anguillarum challenge. Frontiers in Microbiology, 6: 1 093. DOI:10.3389/fmicb.2015.01093
Crowther T W, Boddy L, Hefin Jones T. 2012. Functional and ecological consequences of saprotrophic fungus-grazer interactions. The ISME Journal, 6: 1 992-2 001. DOI:10.1038/ismej.2012.53
Cui L J, Morris A, Ghedin E. 2013. The human mycobiome in health and disease. Genome Medicine, 5: 63. DOI:10.1186/gm467
Denstadli V, Vegusdal A, Krogdahl Å, Bakke-McKellep A M, Berge G M, Holm H, Hillestad M, Ruyter B. 2004. Lipid absorption in different segments of the gastrointestinal tract of Atlantic salmon (Salmo salar L.). Aquaculture, 240(1-4): 385-398. DOI:10.1016/j.aquaculture.2004.06.030
Edgar R C. 2013. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nature Methods, 10: 996-998. DOI:10.1038/nmeth.2604
Eichmiller J J, Hamilton M J, Staley C, Sadowsky M J, Sorensen P W. 2016. Environment shapes the fecal microbiome of invasive carp species. Microbiome, 4(1): 44. DOI:10.1186/s40168-016-0190-1
Ellison M D, Hung R T, Harris K, Campbell B H. 1998. Report of the first case of invasive fungal sinusitis caused by Scopulariopsis acremonium: review of scopulariopsis infections. Archives of Otolaryngology-Head & Neck Surgery, 124(9): 1 014-1 016. DOI:10.1001/archotol.124.9.1014
Gadd G M. 2007. Geomycology: biogeochemical transformations of rocks, minerals, metals and radionuclides by fungi, bioweathering and bioremediation. Mycological Research, 111(1): 3-49. DOI:10.1016/j.mycres.2006.12.001
Gatesoupe F J. 2007. Live yeasts in the gut: natural occurrence, dietary introduction, and their effects on fish health and development. Aquaculture, 267(1-4): 20-30. DOI:10.1016/j.aquaculture.2007.01.005
Hassaan M S, Soltan M A, Jarmolowicz S, Abdo H S. 2018. Combined effects of dietary malic acid and Bacillus subtilis on growth, gut microbiota and blood parameters of Nile tilapia (Oreochromis niloticus). Aquaculture Nutrition, 24(1): 83-93. DOI:10.1111/anu.12536
Kashinskaya E N, Simonov E P, Kabilov M R, Izvekova G I, Andree K B, Solovyev M M. 2018. Diet and other environmental factors shape the bacterial communities of fish gut in an eutrophic lake. Journal of Applied Microbiology, 125(6): 1 626-1 641. DOI:10.1111/jam.14064
Le H T M D, Shao X T, Krogdahl Å, Kortner T M, Lein I, Kousoulaki K, Lie K K, Sæle Ø. 2019. Intestinal function of the stomachless fish, Ballan wrasse (Labrus bergylta). Frontiers in Marine Science, 6: 140. DOI:10.3389/fmars.2019.00140
Li J W, Liu G, Li C W, Deng Y L, Tadda M A, Lan L H, Zhu S M, Liu D Z. 2018a. Effects of different solid carbon sources on water quality, biofloc quality and gut microbiota of Nile tilapia (Oreochromis niloticus) larvae. Aquaculture, 495: 919-931. DOI:10.1016/j.aquaculture.2018.06.078
Li X H, Yu Y H, Li C, Yan Q Y. 2018b. Comparative study on the gut microbiotas of four economically important Asian carp species. Science China Life Sciences, 61: 696-705. DOI:10.1007/s11427-016-9296-5
Li X M, Zhu Y J, Ringø E, Wang X G, Gong J L, Yang D G. 2018c. Intestinal microbiome and its potential functions in bighead carp (Aristichthys nobilis) under different feeding strategies. PeerJ, 6: e6000. DOI:10.7717/peerj.6000
Lin K T, Wang W X, Ruan H T, Dai J G, Sun J J, Liu L, Huang X D. 2019. Transcriptome analysis of differentially expressed genes in the fore- and hind-intestine of ovate pompano Trachinotus ovatus. Aquaculture, 508: 76-82. DOI:10.1016/j.aquaculture.2019.04.078
Lynch S V, Pedersen O. 2016. The human intestinal microbiome in health and disease. The New England Journal of Medicine, 375: 2 369-2 379. DOI:10.1056/NEJMra1600266
Magoč T, Salzberg S L. 2011. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics, 27(21): 2 957-2 963. DOI:10.1093/bioinformatics/btr507
Marden C L, McDonald R, Schreier H J, Watts J E M. 2017. Investigation into the fungal diversity within different regions of the gastrointestinal tract of Panaque nigrolineatus, a wood-eating fish. AIMS Microbiology, 3(4): 749-761. DOI:10.3934/microbiol.2017.4.749
Nash A K, Auchtung T A, Wong M C, Smith D P, Gesell J R, Ross M C, Stewart C J, Metcalf G A, Muzny D M, Gibbs R A, Ajami N J, Petrosino J F. 2017. The gut mycobiome of the Human Microbiome Project healthy cohort. Microbiome, 5: 153. DOI:10.1186/s40168-017-0373-4
Nayak S K. 2010. Role of gastrointestinal microbiota in fish. Aquaculture Research, 41(11): 1 553-1 573. DOI:10.1111/j.1365-2109.2010.02546.x
Nazir R, Mazurier S, Yang P, Lemanceau P, van Elsas J D. 2017. The ecological role of type three secretion systems in the interaction of bacteria with fungi in soil and related habitats is diverse and context-dependent. Frontiers in Microbiology, 8: 38,
Nguyen N H, Song Z W, Bates S T, Branco S, Tedersoo L, Menke J, Schilling J S, Kennedy P G. 2016. FUNGuild: An open annotation tool for parsing fungal community datasets by ecological guild. Fungal Ecology, 20: 241-248. DOI:10.1016/j.funeco.2015.06.006
Nguyen T A, Cissé O H, Wong J Y, Zheng P, Hewitt D, Nowrousian M, Stajich J E, Jedd G. 2017. Innovation and constraint leading to complex multicellularity in the Ascomycota. Nature Communications, 8: 14 444. DOI:10.1038/ncomms14444
Øverland M, Karlsson A, Mydland L T, Romarheim O H, Skrede A. 2013. Evaluation of Candida utilis, Kluyveromyces marxianus and Saccharomyces cerevisiae yeasts as protein sources in diets for Atlantic salmon (Salmo salar). Aquaculture, 402-403: 1-7. DOI:10.1016/j.aquaculture.2013.03.016
Peay K G, Kennedy P G, Talbot J M. 2016. Dimensions of biodiversity in the Earth mycobiome. Nature Reviews Microbiology, 14: 434-447. DOI:10.1038/nrmicro.2016.59
Qin J J, Li R Q, Raes J, Arumugam M, Burgdorf K S, Manichanh C, Nielsen T, Pons N, Levenez F, Yamada T, Mende D R, Li J H, Xu J M, Li S C, Li D F, Cao J J, Wang B, Liang H Q, Zheng H S, Xie Y L, Tap J, Lepage P, Bertalan M, Batto J M, Hansen T, Le Paslier D, Linneberg A, Nielsen H B, Pelletier E, Renault P, Sicheritz-Ponten T, Turner K, Zhu H M, Yu C, Li S T, Jian M, Zhou Y, Li Y R, Zhang X Q, Li S G, Qin N, Yang H M, Wang J, Brunak S, Doré J, Guarner F, Kristiansen K, Pedersen O, Parkhill J, Weissenbach J, MetaHIT Consortium, Bork P, Ehrlich S D, Wang J. 2010. A human gut microbial gene catalogue established by metagenomic sequencing. Nature, 464(7285): 59-65. DOI:10.1038/nature08821
Rooks M G, Garrett W S. 2016. Gut microbiota, metabolites and host immunity. Nature Reviews Immunology, 16: 341-352. DOI:10.1038/nri.2016.42
Sarlin P J, Philip R. 2011. Efficacy of marine yeasts and baker's yeast as immunostimulants in Fenneropenaeus indicus: a comparative study. Aquaculture, 321(3-4): 173-178. DOI:10.1016/j.aquaculture.2011.08.039
Schoch C L, Sung G H, López-Giráldez F, Townsend J P, Miadlikowska J, Hofstetter V, Robbertse B, Matheny P B, Kauff F, Wang Z, Gueidan C, Andrie R M, Trippe K, Ciufetti L M, Wynns A, Fraker E, Hodkinson B P, Bonito G, Groenewald J Z, Arzanlou M, de Hoog G S, Crous P W, Hewitt D, Pfister D H, Peterson K, Gryzenhout M, Wingfield M J, Aptroot A, Suh S O, Blackwell M, Hillis D M, Griffith G W, Castlebury L A, Rossman A Y, Lumbsch H T, Lücking R, Büdel B, Rauhut A, Diederich P, Ertz D, Geiser D M, Hosaka K, Inderbitzin P, Kohlmeyer J, Volkmann-Kohlmeyer B, Mostert L, O'Donnell K, Sipman H, Rogers J D, Shoemaker R A, Sugiyama J, Summerbell R C, Untereiner W, Johnston P R, Stenroos S, Zuccaro A, Dyer P S, Crittenden P D, Cole M S, Hansen K, Trappe J M, Yahr R, Lutzoni F, Spatafora J W. 2009. The Ascomycota tree of life: a phylum-wide phylogeny clarifies the origin and evolution of fundamental reproductive and ecological traits. Systematic Biology, 58(2): 224-239. DOI:10.1093/sysbio/syp020
Seed P C. 2014. The human mycobiome. Cold Spring Harbor Perspectives in Medicine, 5: a019810. DOI:10.1101/cshperspect.a019810
Talbot J M, Bruns T D, Smith D P, Branco S, Glassman S I, Erlandson S, Vilgalys R, Peay K G. 2013. Independent roles of ectomycorrhizal and saprotrophic communities in soil organic matter decomposition. Soil Biology and Biochemistry, 57: 282-291. DOI:10.1016/j.soilbio.2012.10.004
Talwar C, Nagar S, Lal R, Negi R K. 2018. Fish gut microbiome: current approaches and future perspectives. Indian Journal of Microbiology, 58: 397-414. DOI:10.1007/s12088-018-0760-y
Teixeira M M, Moreno L F, Stielow B J, Muszewska A, Hainaut M, Gonzaga L, Abouelleil A, Patané J S L, Priest M, Souza R, Young S, Ferreira K S, Zeng Q, da Cunha M M L, Gladki A, Barker B, Vicente V A, de Souza E M, Almeida S, Henrissat B, Vasconcelos A T R, Deng S, Voglmayr H, Moussa T A A, Gorbushina A, Felipe M S S, Cuomo C A, de Hoog G S. 2017. Exploring the genomic diversity of black yeasts and relatives (Chaetothyriales, Ascomycota). Studies in Mycology, 86: 1-28. DOI:10.1016/j.simyco.2017.01.001
Wang A R, Ran C, Ringø E, Zhou Z G. 2018. Progress in fish gastrointestinal microbiota research. Reviews in Aquaculture, 10(3): 626-640. DOI:10.1111/raq.12191
Wang M, Liu G B, Lu M X, Ke X L, Liu Z G, Gao F Y, Cao J M, Zhu H P, Yi M M, Yu D G. 2017. Effect of Bacillus cereus as a water or feed additive on the gut microbiota and immunological parameters of Nile tilapia. Aquaculture Research, 48(6): 3 163-3 173. DOI:10.1111/are.13146
Wijayawardene N N, Hyde K D, McKenzie E H C, Wang Y. 2018. Notes for genera update—Ascomycota: 6822-6917. Mycosphere, 9(6): 1 222-1 234. DOI:10.5943/mycosphere/9/6/11
Yang E C, Xu L L, Yang Y, Zhang X Y, Xiang M C, Wang C S, An Z Q, Liu X Z. 2012. Origin and evolution of carnivorism in the Ascomycota (fungi). Proceedings of the National Academy of Sciences of the United States of America, 109(27): 10 960-10 965. DOI:10.1073/pnas.1120915109
Yu L L, Qiao N Z, Li T Q, Yu R P, Zhai Q X, Tian F W, Zhao J X, Zhang H, Chen W. 2019. Dietary supplementation with probiotics regulates gut microbiota structure and function in Nile tilapia exposed to aluminum. PeerJ, 7: e6963. DOI:10.7717/peerj.6963
Zheng Y, Wu W, Hu G D, Qiu L P, Meng S L, Song C, Fan L M, Zhao Z X, Bing X W, Chen J Z. 2018. Gut microbiota analysis of juvenile genetically improved farmed tilapia (Oreochromis niloticus) by dietary supplementation of different resveratrol concentrations. Fish and Shellfish Immunology, 77: 200-207. DOI:10.1016/j.fsi.2018.03.040
Zhou L, Lin K T, Gan L, Sun J J, Guo C J, Liu L, Huang X D. 2019. Intestinal microbiota of grass carp fed faba beans: a comparative study. Microorganisms, 7(10): 465. DOI:10.3390/microorganisms7100465