Identifying Mobile Genetic Elements in the Ruminal Microbiome of Nellore Cattle: An Initial Investigation

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

Download PDF

Article

Identifying Mobile Genetic Elements in the Ruminal Microbiome of Nellore Cattle: An Initial Investigation

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

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

Metagenomics has made it feasible to elucidate the intricacies of the ruminal microbiome and its role in the differentiation of animal production phenotypes of significance. The search for mobile genetic elements (MGEs) has taken on great importance, as they play a critical role in the transfer of genetic material between organisms. Furthermore, these elements serve a dual purpose by controlling populations through lytic bacteriophages, thereby maintaining ecological equilibrium and driving the evolutionary progress of host microorganisms. In this study, we aimed to identify the association between ruminal bacteria and their MGEs in Nellore cattle using physical chromosomal links through the Hi-C method. Shotgun metagenomic sequencing and the proximity ligation method ProxiMeta™ were used to analyze DNA, getting 1,713,111,307 bp, which gave rise to 107 metagenome-assembled genomes from rumen samples of four Nellore cows maintained on pasture. Taxonomic analysis revealed that most of the bacterial genomes belonged to the families Lachnospiraceae, Bacteroidaceae, Ruminococcaceae, Saccharofermentanaceae, and Treponemataceae and mostly encoded pathways for central carbon and other carbohydrate metabolisms. A total of 31 associations between host bacteria and MGE were identified, including 17 links to viruses and 14 links to plasmids. Additionally, we found 12 antibiotic resistance genes. To our knowledge, this is the first study in Brazilian cattle that connect MGEs with their microbial hosts. It identifies MGEs present in the rumen of pasture-raised Nellore cattle, offering insights that could advance biotechnology for food digestion and improve ruminant performance in production systems.

Biological sciences/Microbiology/Bacteriophages

Biological sciences/Genetics/Microbial genetics

Advancements in metagenomics in rumen environments have enabled us to identify and classify the taxonomy and functional capabilities of bacteria that specialize in digesting plant material, which is the primary food source for ruminants. Given the impact of microorganisms on host cattle performance, microbiota modulation is crucial for developing a sustainable and productive livestock system¹. For instance, feed efficiency is a trait with implications for economic and environmental sustainability² and is strongly associated with host animal microbiota profiles^3–5. Ruminal methane production, which results from microbial fermentation in the rumen, is another example of trait that can be explained by genetic variation and the composition of the host animal microbiota^6–8. Approximately 65% of rumen prokaryotic species have been classified and described in taurine breeds⁹; however, less information is available on Indian cattle, which are the majority in the world, especially in Brazil.

Understanding the role of specific phylotypes in the rumen during fermentation and identifying new enzymes through metagenomic analysis can improve biotechnology for food digestion and ruminant performance in production systems¹⁰. Underutilized components of the rumen, known as Mobile Genetic Elements (MGEs), possess great potential for manipulation and biotechnological applications owing to their ability to confer benefits to their hosts and participate in crucial functional modifications of their hosts, carrying auxiliary metabolic genes such as glycosidic hydrolases (GH), which contribute to the breakdown of complex carbohydrates.

In addition, viruses carry genes that contribute positively their hosts in the metabolism of nutrients and stages of the biogeochemical cycles carried out by the hosts ^11,12. They are classified as mobile DNA elements that move within and between bacterial cells through agents, such as plasmids, bacteriophages, and transposons, resulting in the introduction of accessory genes that contribute to the evolution of host cells^13,14. Solden et al.¹⁵ reported that viruses play a crucial role in controlling the ecosystem functions in the rumen. These viruses have metabolic interdependence with their microbial hosts, carry genes that can redirect carbon metabolism, and infect dominant microorganisms that degrade carbohydrates.

Early studies on bacteriophages in the rumen revealed that they can eliminate specific bacteria¹⁶. Recent research has indicated the potential use of lytic phages to regulate microorganisms in the rumen responsible for an increase in acidosis and methane production^17,18. This suggests that the components in the rumen are interconnected with viruses that infect ruminal bacteria, the microbiome, and the host animal interacting and affecting each other. This can change the ruminal function and the nutritional capacity of the host animal¹⁹.

Recent studies have employed the Hi-C approach, which captures intimate interactions between bacterial chromosomes within living cells, to investigate MGEs in microbial communities. This method has been used in conjunction with metagenomic sequencing and has proven effective in assembling genomes and uncovering new relationships between hosts and MGEs^20–23. The objective of this study was to determine the existence of MGEs in the rumen content of Nellore cattle in Brazil, using the Hi-C method to establish a link between these elements and their microbial hosts.

Metagenome of bovine rumen content samples

In this study, the metagenomic sequencing was performed with DNA extracted from a pool of rumen content samples of four 4 Nellore cows, applying ProxiMeta™, a proximity linkage method based on Hi-C. This method used contigs derived from shotgun sequencing and reads from a library Hi-C, and generated an assembly length of 1,713,111,307 bp, with 834,164 contigs total and 16,035 contigs in clusters allowing the assembly of 107 genomes (Fig. 1A). The most completely assembled genome was bin_4 (99.30%), which corresponded to Treponema D. sp016288205. Other genomes classified as complete, with completeness greater than 95%, were UBA2868 sp003535955 (bin_1), Treponema D. sp017421145 (bin_9), Treponema D. sp902783155 (bin_6) and UBA2868 (bin _12). The clusters classified as excellent (> 90%) were dominated by species from the Lachnospiraceae family, namely UBA3766 sp902803015 (bin_22), UBA1712 (bin_14), UBA2868 (bin_3), UBA2868 sp003535955 (bin_17), NK4A144 sp902783395 (bin_7), Acetatifactor sp900066565 (bin_21), UBA1711 sp900317125 (bin_11), and UBA1711 sp001543385 (bin_18). Clusters belonging to the families Saccharofermentanaceae with Saccharofermentans (bin_19), and Bacteroidaceae with Prevotella sp000702825 (bin_5) also had an excellent completeness. The clusters considered good are presented in the Supplementary Table 1, and the 107 clusters subjected to taxonomic classification are represented in Fig. 1B.

