1. Introduction
The gut microbiota refers to the diverse microorganisms present in the digestive system of animals, which plays a vital role in animal metabolism, immunity, and reproduction [1,2,3]. Sometimes it is called a “forgotten organ” [4]. For the past few years, depending on the rapid development of bioinformatics, especially the rise of metagenomics, the functions of gut microbiota are being gradually understood [5,6]. The most important function of the intestinal microbiota is the nutritional function, providing energy to the host. Up to 35% of the digestive and metabolic enzymes in mammals are secreted by gut microbiota [7]. The gut microbiota is not static. It is diverse and unstable and is highly susceptible to external environmental influences, such as food [8,9], age [10], disease [11], and living areas [12,13]. The gut microbiota of the same species can vary greatly at different times and in different environments, which can help the host to adapt to its surroundings by influencing host energy metabolism or other aspects.
Recently, numerous studies have revealed that animal reproduction is closely linked to gut microbiota. Clostridium scindens American Type Culture Collection 35,704 can convert glucocorticoids to androgens via side-chain cleavage [14]. The gut microbiota can also be involved in gut metabolism and deglucuronidation of dihydrotestosterone (DHT) and testosterone (T) and results in higher DHT levels in the colon of young healthy mice than in germ-free mice [15]. Accordingly, studies for the composition and function of the gut microbiota might be essential for further research on animal reproduction.
Seasonal breeding is a phenomenon in which some animals mate only at certain times of the year. Seasonal breeding activity may be influenced mainly by photoperiodism [16]. Photoperiodic changes are sensed by the pineal gland in the brain. It secretes melatonin at night to regulate the secretion of the Gonadotropin-releasing hormone and the Luteinizing hormone, leading to seasonal changes in reproductive activity [17]. Meanwhile, several studies have confirmed that many other factors, such as estrogen, thyroid hormone, kisspeptin, and living environment, affect seasonal breeding animals [18,19,20,21]. According to a recent study, the gut microbiota may regulate seasonal breeding and behavior based on photoperiodic timing in rodents via the Hypothalamic Pituitary Gonadal (HPG) axis, melatonin, and the Kisspeptin/G-protein coupled receptor 54 (GPR54) system in the hypothalamus [22]. This report has aroused our interest in the further investigation of the association between gut microbiota and seasonal breeding.
The muskrat (Ondatra zibethicus) is a seasonal breeding animal. It has a breeding season from March to October and a non-breeding season from November to February [23]. During the breeding season, muskrats secrete a substance called musk, which is highly valued for its medicinal properties [24,25]. Most of the previous studies on muskrats have focused on the differences between the reproductive organs or secretory glands of muskrats during the breeding and non-breeding seasons [23,26,27]. Nevertheless, the variation in the composition and function of the gut microbiota of muskrats in different seasons is unknown, and the correlation between gut microbiota and seasonal reproductive phenomena in muskrats has not been explored.
The present study utilized high-throughput sequencing by metagenome sequencing to analyze the variation in the structure and function of the muskrat gut microbiota during different seasons. α and β diversity analysis, LEfSe analysis, Metastats analysis, and functional difference analysis on sequencing results were performed to explain differences in gut microbiota. This study aimed to investigate the seasonal variation of composition and function in the gut microbiota of muskrats, further study the secrets of seasonal breeding in muskrats and provide a new theoretical basis for the progress of captive breeding of muskrats.
2. Materials and Methods
2.1. Sample Collection
Adult male muskrats were obtained in January 2022 (the non-breeding season, N = 5) and May 2022 (the breeding season, N = 5) from a muskrat breeding base in Xinji, Hebei, China. Based on the references given by the breeding base, the nutritional compositions of the diets of muskrats for each season are as follows. Breeding season: total energy 19.84 MJ/kg, dry matter 36 g, crude protein 7.2 g, crude fiber 5.4 g, and crude fat 1.4 g. Non-breeding season: total energy 15.95 MJ/kg, dry matter 30 g, crude protein 3.6 g, crude fiber 9.0 g, and crude fat 1.1 g. Additionally, vitamins A, D, and E were added to the diet during the breeding season for better reproduction of muskrats [28,29]. Animals were executed after being anesthetized (CO2 inhalation) [30], and the animals were dissected using antiseptic equipment to collect cecum contents. As soon as the cecum contents were collected, they were frozen in liquid nitrogen. All animal experimental procedures were approved by the Institutional Review Board (or Ethics Committee) of Beijing Forestry University (protocol code EAWC_BJFU_2021004).
2.2. Metagenomic Sequencing and Data Processing
Using the TIANamp Stool DNA Kit (TIANGEN BIOTECH, Beijing, China), the total genomic DNA was extracted from cecum contents. The libraries were constructed through the Enzymic Universal DNAseq Library Prep Kit (Kaitai-Bio, Zhejiang, China), and then the libraries were tested for quality control. The qualifying libraries were sequenced on the Illumina NovaSeq PE150 platforms (Illumina, San Diego, CA, USA). To obtain clean data for further analysis, the raw data were pre-processed to remove low-quality and host sequence contamination. The clean metagenome data assembly was performed using Megahit 1.2.9 [31]. In the gene assembly, we keep the sequences (contigs) with lengths longer than 500 bp [32,33,34] for subsequent analysis. Gene prediction was performed on contigs using MetaGeneMark v.3.38 (Georgia Institute of Technology, Atlanta, GA, USA) [35,36,37], clustering with identity 95%, converge 90% [38,39], and statistical Unigene abundance information by Salmon 1.8.0 [40].
2.3. Analysis of Species Composition and Function
Gene prediction information was compared with the Non-Redundant Protein Sequence Database (NR), evolutionary genealogy of genes: Non-supervised Orthologous (eggNOG) database [41,42] and Kyoto Encyclopedia of Genes and Genomes (KEGG) PATHWAY database [43,44] using DIAMOND 2.0.7 (Max Planck Institute for Biology, Tübingen, Germany)(blastp, evalue ≤ 1 × 10−5) [45,46]. Species-specific taxonomic information at all levels (phylum, class, order, family, genus, and species) was obtained by MEGAN 6.21.2 (Tübingen University, Tübingen, Germany) [47,48]. The Bray–Curtis distance used for β-diversity was calculated and visualized with principal coordinate analysis (PCoA). Biomarkers with significant variation between groups were ascertained using linear discriminant analysis effect size (LEfSe) [49]. The images were drawn through R 4.1.1 or the online site ImageGP [50]. Analysis of functional pathways with significant differences was conducted using STAMP 2.1.3 [51] and the Metastats method [52].
3. Results
3.1. Metagenomic Sequencing and Gene Prediction
We performed sequencing on muskrat cecum contents samples and obtained a total of 117.50 gigabases (Gbps) of raw data. The raw data were quality controlled to obtain 115.52 Gbps of clean data. The quality control efficiency of the clean data was 98.29% (Table S1). A total of 1,944,649 contigs with the longest length of 585,984 bp were obtained (Table S2).
A total of 7,414,843 open reading frames and a total length of 4411.7 Mbp Unigene were obtained for gene prediction (Table S3). The gene information was used to construct a core-pan gene rarefaction curve to evaluate the sample sequencing depth (Figure 1A,B). The curve tends to increase and flatten out with increasing sample size, proving that our sequencing results have largely covered all species. The sample size was reasonable for our study. Meanwhile, the correlation heat map was constructed, and the richness varied greatly between groups, which further proved the reliability of the experiment (Figure 1C).
3.2. Gut Microbiota Composition and Difference
The following analyses were performed using the abundance information obtained from the species annotations. The analysis of the top 10 microbiota in relative abundance at the phylum level and genus level were plotted separately (Figure 2). At the phylum level, the dominant phylum was Bacteroidetes and Firmicutes in both the breeding and non-breeding seasons, and the mean abundance of Firmicutes was higher in the non-breeding season (41.30%) than in the breeding season (23.94%), while the mean abundance of Bacteroidetes was lower (21.71%) than in the breeding season (36.19%). From the genus perspective, only 24.76% of sequences were unclassified. Prevotella (6.63%), Muribaculaceae_noname (6.15%), Bacteroides (5.20%), and Bacteria_noname (5.19%) had high relative abundance (>5%) during both the breeding and non-breeding seasons. The detailed abundance table was shown in Table S4.
α-Diversity analysis showed significantly higher gut microbiota richness in the non-breeding season muskrats than in the breeding season (Figure 3A). PCoA at the genus level can be used to reveal the variation in the gut microbiota of muskrats between the breeding and non-breeding seasons (Figure 3B). Samples from different seasons were well clustered together, with PCo1 explaining 68.34% of the total variation in samples and PCo2 explaining 11.03% of the total variation in samples, indicating that the gut microbiota of muskrats varied at the genus level in different seasons. Prevotella, Bacteroides, and Alistipes were considerably higher in the breeding season than in the non-breeding season, whilst Flavonifractor, Oscillibacter, and Colidextribacter were significantly higher in the non-breeding season than in the breeding season (Figure 3C). LEfSe demonstrated significantly different biomarkers between groups (Figure 3D,E). The non-breeding season microbiota substantially enriched in phylum Firmicutes and phylum Actinobacteria. In contrast, phylum Tenericutes and order family Acetobacteraceae enriched in the breeding season. These results further validated the difference in gut microbiota between the breeding season and the non-breeding season.
3.3. Functional Analysis of Gut Microbiota in Muskrat
To investigate the metabolic-related changes in the gut microbiota during the breeding and non-breeding seasons, the KEGG pathway (Figure 4A) and eggNOG (Figure 5A) analyses were performed using DIAMOND. KEGG pathways were mainly enriched in metabolism (Figure 4A), including carbohydrate metabolism (12.65%), amino acid metabolism (9.26%), metabolism of cofactors and vitamins (6.90%), energy metabolism (6.40%), and glycan biosynthesis and metabolism (4.88%). Other pathways were genetic information processing (18.95%), environmental information processing (13.24%), cellular processes (9.66%), human diseases (6.70%), and organismal systems (3.12%) (Table S5). Difference analysis of the KEGG level2 pathways (p < 0.05) showed that there are 10 differentially significant pathways enriched in the non-breeding season, such as the carbohydrate metabolism, cell motility, signal transduction, drug resistance: antimicrobial, endocrine and metabolic disease, infectious disease: parasitic, cellular community-prokaryotes, infectious disease: viral, membrane transport, and signaling molecules and interaction. There are 12 differentially significant pathways enriched in the breeding season, such as the circulatory system, development and regeneration, environmental adaptation, folding, sorting and degradation, neurodegenerative disease, immune disease, glycan biosynthesis and metabolism, lipid metabolism, metabolism of cofactors and vitamins, metabolism of terpenoids and polyketides, translation and transport, and catabolism (Figure 4B).
PCoA of specific pathways showed a clear separation between the breeding and non-breeding samples (Figure 4C). The heat map showed the clustering of metabolic pathways with significant differences in relative abundance in the top 30 (p < 0.01). The breeding season was enriched with many functional pathways related to the amino acid, cofactor, and vitamin metabolism, such as biotin metabolism, D-Glutamine and D-glutamate metabolism, folate biosynthesis, and nicotinate and nicotinamide metabolism (Figure 4D). The plot of all significantly different KEGG pathway analyses was shown in Figure S1.
The 25 functions of the eggNOG database were annotated (Figure 5A). The major functions (annotated Unigene > 100,000) were translation, ribosome structure, and biogenesis (J), carbohydrate transport and metabolism (G), cell wall/membrane/envelope biogenesis (M), replication, recombination, and repair (L), amino acid transport and metabolism (E), transcription (K), general function prediction only (R), signal transduction mechanism (T), energy production and conversion (C), coenzyme transport and metabolism (H), defense mechanism (V), posttranslational modification, protein turnover, and chaperones (O), inorganic ion transport and metabolism (P), nucleotide transport and metabolism (F), cell cycle control, cell division, and chromosome partitioning (D), lipid transport and metabolism (L), and function unknown (S). The top 10 significantly different functions obtained using the Metastats method were shown in Figure 5B. During the breeding season, several metabolism-related functions were enriched: coenzyme transport and metabolism (H), lipid transport and metabolism (I), and inorganic ion transport and metabolism (P). The functional analysis illustrated that the gut microbiota of breeding season and non-breeding season muskrats differed significantly in the roles that they play in their hosts.
4. Discussion
As a typical seasonal breeding animal, changes in the gut microbiota of the muskrat with the seasons are of great curiosity. There are many studies describing the phenomenon of seasonal changes in gut microbiota, such as rats (Rattus norregicus) [53], Tibetan macaques (Macaca thibetana) [54], and Siberian hamsters (Phodopus sungorus) [55]. Seasonal differences were also found in the gut microbiota of forest musk deer (Moschus spp.), an animal capable of musk secretion similar to the muskrat [56]. The present study analyzed the changes in the gut microbiota of muskrats during different breeding seasons. The results showed that there was significant variation in the structure and function of the gut microbiota between the breeding and non-breeding muskrats. These results suggested the possible involvement of gut microbes in meeting the various demands of seasonal breeding in muskrats.
Seasonal changes may influence the structure of the muskrats’ gut microbiota. Results of alpha diversity analysis showed that the breeding and non-breeding seasons differed significantly and that the gut microbial richness was higher in the non-breeding season. Previous studies have shown that a higher abundance of gut microbiota helps hosts adapt to the external environment and enhances their resistance to adverse external conditions [57,58]. Therefore, we hypothesized that the increased abundance of gut microbiota in the non-breeding muskrats during winter may enhance resistance to factors such as the cold. The dominant phyla of the gut microbiota of the muskrat were Bacteroidetes and Firmicutes, which was consistent with our previous findings on another rodent seasonal breeder, the wild ground squirrel (Spermophilus dauricus) [59]. The dominant phylum of other rodents such as arctic ground squirrels (Urocitellus parryii) was also Bacteroidetes and Firmicutes [60]. Bacteroidetes and Firmicutes are widespread in animals and are a major component of healthy gut microbiota [61]. Their main function is to break down carbohydrates in the gut and synthesize short-chain fatty acids (SCFAs) to provide energy to the host [62]. Many studies have shown that changes in the ratio of Bacteroidetes and Firmicutes may be associated with obesity and that an increased Firmicutes/Bacteroidetes ratio may lead to obesity and other diseases [63,64]. The abundance of Bacteroidetes increased in the breeding season, whilst Firmicutes significantly increased in the non-breeding season. This change may lead to more stable energy metabolism during the breeding season and safeguard the breeding activity of muskrats.
Muskrats in the breeding season have a higher protein requirement, about 7.2 g of protein per 36 g of dry matter, which is greater than in the non-breeding season (winter). The winter diet should not be high in protein and fat to avoid over-fattening the animal [65]. Genus Prevotella and genus Alistipes were enriched in the breeding season, genus Oscillibacter, order Lactobacillales, family Lachnospiraceae, and family Ruminococcaceae were enriched in the non-breeding season. In a previous study, Prevotella and Alistipes were enriched in Brandt’s voles (Lasiopodomys brandtii) with long photoperiods [22], which were reflected in the present study for the breeding season. Alistipes were positively correlated with cholesterol metabolism [66], which may be related to the high-fat, high-protein diet structure of breeding muskrats. Prevotella had a higher diversity, and it has been suggested that it may promote increased glycogen stores [67], while the study by Chen et al. also demonstrated that Prevotella was positively associated with fat accumulation in pigs [68].
Genus Oscillibacter has been reported to be positively associated with a reduction in body weight, which is an important reference for weight control in non-breeding muskrats [69]. The enrichment of some bacteria of order Lactobacillales has been found to cause a decrease in sperm viability, which may be the reason why they were enriched in the non-breeding season rather than the breeding season [70]. Additionally, the family Lachnospiraceae, order Lactobacillales, and family Ruminococcaceae can hydrolyze starch and other sugars to produce butyrate and other SCFAs [71,72]. The non-breeding muskrats with increased carbohydrate content in their diet and reduced total food intake are likely to encounter energy-deficient conditions [65], and these gut microorganisms provide effective energy for muskrats over the winter.
The present study also predicted and analyzed the function of the muskrat gut microbiota based on the KEGG and eggNOG databases. The results of the analysis showed that most of the functions of muskrat gut microbiota were enriched in metabolism-related pathways, especially carbohydrate metabolism. Sixty percent of the muskrat’s diet is composed of carbohydrates, which are the main source of energy for the muskrat. Other major pathways involved are amino acid metabolism, metabolism of cofactors and vitamins, energy metabolism, lipid metabolism, etc. There are significant seasonal differences in all these functions. Muskrats, similar to other monogastric animals, can automatically regulate the amount of food that they eat to meet their energy requirements. During the breeding season, the muskrat’s diet contains more animal feed and higher levels of proteins, lipids, and other nutrients, so the metabolism pathways associated with amino acids, lipids, vitamins, and cofactors increase significantly during the breeding season to meet the complex food structure of the muskrat during that season. In contrast, during the non-breeding season, muskrats eat less and have a high proportion of plant-based feeds, while the lower ambient temperature was not conducive to their energy supply, and the enriched carbohydrate metabolic pathways at this time can help muskrats better survive in this period.
This study investigated seasonal differences in the gut microbiota of muskrats using a metagenomic approach. The present study initially explored the correlation between gut microbiota and seasonal breeding and speculated on the possible regulatory mechanisms of gut microbiota on seasonal breeding in muskrats. This research can provide a new theoretical basis for muskrat captive breeding and help breeders understand the seasonal changes in muskrat metabolism. They can use this as a basis to better regulate the diet structure of muskrats in different seasons. Furthermore, more work needs to be conducted in the future to clarify the composition and function of the gut microbes in muskrats and to further explore the important role of the gut microbiota in seasonal breeding animals.
5. Conclusions
There were notable seasonal differences in the species composition and structure of the gut microbiota of muskrats. Gut microbiota richness increased in muskrats during the non-breeding season, which may help muskrats to increase their resistance to the external environment. At the phylum level, the relative abundance of Bacteroidetes was increased and the relative abundance of Firmicutes was decreased in the gut microbiota of breeding muskrats. The genus Prevotella and genus Alistipes enriched in the breeding season, and the genus Oscillibacter, order Lactobacillales, family Lachnospiraceae, and family Ruminococcaceae enriched in the non-breeding season. The muskrat’s gut microbiota was highly enriched in nutrient metabolic pathways. Differential microbiotas during the breeding season were mainly enriched in amino acid, lipid, vitamin, and cofactor metabolism pathways to help muskrats better reproductive activities. This study may provide a new theoretical basis for managing muskrat captive breeding and lay the foundation for further ensuring the efficient breeding of muskrats.
Conceptualization, F.S., H.Z., Y.H., Q.W. and Z.Y.; investigation, F.S.; resources, F.S., Y.X., P.P., H.L. and R.Z.; supervision, Z.Y.; writing–original draft, F.S.; writing—review and editing, F.S., Y.X., P.P., H.L., R.Z., H.Z., Y.H., Q.W. and Z.Y. All authors have read and agreed to the published version of the manuscript.
The animal study protocol was approved by the Institutional Review Board (or Ethics Committee) of Beijing Forestry University (protocol code EAWC_BJFU_2021004). Approval Date: 11 March 2021.
Not applicable.
Data is contained within the article and
The authors declare no conflict of interest.
Footnotes
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Enlarge this image.Figure 1. Statistical analysis of gene information. The rarefaction curve for pan gene (A) and core gene (B) (“+” and “○” represent outliers). (C) The heat map of correlation among samples.
Enlarge this image.Figure 2. Species composition between the breeding (B) and non-breeding seasons (NB). The top ten species in relative abundance for each taxon were selected for plotting, and the remaining species were categorized as others. (A) The stack bar plot of phylum-level species composition. (B) The stack bar plot of genus-level species composition.
Enlarge this image.Figure 3. Bioinformatics analysis of the differences in gut microbiota composition between the breeding (B) and non-breeding (NB) seasons. (A) The Shannon’s index of α-diversity analysis (p < 0.05). (B) The Principal Coordinate Analysis (PCoA) is based on the Bray–Curtis distance with significant differences in genus levels between the two groups. (C) The genus-level abundance heat map shows the top 20 genera in relative abundance. (D,E) The LEfSe plots show taxa differing between groups. (D) The bar plot shows the distribution of LDA values for the differential biomarkers, and the histogram of LDA value distribution shows the species with significant differences (LDAScore > 4). (E) The cladogram shows the evolutionary branching of divergent species.
Enlarge this image.Figure 4. Functional prediction analysis of the muskrat gut microbiota based on the KEGG database. (A) The statistical plot of KEGG annotations of the muskrat gut microbiota. (B) The top 20 significant functional difference analyses of muskrat gut microbiota in the KEGG level 2 pathway (p < 0.05). (C) The functional prediction of muskrat gut microbiota in KEGG pathway by Principal coordinate analysis (PCoA). (D) The heat map of significantly different top 30 KEGG pathways (p < 0.01).
Enlarge this image.Figure 5. Functional prediction analysis of the muskrat gut microbiota based on the eggNOG database. (A) The bar plot of eggNOG annotations of the muskrat gut microbiota. (B) The box plot based on Metastats analysis data shows the top 10 functions with significant differences.
Supplementary Materials
The following supporting information can be downloaded at:
References
1. Hsiao, W.W.; Metz, C.; Singh, D.P.; Roth, J. The microbes of the intestine: An introduction to their metabolic and signaling capabilities. Endocrinol. Metab. Clin. North Am.; 2008; 37, pp. 857-871. [DOI: https://dx.doi.org/10.1016/j.ecl.2008.08.006] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19026936]
2. Ley, R.E.; Peterson, D.A.; Gordon, J.I. Ecological and evolutionary forces shaping microbial diversity in the human intestine. Cell; 2006; 124, pp. 837-848. [DOI: https://dx.doi.org/10.1016/j.cell.2006.02.017] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16497592]
3. Morrison, D.J.; Preston, T. Formation of short chain fatty acids by the gut microbiota and their impact on human metabolism. Gut Microbes; 2016; 7, pp. 189-200. [DOI: https://dx.doi.org/10.1080/19490976.2015.1134082] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26963409]
4. O’Hara, A.M.; Shanahan, F. The gut flora as a forgotten organ. EMBO Rep.; 2006; 7, pp. 688-693. [DOI: https://dx.doi.org/10.1038/sj.embor.7400731]
5. Chen, K.; Pachter, L. Bioinformatics for whole-genome shotgun sequencing of microbial communities. PLoS Comput. Biol.; 2005; 1, pp. 106-112. [DOI: https://dx.doi.org/10.1371/journal.pcbi.0010024]
6. Handelsman, J.; Rondon, M.R.; Brady, S.F.; Clardy, J.; Goodman, R.M. Molecular biological access to the chemistry of unknown soil microbes: A new frontier for natural products. Chem Biol.; 1998; 5, pp. R245-R249. [DOI: https://dx.doi.org/10.1016/S1074-5521(98)90108-9]
7. Gill, S.R.; Pop, M.; Deboy, R.T.; Eckburg, P.B.; Turnbaugh, P.J.; Samuel, B.S.; Gordon, J.I.; Relman, D.A.; Fraser-Liggett, C.M.; Nelson, K.E. Metagenomic analysis of the human distal gut microbiome. Science; 2006; 312, pp. 1355-1359. [DOI: https://dx.doi.org/10.1126/science.1124234]
8. Fan, C.; Zhang, L.; Jia, S.; Tang, X.; Fu, H.; Li, W.; Liu, C.; Zhang, H.; Cheng, Q.; Zhang, Y. Seasonal variations in the composition and functional profiles of gut microbiota reflect dietary changes in plateau pikas. Integr. Zool; 2022; 17, pp. 379-395. [DOI: https://dx.doi.org/10.1111/1749-4877.12630]
9. Guo, N.; Wu, Q.; Shi, F.; Niu, J.; Zhang, T.; Degen, A.A.; Fang, Q.; Ding, L.; Shang, Z.; Zhang, Z. et al. Seasonal dynamics of diet-gut microbiota interaction in adaptation of yaks to life at high altitude. NPJ Biofilms Microbiomes; 2021; 7, 38. [DOI: https://dx.doi.org/10.1038/s41522-021-00207-6]
10. O’Toole, P.W. Changes in the intestinal microbiota from adulthood through to old age. Clin. Microbiol. Infect.; 2012; 18, pp. 44-46. [DOI: https://dx.doi.org/10.1111/j.1469-0691.2012.03867.x]
11. Gomaa, E.Z. Human gut microbiota/microbiome in health and diseases: A review. Antonie Van Leeuwenhoek; 2020; 113, pp. 2019-2040. [DOI: https://dx.doi.org/10.1007/s10482-020-01474-7]
12. Sun, Y.; Sun, Y.; Shi, Z.; Liu, Z.; Zhao, C.; Lu, T.; Gao, H.; Zhu, F.; Chen, R.; Zhang, J. et al. Gut Microbiota of Wild and Captive Alpine Musk Deer (Moschus chrysogaster). Front. MicroBiol.; 2019; 10, 3156. [DOI: https://dx.doi.org/10.3389/fmicb.2019.03156]
13. Wong, S.; Rawls, J.F. Intestinal microbiota composition in fishes is influenced by host ecology and environment. Mol. Ecol.; 2012; 21, pp. 3100-3102. [DOI: https://dx.doi.org/10.1111/j.1365-294X.2012.05646.x]
14. Ridlon, J.M.; Ikegawa, S.; Alves, J.M.P.; Zhou, B.; Kobayashi, A.; Iida, T.; Mitamura, K.; Tanabe, G.; Serrano, M.; De Guzman, A. et al. Clostridium scindens: A human gut microbe with a high potential to convert glucocorticoids into androgens. J. Lipid Res.; 2013; 54, pp. 2437-2449. [DOI: https://dx.doi.org/10.1194/jlr.M038869]
15. Colldén, H.; Landin, A.; Wallenius, V.; Elebring, E.; Fändriks, L.; Nilsson, M.E.; Ryberg, H.; Poutanen, M.; Sjögren, K.; Vandenput, L. et al. The gut microbiota is a major regulator of androgen metabolism in intestinal contents. Am. J. Physiol Endocrinol. Metab.; 2019; 317, pp. E1182-E1192. [DOI: https://dx.doi.org/10.1152/ajpendo.00338.2019]
16. Matthews, C.D.; Guerin, M.V.; Deed, J.R. Melatonin and photoperiodic time measurement: Seasonal breeding in the sheep. J. Pineal Res.; 1993; 14, pp. 105-116. [DOI: https://dx.doi.org/10.1111/j.1600-079X.1993.tb00493.x]
17. Malpaux, B.; Tricoire, H.; Mailliet, F.; Daveau, A.; Migaud, M.; Skinner, D.C.; Pelletier, J.; Chemineau, P. Melatonin and seasonal reproduction: Understanding the neuroendocrine mechanisms using the sheep as a model. Reprod Suppl.; 2002; 59, pp. 167-179.
18. Dardente, H.; Lomet, D.; Robert, V.; Decourt, C.; Beltramo, M.; Pellicer-Rubio, M.T. Seasonal breeding in mammals: From basic science to applications and back. Theriogenology; 2016; 86, pp. 324-332. [DOI: https://dx.doi.org/10.1016/j.theriogenology.2016.04.045]
19. Nakao, N.; Ono, H.; Yoshimura, T. Thyroid hormones and seasonal reproductive neuroendocrine interactions. Reproduction; 2008; 136, pp. 1-8. [DOI: https://dx.doi.org/10.1530/REP-08-0041]
20. Smith, J.T. The role of kisspeptin and gonadotropin inhibitory hormone in the seasonal regulation of reproduction in sheep. Domest. Anim. Endocrinol.; 2012; 43, pp. 75-84. [DOI: https://dx.doi.org/10.1016/j.domaniend.2011.11.003]
21. Smith, J.T.; Clarke, I.J. Seasonal breeding as a neuroendocrine model for puberty in sheep. Mol. Cell Endocrinol.; 2010; 324, pp. 102-109. [DOI: https://dx.doi.org/10.1016/j.mce.2010.03.007] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20298744]
22. Zhu, H.; Li, G.; Liu, J.; Xu, X.; Zhang, Z. Gut microbiota is associated with the effect of photoperiod on seasonal breeding in male Brandt’s voles (Lasiopodomys brandtii). Microbiome; 2022; 10, 194. [DOI: https://dx.doi.org/10.1186/s40168-022-01381-1] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/36376894]
23. Zhang, F.; Liu, Q.; Wang, Z.; Xie, W.; Sheng, X.; Zhang, H.; Yuan, Z.; Han, Y.; Weng, Q. Seasonal Expression of Oxytocin and Oxytocin Receptor in the Scented Gland of Male Muskrat (Ondatra zibethicus). Sci. Rep.; 2017; 7, 16627. [DOI: https://dx.doi.org/10.1038/s41598-017-16973-3] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29192229]
24. Chen, Y.S.; Jin, S.D.; Dong, W.C.; Tong, Y.R.; Hua, S.F.; Zhao, H.Y. Effect of muskrat on the cardiovascular system in anesthetized dogs. Zhongguo Zhong Yao Za Zhi; 1989; 14, pp. 241–243, 256. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/2505807]
25. Li, Y.; Zhang, T.; Fan, M.; Zhou, J.; Yang, S.; Zhang, M.; Qi, L.; Lin, S.; Hu, D.; Liu, S. Comparison of amino acid profiles and metabolic gene expression in muskrat scented glands in secretion and non-secretion season. Sci. Rep.; 2017; 7, 41158. [DOI: https://dx.doi.org/10.1038/srep41158]
26. Li, Y.; Zhang, T.; Zhou, J.; Yang, S.; Fan, M.; Sun, X.; Zhang, M.; Xu, S.; Cha, M.; Hu, X. et al. Transcriptome analysis of muskrat scented glands degeneration mechanism. PLoS ONE; 2017; 12, e0176935. [DOI: https://dx.doi.org/10.1371/journal.pone.0176935]
27. Cao, H.; Wang, L.; Zhang, S.; Lu, L.; Sheng, X.; Han, Y.; Yuan, Z.; Weng, Q. Seasonal Expression of Prolactin Receptor in the Scented Gland of Male Muskrat (Ondatra zibethicus). Sci. Rep.; 2015; 5, 15036. [DOI: https://dx.doi.org/10.1038/srep15036]
28. Hua, S.F. Muskrat Breeding Technology; China Agricultural Press: Beijing, China, 1989; pp. 52-53.
29. Wang, Z.Y. Practical Breeding Techniques for Muskrat; Xinji Jinmu Technology Co.: Xinji, China, 2010; pp. 21-30.
30. Shomer, N.H.; Allen-Worthington, K.H.; Hickman, D.L.; Jonnalagadda, M.; Newsome, J.T.; Slate, A.R.; Valentine, H.; Williams, A.M.; Wilkinson, M. Review of Rodent Euthanasia Methods. J. Am. Assoc. Lab Anim. Sci.; 2020; 59, pp. 242-253. [DOI: https://dx.doi.org/10.30802/AALAS-JAALAS-19-000084]
31. Li, D.; Liu, C.-M.; Luo, R.; Sadakane, K.; Lam, T.-W. MEGAHIT: An ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics; 2015; 31, pp. 1674-1676. [DOI: https://dx.doi.org/10.1093/bioinformatics/btv033]
32. Zeller, G.; Tap, J.; Voigt, A.Y.; Sunagawa, S.; Kultima, J.R.; Costea, P.I.; Amiot, A.; Böhm, J.; Brunetti, F.; Habermann, N. et al. Potential of fecal microbiota for early-stage detection of colorectal cancer. Molecular Systems Biology; 2014; 10, 766. [DOI: https://dx.doi.org/10.15252/msb.20145645]
33. Qin, J.; Li, R.; Raes, J.; Arumugam, M.; Burgdorf, K.S.; Manichanh, C.; Nielsen, T.; Pons, N.; Levenez, F.; Yamada, T. et al. A human gut microbial gene catalogue established by metagenomic sequencing. Nature; 2010; 464, pp. 59-65. [DOI: https://dx.doi.org/10.1038/nature08821]
34. Nielsen, H.B.; Almeida, M.; Juncker, A.S.; Rasmussen, S.; Li, J.; Sunagawa, S.; Plichta, D.R.; Gautier, L.; Pedersen, A.G.; Le Chatelier, E. et al. Identification and assembly of genomes and genetic elements in complex metagenomic samples without using reference genomes. Nat. Biotechnol.; 2014; 32, pp. 822-828. [DOI: https://dx.doi.org/10.1038/nbt.2939]
35. Daniel, R.M.; Alison, S.W.; Shinichi, S.; Aino, I.J.; Aino, I.J.; Michelle, M.Y.C.; Manimozhiyan, A.; Jeroen, R.; Peer, B. Assessment of metagenomic assembly using simulated next generation sequencing data. PLoS ONE; 2012; [DOI: https://dx.doi.org/10.1371/journal.pone.0031386]
36. Julia, O.; Allyson, L.B.; Clay, D.; Sean, C.; Nisc Comparative Sequencing, P.; Heidi, H.K.; Julia, A.S. Biogeography and individuality shape function in the human skin metagenome. Nature; 2014; 514, pp. 59-64. [DOI: https://dx.doi.org/10.1038/nature13786]
37. Jun, L.; Junhua, L.; Huijue, J.; Xianghang, C.; Huanzi, Z.; Qiang, F.; Shinichi, S.; Manimozhiyan, A.; Jens Roat, K.; Edi, P. et al. An integrated catalog of reference genes in the human gut microbiome. Nat. Biotechnol.; 2014; 32, pp. 834-841. [DOI: https://dx.doi.org/10.1038/nbt.2942]
38. Georg, Z.; Julien, T.; Anita, Y.V.; Shinichi, S.; Jens Roat, K.; Paul, I.C.; Aurelien, A.; Jürgen, B.; Francesco, B.; Nina, H. et al. Potential of fecal microbiota for early-stage detection of colorectal cancer. Mol. Syst. Biol.; 2014; 10, 766. [DOI: https://dx.doi.org/10.15252/msb.20145645]
39. Nan, Q.; Fengling, Y.; Ang, L.; Edi, P.; Yanfei, C.; Li, S.; Jing, G.; Jian, Y.; Lingjiao, W.; Jiawei, Z. et al. Alterations of the human gut microbiome in liver cirrhosis. Nature; 2014; 513, pp. 59-64. [DOI: https://dx.doi.org/10.1038/nature13568]
40. Rob, P.; Geet, D.; Michael, I.L.; Rafael, A.I.; Carl, K. Salmon provides fast and bias-aware quantification of transcript expression. Nat. Methods; 2017; 14, pp. 417-419. [DOI: https://dx.doi.org/10.1038/nmeth.4197]
41. Powell, S.; Forslund, K.; Szklarczyk, D.; Trachana, K.; Roth, A.; Huerta-Cepas, J.; Gabaldón, T.; Rattei, T.; Creevey, C.; Kuhn, M. et al. eggNOG v4.0: Nested orthology inference across 3686 organisms. Nucleic Acids Res.; 2014; 42, pp. D231-D239. [DOI: https://dx.doi.org/10.1093/nar/gkt1253]
42. Jensen, L.J.; Julien, P.; Kuhn, M.; von Mering, C.; Muller, J.; Doerks, T.; Bork, P. eggNOG: Automated construction and annotation of orthologous groups of genes. Nucleic Acids Res.; 2008; 36, pp. D250-D254. [DOI: https://dx.doi.org/10.1093/nar/gkm796]
43. Kanehisa, M.; Goto, S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res.; 2000; 28, pp. 27-30. [DOI: https://dx.doi.org/10.1093/nar/28.1.27] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/10592173]
44. Kanehisa, M.; Goto, S.; Sato, Y.; Kawashima, M.; Furumichi, M.; Tanabe, M. Data, information, knowledge and principle: Back to metabolism in KEGG. Nucleic Acids Res.; 2014; 42, pp. D199-D205. [DOI: https://dx.doi.org/10.1093/nar/gkt1076]
45. Buchfink, B.; Xie, C.; Huson, D.H. Fast and sensitive protein alignment using DIAMOND. Nat. Methods; 2015; 12, pp. 59-60. [DOI: https://dx.doi.org/10.1038/nmeth.3176] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25402007]
46. Buchfink, B.; Reuter, K.; Drost, H.-G. Sensitive protein alignments at tree-of-life scale using DIAMOND. Nat. Methods; 2021; 18, pp. 366-368. [DOI: https://dx.doi.org/10.1038/s41592-021-01101-x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33828273]
47. Karlsson, F.H.; Tremaroli, V.; Nookaew, I.; Bergström, G.; Behre, C.J.; Fagerberg, B.; Nielsen, J.; Bäckhed, F. Gut metagenome in European women with normal, impaired and diabetic glucose control. Nature; 2013; 498, pp. 99-103. [DOI: https://dx.doi.org/10.1038/nature12198]
48. Bağcı, C.; Patz, S.; Huson, D.H. DIAMOND+MEGAN: Fast and Easy Taxonomic and Functional Analysis of Short and Long Microbiome Sequences. Curr. Protoc.; 2021; 1, e59. [DOI: https://dx.doi.org/10.1002/cpz1.59]
49. Segata, N.; Izard, J.; Waldron, L.; Gevers, D.; Miropolsky, L.; Garrett, W.S.; Huttenhower, C. Metagenomic biomarker discovery and explanation. Genome Biol.; 2011; 12, R60. [DOI: https://dx.doi.org/10.1186/gb-2011-12-6-r60]
50. Chen, T.; Liu, Y.-X.; Huang, L. ImageGP: An easy-to-use data visualization web server for scientific researchers. iMeta; 2022; 1, e5. [DOI: https://dx.doi.org/10.1002/imt2.5]
51. Parks, D.H.; Tyson, G.W.; Hugenholtz, P.; Beiko, R.G. STAMP: Statistical analysis of taxonomic and functional profiles. Bioinformatics; 2014; 30, pp. 3123-3124. [DOI: https://dx.doi.org/10.1093/bioinformatics/btu494]
52. White, J.R.; Nagarajan, N.; Pop, M. Statistical methods for detecting differentially abundant featuRes. in clinical metagenomic samples. PLoS Comput. Biol.; 2009; 5, e1000352. [DOI: https://dx.doi.org/10.1371/journal.pcbi.1000352]
53. Liu, J.; Song, Y.; Lu, X.; Chen, T.; Guo, W.; Fan, Z.; Kang, X.; Wang, Y.; Wang, Y. Seasonal Variation Influences on Intestinal Microbiota in Rats. Curr. MicroBiol.; 2018; 75, pp. 1006-1010. [DOI: https://dx.doi.org/10.1007/s00284-018-1480-6]
54. Sun, B.; Wang, X.; Bernstein, S.; Huffman, M.A.; Xia, D.P.; Gu, Z.; Chen, R.; Sheeran, L.K.; Wagner, R.S.; Li, J. Marked variation between winter and spring gut microbiota in free-ranging Tibetan Macaques (Macaca thibetana). Sci. Rep.; 2016; 6, 26035. [DOI: https://dx.doi.org/10.1038/srep26035]
55. Shor, E.K.; Brown, S.P.; Freeman, D.A. A novel role for the pineal gland: Regulating seasonal shifts in the gut microbiota of Siberian hamsters. J. Pineal Res.; 2020; 69, e12696. [DOI: https://dx.doi.org/10.1111/jpi.12696]
56. Jiang, F.; Gao, H.; Qin, W.; Song, P.; Wang, H.; Zhang, J.; Liu, D.; Wang, D.; Zhang, T. Marked Seasonal Variation in Structure and Function of Gut Microbiota in Forest and Alpine Musk Deer. Front MicroBiol.; 2021; 12, 699797. [DOI: https://dx.doi.org/10.3389/fmicb.2021.699797]
57. Stoffel, M.A.; Acevedo-Whitehouse, K.; Morales-Durán, N.; Grosser, S.; Chakarov, N.; Krüger, O.; Nichols, H.J.; Elorriaga-Verplancken, F.R.; Hoffman, J.I. Early sexual dimorphism in the developing gut microbiome of northern elephant seals. Mol. Ecol.; 2020; 29, pp. 2109-2122. [DOI: https://dx.doi.org/10.1111/mec.15385]
58. Videvall, E. Microbiome maturation during a unique developmental window. Mol. Ecol.; 2020; 29, pp. 1941-1943. [DOI: https://dx.doi.org/10.1111/mec.15436]
59. Yang, X.; Yao, Y.; Zhang, X.; Zhong, J.; Gao, F.; Zhang, H.; Han, Y.; Weng, Q.; Yuan, Z. Seasonal Changes in the Distinct Taxonomy and Function of the Gut Microbiota in the Wild Ground Squirrel (Spermophilus dauricus). Animals; 2021; 11, 2685. [DOI: https://dx.doi.org/10.3390/ani11092685]
60. Stevenson, T.J.; Duddleston, K.N.; Buck, C.L. Effects of season and host physiological state on the diversity, density, and activity of the arctic ground squirrel cecal microbiota. Appl. Environ. MicroBiol.; 2014; 80, pp. 5611-5622. [DOI: https://dx.doi.org/10.1128/AEM.01537-14]
61. Jandhyala, S.M.; Talukdar, R.; Subramanyam, C.; Vuyyuru, H.; Sasikala, M.; Nageshwar Reddy, D. Role of the normal gut microbiota. World J. Gastroenterol; 2015; 21, pp. 8787-8803. [DOI: https://dx.doi.org/10.3748/wjg.v21.i29.8787]
62. Macfarlane, S.; Macfarlane, G.T. Regulation of short-chain fatty acid production. Proc. Nutr. Soc.; 2003; 62, pp. 67-72. [DOI: https://dx.doi.org/10.1079/PNS2002207]
63. Indiani, C.; Rizzardi, K.F.; Castelo, P.M.; Ferraz, L.F.C.; Darrieux, M.; Parisotto, T.M. Childhood Obesity and Firmicutes/Bacteroidetes Ratio in the Gut Microbiota: A Systematic Review. Child Obes.; 2018; 14, pp. 501-509. [DOI: https://dx.doi.org/10.1089/chi.2018.0040] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30183336]
64. Magne, F.; Gotteland, M.; Gauthier, L.; Zazueta, A.; Pesoa, S.; Navarrete, P.; Balamurugan, R. The Firmicutes/Bacteroidetes Ratio: A Relevant Marker of Gut Dysbiosis in Obese Patients?. Nutrients; 2020; 12, 1474. [DOI: https://dx.doi.org/10.3390/nu12051474] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32438689]
65. Feng, X.L. Characteristics of nutritional metabolism and energy balance of muskrats. Mod. Livest. Poult. Breed. Ind.; 2012; 10, pp. 48-49. (In Chinese) [DOI: https://dx.doi.org/10.14070/j.cnki.15-1150.2012.10.050]
66. Le Roy, T.; Lécuyer, E.; Chassaing, B.; Rhimi, M.; Lhomme, M.; Boudebbouze, S.; Ichou, F.; Haro Barceló, J.; Huby, T.; Guerin, M. et al. The intestinal microbiota regulates host cholesterol homeostasis. BMC Biol.; 2019; 17, 94. [DOI: https://dx.doi.org/10.1186/s12915-019-0715-8]
67. Kovatcheva-Datchary, P.; Nilsson, A.; Akrami, R.; Lee, Y.S.; De Vadder, F.; Arora, T.; Hallen, A.; Martens, E.; Björck, I.; Bäckhed, F. Dietary Fiber-Induced Improvement in Glucose Metabolism Is Associated with Increased Abundance of Prevotella. Cell Metab.; 2015; 22, pp. 971-982. [DOI: https://dx.doi.org/10.1016/j.cmet.2015.10.001]
68. Chen, C.; Fang, S.; Wei, H.; He, M.; Fu, H.; Xiong, X.; Zhou, Y.; Wu, J.; Gao, J.; Yang, H. et al. Prevotella copri increases fat accumulation in pigs fed with formula diets. Microbiome; 2021; 9, 175. [DOI: https://dx.doi.org/10.1186/s40168-021-01110-0]
69. Konikoff, T.; Gophna, U. Oscillospira: A Central, Enigmatic Component of the Human Gut Microbiota. Trends MicroBiol.; 2016; 24, pp. 523-524. [DOI: https://dx.doi.org/10.1016/j.tim.2016.02.015]
70. Wang, H.; Chen, T.; Chen, Y.; Luo, T.; Tan, B.; Chen, H.; Xin, H. Evaluation of the inhibitory effects of vaginal microorganisms on sperm motility in vitro. Exp. Ther. Med.; 2020; 19, pp. 535-544. [DOI: https://dx.doi.org/10.3892/etm.2019.8237]
71. Vacca, M.; Celano, G.; Calabrese, F.M.; Portincasa, P.; Gobbetti, M.; De Angelis, M. The Controversial Role of Human Gut Lachnospiraceae. Microorganisms; 2020; 8, 573. [DOI: https://dx.doi.org/10.3390/microorganisms8040573]
72. Zhuang, M.; Shang, W.; Ma, Q.; Strappe, P.; Zhou, Z. Abundance of Probiotics and Butyrate-Production Microbiome Manages Constipation via Short-Chain Fatty Acids Production and Hormones Secretion. Mol. Nutr. Food Res.; 2019; 63, e1801187. [DOI: https://dx.doi.org/10.1002/mnfr.201801187]
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Details
; Han, Yingying; Weng, Qiang; Yuan, Zhengrong