Next Article in Journal
Do You Think I Am Living Well? A Four-Season Hair Cortisol Analysis on Leisure Horses in Different Housing and Management Conditions
Next Article in Special Issue
Effect of Inulin Source and a Probiotic Supplement in Pig Diets on Carcass Traits, Meat Quality and Fatty Acid Composition in Finishing Pigs
Previous Article in Journal
The Effect of Demonstrator Social Rank on the Attentiveness and Motivation of Pigs to Positively Interact with Their Human Caretakers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Diversity, Composition, and Metabolic Pathways of Archaea in Pigs

1
Guangdong Provincial Key Laboratory of Animal Molecular Design and Precise Breeding, College of Life Science and Engineering, Foshan University, Foshan 528225, China
2
School of Life Science and Engineering, Foshan University, Foshan 528225, China
3
Department of Animal Science, Division of Agriculture, University of Arkansas, Fayetteville, AR 72701, USA
4
State Key Laboratory for Managing Biotic and Chemical Threats to the Quality and Safety of Agro-Products, Institute of Agro-Product Safety and Nutrition, Zhejiang Academy of Agricultural Sciences, Hangzhou 310021, China
*
Author to whom correspondence should be addressed.
Submission received: 8 June 2021 / Revised: 14 July 2021 / Accepted: 15 July 2021 / Published: 20 July 2021

Abstract

:

Simple Summary

Archaea is identified as the key link in the interaction between gut microbiota and host metabo-lism. Studies on human and mice have reported archaea, especially methanogenic archaea, makes an important impact on the energy harvesting capacity of the host by improving fermentation. But, in pigs, the metabolic potential of archaea at different production stages are still largely unknown. Herein, we re-analyzed 276 metagenomic samples to explore the diversity, composi-tion, and potential functions of archaea in pigs. The results showed significant regional variations in archaeal composition. Furthermore, the Metacyc pathway related to hydrogen consumption (METHANOGENESIS-PWY) was only observed in archaeal reads, and archaea may be involved in carbohydrate metabolism and de novo synthesis of some kinds of essential amino acid. Overall, metagenomic re-analysis revealed that the composition and functional potential of archaea in the swine gut and suggested that archaea may make an important function in pigs.

Abstract

Archaea are an essential class of gut microorganisms in humans and animals. Despite the substantial progress in gut microbiome research in the last decade, most studies have focused on bacteria, and little is known about archaea in mammals. In this study, we investigated the composition, diversity, and functional potential of gut archaeal communities in pigs by re-analyzing a published metagenomic dataset including a total of 276 fecal samples from three countries: China (n = 76), Denmark (n = 100), and France (n = 100). For alpha diversity (Shannon Index) of the archaeal communities, Chinese pigs were less diverse than Danish and French pigs (p < 0.001). Consistently, Chinese pigs also possessed different archaeal community structures from the other two groups based on the Bray–Curtis distance matrix. Methanobrevibacter was the most dominant archaeal genus in Chinese pigs (44.94%) and French pigs (15.41%), while Candidatus methanomethylophilus was the most predominant in Danish pigs (15.71%). At the species level, the relative abundance of Candidatus methanomethylophilus alvus, Natrialbaceae archaeon XQ INN 246, and Methanobrevibacter gottschalkii were greatest in Danish, French, and Chinese pigs with a relative abundance of 14.32, 11.67, and 16.28%, respectively. In terms of metabolic potential, the top three pathways in the archaeal communities included the MetaCyc pathway related to the biosynthesis of L-valine, L-isoleucine, and isobutanol. Interestingly, the pathway related to hydrogen consumption (METHANOGENESIS-PWY) was only observed in archaeal reads, while the pathways participating in hydrogen production (FERMENTATION-PWY and PWY4LZ-257) were only detected in bacterial reads. Archaeal communities also possessed CAZyme gene families, with the top five being AA3, GH43, GT2, AA6, and CE9. In terms of antibiotic resistance genes (ARGs), the class of multidrug resistance was the most abundant ARG, accounting for 87.41% of archaeal ARG hits. Our study reveals the diverse composition and metabolic functions of archaea in pigs, suggesting that archaea might play important roles in swine nutrition and metabolism.

1. Introduction