Functional and carbohydrate-active enzymes (CAZYmes) characterization of bins

Before analyzing the MGEs, we conducted a functional analysis to understand the functional characteristics encoded within metagenome-assembled genomes (MAGs) from the ruminal microbiome of Nellore cattle. For this, MAGs with completeness > 50% and those associated with MGEs were used, totaling 43 genomes (Supplementary Tables 8 and 9). This step provided essential insights into the interplay between microbial communities and the functional potential associated with MGEs in the rumen ecosystem. Our study evaluated the involvement of rumen microorganisms in these biogeochemical pathways, and we found that the majority were associated with the stages of the carbon, sulfur, nitrogen and iron cycles (Fig. 2). Regarding the carbon cycle, most of the genomes identified demonstrated involvement in the acetogenesis process and carbon degradation complex. More specifically, we identified two genomes belonging to Lachnospiraceae and Bacteroidales encoding the enzyme methane monooxygenase (mmoBD), which act in methane oxidation, and members of these clades also participated in the oxidation step of the sulfur cycle. In the nitrogen cycle, phylotypes belonging to Saccharofermentans and Treponema acted in the nitrogen fixation stage and members of the genus Provotella encoded genes for the reduction of nitrite into ammonia. All taxonomic classes presented members active in the iron oxidation stage, which is fundamental for most microorganisms, due to its redox properties, being used as an electron transporter by bacteria ^24,25.

The most frequent processes included steps in the central carbon metabolism pathway, the metabolism of other carbohydrates, and purine metabolism. Among the bins examined, 69% contained genes that encoded the phosphate acetyltransferase-acetate kinase pathway, a reversible process where acetyl-CoA can be converted into acetate or acetate into acetyl-CoA. Additionally, the aromatic amino acid metabolism pathway with tryptophan biosynthesis was present in 64% of the bins (Fig. 3A).

The metagenomic approach allows for the evaluation of the potential of rumen microbial populations to break down plant lignocellulolytic materials, by annotating the genes associated with this process, known as CAZymes¹⁰. CAZymes are classified into several categories based on their degradation capabilities, including glycoside hydrolases (GHs), carbohydrate esterases (CEs), polysaccharide lyases (PLs), and others²⁶. In this study, the bins identified were assessed for their ability to encode CAZymes, and the presence of GHs and PLs was confirmed (Supplementary Table 3). We found that GH2, GH3, GH5, GH13, and GH43 were the most abundant enzymes mapped in the bins (Fig. 3B).

Assignment of mobile genetic elements to their hosts

Clustering with Hi-C allows not only the construction of MAGs but also the sequence grouping of mobile genetic elements with their core genome²⁷. The sequences were obtained with sequencing and the use of the proximity linkage method to assign viruses and plasmids to their microbial hosts, the data obtained is summarized in Fig. 4.

Plasmids and the Identification of microbial resistance genes

The reads derived from shotgun sequencing were used to evaluate the presence of plasmid sequences. Thirty-three DNA sequences were identified as plasmids, which were not integrated into the host genome. Among these, 14 had at least one connection with the prokaryotic host genome through Hi-C reads. Furthermore, two putative plasmids contigs k141_1760976 and k141_5289202 were shared between multiple host clusters (bin_3, bin_5, bin_20, bin_80) and (bin_1; bin_3; bin_11; bin_30; bin_88) (Supplementary Table 4). The bins shared by contig k141_1760976 belong to different genus, including UBA2868, Prevotella and Saccharofermentans, which characterizes horizontal gene transfer mediated by plasmids. The plasmid k141_5289202 demonstrated greater specificity, being shared by three individuals belonging to Lachnospiracea and one belonging to Bacteroidales.

Additionally, we identified 19 plasmid sequences integrated into the host genome sequences (Fig. 4A). Within this set, it was observed that 50% of the total antibiotic resistance genes (ARGs) validated in our database were present (Fig. 5). All the identified resistance genes had an identity of more than 90% with the reference sequence (Supplementary Table 5). The tet32 gene sequence was associated with bin_14, which belongs to the Lachnospiraceae family.

Virus found in the rumen

A total of 11 high-quality metagenome-assembled viral genomes (vMAGs) were identified and characterized as viral sequences (as shown in Fig. 4B). They had completeness rates of 92–100% and a length of 35 to 126 kbp. Another 22 MAGs were obtained with medium quality, of which 13 had completeness rates of 71–89% and a length of 25 to 90 kbp, and nine had completeness rates of 68–50% and lengths of 20 and 130 kbp. Any genomes with less than 50% completeness were considered low quality. All of the assembled genomes were determined to be free of contamination (Supplementary Table 6).

The majority of vMAGs had their taxonomy associated with genomes that have not yet been assigned to the ICTV phage database²⁸. There were vMAGs associated with the same viral cluster (VC); however, these do not yet have taxonomic assignment in the database: vMAG_4 and vMAG_35 were associated with the viral cluster VC_146, as well as vMAG_5 and vMAGs_18, with VC_147 and vMAG_25, and vMAG_56, with VC_200 (Supplementary Table 7).

