Mastitis, the inflammation of mammary gland, is one of the most prevalent diseases affecting dairy cows in worldwide [1] . Mastitis can be classified into three major types: Clinical Mastitis (CM), Sub-Clinical Mastitis (SCM) and Chronic Mastitis [2] [3] [4] . Among them, there is one common method for diagnosing CM by abnormalities including flakes, pus and the altering in color in the milk [5] . Because fungi can cause infections in animals [6] , new methods of diagnosing CM have also been found from the perspective of observing the fungal community.

Fungi may sporadically infect cows through the injuries of the mammary gland caused by milking machine [7] . Extensive and indiscriminate use of antibiotics for treatment of mastitis can also make it easier for cows to encounter fungi invasion [8] . According to previous literature data, cows which are infected by fungi account for 2% - 13% of all the cows with mastitis [9] [10] . The incidence of mastitis due to fungi has also been reported by various workers [11] [12] .

Many relevant studies have been focused on the fungal community in the cows with CM. Spanamberg [6] has found that Candida species are the most frequent organisms among the fungal community isolated from infected glands in cows with CM. Among them, the occurrence of Candida tropicalis is most frequent in clinical cases of mastitis [6] [11] . Besides, some researchers study the cows with CM and come to the conclusion that Aspergillus spp. was one of the causative fungi pathogens [13] [14] . Geotrichum candidum, Trichosporon cutaneum, Saccharomyces cerevisiae and Rhodotorula rubra have also been detected in the cows with CM [15] .

In the literature published previously, the fungal community of the raw milk [15] and the rumen fluid of cows with CM have already been studied [16] . Despite researchers are aware of the severity of the Clinical Mastitis caused by fungi, few studies have been conducted in the cows to better understand how CM can affect the gut fungal community. It is well known that the gastrointestinal community plays an important role in maintaining host health [17] . And the dysbiosis of the gut fungal community can lead to high incidence of CM in cows and cause huge economic loss in milk yield [18] [19] [20] . Therefore, the fungal community of gut in cows with CM should be described in detail for further diagnosis. The present study is undertaken to detect the difference of gut fungal community between healthy cows and individuals with CM.

The Shannon Wiener curves generated by MOTHUR plotted the number of qualified sequencing reads. It reflected the gut fungal diversity of samples. The curves tend to approach a horizontal asymptote, which indicated the sequencing effort was saturated. We classified cows named CM1, CM2 into the case group. Besides, the control group included four cows named H1, H2, H3, and H4. As was depicted in Figure 1, the diversity of the gut fungal community in the case group was lower than that in the control group.

Detailed information of six cows from which fecal samples were collected was presented (Table 1). We could get the times of calving of six cows from the row

of parity clearly. The total number of OTUs, sequences, coverage and statistical estimates of species richness and diversity of gut fugal community were presented in Table 2. The Chao and Ace indices showed that the richness of gut fungal community in the case group was lower than the control group. What’s more, Simpson and Shannon indices confirmed that there was lower level of gut fungal community diversity in the case group.

The rank-abundance distribution curve could be used to visualized species richness and species evenness. The curves of the case group were kept in a steep gradient, indicating low evenness as the high-ranking species have much higher abundance than the low-ranking species. Conversely, it was also clear that the OTU rank of CM1 and CM2 were lower, the curves of the control group was kept in a shallow gradient, indicating high evenness as the abundance of different species were similar (Figure 2).

All sequences were classified in phylum and genus using the program Mothur with the default setting. At the phylum level, the relative abundance of Ascomycota and Zygomycota in the case group was similar to the control group. Ascomycota was dominant both in the case group and the control group. However,

the relative abundance of Basidiomycota was lower and the relative abundance of unclassified fungi was higher in the case group compared with the control group (Figure 3(a), Figure 3(b)).

Main fungi genus have been detected and heatmap was generated to graphically represent the relative abundance of gut fungal genera. Pseudallescheria was dominated followed by Microascaceae-unclassified in the control group. In the case group, the proportion of many unclassified species which belonged to Saccharomycetales were in highest in the case group. It may suggest that there may be a wide range of Saccharomycetales disorders in the gut fungal community of cows with CM. What’s more, the relative abundance of Pseudallescheria, Microascaceae-unclassified, Scedosporium and Trichosporon were lower than the control group (Figure 4(a), Figure 4(b)). From the figure presented above, the diversity of gut fungal community in the case group was lower than the control group at the genus level.

