A Bacterial Genome And Culture Collection of Gut Microbial in Weanling Piglet

doi:10.21203/rs.3.rs-988012/v1

Download PDF

Research

A Bacterial Genome And Culture Collection of Gut Microbial in Weanling Piglet

https://doi.org/10.21203/rs.3.rs-988012/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 16 Feb, 2022

Read the published version in Microbiology Spectrum →

Version 1

posted

You are reading this latest preprint version

Background: The microbiota hosted in the pig gastrointestinal tract are important to health of this biomedical model. However, the individual species and functional repertoires that make up the pig gut microbiome remain largely undefined.

Results: Here we comprehensively investigated the genomes and functions of the piglet gut microbiome using culture-based and metagenomics approaches. A collection included 266 cultured genomes and 482 metagenome-assembled genomes (MAGs) that were clustered to 428 species across 10 phyla was established. Among these clustered species, 333 genomes represent potential new species. Less matches between cultured genomes and MAGs revealed a substantial bias for the acquisition of reference genomes by the two strategies. Glycoside hydrolases was the dominant category of carbohydrate-active enzymes. 445 secondary metabolite biosynthetic genes were predicted from 292 genomes with bacteriocin being the most. Pan genome analysis of Limosilactobacillus reuteri uncover the biosynthesis of reuterin was strain-specific and the production was experimentally determined.

Conclusions: A total of 266 isolated bacterial genomes and 482 MAGs were obtained and investigated their functional repertoires. This study provides a comprehensive view of the microbiome composition and the function landscape of the gut of weanling piglets and a valuable bacterial resource for further experimentations.

General Microbiology

genome collection

metagenome-assembled genomes

weanling piglet

functional repertoires

limosilactobacillus reuteri

Pigs are economically important livestock, widely used as monogastric animal model for enteric microbiological studies, and serve as most consumed meat for human worldwide^{1, 2}. A large amount of microorganisms harbored in the gut of pig. The research of gut microbiome of pig has increasingly advanced our understanding of the role of the microbiota in feed conversion efficiency, pig health, and production in recent years^{3, 4}.

Metagenomic studies have uncovered the genes, species, and functional diversity of bacteria in mammalian intestines with the advantage in avoiding the time-consuming of traditional culture-based approach and the inaccuracy of 16S rRNA sequences analysis. A gene catalog consisting of 7.7 million non-redundant genes has been constructed that draw a basic function map of the pig gut microbiota⁵. PIGC, the pig integrated gene catalog, a recent expanded gene catalog comprised 17 million complete genes that provide an expanded resource for pig gut microbiome research⁶. Metagenomic binning is an essential tool for acquiring reference genomes readily, particularly for uncultured microbiota, that has been widely used in study of human gut microbiota, but rarely applied for pigs. However, misassembles and chimeric contigs from MAGs result in substantial biases for the analysis of the gut microbiome^{7, 8}. These metagenomic studies have greatly increased our understanding of pig gut microbiome^{6, 9}. Nevertheless, the lack of bacterial isolates and high-quality genomes limited our comprehensively understanding of the structure and function of the gut microbiota of pig and the development of probiotics for pig farming.

The culture-based approach has been well used for establishment of bacterial collection for the human gut microbiota¹⁰. Over 1,500 species have been successfully isolated from the human gut by using culturomics, introduced by Lagier et al. and 247 new species have been unveiled¹¹. Culture-based studies enable the acquisition of both live bacteria and high-quality reference genomes and provide experimental access for function exploration and intervention trials with probiotics. PiBAC, a pig intestinal bacterial collection, was constructed by cultivation of 110 species and described 38 novel species from the gut of 19 pigs in Germany, United States, and Canada¹². This study highlighted the importance and necessity of continuously isolating and characterizing microbial taxa from pig gut.

Here we performed a study of bacterial cultivation and genome sequencing along with deep metagenome sequencing of 14 samples, including ileum and colon contents, and feaces, from weanling piglets. A total of 266 isolated bacterial genomes and 482 MAGs were obtained and investigated their functional repertoires. This study provides a comprehensive view of the microbiome composition and the function landscape of the gut of weanling piglets and a valuable bacterial resource for further experimentations.

The metagenome sequencing described community composition of microbiota in different gut regions of weanling piglets

A total of 42 samples, including 15 ileum contents, 18 colon contents and 9 faeces, were collected from 42 weanling piglets from Guangdong Academy of Agricultural Sciences. Every 3 samples from one line were combined into one sample which ultimately resulted in 14 samples. To investigate the microbiota composition in each sample, we first performed deep metagenomic sequencing for these samples and generated 20 Gb of data for each sample on average. Notably, over 12% of reads mapped to bacteria that could not be classified at the species level, which represent novel species in the gut of weanling piglets (Supplementary Fig. 1). 99.8% of the bacterial reads were assigned into 10 phyla. The phylum Firmicutes dominated the gut microbiota in the ileum, colon, and faeces, accounting for 98.9%, 67%, and 76.8% of the bacterial sequences, respectively (Figs. 1a, b). We also discovered Lactobacillus were common present in 3 regions and accounted for the highest proportion in samples from faeces and ileum and some samples from colon (Supplementary Fig. 2a, b). The higher microbial diversity, including the Shannon index, was observed in the colon as compared with the faeces and ileum (Fig. 1c), consistent with previous studies^{13, 14}. The principal coordinates analysis showed that the microbiota composition of the faeces and colon were similar and were different from that of the ileum (Fig. 1d). These results indicated the regional similarities and differences in microbiota composition and the beneficial bacteria, like Lactobacillus dominated the gut microbiota of weanling piglets, were valuable resource for the furtherly cultured-based study.