Among the viruses that could be assigned to known genomes, vMAG_1 and vMAG_48 were assigned to the Myoviridae family VC_155, vMAG_40, vMAG_46 and vMAG_57 were assigned to the Siphoviridae family, VC_102, vMAG_7, vMAG_16, 29 and vMAG_52 were also assigned to the Siphoviridae family, but with genomes associated with Streptomyces phage, VC_87 (Fig. 6).

Virus-host association

The viral clusters associated with the bacterial genomes determined by Hi-C mostly exhibited a lytic cycle. The vMAG_54, associated with the host UBA2868 sp003535955 (Lachnospiraceae), presented high quality, achieved 100% completeness in its genome, and was characterized by a lytic lifestyle (Table 1).

Two co-infection events were identified, in which bacterial cells presented the assignment of more than one bacteriophage simultaneously, that is, different bacteriophages infecting the same bacterial host. Bacteriophages corresponding to vMAGs 18 and 45 infected the same bacterial cell corresponding to bin_28, and the same occurred with vMAGs 12 and 54, which infected the same bacterial cell corresponding to bin_30 (Fig. 7).

Table 1

Viruses associated with bacterial hosts by Hi-C
Name cluster host	Host	Name Cluster (phage)	Life Cycle	CheckV quality	Completeness (phages)
bin_30	UBA2868 sp003535955 (Lachnospiraceae)	vMAG_54	Lytic	High-quality	100.0
bin_41	g_UBA2868 (Lachnospiraceae)	vMAG_13	-	Medium-quality	87.8
bin_12	UBA2868 sp002368575 (Lachnospiraceae)	vMAG_4	Lytic	Medium-quality	85.14
bin_28	Saccharofermentans sp003543635	vMAG_45	Lytic	Medium-quality	72.67
bin_22	UBA3766 sp902803015 (Lachnospiraceae)	vMAG_35	Lytic	Medium-quality	60.16
bin_30	UBA2868 sp003535955 (Lachnospiraceae)	vMAG_12	Lytic	Medium-quality	54.18
bin_28	Saccharofermentans sp003543635	vMAG_18	Lytic	Medium-quality	50.67
bin_16	RUG12461 sp016286115 (Lachnospiraceae)	vMAG_22	Lytic	Low-quality	46.05
bin_102	Non-classified	vMAG_43	Lytic	Low-quality	39.81
bin_27	Butyrivibrio sp017620235 (Lachnospiraceae)	vMAG_17	Lytic	Low-quality	39.28
bin_58	g_F23-D06	vMAG_51	-	Low-quality	37.32
bin_74	g_UBA1711 (P3)	vMAG_42	Lytic	Low-quality	34.54
bin_73	Non-classified	vMAG_27	-	Low-quality	25.64
bin_49	Non-classified	vMAG_6	-	Low-quality	24.02
bin_59	RUG191 sp002373675 (Lachnospiraceae)	vMAG_44	Lytic	Low-quality	23.57
bin_7	NK4A144 sp902783395 (Lachnospiraceae)	vMAG_24	Lytic	Low-quality	8.71
bin_26	Saccharofermentans sp902782235	vMAG_49	Lytic	Not-determined	NA

g_(genus).

Of the 17 bins identified to be associated with viruses, 8 belonged to the Lachnospiraceae family. The bin_30 and bin_41 belong to the same genus UBA2868, and the other bins assigned to this family were: (bin_7) NK4A144 sp902783395, (bin_12) UBA2868 sp002368575, (bin_16) RUG12461 sp016286115, (bin_22) UBA3766 sp902803015, (bin_27) Butyrivibrio sp017620235, and (bin_59) RUG191 sp002373675. The bin_26 and bin_28 were assigned to Saccharofermentans sp902782235 and sp003543635, respectively, and bin_58 was assigned to genus F23-D06 and bin_74 to genus UBA1711, both belonging to the order Bacteriodales.

Metagenomics has enabled the efficient identification of microorganisms, enzymes, and metabolic pathways that play a role in plant breakdown in the bovine rumen²⁹. This tool has also proven to be crucial for identifying the microbiota and characterizing important traits such as methane production and feed efficiency^3,8,30. This study is the first to explore rumen mobile genetic elements of Nellore cattle fed a grass diet by associating viruses and plasmids with their hosts through physical links between DNA molecules from the same genome using the Hi-C method²¹. We identified 31 links between bacteria and mobile genetic elements, with 14 of these links being plasmids and 17 links being viruses. Furthermore, we identified 12 genes conferring antibiotic resistance to tetracycline (6 genes), nitroimidazole (1 gene), lincosamide (1 gene), beta-lactam (1 gene), aminoglycoside (1 gene) and macrolide (1 gene). Additionally, we assembled a collection of 107 bacterial genomes, which mostly encoded pathways for central carbon and other carbohydrate metabolisms.