Based on the principal component analysis, the fungal communities of healthy cows and cows with CM showed distinction (Figure 5(a)). The fungal communities of CM1, CM2 clustered closely and didn’t cluster with H1, H2, H3, H4. As the results revealed, H1 wasn’t close to the other three healthy cows H2, H3, H4. We attempted to make Hierarchical clustering analysis according to beta diversity distance matrix. And then use unweighted pair group method with arithmetic mean to build tree structure. The length of branch represented the

distance between samples and vertical bar in the end represent clustering [21] (Figure 5(b)). In conclusion, Principal coordinates analysis and the multiple

samples similarity tree showed that the gut fungal community of the case group was clearly distinguishable from that of the control group at the level of the OUT (Figure 5(a), Figure 5(b)).

Product safety and animal welfare in dairy management system can be improved by means of early detection of diseases such as mastitis [22] . Since the heritability for CM was extremely low [23] , hardly can we predict the occurrence of CM according to genetic heritability. In addition, somatic cell counts (SCC) becomes the most common measurement of milk quality and udder health [24] . SCC is generally both an indicator of CM. It has been indeed demonstrated that mastitis causes an increase in SCC of cattle [25] . Though it is easy to identify Clinical Mastitis, more precise detection tests need to be explored.

In our study, firstly we found the difference of intestinal fungal community between the case group and the control group without limitation of observing bovine mammary gland tissue. Compared with the control group, the gut fungal community showed less diversity and richness in the case group. The result can be explained by the different relative abundance both in phylum and genera.

In phylum, the relative abundance of Basidiomycota in the case group was lower than the control group. The relative abundance of Ascomycota and Zygomycota in the case group was similar to the control group, which showed that the relative abundance of Ascomycota and Zygomycota were stable.

In genus, in comparison to the control group, the relative abundance of Saccharomycetales-unclassified was higher. On the contrary, the relative abundance of Pseudallescheria, Microascaceae-unclassified, Scedosporium and Trichosporon were lower than the control group. In the present work, it is found that a large number of unclassified genera and phyla in the feces of cows belong to unclassified, since little work on gut fungal community has been done. Although the incidence of CM caused by Candida spp. was high [6] , the relative abundance was lower in the case group than the control group. It could be concluded that another fungi was the main pathogeny instead of Candida. Among them, unclassified species which belonged to Saccharomycetales was dominant, showing that it may be the core genus in the case group.

Nowadays, papers about mastitis and its optimal technical management and treatment are available [26] [27] . The prevalence of mastitis in organized dairy farm is decreased by improving preventive measures and using antibiotics. However, it is still of a challenge as many of these fungi do not respond to the antibiotics rather they use some of the antibiotic like tetracyclineas their source of energy [28] . Hence, it is important to study out new treatment according to the specific fungal community in cows with CM mentioned above.

This study used high-throughput sequencing technology is more comprehensive and systematic to measure the fungi in the feces of Chinese Holstein cows. The results suggest that Clinical Mastitis can cause the shifts of gut fungal community, showing that fungi diversity and abundance in the feces of Chinese Holstein cows with CM are lower than healthy cows with CM under the same diet.

The diversity and richness of gut fungal community in the case group are lower than those in the control group. It is found that Clinical Mastitis gives rise to gut fungal community disorder, which allows us to detect CM according to the microbial difference in phyla and genera. As a consequence, fungal detection of CM can be helpful to recover benefits of dairy cooperation and reduce the mortality of cows caused by CM. However, the relative abundance of unclassified fungi is predominant in the case group both in phylum and in genus. Thus, further study on fungal mastitis pathogens and relevant methods for treatment needs to be continued.

The farmers treating cows suffering from clinical mastitis are particularly thanked for giving me access to data and helping us collect samples. Ilumina MiSeq platform is gratefully acknowledged for providing competent contributions to sequencing. This work was supported by undergraduate Training Programs for Innovation and Entrepreneurship of Jiangsu Province (201710320106X), to which we also show gratitude.

Wen, H.Y., Lu, C.X., Yuan, Z.Y., Wang, X.Y. and Su, S.T. (2018) Analysis of Gut Fungal Community of Cows with Clinical Mastitis. Advances in Microbiology, 8, 366-377. https://doi.org/10.4236/aim.2018.85024