Bacteria culturing and genome sequencing of gut of weanling piglets

To advance our understanding of the diversity of gut microbiota and acquire bacterial isolates from the weanling pig gut, we used a variety of culture methods, including nutrient-rich medium (Peptone Yeast Extract Glucose Medium, Modified, abbreviated as MPYG, with or without sheep blood), oligotrophic medium (R2A), Columbia Agar, spore medium, etc. In total, 25 culture conditions were used (Supplementary method). Most of the bacteria were cultured in sheep blood-MPYG, MPYG, and sheep blood-spore medium. A total of 1,476 strains were obtained and identified using 16S rRNA gene sequence, the detail information is provided in Supplementary Table 2. These isolates belonged to 5 different phyla, including Firmicutes, Bacteroidetes, Actinobacteria, Proteobacteria, and Fusobacteria (Supplementary Fig. 3). Limosilactobacillus reuteri was the most abundant, which represented 18% of the taxa isolated and half of these 1,476 bacteria were Lactobacillus (recently divided as Lactobacillus, Ligilactobacillus, and Limosilactobacillus), representing the beneficial bacteria in gut of piglets (Supplementary Fig. 4). Most of these bacteria were the first time isolated from the pig gut microbiota, which expanded the bacteria diversity of the gut microbiota of pig. Notably, the isolated bacteria from the ileum, colon, and faeces were much different with highest species diversity of bacteria from faeces sample, that consistent with the analysis of metagenome (Supplementary Fig. 4). All bacteria have been deposited in the China National Genebank (CNGB) Shenzhen for public accessibility.

Comparison of cultured genomes and MAGs

We selected 266 representative strains that covered the bacterial diversity of the 1,476 isolates for genome sequencing. After assembly, 266 high-quality genomes with completeness more than 90% and contamination less than 5% were obtained (Supplementary Table 3a), of which 195 genomes (73.31%) were more than 99% completeness (Supplementary Fig. 5a), indicating that the great majority of the genomes were relatively complete. Based on the 95% average nucleotide identity (ANI), we clustered the 266 genomes into 59 species-level clusters. The vast majority of clusters were Firmicutes (242 genomes, 52 clusters), followed by Proteobacteria (18 genomes, 3 clusters), Actinobacteriota (4 genomes, 2 clusters), Bacteroidota (1 genome, 1 cluster) and Fusobacteriota (1 genome, 1 cluster) according to the taxonomic annotation of these genomes (Fig. 2a). Notably, 10 out of 59 clusters (23 genomes) lacked a species-level match with the GTDB-Tk (The Genome Taxonomy Database Toolkit) database (Fig. 2a and Supplementary Table 3a).

By Assembling and binning metagenomic sequence data from 14 samples, we reconstructed 482 non-redundant MAGs with completeness more than 70% and contamination less than 5% (Supplementary Table 3b). Only 223 (46.27%) MAGs were high-quality genomes (completeness more than 90%), 11 of which with completeness more than 99% (Supplementary Fig. 5b), which indicated that the quality of the genomes generated by culture-based methods was generally higher than metagenomic assembly. To evaluate the species bias caused by the two methods, we clustered 489 high-quality genomes. The results showed that MAGs and isolated genomes had 74.63% and 25.42% of unique species-clusters, respectively (Supplementary Fig. 6a, b), illustrating that the combination of two methods is needed to more comprehensively represent the species of the pig gut microbiota. We thus clustered the 748 genomes into 428 species-clusters to integrate the reference genomes. According to the GTDB classification annotations, 5 clusters (5 MAGs) were Archaea, and the remaining clusters were bacteria from 10 phyla (Fig. 2a and Supplementary Table 4). It is worth noting that 77.80% of clusters (333 clusters) could not be matched with any existing species, which represent novel species, and 96.70% (322 clusters) of the novel species without isolate genome representative in this study (Fig. 2a).

Previous studies have used metagenomics or culture-based method to study the intestinal microbes of pigs. Wang et al.⁹ generated 360 high-quality assembled genomes for pig fecal microbiome. And Wylensek et al.¹² established the pig intestinal bacterial collection (PiBAC), a resource of cultured bacteria from the pig intestine. We evaluated the novelty of our genomes by mapping 748 genomes against these two reference datasets. The result showed that our MAGs and isolated genomes had 65.89% and 50.85% unique clusters (novel species), respectively (Fig. 2b, c and Supplementary Fig. 7a, b), which contribute new resources for the research of pig gut microbiota. In addition, most of the clusters were only detected in the MAG dataset, but not in any culture studies, reflecting the lack of culture-based studies of pig intestines.

The assembled genes contributed to the present gene catalog

We next predicted genes from the 14 metagenomic samples and generated a non-redundant gene catalog with a number of 5,283,405, which covered 22.01% of the reference gene catalog established by Xiao et al.⁵ (Fig. 3a, b). The samples from faeces contained the most genes, followed by the colon and ileum. Each sample contributed more than 50% of novel genes for the gene catalog. Among them, samples from ileum have the largest proportion of novel genes, followed by colon and faeces (Fig. 3a). Considering the samples of the reference gene catalog by Xiao et al. were derived from pig faeces, we thought that construct a more complete gene catalog of pig intestinal needs to include samples from multiple parts. Altogether, we expanded the present gene catalog to 10.91 Mb (Fig. 3b). In addition, the coverage of genes obtained by MAGs and isolated genomes were 12.97% and 1.29%, respectively, whereas 62.97% of the genes in the isolated genomes could not be detected by metagenomics (Fig. 3c). This indicates that the missed genes from metagenomic analyses can be detected by culture-dependent methods.

Functional insight of gut microbial in weanling piglets

The comprehensive gene catalog of pig intestines enables a higher-resolution functional analysis to better understand the interaction between gut microbes and pigs. We subsequently performed a KEGG pathway annotation of total protein sequences and found that 213 pathways were annotated in at least one genome. The most general pathway was Metabolic pathways, followed by Biosynthesis of secondary metabolites, Biosynthesis of antibiotics, Microbial metabolism in diverse environments, and Biosynthesis of amino acids, which were extensively annotated in all genomes (Supplementary Fig. 8 and Supplementary Table 5). In addition, we predicted 229 carbohydrate-active enzymes (CAZymes) in all genomes to explore the ability to metabolize carbohydrates (Fig. 4a and Supplementary Table 5). The result showed that the CAZymes of bacteria are mainly glycoside hydrolases (GHs), whereas archaea are glycosyl transferases (GTs) (Fig. 4a). This indicates that bacteria obtained energy for growth by metabolizing carbohydrates, while archaea obtain energy through other pathways and synthesize sugars substance.

Antibiotics are widely used for the treatment and prevention of infection during the suckling period. Microbes may accumulate antibiotic resistance genes (ARGs) under the exposure of antibiotics, and can eventually evolve into drug resistance. We have predicted a total of 2,113 ARGs that can be classified into 38 drug classes from metagenomics, isolate genomes, and MAGs (Fig. 4b). Comparison of the ARGs predicted by the three approaches showed that metagenomics can detect the largest number and types of ARGs, followed by culture-dependent, while a large number of ARGs will be lost in the process of reconstructed genomes from metagenomic bins (Fig. 4b). For analysis of virulence factors (VFs), we used the predicted genes for blast to VFDB⁴⁵. A total of 1,866 VFs were predicted in 748 genomes (Supplementary Fig. 9). Similarly, a large amount of VFs cannot be detected in MAG. We note that Escherichia flexneri not only had most various ARGs, but also contained a large number of VFs.

Discovery of novel SMBGs in the gut microbiome of weanling piglets

Microbes produce a series of secondary metabolites that are not necessary for life activities but have biological activity, which usually mediate important interactions between microbe and microbe-hosts. We explored secondary metabolite biosynthetic gene clusters (SMBGs) in the 748 genomes, and identified 445 SMBGs that could be classified into 19 types from 292 genomes (Supplementary Table 6). Most of these SMBGs are responsible for the synthesis of bacteriocin, followed by sactipetide, arylpolyene, and lantipeptide (Fig. 4c). These SMBGs are predicted from the genomes from 6 phyla, among which Firmicutes contained the largest number and widest variety of SMBGs (Fig. 4c).

For clustering the 445 SMBGs, we generated a total of 231 families, which were distributed in 5 classes of ribosomally synthesized and post-translationally modified peptides (RiPPs), terpene, non-ribosomal peptide synthetases (NRPs), polyketide synthases (PKS) and Others. RiPPs make up the largest class, but only 4 families of RiPPs that are included the reference with known functions in the MIBiG database. Salivaricin CRL1328 is a bacteriocin produced by Ligilactobacillus salivarius^{15, 16}, which has activity against pathogenic bacteria such as Enterococcus faecalis and Enterococcus faecium. We have predicted SMBGs related to the biosynthesis of Salivicin CRL1328 in 22 isolated L. salivarius genomes (Fig. 4d). We speculating that these strains have the potential to inhibit infections. Gassericin, derived from Lactobacillus gasseri, is another important bacteriocin. Studies have shown that Gassericin A can confer diarrhea resistance in pigs¹⁷. We discovered SMBGs that synthesize Gassericin E and Gassericin T in isolated genomes of L. johnsonii (Fig. 4d).

The pangenome analysis of representative species

To extend the phylogenetic analysis of representative species in the gut of weanling piglets at a genome-wide level, we conducted a pangenome analysis of 8 species with the genome number more than 10, including Escherichia flexneri (n=16), Enterococcus faecalis (n=22), Enterococcus faecium (n=21), which represent opportunistic pathogens, and Lactobacillus amylovorus (n=16), Ligilactobacillus salivarius (n=43), Limosilactobacillus mucosae (n=13), Limosilactobacillus reuteri (n=21), Lactobacillus johnsonii (n=17), which represent probiotics. As expected, the resulting accumulation curves showed that the gene repertoires of pan-genome of all the representatives were increased on addition of a new genome and the core genome decreased in contrast (Fig 5a and Supplementary Fig. 10). However, the number of gene families did not show a rapid increase in the pan genome. E. flexneri contained the largest pan genome size of 8,143 genes, 3,111 of which formed the core genome (Supplementary Fig. 11 and Supplementary Table 7). Over half of core genes present in the genomes of E. faecalis and E. faecium, indicating that gene loss and acquisition happened less frequently in these two species. The size of the species-specific pan genomes of the 5 lactobacilli varied from 3,194 to 4,187 genes, respectively (Supplementary Fig. 11 and Supplementary Table 7).

We next analyzed the function details of the core genomes and pan genomes, including the Clusters of Orthologous Groups (COGs) and Kyoto Encyclopedia of Genes and Genomes (KEGG) classification, biosynthesis of bacteriocin, and ARGs. It is obvious that different COG functional classes were enriched in both the core and pan genomes between the opportunistic pathogens and probiotics. COG0438 (Glycosyltransferase involved in cell wall biosynthesis), COG4690 (Dipeptidase), and COG1307 (Fatty acid-binding protein DegV) were enriched in the core genomes compared to their respective accessory and unique genomes of the 5 lactobacilli (Supplementary Table 8), indicating that these COG categorizations represented housekeeping functions involved in implementing the basic growth and metabolism of these probiotics. However, in contrast to lactobacilli, E. flexneri, E. faecalis, and E. faecium possessed more abundant COG1609 (DNA-binding transcriptional regulator, LacI/PurR family) and COG1132 (ABC-type multidrug transport system, ATPase and permease component) in their core genomes (Supplementary Table 8).

For determining the ARGs distribution in the core and pan genome of the 8 species, we discovered tetracycline and macrolide were distributed in the core genome of E. flexneri, E. faecalis, and E. faecium, respectively (Fig. 5b), indicating a wide prevalence of tetracycline and macrolide resistance among these species. Their presence may have resulted from the overuse of antibiotics in farmed pigs for disease prevention in past years. Nevertheless, rare ARGs present both in the core and pan genome of other lactobacilli (Fig. 5b and Supplementary Table 9), suggests the safety of these species for using as probiotics in feeding of weanling piglet.

Reuterin, produced by some strains of L. reuteri, has antimicrobial properties¹⁸. We investigated the prevalence of biosynthesis genes (pdu-cbi-cob-hem) related to reuterin in the core and dispensable genomes and detected the production of reuterin experimentally. In total 41 genes related to the synthesis of reuterin were predicted in the pan genome of L. reuteri, including cbiA~cbiT (15 genes), cobB~cobU (10 genes), hemA~hemL (7 genes), and pduA~pduV (9 genes) (Fig. 5c and Supplementary Table 10). Notably, only cobB, cobC_1, and cobC_2 were found in the core genomes for all the 21 genomes. The rest were present in the accessory or unique genomes. 9 of 21 genomes encoded completed biosynthesis gene clusters (pdu-cbi-cob-hem), indicating these strains were most likely to be reuterin producers. We next detected the production of Reuterin in the fermentation supernatant of the culture of 19 isolated L. reuteri. As expected, 4 strains with completed biosynthesis gene clusters (pdu-cbi-cob-hem) yield reuterin ranged from 0.6 to 2.54 nm/L (Supplementary Table 10).

Previous culture-independent methods have uncovered that the gut microbiota is associated with the survival rate of piglet after weaning¹⁹. In this study, we investigated the gut microbiota of piglets combining with culture-based and metagenomic strategies. A collection including 1,476 cultured bacteria, 748 reference genomes catalog, comprising 266 cultured genomes and 482 MAGs, which represent the gut microbiome of piglet was constructed. 70.80% of these were represent novel species according to the present database. We also found 5 MAGs were archaea. We contributed 3,630,043 novel genes to the present pig gene catalog and expanded to 10.91 Mb. In our collection, over half of the taxa were lactobacilli, that will be a valuable resource for probiotics development and diversity study for different niche of these species.

We discovered 84.20% of the clusters from MAGs could not be assigned to any known species, only 3 clusters matched 11 cultured genomes with novelty. Prevotella, a dominant genus in the gut of piglet, can be detected using metagenomic analysis and a total of 26 genomes from the metagenome bins were annotated as Prevotella, but none were included in our culture collection constructed in this study. This emphasizes the necessity of continuing effort for developing more culture methods that are appropriate for more bacteria.

Our genome-based functional analysis provides a better understanding of the gut microbiota of weanling piglets. Genes related to the metabolic pathways, biosynthesis of secondary metabolites, and Biosynthesis of antibiotics were highly enriched in all genomes. We found that the gut microbiota of weaning piglets contained a large number of glycoside hydrolase enzymes, which can be used to metabolize major dietary carbohydrates such as corn, soybean meal, wheat, and rice bran to provide most of the energy and nutrients needed for daily life²⁰. The excavation of SMBGs found antibiotics that inhibit pathogenic bacteria such as E. faecalis and E. faecium, emphasized the effect of gut microbiomes in suppressing infections in piglets. Finally, we annotated the ARGs and VFs from metagenomes and genomes. It is worth noting that the metagenomes contained abundant ARGs and VFs, which were distributed in specific isolated genomes, whereas the MAGs had a very small number of annotated ARGs and VFs. We highlighted that both techniques are complement for the detection of harmful genes.

Pan genome analyses provide comprehensive genetic landscapes of the representative species and identify the characteristics of core and strain-specific genes. Lactobacillus, which represent the probiotics in the pig gut, implement functions with COG0438 (Glycosyltransferase involved in cell wall biosynthesis), COG4690 (Dipeptidase), and COG1307 (Fatty acid-binding protein DegV). Bacteriocin, like reuterin, can be used as a broad-spectrum antimicrobial agent to prevent piglet diarrhea¹⁷. We discovered the completeness of biosynthesis pathway for reuterin is strain-specific, the reuterin-related genes were present in the dispensable genome, indicating when exploitation and application of bacteriocin, the strain level diversity should be considered. In this collection, we identified 4/21 L. reuteri stains that had the capacity to produce reuterin, indicating those strains have potential in prevention of pathogen infections and colonization of the piglet gastrointestinal tract.