In this study, we identified 52 plasmid sequences, 14 of which were linked to their hosts using Hi-C. Similarly, in their evaluation of the canine fecal metagenome, Cuscó et al.²³ found that six Hi-C plasmids were ligated into their bacterial hosts, including five circular plasmids and a plasmid carrying the linA resistance gene. In our investigation, we uncovered two plasmids, k141_1760976 and k141_5289202, which are associated with multiple host clusters. As previously reported by Mo et al.³¹ and Stewart et al.²⁷, sequences may be shared among various bacterial species found in the rumen. Through the Hi-C approach, Stalder et al.³² found that plasmids are efficient vectors for horizontal gene transfer and that this method is useful for tracing microorganisms that harbor antibiotic-resistant genes. We identified 12 genes conferring antibiotic resistance, and among these, six genes were found to confer resistance to tetracycline: tet32, tet40, tet44, tetO, tetQ, and tetW. These findings are consistent with a previous study by Jing and Yan³³, in which tetracycline resistance genes (tet44, tetQ, and tetW) were the most prevalent in rumen content samples and were the most abundant in the genome of Prevotellaceae, but Lachnospiraceae were also present. Furthermore, tetQ, tetW, and tet40 were also detected in most samples derived from the rumen of 48 beef cattle from three different taurine breeds³⁴.

In addition to the genes that confer resistance to tetracycline, other classes of antibiotics, such as nitroimidazole nimJ, macrolide mefA, lincosamide lnuC, beta-lactam blaACI-1, and aminoglycoside aadE, were also identified. Genes that confer resistance to macrolides, lincosamide, and aminoglycosides were the most prevalent among ARGs in the study conducted by Ma et al.³⁴, with the mefA gene being one of the most common. Genes that confer resistance to macrolides, lincosamide, and aminoglycosides were the most prevalent among ARGs in the study conducted by Ma et al.³⁴, with the mefA gene being one of the most common. Additionally, ARGs associated with plasmids were identified, indicating that plasmids should be the focus of future investigations on antimicrobial resistance in livestock. Auffret et al.³⁵ also found an increase in ARGs in rumen content samples from cattle fed diets with higher concentrate content, with an abundance of genes related to resistance to macrolides and beta-lactams, suggesting that this increase is associated with dysbiosis caused by more concentrated diets. The findings of this study are in agreement with those of an earlier in silico study conducted by Sabino et al.³⁶, who analyzed the rumen resistome in 435 microbial genomes and found a high prevalence of antibiotic resistance genes (ARGs) that provide resistance to β-lactams, glycopeptides, tetracyclines, and aminoglycosides. Additionally, recent metagenomic studies have indicated that bovine rumen is a significant contributor to antibiotic resistance³³.

However, advancements in bioinformatics have contributed to a deeper understanding and better characterization of the rumen virome. Yan et al.³⁷ identified 397,180 viral operational taxonomic units (vOTUs) from 975 metagenomes of 13 ruminant species, with the majority of metagenomes derived from taurine animals. Only 23 metagenomes were from Bos indicus raised in Kenya, and none were from Brazil, supporting the significance of this research. The viral metagenomic assemblages discovered in this study were mostly assigned to unclassified genomes, which is probably because of the limited rumen virome of Bos indicus data in the ICTV phage database. This finding is consistent with those of Sato et al.³⁸, who observed low numbers of viral operational taxonomic units shared with the RefSeq database for samples from the rumen of Japanese cattle. The assigned vMAGs belong to the families Myoviridae and Siphoviridae, which are commonly found in the rumen environment and constitute the most abundant group, along with Mimiviridae and Podoviridae ^12,15,38,39.

A group of bacteria classified as UBA2868 sp003535955 (unclassified Lachnospiraceae bacterium, NCBI) was found to have the highest percentage of complete vMAGs (100%) and was characterized as high and medium quality. This genus was previously identified in a community of uncultured microorganisms from the intestines of pigs through sequencing of fecal samples^40,41. Bacteria from the unclassified family Lachnospiraceae are prevalent in the ruminal environment and serve as the central microbiome in cattle⁹. This allowed us to infer that population-associated viruses play a crucial role in the bovine rumen, rather than populations present in low abundance, which remains a challenge²⁰.

Friedersdorff et al.⁴² successfully isolated and sequenced active lytic phages belonging to the Siphoviridae family that infect Butyrivibrio fibrisolvens, which was found in both the rumen and feces of cattle and sheep. We discovered a lytic phage in Butyrivibrio sp017620235, which is commonly found in the rumen and promotes the degradation of lignocellulose and fermented carbohydrates into butyrate, formate, lactate, and acetate^43–45. Two bacterial species classified as Rumen Uncultured Genomes (RUGs), (bin_16) RUG12461 sp016286115 and (bin_59) RUG191 sp002373675, were also associated with viruses having a lytic life cycle.

The study's assembly of 107 bacterial genomes achieved a high completeness (> 95%) for members of the Treponema genus, which are known as spirochetes and have been positively correlated with feed conversion in cattle³. These species are considered to be one of the main fiber degraders in the rumen of Gir cattle, along with other bacteria, such as Clostridium, Ruminococcus, Eubacterium, Butyrivibrio, Roseburia, Caldicellulosiruptor, and Rhodospirillum⁴⁶. The bins identified were primarily assigned to families, such as Lachnospiraceae, Bacteroidaceae, P3, Ruminococcaceae, Saccharofermentanaceae, and Treponemataceae. These bins mainly encoded pathways related to central carbon metabolism and metabolism of various carbohydrates. One notable pathway found between these bins is the phosphate acetyltransferase-acetate kinase pathway, which is activated at high acetate concentrations, requires less ATP, and plays a role in energy metabolism ⁴⁷. Additionally, the aromatic amino acid metabolism pathway with tryptophan biosynthesis was abundant among the bins in the current study. This pathway occurs only in microorganisms and plants, and tryptophan is obtained by animals via symbiosis for protein synthesis⁴⁸.