Before Dr. Carl Woese and his colleagues separated archaea (named “Archaebacteria” at that time) from the bacteria domain based on phylogenetic analysis of ribosomal RNA sequences, archaea were considered a subgroup of bacteria [1]. To highlight the differences between archaea and bacteria, Dr. Carl Woese formally proposed the name of “Archaea” in his publication [2]. With the advancement of research technologies and tools, a large number of novel archaea and their characteristics have been discovered.
Archaea have highly diverse energy sources and unique metabolic characteristics and cell physiology, making it possible for them to live in various extreme environments, such as extreme temperature [3], high-salt [4], extreme alkaline [5], or acidic environments [6]. At the same time, archaea were also detected in various non-extreme ecosystems, including soil, ocean, lake water, and habitats associated with human and animal hosts, and were confirmed to make significant contributions to the ecological cycle by previous studies [7,8]. As a “young” group with unique natural capabilities and biological characteristics, archaea have attracted increasing research interest in the past decade.
Unlike bacteria, archaea are one of the least studied and least understood members of the human microbiota. One of the reasons is that archaea only account for a very small proportion of the microbiota. A recent study among East Asians living in South Korea showed that archaea were detected in 42.47% (381/897) and constituted 10.24 ± 4.58% of total bacteria and archaea among positive subjects [9]. However, the relatively small proportion of archaea does not mean they are not important. For example, some species of archaea have been reported to reduce Trimethylamine-N-oxide (TMAO), a harmful product associated with several diseases [10], through methanogenesis. In addition, methanogenic archaea can remove H2, the end product of bacterial fermentation, to increase the efficiency of fermentation and improve the energy harvesting capacity of the host [11]. Two studies performed on adults and children report that hydrogen-utilizing methanogenic archaea are co-enriched with H2-producing bacteria, Prevotella species, in the gut of obese humans [12]. Furthermore, human gut colonization by archaea has been associated with several other gastrointestinal and metabolic diseases, such as chronic constipation [13], inflammatory bowel disease [14], and colorectal cancer [15].
In pigs, the major ecological niches harboring archaea along the gastrointestinal tract (GIT) include the cecum, colon, and rectum, according to Gresse et al. [16]. The genus Methanobrevibacter is predominant with a relative abundance of 57–100% of the archaeal communities in pigs and sows [16,17,18,19]. Additionally, the archaeal composition in the swine GIT is dynamic, according to Su et al. [17] and Federici et al. [20], and diet-driven [21]. These exploratory studies provided important information on archaea in pigs, however, their sample sizes were small, and the composition and metabolic potential of archaea at different production stages in pigs are still largely unknown.
Xiao et al. [22] established a swine gut microbiome reference gene catalog with metagenome deep sequencing of fecal samples from 287 pigs. In this study, we re-analyzed their data and investigated the diversity, composition, and potential functions of archaea in pigs.

2. Materials and Methods

2.1. Obtaining Data (Metadata (from Nature Microbiology))