Metagenomic binning has been widely used to recover genomes from the fecal samples. But the cultured bacteria are of great importance for the experimental validation of their functions, especially for the development of probiotics. As a result of this study, we have provided both the cultured bacteria and reference genome data, that are publicly available at CNGB. It will be a useful resource for the future studies of microbiota–host interactions and the development as probiotics for disease prevention of weanling piglets.

In this study, we constructed a collection including 1,476 cultured bacteria, 748 reference genomes that comprising of 266 cultured genomes and 482 MAGs, which represent the gut microbiome of piglet. These genomes represent 428 species belonging to10 phyla. This collection expands the bacterial and genomic resources of gut microbiota of weanling piglets by adding 333 new genomic species. Functional exploration revealed these bacteria harbored large amount of glycoside hydrolase enzymes and secondary metabolite biosynthetic genes. More abundant ARGs and VFs were annotated from metagenomes and specific isolated genomes. Pan genome analysis of Limosilactobacillus reuteri showed that the biosynthesis of reuterin was strain-specific and the production was experimentally determined. We believe that our collection provided more valuable resource for the further study of the gut microbiota of weanling piglets and enable explore the probiotics for prevention of disease.

Sample collection and bacteria culturing

42 healthy piglets (Duroc♂×(Landrace×Yorkshire)♀) with good mental state were selected for sampling, of which 12 piglets of lactation were sampled on 18 days for ileum and colon content, respectively, 21 weaned piglets were sampled on 35 days for ileum and colon content, and 9 weaned piglets were sampled on 60 days for faeces. During sampling, the intestinal contents were obtained by dissecting the colon and ileum with sterile instruments. Every 3 samples from one line with same time point and region were pooled together. All samples were stored at -80 ℃ and prepared for bacterial culturing and metagenome sequencing.

For bacteria culture, samples were processed according to our previous bacteria culture-based study²³. To acquire more different bacterial taxa, we applied a batch of rich and selective culture medium which are provided in the Supplementary method. The purified bacteria were maintained with 20% glycerine at -80 ℃. Detailed information about the origin of all strains, including the age of the donors and gut region are presented in Supplementary Data 1.

The acquired bacteria were subject to 16S rRNA gene sequencing. The amplification of 16S rRNA gene was using 27f and 1492r primers and the obtained sequences were trimmed as previously described²⁴. The taxonomy of bacteria was recognized by blast the 16S rRNA gene sequences against reference sequences within the NCBI 16S ribosomal RNA sequence database. We used the recommended threshold of 98.7% and 94.5% as the species and genus boundaries²⁵. 16S rRNA gene sequences were aligned by MAFFT v7.310²⁶ and trimmed by trimAl v1.4. rev 22²⁷ with auto option. Phylogenetic tree based on 16S rRNA gene sequences was reconstructed by using the maximum-likelihood method with FastTree Version 2.1.3 SSE3²⁸.

Genome sequencing and assembly

Genomic DNA of the 226 representative strains was extracted using the method as previously described²³. The draft genomes were sequenced for 100 bp paired-end on the DNBSEQ-T7 platform. The raw reads were filtered using SOAPnuke26²⁹ (v1.5.0; -l 20 -q 0.4 -n 0.1 -d -M 3 --seqType 0 -Q 2 -c 2.66666666666667). After preliminary assessment of the genome size based on kmer, SPAdes (v3.11.1; -T 4 -m 100) was used to assembly.

Metagenome sequencing, assembly, and binning

Metagenomic DNA from 14 samples was extracted and sequenced according to the method described above. Raw reads of each sample were preprocessed by SOAPnuke (v1.5.2)²⁹ using the ‘filter’ module (option -l 20 -q 0.2 -n 0.05 -Q 2 -d -c 0 -5 0 -7 1), and host reads were removed by SOAPaligner³⁰ (v2.22, option -m 4 -s 30 -r 1 -n 50 -x 1000 -v 4). Kraken2³¹ was used to assign taxonomic labels to the metagenome reads and calculate bacterial abundance and the profile was provided in Supplementary Table 1. Shannon alpha diversity indexes (‘diversity’ function) and Bray-Curtis dissimilarity indexes (‘vegdist’ function) for the principal coordinate analysis (PCoA) were identified in R software with species abundances.

Thereafter, IDBA-UD (v 1.1.3)³² was used to assemble and merge the optimal contig. MAGs were generated using three different tools configured in the 'binning' module of metaWRAP (v1.1.5)³³ (integrated with MaxBin2, MetaBAT2, and CONCOCT) and were finally consolidated and optimized using 'Bin_refinement'. dRep³⁴ (v2.5.4, option -p 8 -comp 70 -con 5) was used to de-redundancy and preliminary quality control of the MAGs. The whole-genome assembly was as described by Zou et al²³. The quality of the MAGs and isolated genomes was estimated by CheckM (v1.0.12)³⁵. The classification criteria for high-quality genomes are based on >90% completeness and <5% contamination, while genomes with <70% completeness or >5% contamination were not included in the following analysis.

Species-level clustering, phylogenetic and taxonomic analyses

FastANI (v1.32)³⁶ was used to calculate the pairwise ANI, and genomes were clustered into species-level clusters with a 95% ANI cut-off using the R package hclust. dRep was used to select the best genome from each cluster as the representative. All genomes were taxonomically annotated using GTDB-Tk³⁷ (v1.3.0, database release95³⁸). Any lineage without an annotated species or invalid names with reserved suffixes was considered to represent a potential new species. The maximum-likelihood phylogenetic tree of the representative genomes was constructed with GTDB-Tk (‘infer’ followed by ‘classify_wf’). The phylogenetic tree was viewed using the online tool EVOLVIEW v2³⁹.