The identification of CAZymes in the bacterial genomes revealed the presence of genes encoding glycoside hydrolases and polysaccharide lyase groups, with GH2, GH3, GH5, GH13, and GH43 being the most prevalent. These groups degrade xylan polysaccharides, which are heterogeneous and require different catalytic enzymes¹⁰. In a study by Wang et al.⁴⁹, the same GH families, GH2, GH3, GH13, and GH43, were found to comprise the majority of CAZymes present in the rumen of dairy cows fed different proportions of roughage and concentrate. The GH3 group was particularly more abundant in the high-roughage diet, and these observations are consistent with the current study's findings for pasture-fed Nellore cows, where glycoside hydrolase groups were the most prevalent.

In summary, this study provides a thorough examination of the MGEs in the rumen of Nellore cattle, shedding new light on the connection between these elements and their microbial inhabitants. The impact and function of these MGEs in beef cattle are currently being investigated by our research group.

Animal characterization and sample collection

Samples of ruminal content were collected from four Nellore (Bos taurus indicus) nulliparous cows cannulated in the rumen, with an average weight of (600 ± 50kg) and age between 9 and 6 years of the Ruminant Nutrition Laboratory, Faculty of Veterinary Medicine and Animal Science, FMVZ, Pirassununga campus. The experimental protocols were approved by Ethics Committee on the Use of Animals of Faculty of Veterinary Medicine and Animal Science (CEUA/FMVZ) under number 7333211118. These protocols adhere to both their guidelines and the ARRIVE guidelines. The animals received a diet exclusively based on pasture and mineral salts. Samples were kept on ice during collection until processing and were subsequently stored at -80°C.

DNA extraction and sample preparation for Hi-C

DNA extraction was performed using a QIAmp® Fast DNA Stool Mini Kit (Qiagen, USA). The crosslinking procedure was performed, in which the samples were processed according to the protocol recommended by the company Phase Genomics (Seattle, WA, USA) adapted by Burton et al.⁵⁰. Briefly, the sample containing liquid and fiber (5 mL) were suspended in 1% formaldehyde (45 ml) and incubated (20 min) at room temperature with periodic shaking. Subsequently, glycine (1 g/100 mL) was added, which was then incubated (15 min) at room temperature with periodic shaking. The sample was centrifuged (1000 × g for 1 min) and washed with PBS. It was then centrifuged again and the supernatant was removed. The pellets obtained formed a pool and were transferred to a cryotube (2 mL), stored at -80°C until sent along with the DNA extracted from the same sample, to be processed and sequenced by Phase Genomics.

Hi-C library preparation and short-read sequencing

To prepare Hi-C libraries, the Phase Genomics ProxiMeta Hi-C kit v4.0 was used as described by the manufacturer⁵¹. After going through the crosslinking process, the pellets underwent lysis and their DNA-protein complexes were digested by the restriction enzymes (Sau3AI and MlucI)³¹, creating free ends on the DNA strands that had labeled nucleotides with biotin to create chimeric molecules composed of fragments from different physically close genome regions in vivo. Proximity-linked DNA molecules were drawn down using streptavidin beads and processed into an Illumina-compatible sequencing library.

Separately, DNA was extracted from an aliquot of the original sample, and a shotgun metagenomic library was prepared using ProxiMeta library preparation reagents. Sequencing was performed using Illumina NovaSeq, generating pairs of PE150 reads for the Hi-C and shotgun libraries. Hi-C and shotgun metagenomic sequencing files were uploaded to the Phase Genomics cloud-based bioinformatics portal for subsequent analyses.

Clustering, quality assessment, taxonomic and functional assignment of assembled microorganisms

Shotgun-derived reads were filtered and adjusted for quality, normalized using fastp⁵², and assembled with MEGAHIT⁵³ using default options. Pairs of Hi-C reads were aligned to each assembly using BWA-MEM⁵⁴ with the − 5SP options specified and all other default options to disable attempts to pair the reads according to the normal Illumina settings. The SAMBLASTER program⁵⁵ was used to flag PCR duplicates that were subsequently excluded from the analysis. The alignments were then filtered with SAMtools⁵⁶ using the -F 2304 filter flag to remove nonprimary and secondary alignments.

Deconvolution, a method that uses the intracellular proximity signal captured by Hi-C as an indicator of the cellular origin of metagenome sequences, was performed using ProxiMeta^27,31 to generate putative genomes and clusters of genome fragments. Prokaryotic clusters were assessed for quality using CheckM⁵⁷ and received preliminary taxonomic classifications using Mash, a generalized and automated method that identifies reference genomes compatible with the NCBI RefSeq database⁵⁸ (Supplementary Table 1). Subsequently, the bins were submitted to the Genome Taxonomy Database (GTDB)⁵⁹ for taxonomic assignment (Supplementary Table 2).

The functional potential of the bins was evaluated by METABOLIC (METabolic and BiogeOchemistry analyses in microCrobes)⁶⁰, which detected the presence or absence of KEGG modules for the most abundant central pathways present in the genome of the microorganisms evaluated, as well as how they attributed the CAZymes and the main biogeochemical cycles present.

Mobile genetic elements

Viral contigs were annotated and used for lifestyle and viral protein annotation per contig using VIBRANT⁶¹. The contigs were grouped into assembled genomes using ProxiPhage viral binning. Sequence integrity was validated with CheckV, viral genomes were classified into quality levels based on AAI or HMM estimation, high confidence (≥ 90% complete), medium confidence (80–90% complete) or low confidence (< 80% complete)^22,62. For taxonomic assignment, the protein sequences of the viral clusters (vMAGs) were predicted using Prodigal⁶³, and these sequences were used to build protein cluster profiles that were compared with sequences present in the ICTV phage database, generating networks of similarities between viral clusters with the using vCONTACT⁶⁴.

To identify sequences with plasmid characteristics, a BLAST (Basic Local Alignment Search Tool) alignment was performed against the NCBI plasmid database³¹. Sequences annotated as plasmids were associated with the host using a host-finding algorithm²². Microbial resistance genes, in turn, were identified with AMRFinderPlus, which allows the assessment of Antimicrobial Resistance gene content in the NCBI database⁶⁵.

Assignment of mobile genetic elements to their hosts

The step of assigning mobile genetic elements to their hosts was composed of several filtering criteria, with the first phase of filtering restricting connections with less than 2 Hi-C read links. To remove false positives, a connectivity rate (≥ 0.1) was used for the number of copies of mobile elements per cell (for repeated vMAGs the highest number of counts per cell was considered) and host intra-MAG connectivity (≥ 10links). Finally, a receiver operating characteristic (ROC) curve was performed, which determines the ideal copy count cutoff value (≥ 0.1). This value is expected to be close to 1, as values close to or above 1 are considered reliable for attributing a real host²². In cases where the optimal copy count threshold is excessively low (< 0.1), it is replaced with 0.1. This occurred in the present study, where the calculated value was 0.01.

Data Availability

The data generated from the present study will be available in the Sequence Read Archive (SRA) database https://dataview.ncbi.nlm.nih.gov/object/PRJNA1054691?reviewer=vd7njv2i8dqvoj4d6tmt82kr1t

Acknowledgements

The authors express their gratitude to the funding agencies that facilitated the realization of this research. The Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - Brasil (CAPES) - Finance Code 001, Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) grant 2022/05541-9 and the Fundação de Estudos Agrários Luiz de Queiroz (FEALQ).

Contributions

C.A.F collected the samples, performed the laboratory procedures until sending them to Phase Genomics for sequencing, analyzed the data obtained, compiled the results and wrote the initial manuscript. A.T.N. collaborated with sampling and writing the manuscript. M.D.P. contributed to data analysis and preparation of figures. E.C.M.O., O.S.G. and P.A.A. contributed with bioinformatics analyses. F.P.J. and P.H.M.R collected the samples and provided the animal data. HF was the project supervisor and reviewed the initial manuscript. All authors reviewed the final version of the manuscript.

Ethics declarations

Competing Interests

The authors declare no competing interests.

