Next Article in Journal
Association of Plasma Lipopolysaccharide-Binding Protein Concentration with Dietary Factors, Gut Microbiota, and Health Status in the Japanese General Adult Population: A Cross-Sectional Study
Next Article in Special Issue
Inter-Day Variation in the Fasting Plasma Lipopolysaccharide Concentration in the Morning Is Associated with Inter-Day Variation in Appetite in Japanese Males: A Short-Term Cohort Study
Previous Article in Journal
Anti-Inflammatory Activity of Compounds Derived from Vitex rotundifolia
Previous Article in Special Issue
Target Metabolites to Slow Down Progression of Amyotrophic Lateral Sclerosis in Mice
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Metabolites
Volume 13
Issue 2
10.3390/metabo13020248
metabolites-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Seasonal Changes in the Structure and Function of Gut Microbiota in the Muskrat (Ondatra zibethicus)

by
Fengcheng Song
,
Yishu Xu
,
Peng Peng
,
Hongxu Li
,
Ranxi Zheng
,
Haolin Zhang
,
Yingying Han
,
Qiang Weng
and
Zhengrong Yuan
*
College of Biological Sciences and Technology, Beijing Forestry University, Beijing 100083, China
*
Author to whom correspondence should be addressed.
Metabolites 2023, 13(2), 248; https://doi.org/10.3390/metabo13020248
Submission received: 4 December 2022 / Revised: 19 January 2023 / Accepted: 29 January 2023 / Published: 9 February 2023
(This article belongs to the Special Issue The Role of Gut Microbes in Metabolism Regulation)
Browse Figures
Versions Notes

Abstract

:
The gut microbiota plays a crucial role in the nutrition, metabolism, and immune function of the host animal. The muskrat (Ondatra zibethicus) is a typical seasonal breeding animal. The present study performed a metagenomic analysis of cecum contents from muskrats in the breeding and non-breeding seasons. The results indicated that the breeding muskrats and non-breeding muskrats differed in gut microbiota structure and function. During the breeding season, the relative abundance of phylum Bacteroidetes, genus Prevotella, and genus Alistipes increased, while the relative abundance of phylum Firmicutes and phylum Actinobacteria decreased. The muskrat gut microbiota was enriched in the metabolism-related pathways, especially amino acid and vitamin metabolism, and genetically related metabolites in the breeding season. We presumed that the muskrat gut microbiota might seasonally change to secure reproductive activity and satisfy the metabolic demands of different seasons. This study could explore potential mechanisms by which gut microbiota affects reproduction. Moreover, this study may provide a new theoretical basis for the management of muskrat captive breeding.

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.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/metabo13020248/s1, Figure S1: Analysis of all p < 0.01 differential metabolic pathways. Table S1: Clean data information after quality control of raw data from sequencing. Table S2: Gene assembly result evaluation. Table S3: Gene prediction and abundance statistics. Table S4: Abundance statistics for all taxonomic levels (phylum, order, family, genus, and species). Table S5: Abundance statistics for KEGG pathways.

Author Contributions

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.

Funding

This research work is supported by the National Natural Science Foundation of China (32272859), and the National Training Program of Innovation and Entrepreneurship for Undergraduates (S202210022048, X202210022057).

Institutional Review Board Statement

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.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article and Supplementary Material.

Conflicts of Interest