A total of 287 fecal samples collected from China (n = 87), France (n = 100), and Denmark (n = 100) were sequenced using shotgun metagenomic sequencing technology by Xiao et al. [22], of them, 10 samples from China were omitted because of low sequencing depth, and we downloaded 276 of these sequences from the European Nucleotide Archive (ENA) under accession code PRJEB11755. We filtered 11 samples due to their low sequencing depths. In addition, we downloaded the predicted gene catalog generated by Liang et al., from GiGaDB (http://gigadb.org/dataset/view/id/100187/token/F4CDHYruxobOKmsE, accessed on 11 November 2020) to annotate the CAZy genes.

2.2. Raw Reads Pre-Processing

Downloaded raw reads were filtered using Kneaddata v0.7.2 (https://huttenhower.sph.harvard.edu/kneaddata/), an integrated pipeline used to perform quality control and remove host contamination from shotgun metagenomic data. In brief, reads were trimmed using Trimmomatic v0.39 [23] with parameter SLIDINGWINDOW:4:20 and MINLEN:60. Potential host contamination was identified and removed by mapping all reads to the swine reference genome (Scrofa 11.1, ftp://ftp.ncbi.nlm.nirefseq/vertebrate_mammalian/h.gov/genomes/Sus_scrofa/representative/GCF_000003025.6_Sscrofa11.1) using bmtagger v.3.102.4 [24] with default parameters. The number of clean metagenomic reads for fecal samples was rarefied to 16,230,988 paired reads for each sample for a consistent sequencing depth by random subsampling.

2.3. Taxonomic Profiling and Diversity

Kraken2 [25], a software used for short reads classification with the lowest common ancestor (LCA) algorithms, was used to classify clean short reads against the bacterial genome and archaeal genome databases. The bacterial genome database was downloaded using the Kraken2 standard command, which downloads NCBI taxonomic information and complete bacterial genomes in RefSeq (1 July 2020) [26]. Then, clean reads were classified as bacterial taxa using Kraken2, and reads identified as bacteria were retained for downstream analysis. The method used for archaeal classification was slightly different compared with that of bacteria. The archaeal genome database was built using Kraken2 from two sources: (1) RefSeq (1 July 2020) and (2) a curated archaeal reference dataset [27]. Similarly, Kraken2 was used to assign clean reads into proper archaeal taxa based on this archaeal database, and archaeal reads were separated and retained.
Reads count table of archaea and bacteria at the species level were entered into QIIME2 [28] platform to rarefy and calculate alpha diversity (Shannon Index) and beta diversity (Bray–Curtis).

2.4. CAZyme and ARG Analysis

HUMAnN3 [29] pipeline was used to characterize the metabolic pathways of archaea and bacteria using the MetaCyc database.
CAZyme family annotation was performed on the downloaded gene catalog generated by Xiao et al. via alignment with the dbCAN database [30] by HMMER3 software [31]. The genes identified as CAZy genes were extracted as target genes for further quantitative analysis. Salmon (quant) [32] was run on separated archaeal reads to quantify the archaeal CAZy genes in each sample, and a similar analysis was run on the separated bacterial reads.
Deeparg-SS [33] pipeline was adapted to predict antibiotic resistance genes in both bacterial and archaeal reads with the following critical parameters (--min-prob 0.8 --arg-alignment-identity 30 --arg-alignment-evalue 1e-10).

2.5. Statistical Analysis

For all analyses, statistical significance was determined at p < 0.05 and p < 0.001. Kruskal–Wallis (pairwise) was performed to explore differences in alpha diversities among different countries. Analysis of similarities (ANOSIM) was carried out to test the statistical significance of beta diversity among different countries. Kruskal–Wallis (pairwise) and ANOSIM were run in QIIME2.

3. Results

3.1. The Diversity of Archaea in Pigs

Significant differences in archaeal alpha diversity between different countries were observed. Pairwise comparisons of countries showed that archaeal diversity (Shannon Index) of Chinese pigs were significantly lower than that of the French (Figure 1A: Kruskal–Wallis test, p < 0.001) and Danish groups (p < 0.001). No differences in the swine archaeal alpha diversity were observed between French and Danish (p = 0.918) pigs. Beta diversity of the archaeal communities was visualized on a principal coordinate analysis (PCoA) plot based on the Bray–Curtis distance matrix. Chinese pigs had significantly different archaeal communities from Danish (ANOSIM, p < 0.001, R = 0.808) and French (ANOSIM, p < 0.001, R = 0.634) pigs as demonstrated on the PCoA plot (Figure 1B). However, the archaeal community structures of the Danish and French pigs could not be distinguished (Figure 1B; ANOSIM, p < 0.001, R = 0.109).

3.2. The Archaeal Composition in Pigs

At the genus level, Methanobrevibacter, with an overall average relative abundance of 23.02%, was the most abundant archaeal genus (Figure 2), followed by Candidatus methanomethylophilus, which had a relative abundance of 10.52%. Of note, although Methanobrevibacter is considered the most dominant archaeal genus of pigs in general, it was ranked the second most abundant in Danish pigs. Additionally, Methanobrevibacter was significantly more abundant in Chinese pigs with the relative abundance of 44.94% than in both Danish (13.97%) and French pigs (15.41%), while Candidatus methanomethylophilus was decreased in Chinese pigs (0.98%) compared to Danish (15.71%) and French (12.59%) pigs.
At the species level, the mean relative abundance of the top 15 species is shown in Figure 3. The gut microbiome of pigs is dominated by Natrialbaceae archaeon XQ INN 246 (mean relative abundance, 10.27%), followed by Candidatus methanomethylophilus alvus (9.57%), Methanobrevibacter gottschalkii (6.40%), and Methanobrevibacter smithii (3.42%). The most abundant species in different countries seem to be inconsistent, as well. The relative abundance of Candidatus methanomethylophilus alvus, Natrialbaceae archaeon XQ INN 246, and Methanobrevibacter gottschalkii were the highest in Danish, French, and Chinese pigs with an average relative abundance of 14.32, 11.67, and 16.28%, respectively.

3.3. The metabolic Pathways of the Swine Archaea

To investigate the possible roles of gut archaea in swine, MetaCyc metabolic pathway profiles based on archaeal reads were investigated using HUMAnN3. A total of 77 MetaCyc pathways were observed, with 12 of these MetaCyc pathways present in more than 60% of fecal samples (Table S1). Figure 4A shows the top six most abundant MetaCyc pathways. The most abundant MetaCyc pathway was involved in pyruvate fermentation to isobutanol (PWY-7111). Pathways of VALSYN-PWY and ILEUSYN-PWY ranked as second and third and were related to the biosynthesis of two essential amino acids, L-valine and L-isoleucine, respectively. Interestingly, these pathways, PWY-7111, VALSYN-PWY, and ILEUSYN-PWY, were also among the top 10 most abundant MetaCyc pathways in bacteria (Figure 4B, Table S2). Furthermore, the pathway related to hydrogen consumption (METHANOGENESIS-PWY) was only observed in archaea and not in bacteria (Tables S1 and S2). In contrast, the pathways involved in hydrogen production (FERMENTATION-PWY and PWY4LZ-257) were only observed in bacteria and not in archaea (Figure 5).
We next profiled the archaea-annotated CAZyme gene families in the swine gut. A total of 297 CAZyme gene families were identified from archaeal reads, out of which 99 were detected with a sample coverage of more than 60%. The three major CAZyme families were glycoside hydrolases (GH, 54/99), glycoside transferases (GT, 24/99), and carboxylesterase (CE, 11/99). The most abundant CAZyme category was AA3 (auxiliary activity family AA3), followed by GH43 and GT2 (Figure 6A). For comparison, we also investigated bacterial CAZyme profiles (Figure 6B) and found that six of the top 10 CAZyme families of archaea and bacteria were shared, including AA3, GT2, AA6, GT41, GT5, and GH13.
Finally, we examined the ARGs in the swine gut archaea. In total, 13 classes of ARGs were detected from archaeal reads in all samples. On average, 0.035% of ARG hits were obtained from archaeal reads. The ARG class of multidrug resistance represented the vast majority (87.41%) of ARG hits; of the ARG class of multidrug resistance, the proportion of rpoB2 gene ranged from 87.18 to 100%. The other ARG classes with more than 1% relative abundance included peptide (5.07%), tetracycline (2.29%), mupirocin (1.69%), and unclassified (1.15%). The ARG type of multidrug class was significant enriched in Chinese pigs (Kruskal-Wallis test, p < 0.01), and the relative abundance of the multidrug class in Danish pigs is significantly higher than that in French pigs (Kruskal–Wallis test, p < 0.01). The boxplots of the top six ARG genes are shown in Figure 6C. In contrast, 24 ARGs were identified in bacterial reads, and 0.144% of bacterial reads were identified as ARG hits. The most abundant ARG class in bacteria was tetracycline, accounting for 40.59% of the total ARG hits (Figure 6D), while the class of multidrug ranked second with 19.39%, followed by MLS (17.73%), beta-lactam (7.46%), peptide (5.33%), and aminoglycoside (4.63%).

4. Discussion

Although archaea are an essential class of gut microorganisms in humans and animals, they have been neglected for many years. A few studies have examined archaeal composition by 16S rRNA sequencing methods; however, most studies only focused on the composition of methanogenic archaea of gut archaea. In fact, there are many non-methanogenic archaea such as Thermococcus, Halorubrum, and Halococcus in the swine gut. Furthermore, only a small number of gut archaeal strains were isolated from pig feces [34]. Therefore, compared to bacteria, very limited information is available for the archaeal composition in pigs, and even less is known about the archaeal metabolic pathways in pigs. In this study, by re-analyzing the swine gut metagenome of a large number of pigs from three countries, we characterized the swine gut archaeal diversity and metabolic pathways, which improves our understanding of the potential roles of archaea in mammals. However, it is also worth noting that the lower abundant archaeal genes detected by metagenome analysis could be highly expressed when examined using metatranscriptomic approaches.
Regarding archaeal diversity, Chinese pigs had archaeal communities significantly different from Danish and French pigs in both alpha and beta diversity, consistent with the bacterial communities described by Xiao and colleagues [22]. Chinese pigs were fed two types of antibiotics, including fluorine benzene (400 g/ton) and penicillin (600g/ton), which may cause a significant reduction in microbial richness. However, most archaea were not sensitive to penicillin and fluorine benzene because of the absence of drug targets [35], which means that the administration of antibiotics could not fully explain the decrease of archaeal richness in Chinese pigs. However, it is well known that antibiotics affect bacterial community diversities, which might in turn have influenced the archaeal community due to bacterial–archaeal interactions. In addition, geography, environment, diet, sex, and breed are confounding factors that might contribute to the difference in alpha and beta diversity of the swine gut archaeal diversity between the three countries for archaeal composition, Methanobrevibacter and Candidatus methanomethylophilus were the top two most abundant genera. Methanobrevibacter is a hydrogen-oxidizing genus that produces methane from carbon dioxide and hydrogen [36], while Candidatus methanomethylophilus synthesizes methane through the reduction of trimethylamine and methanol [37]. Methanobrevibacter was the most abundant archaeal genus in general, which was in line with previous studies [16]. It should be noted that variations in the abundance of Methanobrevibacter between countries were observed. It was ranked the second most abundant genus in the Danish group with a relative abundance of 13.97%, following Candidatus methanomethylophilus (relative abundance of 15.71%). The archaeal composition based on metagenomics in this study was different from that revealed by 16S rRNA or mcrA genes in other studies. Mi and colleagues found that the genus Methanobrevibacter was the dominant archaeal taxon with an average relative abundance of 57% in the large intestine of finishing pigs [19]. Furthermore, Methanobrevibacter was also reported as the overwhelmingly dominant archaeal genus and accounted for 95.01–100% of all archaeal reads in pigs in other studies [17,18]. Many factors, such as environment, diet, breed, and age, might contribute to the discrepancies in the archaeal composition between these studies. It is worth noting that technology (i.e., shotgun metagenomics sequencing vs. 16S rRNA sequencing) is likely another important confounder. Amplification bias is an inevitable disadvantage to accurately quantify bacterial and archaeal cells [38]. Of particular note, compared with bacteria, the domain archaea remains poorly understood, and there is a lack of sufficient reference sequences to design primers properly, which causes many archaeal taxa to remain undetected [39]. In addition, the lack of archaeal reference sequences not only affects proper primer design, but also hinders the annotation of archaeal genes, resulting in the elusive view of archaeal functions in pigs.
Previous studies on animals and humans have revealed that archaea were involved in energy metabolism through the consumption of hydrogen and carbon dioxide to produce methane [11,40]. Consistent with those studies, the metabolic pathway related to methanogenesis from H2 and CO2 (METHANOGENESIS-PWY) was only found in archaeal reads. In contrast, hydrogen production pathways were only found in bacterial reads and not in archaeal reads. These results suggest that there is a complementary relationship between archaea and bacteria, and archaea participated in energy metabolism in pigs as well. Samuel and colleagues showed that methanogenic archaea were able to enhance caloric harvest and weight gain by methods other than hydrogen utilization in humans, but the detailed mechanism is unclear [40]. We inferred that archaea have the potential to improve the growth performance of pigs through the synthesis of essential amino acids. In our study, L-valine and L-isoleucine-related pathways were ranked as the top MetaCyc pathways of archaea. L-valine and L-isoleucine are essential amino acids and could improve the growth performance of pigs at multiple stages of growth [41,42]. Of note, the two pathways were also abundant in bacteria. Further studies are desired to elucidate the proportion and magnitude of archaeal contribution to amino acid biosynthesis in pigs.
Interestingly, a large number of CAZyme gene families were also detected in archaeal reads. CAZyme genes are involved in carbohydrate metabolism and are important in the absorption and utilization efficiency of energy for the host. CAZyme gene families abundance in the host gut were confirmed association with growth performance in several domestic animals [43,44,45]. In general, bacteria in the mammal gut, especially species belonging to Firmicutes and Bacteroidetes, were considered the main providers of CAZyme genes for encoding enzymes to digest complex carbohydrates [46]. Nevertheless, recent publications revealed that archaea participate in carbohydrate metabolism through the production of various enzymes [47,48]. In our study, 297 CAZyme gene families were detected in archaea reads, and six of the top 10 CAZyme families of archaea and bacteria were shared. These CAZyme gene families of archaea may be useful supplements for carbohydrate digestion in swine. Therefore, our study suggests that, in addition to H2 depletion, archaea might be directly involved in energy metabolism in the pig gut.
Bacterial ARGs have been investigated in many studies in different animals. However, little is known about ARGs in archaea. In this study, we also explored the archaeal ARGs in pigs in addition to metabolic pathways. The ARG classes of tetracycline, multidrug, and MLS were relatively abundant in bacterial reads of the 287 pigs. However, the ARG class of multidrug resistance represented the overwhelming majority at 87.41% of ARG hits in archaeal reads of these pigs. RpoB2 gene accounted for 87.18–100% of the multidrug resistance type. The rpoB2 gene confers resistance to rifamycin through influencing Rif binding, which decreases antibiotic activity of rifamycin. The effects of low-dose antibiotics on archaeal ARGs in pigs, the generation and spread of ARGs among archaeal and between domains need further investigation.
In this study, we profiled archaeal composition and functional potential in pigs. Although amplification bias was avoided by using shotgun metagenomics data, there are two limitations of our study. First, our study only included a limited number of archaeal reference genomes. Therefore, numerous archaea with no reference genomes could not be identified, which can underestimate the richness of archaea and its functional diversity. Second, compared with bacteria, the relative abundance of archaea is lower by about 1%. Therefore, higher sequencing depth was demanded to quantify archaea accurately.
Although we demonstrated archaea had functional potentials, further studies are needed to link archaeal species or strains with the growth performance of pigs. In addition, methane, the production of methanogenic archaea, is a greenhouse gas with environmental concerns. Thus, the modulation of archaeal communities to maintain a balance between animal production (e.g., growth-promoting functions) and environmental protection (e.g., methane emission) is critical.

5. Conclusions

In conclusion, our study revealed the composition and functional potential of archaea in the swine gut based on the re-analysis of metagenomic data. We have demonstrated different predominant archaeal genera and species in different swine populations. In terms of functional potential, METHANOGENESIS-PWY related to methanogenesis from H2 and CO2 was only found in archaeal reads, which indicated that archaea were involved in energy metabolism by depleting H2 in the swine gut. Moreover, archaea shared several amino acid synthesis pathways and the main CAZyme gene families with bacteria, which suggested archaea are potentially involved in amino acid synthesis and energy metabolism.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/ani11072139/s1, Table S1. MetaCyc pathway detected in archaea reads. Table S2. MetaCyc pathway detected in bacteria reads

Author Contributions

Conceptualization, Y.L. (Ying Li); data curation, F.D., Y.L. (Yushan Li) and X.W. (Xiaoyuan Wei); formal analysis, Y.P., X.W. (Xiaoyuan Wei), X.W. (Xiaofan Wang), H.Y., Y.X. and J.Z.; methodology, F.D. and J.Z.; software, H.L.; writing—original draft, F.D.; writing—review and editing, S.H., J.Z. and Y.L. (Ying Li). All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Guangdong Provincial Key Laboratory of Animal Molecular Design and Precise Breeding (2019B030301010), and the Key Laboratory of Animal Molecular Design and Precise Breeding of Guangdong Higher Education Institutes (2019KSYS011). Computational support was provided by the Extreme Science and Engineering Discovery Environment (XSEDE) [49], which is supported by the National Science Foundation grant number ACI-1548562. Specifically, it used the Bridges-2 system, which is supported by NSF award number ACI-1445606, at the Pittsburgh Supercomputing Center (PSC).

Institutional Review Board Statement

Not applicable. All data used in this manuscript were publicly published by other groups. We re-analyzed those data, so do not involve in any experiment that need ethic approval.

Data Availability Statement

The data generated and analyzed during this study are included in this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Woese, C.R.; Fox, G.E. Phylogenetic structure of the prokaryotic domain: The primary kingdoms. Proc. Natl. Acad. Sci. USA 1977, 74, 5088–5090. [Google Scholar] [CrossRef] [Green Version]
  2. Woese, C.R.; Kandler, O.; Wheelis, M.L. Towards a natural system of organisms: Proposal for the domains Archaea, Bacteria, and Eucarya. Proc. Natl. Acad. Sci. USA 1990, 87, 4576–4579. [Google Scholar] [CrossRef] [Green Version]
  3. Cavicchioli, R. Cold-adapted archaea. Nat. Rev. Microbiol. 2006, 4, 331–343. [Google Scholar] [CrossRef]
  4. Andrei, A.-Ş.; Banciu, H.L.; Oren, A. Living with salt: Metabolic and phylogenetic diversity of archaea inhabiting saline ecosystems. FEMS Microbiol. Lett. 2012, 330, 1–9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Shen, J.p.; Zhang, L.m.; Zhu, Y.g.; Zhang, J.b.; He, J.z. Abundance and composition of ammonia-oxidizing bacteria and ammonia-oxidizing archaea communities of an alkaline sandy loam. Environ. Microbiol. 2008, 10, 1601–1611. [Google Scholar] [CrossRef]
  6. Herbold, C.W.; Lehtovirta-Morley, L.E.; Jung, M.Y.; Jehmlich, N.; Hausmann, B.; Han, P.; Loy, A.; Pester, M.; Sayavedra-Soto, L.A.; Rhee, S.K. Ammonia-oxidising archaea living at low pH: Insights from comparative genomics. Environ. Microbiol. 2017, 19, 4939–4952. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Moissl-Eichinger, C.; Pausan, M.; Taffner, J.; Berg, G.; Bang, C.; Schmitz, R.A. Archaea are interactive components of complex microbiomes. Trends Microbiol. 2018, 26, 70–85. [Google Scholar] [CrossRef] [PubMed]
  8. Borrel, G.; Brugere, J.-F.; Gribaldo, S.; Schmitz, R.A.; Moissl-Eichinger, C. The host-associated archaeome. Nat. Rev. Microbiol. 2020, 18, 622–636. [Google Scholar] [CrossRef]
  9. Kim, J.Y.; Whon, T.W.; Lim, M.Y.; Kim, Y.B.; Kim, N.; Kwon, M.-S.; Kim, J.; Lee, S.H.; Choi, H.-J.; Nam, I.-H. The human gut archaeome: Identification of diverse haloarchaea in Korean subjects. Microbiome 2020, 8, 114. [Google Scholar] [CrossRef] [PubMed]
  10. Bain, M.A.; Faull, R.; Fornasini, G.; Milne, R.W.; Evans, A.M. Accumulation of trimethylamine and trimethylamine-N-oxide in end-stage renal disease patients undergoing haemodialysis. Nephrol. Dial. Transplant. 2006, 21, 1300–1304. [Google Scholar] [CrossRef] [Green Version]
  11. Turnbaugh, P.J.; Ley, R.E.; Mahowald, M.A.; Magrini, V.; Mardis, E.R.; Gordon, J.I. An obesity-associated gut microbiome with increased capacity for energy harvest. Nature 2006, 444, 1027. [Google Scholar] [CrossRef] [PubMed]
  12. Zhang, H.; DiBaise, J.K.; Zuccolo, A.; Kudrna, D.; Braidotti, M.; Yu, Y.; Parameswaran, P.; Crowell, M.D.; Wing, R.; Rittmann, B.E. Human gut microbiota in obesity and after gastric bypass. Proc. Natl. Acad. Sci. USA 2009, 106, 2365–2370. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Pimentel, M.; Gunsalus, R.P.; Rao, S.S.; Zhang, H. Methanogens in human health and disease. Am. J. Gastroenterol. Suppl. 2012, 1, 28. [Google Scholar] [CrossRef]
  14. Khelaifia, S.; Raoult, D. Haloferax massiliensis sp. nov., the first human-associated halophilic archaea. New Microbes New Infect. 2016, 12, 96. [Google Scholar] [CrossRef] [Green Version]
  15. Coker, O.O.; Wu, W.K.K.; Wong, S.H.; Sung, J.J.; Yu, J. Altered gut archaea composition and interaction with bacteria are associated with colorectal cancer. Gastroenterology 2020, 159, 1459–1470.e1455. [Google Scholar] [CrossRef] [PubMed]
  16. Gresse, R.; Chaucheyras Durand, F.; Dunière, L.; Blanquet-Diot, S.; Forano, E. Microbiota composition and functional profiling throughout the gastrointestinal tract of commercial weaning piglets. Microorganisms 2019, 7, 343. [Google Scholar] [CrossRef] [Green Version]
  17. Su, Y.; Bian, G.; Zhu, Z.; Smidt, H.; Zhu, W. Early methanogenic colonisation in the faeces of Meishan and Yorkshire piglets as determined by pyrosequencing analysis. Archaea 2014, 2014, 547908. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Luo, Y.-h.; Su, Y.; Wright, A.-D.G.; Zhang, L.-l.; Smidt, H.; Zhu, W.-Y. Lean breed Landrace pigs harbor fecal methanogens at higher diversity and density than obese breed Erhualian pigs. Archaea 2012, 2012, 605289. [Google Scholar] [CrossRef]
  19. Mi, J.; Peng, H.; Wu, Y.; Wang, Y.; Liao, X. Diversity and community of methanogens in the large intestine of finishing pigs. BMC Microbiol. 2019, 19, 83. [Google Scholar] [CrossRef]
  20. Federici, S.; Miragoli, F.; Pisacane, V.; Rebecchi, A.; Morelli, L.; Callegari, M.L. Archaeal microbiota population in piglet feces shifts in response to weaning: Methanobrevibacter smithii is replaced with Methanobrevibacter boviskoreani. FEMS Microbiol. Lett. 2015, 362, fnv064. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Luo, Y.; Chen, H.; Yu, B.; He, J.; Zheng, P.; Mao, X.; Tian, G.; Yu, J.; Huang, Z.; Luo, J. Dietary pea fiber increases diversity of colonic methanogens of pigs with a shift from Methanobrevibacter to Methanomassiliicoccus-like genus and change in numbers of three hydrogenotrophs. BMC Microbiol. 2017, 17, 17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Xiao, L.; Estellé, J.; Kiilerich, P.; Ramayo-Caldas, Y.; Xia, Z.; Feng, Q.; Liang, S.; Pedersen, A.Ø.; Kjeldsen, N.J.; Liu, C. A reference gene catalogue of the pig gut microbiome. Nat. Microbiol. 2016, 1, 16161. [Google Scholar] [CrossRef]
  23. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. BMTagger: Best Match Tagger for Removing Human Reads from Metagenomics Datasets. Available online: Ftp://ftp.ncbi.nlm.nih.gov/pub/agarwala/bmtagger/ (accessed on 1 July 2020).
  25. Wood, D.E.; Lu, J.; Langmead, B. Improved metagenomic analysis with Kraken 2. Genome Biol. 2019, 20, 257. [Google Scholar] [CrossRef] [Green Version]
  26. Haft, D.H.; DiCuccio, M.; Badretdin, A.; Brover, V.; Chetvernin, V.; O’Neill, K.; Li, W.; Chitsaz, F.; Derbyshire, M.K.; Gonzales, N.R. RefSeq: An update on prokaryotic genome annotation and curation. Nucleic Acids Res. 2018, 46, D851–D860. [Google Scholar] [CrossRef]
  27. Kirstahler, P.; Bjerrum, S.S.; Friis-Møller, A.; La Cour, M.; Aarestrup, F.M.; Westh, H.; Pamp, S.J. Genomics-based identification of microorganisms in human ocular body fluid. Sci. Rep. 2018, 8, 1–14. [Google Scholar] [CrossRef] [Green Version]
  28. Bolyen, E.; Rideout, J.R.; Dillon, M.R.; Bokulich, N.A.; Abnet, C.C.; Al-Ghalith, G.A.; Alexander, H.; Alm, E.J.; Arumugam, M.; Asnicar, F. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 2019, 37, 852–857. [Google Scholar] [CrossRef]
  29. Franzosa, E.A.; McIver, L.J.; Rahnavard, G.; Thompson, L.R.; Schirmer, M.; Weingart, G.; Lipson, K.S.; Knight, R.; Caporaso, J.G.; Segata, N. Species-level functional profiling of metagenomes and metatranscriptomes. Nat. Methods 2018, 15, 962–968. [Google Scholar] [CrossRef]
  30. Zhang, H.; Yohe, T.; Huang, L.; Entwistle, S.; Wu, P.; Yang, Z.; Busk, P.K.; Xu, Y.; Yin, Y. dbCAN2: A meta server for automated carbohydrate-active enzyme annotation. Nucleic Acids Res. 2018, 46, W95–W101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Eddy, S.R. A new generation of homology search tools based on probabilistic inference. Genome Inform. 2009, 23, 205–211. [Google Scholar] [PubMed]
  32. Patro, R.; Duggal, G.; Love, M.I.; Irizarry, R.A.; Kingsford, C. Salmon provides fast and bias-aware quantification of transcript expression. Nat. Methods 2017, 14, 417–419. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Arango-Argoty, G.; Garner, E.; Pruden, A.; Heath, L.S.; Vikesland, P.; Zhang, L. DeepARG: A deep learning approach for predicting antibiotic resistance genes from metagenomic data. Microbiome 2018, 6, 1–15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Miller, T.L.; Wolin, M.; Kusel, E. Isolation and characterization of methanogens from animal feces. Syst. Appl. Microbiol. 1986, 8, 234–238. [Google Scholar] [CrossRef]
  35. Al-Mailem, D.; Sorkhoh, N.; Al-Awadhi, H.; Eliyas, M.; Radwan, S. Biodegradation of crude oil and pure hydrocarbons by extreme halophilic archaea from hypersaline coasts of the Arabian Gulf. Extremophiles 2010, 14, 321–328. [Google Scholar] [CrossRef] [PubMed]
  36. Dong, X.; Plugge, C.M.; Stams, A.J. Anaerobic degradation of propionate by a mesophilic acetogenic bacterium in coculture and triculture with different methanogens. Appl. Environ. Microbiol. 1994, 60, 2834–2838. [Google Scholar] [CrossRef] [Green Version]
  37. Noel, S.J.; Højberg, O.; Urich, T.; Poulsen, M. Draft genome sequence of “Candidatus Methanomethylophilus” sp. 1R26, enriched from bovine rumen, a methanogenic archaeon belonging to the Methanomassiliicoccales order. Genome Announc. 2016, 4, e01734-15. [Google Scholar]
  38. Kennedy, K.; Hall, M.W.; Lynch, M.D.; Moreno-Hagelsieb, G.; Neufeld, J.D. Evaluating bias of Illumina-based bacterial 16S rRNA gene profiles. Appl. Environ. Microbiol. 2014, 80, 5717–5722. [Google Scholar] [CrossRef] [Green Version]
  39. Eloe-Fadrosh, E.A.; Ivanova, N.N.; Woyke, T.; Kyrpides, N.C. Metagenomics uncovers gaps in amplicon-based detection of microbial diversity. Nat. Microbiol. 2016, 1, 1–4. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Samuel, B.S.; Gordon, J.I. A humanized gnotobiotic mouse model of host–archaeal–bacterial mutualism. Proc. Natl. Acad. Sci. USA 2006, 103, 10011–10016. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Mavromichalis, I.; Kerr, B.; Parr, T.; Albin, D.; Gabert, V.; Baker, D. Valine requirement of nursery pigs. Anim. Sci. J. 2001, 79, 1223–1229. [Google Scholar] [CrossRef] [PubMed]
  42. Figueroa, J.; Lewis, A.; Miller, P.; Fischer, R.; Diedrichsen, R. Growth, carcass traits, and plasma amino acid concentrations of gilts fed low-protein diets supplemented with amino acids including histidine, isoleucine, and valine. Anim. Sci. J. 2003, 81, 1529–1537. [Google Scholar] [CrossRef] [Green Version]
  43. Alfaia, C.; Pestana, J.; Rodrigues, M.; Coelho, D.; Aires, M.; Ribeiro, D.; Major, V.; Martins, C.; Santos, H.; Lopes, P. Influence of dietary Chlorella vulgaris and carbohydrate-active enzymes on growth performance, meat quality and lipid composition of broiler chickens. Poult. Sci. 2021, 100, 926–937. [Google Scholar] [CrossRef] [PubMed]
  44. Jiang, B.; Wang, T.; Zhou, Y.; Li, F. Effects of enzyme+ bacteria treatment on growth performance, rumen bacterial diversity, KEGG pathways, and the CAZy spectrum of Tan sheep. Bioengineered 2020, 11, 1221–1232. [Google Scholar] [CrossRef]
  45. Fang, S.; Chen, X.; Ye, X.; Zhou, L.; Xue, S.; Gan, Q. Effects of gut microbiome and short-chain fatty acids (SCFAs) on finishing weight of meat rabbits. Front. Microbiol. 2020, 11, 1835. [Google Scholar] [CrossRef] [PubMed]
  46. White, B.A.; Lamed, R.; Bayer, E.A.; Flint, H.J. Biomass utilization by gut microbiomes. Annu. Rev. Microbiol. 2014, 68, 279–296. [Google Scholar] [CrossRef] [PubMed]
  47. Kırtel, O.; Lescrinier, E.; Van den Ende, W.; Öner, E.T. Discovery of fructans in Archaea. Carbohydr. Polym. 2019, 220, 149–156. [Google Scholar] [CrossRef] [PubMed]
  48. Suleiman, M.; Krüger, A.; Antranikian, G. Biomass-degrading glycoside hydrolases of archaeal origin. Biotechnol. Biofuels 2020, 13, 1–14. [Google Scholar] [CrossRef]
  49. Towns, J.; Cockerill, T.; Dahan, M.; Foster, I.; Gaither, K.; Grimshaw, A.; Hazlewood, V.; Lathrop, S.; Lifka, D.; Peterson, G.D. XSEDE: Accelerating scientific discovery. Comput. Sci. Eng. 2014, 16, 62–74. [Google Scholar] [CrossRef]
Figure 1. Swine archaeal alpha and beta diversities between three countries. Alpha diversity of archaeal communities in pigs was measured using the Shannon Index (A). Beta diversity was illustrated on a principal coordinate analysis (PCoA) plot based on Bray–Curtis distance (B).
Figure 1. Swine archaeal alpha and beta diversities between three countries. Alpha diversity of archaeal communities in pigs was measured using the Shannon Index (A). Beta diversity was illustrated on a principal coordinate analysis (PCoA) plot based on Bray–Curtis distance (B).
Animals 11 02139 g001
Figure 2. Relative abundance of the top 15 archaeal genera in swine fecal samples from different countries. Each column represents one pig, and each color represents one archaeal genus. The y-axis represents relative abundance ranging from 0 to 100%.
Figure 2. Relative abundance of the top 15 archaeal genera in swine fecal samples from different countries. Each column represents one pig, and each color represents one archaeal genus. The y-axis represents relative abundance ranging from 0 to 100%.
Animals 11 02139 g002
Figure 3. Relative abundance of the top 15 archaeal species in swine fecal samples from different countries. Each column represents one pig, and each color represents one archaeal species. The y-axis represents relative abundance ranging from 0 to 100%.
Figure 3. Relative abundance of the top 15 archaeal species in swine fecal samples from different countries. Each column represents one pig, and each color represents one archaeal species. The y-axis represents relative abundance ranging from 0 to 100%.
Animals 11 02139 g003
Figure 4. Summary of the predominant MetaCyc pathways of archaea and bacteria in pigs. The top six most abundant MetaCyc pathways of archaea (A) and bacteria (B) are displayed. In (A,B), the y-axis represents the number of reads mapped to reference genes involved in pathways. Data have been normalized to RPKM (reads per kilobase of transcript per million mapped reads).
Figure 4. Summary of the predominant MetaCyc pathways of archaea and bacteria in pigs. The top six most abundant MetaCyc pathways of archaea (A) and bacteria (B) are displayed. In (A,B), the y-axis represents the number of reads mapped to reference genes involved in pathways. Data have been normalized to RPKM (reads per kilobase of transcript per million mapped reads).
Animals 11 02139 g004
Figure 5. Overview of the hydrogen production pathways of bacteria and consumption pathway of archaea.
Figure 5. Overview of the hydrogen production pathways of bacteria and consumption pathway of archaea.
Animals 11 02139 g005
Figure 6. Summary of predominant CAZyme gene families (A,B) and antibiotic resistance genes (C,D) of archaea and bacteria in pigs. The top six CAZy gene families of archaea (A) and bacteria (B) were normalized by TPM (trans per million) using salmon software. ARG raw counts were identified using deeparg software. The ARG reads count of archaea (C) and bacteria (D) were normalized by CPM (counts per million).
Figure 6. Summary of predominant CAZyme gene families (A,B) and antibiotic resistance genes (C,D) of archaea and bacteria in pigs. The top six CAZy gene families of archaea (A) and bacteria (B) were normalized by TPM (trans per million) using salmon software. ARG raw counts were identified using deeparg software. The ARG reads count of archaea (C) and bacteria (D) were normalized by CPM (counts per million).
Animals 11 02139 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Deng, F.; Li, Y.; Peng, Y.; Wei, X.; Wang, X.; Howe, S.; Yang, H.; Xiao, Y.; Li, H.; Zhao, J.; et al. The Diversity, Composition, and Metabolic Pathways of Archaea in Pigs. Animals 2021, 11, 2139. https://0-doi-org.brum.beds.ac.uk/10.3390/ani11072139

AMA Style

Deng F, Li Y, Peng Y, Wei X, Wang X, Howe S, Yang H, Xiao Y, Li H, Zhao J, et al. The Diversity, Composition, and Metabolic Pathways of Archaea in Pigs. Animals. 2021; 11(7):2139. https://0-doi-org.brum.beds.ac.uk/10.3390/ani11072139

Chicago/Turabian Style

Deng, Feilong, Yushan Li, Yunjuan Peng, Xiaoyuan Wei, Xiaofan Wang, Samantha Howe, Hua Yang, Yingping Xiao, Hua Li, Jiangchao Zhao, and et al. 2021. "The Diversity, Composition, and Metabolic Pathways of Archaea in Pigs" Animals 11, no. 7: 2139. https://0-doi-org.brum.beds.ac.uk/10.3390/ani11072139

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