Alignment with other genome collections

We downloaded the genomes from two studies, which represent the reference dataset for the MAGs and the isolated genome of the pig gut microbiome. Wang et al.⁹ generated a total of 360 substantially complete MAGs (>70% completeness, <5% contamination) as the first metagenomic reference for swine intestinal microbiota. PiBAC¹² comprised 117 high-quality isolated genomes representing 110 species. Both datasets were downloaded from NCBI (BioProject: PRJNA494875 and PRJNA561470, respectively). We performed a pairwise ANI calculation for all genomes, and only those ANI>95% were considered to be matched.

Gene prediction and function annotation

GeneMark.hmm PROKARYOTIC (v3.38)⁴⁰ is a gene prediction program for metagenomes. Protein-coding sequences (CDS) for each genome were predicted and annotated by Prokka v1.14.6⁴¹. The metagenome non-redundancy gene catalog was generated by CD-HIT v4.6.3⁴² (cd-hit-est with option -c 0.95 defining protein identity of 95%). Based on 95% protein similarity, we clustered the gene catalog with Xiao⁵ and calculated the gene overlap of the metagenome and genome.

Functional characterization of all CDS was performed by BLAST analysis with the KEGG database (release 81)⁴³, CAZyDB⁴⁴ (downloaded from http://www.cazy.org/, as of April 2016), and VFDB (setB, 2021-07)⁴⁵. ARGs were identified with rgi (5.2.0) using reference data from the Comprehensive Antibiotic Resistance Database (CARD, v 3.1.2)⁴⁶.

The SMBGs were identified by using anti-SMASH (v4.2.0)⁴⁷ and were clustered by using MIBiG (version 1.4)⁴⁸ with reference BGCs by BiG-SCAPE⁴⁹ with default parameters. The similarity networks of the SMBGs were displayed by using Cytoscape (v3.8.2)⁵⁰.

Pan genome analysis of the representative species

The Bacterial Pan Genome Analysis tool (BPGA) pipeline⁵¹ was used for protein clustering (usearch, 80% identity value), identification of core, accessory, and unique genes of each cluster, and generation of pan-genome accumulation curves for the 171 genomes in eight clusters. Clusters of Orthologous Groups (COGs) of proteins based on Prokka annotation described above were collected to identify the core COGs specific to genome clusters with an in-house script. Total ARGs of these clusters were summarized to determine the presence and absence of them in each genome and to identify the core and dispensable genomes of antibiotic resistance with an in-house script.

Detecting the production of reuterin

Reuterin production from L. reuteri cultures was determined according to the previous described method¹⁸. In brief, cell pellets were harvested by centrifugation from overnight grown cultures of L. reuteri in MRS broth. The pellets were washed twice with 0.1 M potassium phosphate buffer, resuspended in 70% glycerol, and then incubated at 25°C for 2.5 h under anaerobic conditions. Quantification of reuterin in the suspension was determined by the colorimetric method as described¹⁸.

MAG: Metagenome-assembled genome

MPYG: Peptone Yeast Extract Glucose Medium, Modified

CNGB: China National Genebank

GTDB-Tk: Genome Taxonomy Database Toolkit

PiBAC: pig intestinal bacterial collection

CAZyme: carbohydrate-active enzyme

GH: glycoside hydrolase

GT: glycosyl transferase

ARG: antibiotic resistance gene

VF: virulence factor

SMBG: secondary metabolite biosynthetic gene cluster

RiPP: ribosomally synthesized and post-translationally modified peptide

NRP: non-ribosomal peptide synthetase

PKS: polyketide synthases

KEGG: Kyoto Encyclopedia of Genes and Genomes

COG: Clusters of Orthologous Group

Author Contributions

Conceived and designed the experiments: Y.Z. and W.W. Performed the experiments: B.D, Y.Z., X.L, and X.J. Analyzed the data: Y.Z., X.L, T.H., J.Z., J.C, and G.L. Contributed reagents/materials/analysis tools: B.D., Y.Z., X.J., Z.C., J.L., and L.X. Supervised the work: L.X., B.W., Y.H., S.L., and J.L. Wrote the paper: B.D., Y.Z., X.L, and T.H. Revised the paper: B.D., Y.Z., and W.W. All authors commented on the manuscript.

Acknowledgements

We appreciate all efforts made by all authors in GDAAS and BGI research.

Funding

This work was supported by grants from National Natural Science Foundation of China (No. 32002287), Guangzhou Municipal Science and Technology Bureau (No. 201804010338), Science and Technology Program of Guangdong Province (2014B020201002), the Presidential Foundation of Guangdong Academy of Agricultural Sciences (grant number 201918), Science and Technology Planning Project of Guangdong Province (2016A030303034), Grant from Guangzhou Municipal Science and Technology Bureau (No. 202102080151), the Key Laboratory of Modern Biological Seed Industry in South China, Ministry of Agriculture and Rural Affairs (2105-000000-20-03-457451), and Guangdong Provincial Key Laboratory of Genome Read and Write (No. 2017B030301011). We also thank the colleagues at BGI-Shenzhen for sample collection, and discussions, and China National Genebank (CNGB) Shenzhen for DNA extraction, library construction, and sequencing.

Ethics approval and consent to participate

Pigs were managed according to the GDAAS Animal Care and Use Committee approved protocol #1008.

Consent for publication

Not applicable

Disclosure statement

The authors declare no competing financial interests.

Data availability statement

The data that support the findings of this study have been deposited into CNGB Sequence Archive (CNSA)²¹ of China National GeneBank DataBase (CNGBdb)²² with accession number CNP0002075 and CNP0002072 for metagenomic and bacterial genomic data, respectively. All the bacterial strains in our collection have been deposited in China National GeneBank (CNGB), a non-profit, public-service-oriented organization in China.

Competing interests

The authors declare that they have no competing interests.

Author details

¹Guangdong Provincial Key Laboratory for Crop Germplasm Resources Preservation and Utilization, Agro-biological Gene Research Center, Guangdong Academy of Agricultural Sciences, Guangzhou, 510640, China

²Qingdao-Europe Advanced Institute for Life Sciences, BGI-Shenzhen, Qingdao 266555, China

³China National Genebank, BGI-Shenzhen, Shenzhen 518120, China

⁴ BGI-Shenzhen, Shenzhen 518083, China

⁵School of Bioscience and Biotechnology, South China University of Technology, Guangzhou 510006, China

⁶BGI-Qingdao, BGI-Shenzhen, Qingdao 266555, China

⁷Laboratory of Genomics and Molecular Biomedicine, Department of Biology, University of Copenhagen, Universitetsparken 13, 2100 Copenhagen, Denmark

⁸ Shenzhen Engineering Laboratory of Detection and Intervention of human intestinal microbiome, BGI-Shenzhen, Shenzhen, China

⁹ Guangdong Provincial Key Laboratory of Genome Read and Write, BGI-Shenzhen, Shenzhen, 518120, China

¹⁰Guangdong Laboratory for Lingnan Modern Agriculture, Guangzhou, Guangdong 510642, China

Maes DGD, Dewulf J, Pineiro C, Edwards S, Kyriazakis I. A critical reflection on intensive pork production with an emphasis on animal health and welfare. J Anim Sci 98, S15-S26 (2020).
Chen M, Hu E, Kuen LL, Wu L. Study on Consumer Preference for Traceable Pork With Animal Welfare Attribute. Front Psychol 12, 675554 (2021).
McCormack UM, et al. Porcine Feed Efficiency-Associated Intestinal Microbiota and Physiological Traits: Finding Consistent Cross-Locational Biomarkers for Residual Feed Intake. mSystems 4, e00324-18 (2019).
Torres-Pitarch A, et al. Effect of cereal fermentation and carbohydrase supplementation on growth, nutrient digestibility and intestinal microbiota in liquid-fed grow-finishing pigs. Sci Rep 10, 13716 (2020).
Xiao L, et al. A reference gene catalogue of the pig gut microbiome. Nat Microbiol 1, 16161 (2016).
Chen C, et al. Expanded catalog of microbial genes and metagenome-assembled genomes from the pig gut microbiome. Nat Commun 12, 1106 (2021).
Teeling H, Glockner FO. Current opportunities and challenges in microbial metagenome analysis--a bioinformatic perspective. Brief Bioinform 13, 728–742 (2012).
Howe A, Chain PS. Challenges and opportunities in understanding microbial communities with metagenome assembly (accompanied by IPython Notebook tutorial). Front Microbiol 6, 678 (2015).
Wang W, Hu H, Zijlstra RT, Zheng J, Ganzle MG. Metagenomic reconstructions of gut microbial metabolism in weanling pigs. Microbiome 7, 48 (2019).
Lagier JC, et al. Culturing the human microbiota and culturomics. Nat Rev Microbiol 16, 540–550 (2018).
Lagier JC, et al. Culture of previously uncultured members of the human gut microbiota by culturomics. Nat Microbiol 1, 16203 (2016).
Wylensek D, et al. A collection of bacterial isolates from the pig intestine reveals functional and taxonomic diversity. Nat Commun 11, 6389 (2020).
Kim HB, Isaacson RE. The pig gut microbial diversity: Understanding the pig gut microbial ecology through the next generation high throughput sequencing. Vet Microbiol 177, 242–251 (2015).
Li K, Xiao Y, Chen J, Chen J, He X, Yang H. Microbial composition in different gut locations of weaning piglets receiving antibiotics. Asian-Australas J Anim Sci 30, 78–84 (2017).
Ocana VS, Pesce De Ruiz Holgado AA, Nader-Macias ME. Characterization of a bacteriocin-like substance produced by a vaginal Lactobacillus salivarius strain. Appl Environ Microbiol 65, 5631–5635 (1999).
Vera Pingitore E, Hebert EM, Nader-Macias ME, Sesma F. Characterization of salivaricin CRL 1328, a two-peptide bacteriocin produced by Lactobacillus salivarius CRL 1328 isolated from the human vagina. Res Microbiol 160, 401–408 (2009).
Hu J, et al. A Microbiota-Derived Bacteriocin Targets the Host to Confer Diarrhea Resistance in Early-Weaned Piglets. Cell Host Microbe 24, 817-832 e818 (2018).
Zhang Z, et al. A Phylogenetic View on the Role of Glycerol for Growth Enhancement and Reuterin Formation in Limosilactobacillus reuteri. Front Microbiol 11, 601422 (2020).
Gresse R, Chaucheyras-Durand F, Fleury MA, Van de Wiele T, Forano E, Blanquet-Diot S. Gut Microbiota Dysbiosis in Postweaning Piglets: Understanding the Keys to Health. Trends Microbiol 25, 851–873 (2017).
Lindberg JE. Fiber effects in nutrition and gut health in pigs. J Anim Sci Biotechnol 5, 15 (2014).
Guo X, et al. CNSA: a data repository for archiving omics data. Database (Oxford) 2020, (2020).
Chen FZ, et al. CNGBdb: China National GeneBank DataBase. Yi Chuan 42, 799–809 (2020).
Zou Y, et al. 1,520 reference genomes from cultivated human gut bacteria enable functional microbiome analyses. Nat Biotechnol 37, 179–185 (2019).
Zou Y, et al. Lactobacillus shenzhenensis sp. nov., isolated from a fermented dairy beverage. Int J Syst Evol Microbiol 63, 1817–1823 (2013).
Yarza P, et al. Uniting the classification of cultured and uncultured bacteria and archaea using 16S rRNA gene sequences. Nat Rev Microbiol 12, 635–645 (2014).
Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol 30, 772–780 (2013).
Capella-Gutierrez S, Silla-Martinez JM, Gabaldon T. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics 25, 1972–1973 (2009).
Price MN, Dehal PS, Arkin AP. FastTree 2--approximately maximum-likelihood trees for large alignments. PLoS One 5, e9490 (2010).
Chen Y, et al. SOAPnuke: a MapReduce acceleration-supported software for integrated quality control and preprocessing of high-throughput sequencing data. Gigascience 7, 1–6 (2018).
Gu S, Fang L, Xu X. Using SOAPaligner for Short Reads Alignment. Curr Protoc Bioinformatics 44, 11 11 11–17 (2013).
Wood DE, Lu J, Langmead B. Improved metagenomic analysis with Kraken 2. Genome Biol 20, 257 (2019).
Peng Y, Leung HC, Yiu SM, Chin FY. IDBA-UD: a de novo assembler for single-cell and metagenomic sequencing data with highly uneven depth. Bioinformatics 28, 1420–1428 (2012).
Uritskiy GV, DiRuggiero J, Taylor J. MetaWRAP-a flexible pipeline for genome-resolved metagenomic data analysis. Microbiome 6, 158 (2018).
Olm MR, Brown CT, Brooks B, Banfield JF. dRep: a tool for fast and accurate genomic comparisons that enables improved genome recovery from metagenomes through de-replication. ISME J 11, 2864–2868 (2017).
Parks DH, Imelfort M, Skennerton CT, Hugenholtz P, Tyson GW. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res 25, 1043–1055 (2015).
Jain C, Rodriguez RL, Phillippy AM, Konstantinidis KT, Aluru S. High throughput ANI analysis of 90K prokaryotic genomes reveals clear species boundaries. Nat Commun 9, 5114 (2018).
Chaumeil PA, Mussig AJ, Hugenholtz P, Parks DH. GTDB-Tk: a toolkit to classify genomes with the Genome Taxonomy Database. Bioinformatics, (2019).
Parks DH, Chuvochina M, Chaumeil PA, Rinke C, Mussig AJ, Hugenholtz P. A complete domain-to-species taxonomy for Bacteria and Archaea. Nat Biotechnol 38, 1079–1086 (2020).
He Z, Zhang H, Gao S, Lercher MJ, Chen WH, Hu S. Evolview v2: an online visualization and management tool for customized and annotated phylogenetic trees. Nucleic Acids Res 44, W236-241 (2016).
Besemer J, Lomsadze A, Borodovsky M. GeneMarkS: a self-training method for prediction of gene starts in microbial genomes. Implications for finding sequence motifs in regulatory regions. Nucleic Acids Res 29, 2607–2618 (2001).
Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinformatics 30, 2068–2069 (2014).
Fu L, Niu B, Zhu Z, Wu S, Li W. CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinformatics 28, 3150–3152 (2012).
Kanehisa M, Furumichi M, Tanabe M, Sato Y, Morishima K. KEGG: new perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res 45, D353-D361 (2017).
Lombard V, Golaconda Ramulu H, Drula E, Coutinho PM, Henrissat B. The carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res 42, D490-495 (2014).
Chen L, Zheng D, Liu B, Yang J, Jin Q. VFDB 2016: hierarchical and refined dataset for big data analysis--10 years on. Nucleic Acids Res 44, D694-697 (2016).
Alcock BP, et al. CARD 2020: antibiotic resistome surveillance with the comprehensive antibiotic resistance database. Nucleic Acids Res 48, D517-D525 (2020).
Blin K, et al. antiSMASH 4.0-improvements in chemistry prediction and gene cluster boundary identification. Nucleic Acids Res 45, W36-W41 (2017).
Crits-Christoph A, Diamond S, Butterfield CN, Thomas BC, Banfield JF. Novel soil bacteria possess diverse genes for secondary metabolite biosynthesis. Nature 558, 440–444 (2018).
Navarro-Munoz JC, et al. A computational framework to explore large-scale biosynthetic diversity. Nat Chem Biol 16, 60–68 (2020).
Shannon P, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res 13, 2498–2504 (2003).
Chaudhari NM, Gupta VK, Dutta C. BPGA- an ultra-fast pan-genome analysis pipeline. Sci Rep 6, 24373 (2016).

Download PDF

Journal Publication

published 16 Feb, 2022

Read the published version in Microbiology Spectrum →

Version 1

posted

You are reading this latest preprint version

A Bacterial Genome And Culture Collection of Gut Microbial in Weanling Piglet

Status:

Journal Publication

Version 1

Abstract

Figures

Introduction

Results

Discussion

Conclusion

Materials And Methods

Abbreviations

Declarations

References

Supplementary Files

Status:

Journal Publication

Version 1