González-Recio, O. et al. Invited review: Novel methods and perspectives for modulating the rumen microbiome through selective breeding as a means to improve complex traits: Implications for methane emissions in cattle. Livest Sci 269, 105171 (2023).
Nkrumah, J. D. et al. Relationships of feedlot feed efficiency, performance, and feeding behavior with metabolic rate, methane production, and energy partitioning in beef cattle. J Anim Sci 84, 145–153 (2006).
Auffret, M. D. et al. Identification of Microbial Genetic Capacities and Potential Mechanisms Within the Rumen Microbiome Explaining Differences in Beef Cattle Feed Efficiency. Front Microbiol 11, 1229 (2020).
McGovern, E. et al. Investigation into the effect of divergent feed efficiency phenotype on the bovine rumen microbiota across diet and breed. Scientific Reports 2020 10:1 10, 1–11 (2020).
Liu, Y. et al. Rumen Microbiome and Metabolome of High and Low Residual Feed Intake Angus Heifers. Front Vet Sci 9, 812861 (2022).
Roehe, R. et al. Bovine Host Genetic Variation Influences Rumen Microbial Methane Production with Best Selection Criterion for Low Methane Emitting and Efficiently Feed Converting Hosts Based on Metagenomic Gene Abundance. PLoS Genet 12, e1005846 (2016).
Zhang, Q. et al. Bayesian modeling reveals host genetics associated with rumen microbiota jointly influence methane emission in dairy cows. ISME J 14, 2019–2033 (2020).
Martínez-Álvaro, M. et al. Bovine host genome acts on rumen microbiome function linked to methane emissions. Communications Biology 2022 5:1 5, 1–16 (2022).
Stewart, R. D. et al. Compendium of 4,941 rumen metagenome-assembled genomes for rumen microbiome biology and enzyme discovery. Nature Biotechnology 2019 37:8 37, 953–961 (2019).
Gharechahi, J. et al. Lignocellulose degradation by rumen bacterial communities: New insights from metagenome analyses. Environ Res 229, 115925 (2023).
Luo, X. Q. et al. Viral community-wide auxiliary metabolic genes differ by lifestyles, habitats, and hosts. Microbiome 10, 1–18 (2022).
Anderson, C. L., Sullivan, M. B. & Fernando, S. C. Dietary energy drives the dynamic response of bovine rumen viral communities. Microbiome 2017 5:1 5, 1–19 (2017).
Kav, A. B. et al. Insights into the bovine rumen plasmidome. Proc Natl Acad Sci U S A 109, 5452–5457 (2012).
Lobo, R. R. & Faciola, A. P. Ruminal Phages – A Review. Front Microbiol 12, 763416 (2021).
Solden, L. M. et al. Interspecies cross-feeding orchestrates carbon degradation in the rumen ecosystem. Nature Microbiology 2018 3:11 3, 1274–1284 (2018).
Orpin, C. G. & Munn, E. A. The occurrence of bacteriophages in the rumen and their influence on rumen bacterial populations. Experientia 30, 1018–1020 (1974).
Altermann, E., Schofield, L. R., Ronimus, R. S., Beatty, A. K. & Reilly, K. Inhibition of Rumen methanogens by a novel archaeal lytic enzyme displayed on tailored bionanoparticles. Front Microbiol 9, 2378 (2018).
Park, S. Y. et al. Characterization of two lytic bacteriophages, infecting Streptococcus bovis/equinus complex (SBSEC) from Korean ruminant. Scientific Reports 2023 13:1 13, 1–16 (2023).
Ross, E. M., Petrovski, S., Moate, P. J. & Hayes, B. J. Metagenomics of rumen bacteriophage from thirteen lactating dairy cattle. BMC Microbiol 13, 1–11 (2013).
Bickhart, D. M. et al. Assignment of virus and antimicrobial resistance genes to microbial hosts in a complex microbial community by combined long-read assembly and proximity ligation. Genome Biology 2019 20:1 20, 1–18 (2019).
Marbouty, M., Baudry, L., Cournac, A. & Koszul, R. Scaffolding bacterial genomes and probing host-virus interactions in gut microbiome by proximity ligation (chromosome capture) assay. Sci Adv 3, (2017).
Uritskiy, G. et al. Accurate viral genome reconstruction and host assignment with proximity-ligation sequencing. bioRxiv 2021.06.14.448389 (2021) doi:10.1101/2021.06.14.448389.
Cuscó, A., Pérez, D., Viñes, J., Fàbregas, N. & Francino, O. Novel canine high-quality metagenome-assembled genomes, prophages and host-associated plasmids provided by long-read metagenomics together with Hi-C proximity ligation. Microb Genom 8, 802 (2022).
Braun, V. & Killmann, H. Bacterial solutions to the iron-supply problem. Trends Biochem Sci 24, 104–109 (1999).
Canfield, D. E., Kristensen, E. & Thamdrup, B. The Iron and Manganese Cycles. Adv Mar Biol 48, 269–312 (2005).
Bohra, V., Dafale, N. A. & Purohit, H. J. Understanding the alteration in rumen microbiome and CAZymes profile with diet and host through comparative metagenomic approach. Arch Microbiol 201, 1385–1397 (2019).
Stewart, R. D. et al. Assembly of 913 microbial genomes from metagenomic sequencing of the cow rumen. Nat Commun 9, 870–870 (2018).
Krupovic, M. et al. Bacterial Viruses Subcommittee and Archaeal Viruses Subcommittee of the ICTV: update of taxonomy changes in 2021. Arch Virol 166, 3239–3244 (2021).
Hua, D., Hendriks, W. H., Xiong, B. & Pellikaan, W. F. Starch and Cellulose Degradation in the Rumen and Applications of Metagenomics on Ruminal Microorganisms. Animals 2022, Vol. 12, Page 3020 12, 3020 (2022).
Wallace, R. J. et al. The rumen microbial metagenome associated with high methane production in cattle. BMC Genomics 2015 16:1 16, 1–14 (2015).
MO, P. et al. Hi-C deconvolution of a human gut microbiome yields high-quality draft genomes and reveals plasmid-genome interactions. (2017) doi:10.1101/198713.
Stalder, T., Press, M. O., Sullivan, S., Liachko, I. & Top, E. M. Linking the resistome and plasmidome to the microbiome. The ISME Journal 2019 13:10 13, 2437–2446 (2019).
Jing, R. & Yan, Y. Metagenomic analysis reveals antibiotic resistance genes in the bovine rumen. Microb Pathog 149, 104350 (2020).
Ma, T. et al. Expressions of resistome is linked to the key functions and stability of active rumen microbiome. Anim Microbiome 4, 1–17 (2022).
Auffret, M. D. et al. The rumen microbiome as a reservoir of antimicrobial resistance and pathogenicity genes is directly affected by diet in beef cattle. Microbiome 5, 1–11 (2017).
Sabino, Y. N. V. et al. Characterization of antibiotic resistance genes in the species of the rumen microbiota. Nature Communications 2019 10:1 10, 1–11 (2019).
Yan, M. et al. Interrogating the viral dark matter of the rumen ecosystem with a global virome database. Nature Communications 2023 14:1 14, 1–16 (2023).
Sato, Y. et al. Characteristics of the rumen virome in Japanese cattle. bioRxiv 2023.03.20.532305 (2023) doi:10.1101/2023.03.20.532305.
Berg Miller, M. E. et al. Phage-bacteria relationships and CRISPR elements revealed by a metagenomic survey of the rumen microbiome. Environ Microbiol 14, 207–227 (2012).
Zhou, S. et al. Characterization of Metagenome-Assembled Genomes and Carbohydrate-Degrading Genes in the Gut Microbiota of Tibetan Pig. Front Microbiol 11, 595066 (2020).
Crossfield, M. et al. Archaeal and Bacterial Metagenome-Assembled Genome Sequences Derived from Pig Feces. Microbiol Resour Announc 11, (2022).
Friedersdorff, J. C. A. et al. The Isolation and Genome Sequencing of Five Novel Bacteriophages From the Rumen Active Against Butyrivibrio fibrisolvens. Front Microbiol 11, 522243 (2020).
Kopečný, J., Zorec, M., Mrázek, J., Kobayashi, Y. & Marinšek-Logar, R. Butyrivibrio hungatei sp. nov. and Pseudobutyrivibrio xylanivorans sp. nov., butyrate-producing bacteria from the rumen. Int J Syst Evol Microbiol 53, 201–209 (2003).
Moon, C. D. et al. Reclassification of Clostridium proteoclasticum as Butyrivibrio proteoclasticus comb. nov., a butyrate-producing ruminal bacterium. Int J Syst Evol Microbiol 58, 2041–2045 (2008).
Palevich, N. et al. Complete Genome Sequence of the Polysaccharide-Degrading Rumen Bacterium Pseudobutyrivibrio xylanivorans MA3014 Reveals an Incomplete Glycolytic Pathway. Genome Biol. Evol 12, 1566–1572 (2020).
Pandit, R. J. et al. Microbiota composition, gene pool and its expression in Gir cattle (Bos indicus) rumen under different forage diets using metagenomic and metatranscriptomic approaches. Syst Appl Microbiol 41, 374–385 (2018).
Zhang, S. et al. Metabolic engineering for efficient supply of acetyl-CoA from different carbon sources in Escherichia coli. Microb Cell Fact 18, 1–11 (2019).
Parthasarathy, A. et al. A Three-Ring circus: Metabolism of the three proteogenic aromatic amino acids and their role in the health of plants and animals. Front Mol Biosci 5, 342220 (2018).
Wang, L., Zhang, G., Xu, H., Xin, H. & Zhang, Y. Metagenomic analyses of microbial and carbohydrate-active enzymes in the rumen of holstein cows fed different forage-to-concentrate ratios. Front Microbiol 10, 441658 (2019).
Burton, J. N., Liachko, I., Dunham, M. J. & Shendure, J. Species-level deconvolution of metagenome assemblies with Hi-C-based contact probability maps. G3: Genes, Genomes, Genetics 4, 1339–1346 (2014).
Lieberman-Aiden, E. et al. Comprehensive Mapping of Long-Range Interactions Reveals Folding Principles of the Human Genome. Science (1979) 326, 289–293 (2009).
Chen, S., Zhou, Y., Chen, Y. & Gu, J. fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 34, i884–i890 (2018).
Li, D. et al. MEGAHIT v1.0: A fast and scalable metagenome assembler driven by advanced methodologies and community practices. Methods 102, 3–11 (2016).
Li, H. & Durbin, R. Fast and accurate long-read alignment with Burrows–Wheeler transform. Bioinformatics 26, 589–595 (2010).
Faust, G. G. & Hall, I. M. SAMBLASTER: fast duplicate marking and structural variant read extraction. Bioinformatics 30, 2503–2505 (2014).
Li, H. et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics 25, 2078 (2009).
Parks, D. H., Imelfort, M., Skennerton, C. T., Hugenholtz, P. & Tyson, G. W. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res 25, 1043 (2015).
Ondov, B. D. et al. Mash Screen: High-throughput sequence containment estimation for genome discovery. Genome Biol 20, 1–13 (2019).
Parks, D. H. et al. A standardized bacterial taxonomy based on genome phylogeny substantially revises the tree of life. Nature Biotechnology 2018 36:10 36, 996–1004 (2018).
Zhou, Z. et al. METABOLIC: high-throughput profiling of microbial genomes for functional traits, metabolism, biogeochemistry, and community-scale functional networks. Microbiome 10, 1–22 (2022).
Kieft, K., Zhou, Z. & Anantharaman, K. VIBRANT: Automated recovery, annotation and curation of microbial viruses, and evaluation of viral community function from genomic sequences. Microbiome 8, 1–23 (2020).
Nayfach, S. et al. CheckV assesses the quality and completeness of metagenome-assembled viral genomes. Nature Biotechnology 2020 39:5 39, 578–585 (2020).
Hyatt, D. et al. Prodigal: Prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics 11, 1–11 (2010).
Bolduc, B. et al. vConTACT: An iVirus tool to classify double-stranded DNA viruses that infect Archaea and Bacteria. PeerJ 2017, e3243 (2017).
Feldgarden, M. et al. AMRFinderPlus and the Reference Gene Catalog facilitate examination of the genomic links among antimicrobial resistance, stress response, and virulence. Scientific Reports 2021 11:1 11, 1–9 (2021).

No competing interests reported.

SupplementaryMaterial.docx.xlsx

Download PDF

Editorial decision: Revision requested
23 Apr, 2024
Reviews received at journal
19 Apr, 2024
Reviewers agreed at journal
04 Apr, 2024
Reviews received at journal
26 Feb, 2024
Reviewers agreed at journal
08 Feb, 2024
Reviewers invited by journal
06 Feb, 2024
Editor assigned by journal
29 Jan, 2024
Editor invited by journal
28 Dec, 2023
Submission checks completed at journal
28 Dec, 2023
First submitted to journal
13 Dec, 2023

You are reading this latest preprint version

Identifying Mobile Genetic Elements in the Ruminal Microbiome of Nellore Cattle: An Initial Investigation

Status:

Version 1

Abstract

Figures

Introduction

Results

Metagenome of bovine rumen content samples

Functional and carbohydrate-active enzymes (CAZYmes) characterization of bins

Assignment of mobile genetic elements to their hosts

Plasmids and the Identification of microbial resistance genes

Virus found in the rumen

Virus-host association

Discussion

Material and methods

Animal characterization and sample collection

DNA extraction and sample preparation for Hi-C

Hi-C library preparation and short-read sequencing

Clustering, quality assessment, taxonomic and functional assignment of assembled microorganisms

Mobile genetic elements

Assignment of mobile genetic elements to their hosts

Declarations

References

Additional Declarations

Supplementary Files

Status:

Version 1