The authors declare no conflict of interest.

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, 857–871. [Google Scholar] [CrossRef] [PubMed]
  2. Ley, R.E.; Peterson, D.A.; Gordon, J.I. Ecological and evolutionary forces shaping microbial diversity in the human intestine. Cell 2006, 124, 837–848. [Google Scholar] [CrossRef] [PubMed]
  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, 189–200. [Google Scholar] [CrossRef] [PubMed]
  4. O’Hara, A.M.; Shanahan, F. The gut flora as a forgotten organ. EMBO Rep. 2006, 7, 688–693. [Google Scholar] [CrossRef]
  5. Chen, K.; Pachter, L. Bioinformatics for whole-genome shotgun sequencing of microbial communities. PLoS Comput. Biol. 2005, 1, 106–112. [Google Scholar] [CrossRef]
  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, R245–R249. [Google Scholar] [CrossRef]
  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, 1355–1359. [Google Scholar] [CrossRef]
  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, 379–395. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  10. O’Toole, P.W. Changes in the intestinal microbiota from adulthood through to old age. Clin. Microbiol. Infect. 2012, 18, 44–46. [Google Scholar] [CrossRef] [Green Version]
  11. Gomaa, E.Z. Human gut microbiota/microbiome in health and diseases: A review. Antonie Van Leeuwenhoek 2020, 113, 2019–2040. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  13. Wong, S.; Rawls, J.F. Intestinal microbiota composition in fishes is influenced by host ecology and environment. Mol. Ecol. 2012, 21, 3100–3102. [Google Scholar] [CrossRef]
  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, 2437–2449. [Google Scholar] [CrossRef]
  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, E1182–E1192. [Google Scholar] [CrossRef]
  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, 105–116. [Google Scholar] [CrossRef]
  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, 167–179. [Google Scholar]
  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, 324–332. [Google Scholar] [CrossRef]
  19. Nakao, N.; Ono, H.; Yoshimura, T. Thyroid hormones and seasonal reproductive neuroendocrine interactions. Reproduction 2008, 136, 1–8. [Google Scholar] [CrossRef]
  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, 75–84. [Google Scholar] [CrossRef]
  21. Smith, J.T.; Clarke, I.J. Seasonal breeding as a neuroendocrine model for puberty in sheep. Mol. Cell Endocrinol. 2010, 324, 102–109. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef] [PubMed]
  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, 241–243, 256. [Google Scholar] [PubMed]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  28. Hua, S.F. Muskrat Breeding Technology; China Agricultural Press: Beijing, China, 1989; pp. 52–53. [Google Scholar]
  29. Wang, Z.Y. Practical Breeding Techniques for Muskrat; Xinji Jinmu Technology Co.: Xinji, China, 2010; pp. 21–30. [Google Scholar]
  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, 242–253. [Google Scholar] [CrossRef]
  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, 1674–1676. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  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, 59–65. [Google Scholar] [CrossRef] [Green Version]
  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, 822–828. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  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, 59–64. [Google Scholar] [CrossRef]
  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, 834–841. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  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, 59–64. [Google Scholar] [CrossRef]
  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, 417–419. [Google Scholar] [CrossRef]
  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, D231–D239. [Google Scholar] [CrossRef]
  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, D250–D254. [Google Scholar] [CrossRef]
  43. Kanehisa, M.; Goto, S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28, 27–30. [Google Scholar] [CrossRef] [PubMed]
  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, D199–D205. [Google Scholar] [CrossRef]
  45. Buchfink, B.; Xie, C.; Huson, D.H. Fast and sensitive protein alignment using DIAMOND. Nat. Methods 2015, 12, 59–60. [Google Scholar] [CrossRef] [PubMed]
  46. Buchfink, B.; Reuter, K.; Drost, H.-G. Sensitive protein alignments at tree-of-life scale using DIAMOND. Nat. Methods 2021, 18, 366–368. [Google Scholar] [CrossRef] [PubMed]
  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, 99–103. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  50. Chen, T.; Liu, Y.-X.; Huang, L. ImageGP: An easy-to-use data visualization web server for scientific researchers. iMeta 2022, 1, e5. [Google Scholar] [CrossRef]
  51. Parks, D.H.; Tyson, G.W.; Hugenholtz, P.; Beiko, R.G. STAMP: Statistical analysis of taxonomic and functional profiles. Bioinformatics 2014, 30, 3123–3124. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  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, 1006–1010. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  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, 2109–2122. [Google Scholar] [CrossRef]
  58. Videvall, E. Microbiome maturation during a unique developmental window. Mol. Ecol. 2020, 29, 1941–1943. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  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, 5611–5622. [Google Scholar] [CrossRef]
  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, 8787–8803. [Google Scholar] [CrossRef]
  62. Macfarlane, S.; Macfarlane, G.T. Regulation of short-chain fatty acid production. Proc. Nutr. Soc. 2003, 62, 67–72. [Google Scholar] [CrossRef]
  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, 501–509. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef] [PubMed]
  65. Feng, X.L. Characteristics of nutritional metabolism and energy balance of muskrats. Mod. Livest. Poult. Breed. Ind. 2012, 10, 48–49. (In Chinese) [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  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, 971–982. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  69. Konikoff, T.; Gophna, U. Oscillospira: A Central, Enigmatic Component of the Human Gut Microbiota. Trends MicroBiol. 2016, 24, 523–524. [Google Scholar] [CrossRef]
  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, 535–544. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
Metabolites 13 00248 g001 550
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.
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.
Metabolites 13 00248 g001
Metabolites 13 00248 g002 550
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.
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.
Metabolites 13 00248 g002
Metabolites 13 00248 g003 550
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.
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.
Metabolites 13 00248 g003
Metabolites 13 00248 g004 550
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).
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).
Metabolites 13 00248 g004
Metabolites 13 00248 g005 550
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.
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.
Metabolites 13 00248 g005
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.

Share and Cite

MDPI and ACS Style

Song, F.; Xu, Y.; Peng, P.; Li, H.; Zheng, R.; Zhang, H.; Han, Y.; Weng, Q.; Yuan, Z. Seasonal Changes in the Structure and Function of Gut Microbiota in the Muskrat (Ondatra zibethicus). Metabolites 2023, 13, 248. https://doi.org/10.3390/metabo13020248

AMA Style

Song F, Xu Y, Peng P, Li H, Zheng R, Zhang H, Han Y, Weng Q, Yuan Z. Seasonal Changes in the Structure and Function of Gut Microbiota in the Muskrat (Ondatra zibethicus). Metabolites. 2023; 13(2):248. https://doi.org/10.3390/metabo13020248

Chicago/Turabian Style

Song, Fengcheng, Yishu Xu, Peng Peng, Hongxu Li, Ranxi Zheng, Haolin Zhang, Yingying Han, Qiang Weng, and Zhengrong Yuan. 2023. "Seasonal Changes in the Structure and Function of Gut Microbiota in the Muskrat (Ondatra zibethicus)" Metabolites 13, no. 2: 248. https://doi.org/10.3390/metabo13020248

APA Style

Song, F., Xu, Y., Peng, P., Li, H., Zheng, R., Zhang, H., Han, Y., Weng, Q., & Yuan, Z. (2023). Seasonal Changes in the Structure and Function of Gut Microbiota in the Muskrat (Ondatra zibethicus). Metabolites, 13(2), 248. https://doi.org/10.3390/metabo13020248

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop