Search

The rumen represents a dynamic microbial ecosystem where fermentation metabolites and microbial concentrations change over time in response to dietary changes. The integration of microbial genomic knowledge and dynamic modelling can enhance our system-level understanding of rumen ecosystem's function. However, such an integration between dynamic models and rumen microbiota data is lacking. The objective of this work was to integrate rumen microbiota time series determined by 16S rRNA gene amplicon sequencing into a dynamic modelling framework to link microbial data to the dynamics of the volatile fatty acids (VFA) production during fermentation. For that, we used the theory of state observers to develop a model that estimates the dynamics of VFA from the data of microbial functional proxies associated with the specific production of each VFA. We determined the microbial proxies using CowPi to infer the functional potential of the rumen microbiota and extrapolate their functional modules from KEGG (Kyoto Encyclopedia of Genes and Genomes). The approach was challenged using data from an in vitro RUSITEC experiment and from an in vivo experiment with four cows. The model performance was evaluated by the coefficient of variation of the root mean square error (CRMSE). For the in vitro case study, the mean CVRMSE were 9.8% for acetate, 14% for butyrate and 14.5% for propionate. For the in vivo case study, the mean CVRMSE were 16.4% for acetate, 15.8% for butyrate and 19.8% for propionate. The mean CVRMSE for the VFA molar fractions were 3.1% for acetate, 3.8% for butyrate and 8.9% for propionate. Ours results show the promising application of state observers integrated with microbiota time series data for predicting rumen microbial metabolism., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Davoudkhani et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

A deeper understanding of the relationship between the antimicrobial resistance (AMR) gene carriage and phenotype is necessary to develop effective response strategies against this global burden. AMR phenotype is often a result of multi-gene interactions; therefore, we need approaches that go beyond current simple AMR gene identification tools. Machine-learning (ML) methods may meet this challenge and allow the development of rapid computational approaches for AMR phenotype classification. To examine this, we applied multiple ML techniques to 16,950 bacterial genomes across 28 genera, with corresponding MICs for 23 antibiotics with the aim of training models to accurately determine the AMR phenotype from sequenced genomes. This resulted in a >1.5-fold increase in AMR phenotype prediction accuracy over AMR gene identification alone. Furthermore, we revealed 528 unique (often species-specific) genomic routes to antibiotic resistance, including genes not previously linked to the AMR phenotype. Our study demonstrates the utility of ML in predicting AMR phenotypes across diverse clinically relevant organisms and antibiotics. This research proposes a rapid computational method to support laboratory-based identification of the AMR phenotype in pathogens., (© 2024 Dillon et al.)

The rumen ecosystem harbours a galaxy of microbes working in syntrophy to carry out a metabolic cascade of hydrolytic and fermentative reactions. This fermentation process allows ruminants to harvest nutrients from a wide range of feedstuff otherwise inaccessible to the host. The interconnection between the ruminant and its rumen microbiota shapes key animal phenotypes such as feed efficiency and methane emissions and suggests the potential of reducing methane emissions and enhancing feed conversion into animal products by manipulating the rumen microbiota. Whilst significant technological progress in omics techniques has increased our knowledge of the rumen microbiota and its genome (microbiome), translating omics knowledge into effective microbial manipulation strategies remains a great challenge. This challenge can be addressed by modelling approaches integrating causality principles and thus going beyond current correlation-based approaches applied to analyse rumen microbial genomic data. However, existing rumen models are not yet adapted to capitalise on microbial genomic information. This gap between the rumen microbiota available omics data and the way microbial metabolism is represented in the existing rumen models needs to be filled to enhance rumen understanding and produce better predictive models with capabilities for guiding nutritional strategies. To fill this gap, the integration of computational biology tools and mathematical modelling frameworks is needed to translate the information of the metabolic potential of the rumen microbes (inferred from their genomes) into a mathematical object. In this paper, we aim to discuss the potential use of two modelling approaches for the integration of microbial genomic information into dynamic models. The first modelling approach explores the theory of state observers to integrate microbial time series data into rumen fermentation models. The second approach is based on the genome-scale network reconstructions of rumen microbes. For a given microorganism, the network reconstruction produces a stoichiometry matrix of the metabolism. This matrix is the core of the so-called genome-scale metabolic models which can be exploited by a plethora of methods comprised within the constraint-based reconstruction and analysis approaches. We will discuss how these methods can be used to produce the next-generation models of the rumen microbiome., (Copyright © 2023 The Author(s). Published by Elsevier B.V. All rights reserved.)

Large regions of prokaryotic genomes are currently without any annotation, in part due to well-established limitations of annotation tools. For example, it is routine for genes using alternative start codons to be misreported or completely omitted. Therefore, we present StORF-Reporter, a tool that takes an annotated genome and returns regions that may contain missing CDS genes from unannotated regions. StORF-Reporter consists of two parts. The first begins with the extraction of unannotated regions from an annotated genome. Next, Stop-ORFs (StORFs) are identified in these unannotated regions. StORFs are open reading frames that are delimited by stop codons and thus can capture those genes most often missing in genome annotations. We show this methodology recovers genes missing from canonical genome annotations. We inspect the results of the genomes of model organisms, the pangenome of Escherichia coli, and a set of 5109 prokaryotic genomes of 247 genera from the Ensembl Bacteria database. StORF-Reporter extended the core, soft-core and accessory gene collections, identified novel gene families and extended families into additional genera. The high levels of sequence conservation observed between genera suggest that many of these StORFs are likely to be functional genes that should now be considered for inclusion in canonical annotations., (© The Author(s) 2023. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Background: Manipulating the rhizosphere microbial community through beneficial microorganism inoculation has gained interest in improving crop productivity and stress resistance. Synthetic microbial communities, known as SynComs, mimic natural microbial compositions while reducing the number of components. However, achieving this goal requires a comprehensive understanding of natural microbial communities and carefully selecting compatible microorganisms with colonization traits, which still pose challenges. In this study, we employed multi-genome metabolic modeling of 270 previously described metagenome-assembled genomes from Campos rupestres to design a synthetic microbial community to improve the yield of important crop plants., Results: We used a targeted approach to select a minimal community (MinCom) encompassing essential compounds for microbial metabolism and compounds relevant to plant interactions. This resulted in a reduction of the initial community size by approximately 4.5-fold. Notably, the MinCom retained crucial genes associated with essential plant growth-promoting traits, such as iron acquisition, exopolysaccharide production, potassium solubilization, nitrogen fixation, GABA production, and IAA-related tryptophan metabolism. Furthermore, our in-silico selection for the SymComs, based on a comprehensive understanding of microbe-microbe-plant interactions, yielded a set of six hub species that displayed notable taxonomic novelty, including members of the Eremiobacterota and Verrucomicrobiota phyla., Conclusion: Overall, the study contributes to the growing body of research on synthetic microbial communities and their potential to enhance agricultural practices. The insights gained from our in-silico approach and the selection of hub species pave the way for further investigations into the development of tailored microbial communities that can optimize crop productivity and improve stress resilience in agricultural systems., (© 2023. The Author(s).)

Background: Honeybees use plant material to manufacture their own food. These insect pollinators visit flowers repeatedly to collect nectar and pollen, which are shared with other hive bees to produce honey and beebread. While producing these products, beehives accumulate a considerable number of microbes, including bacteria that derive from plants and different parts of the honeybees' body. Whether bacteria form similar communities amongst beehives, even if located in close proximity, is an ecologically important question that has been addressed in this study. Specific ecological factors such as the surrounding environment and the beekeeping methods used can shape the microbiome of the beehive as a whole, and eventually influence the health of the honeybees and their ecosystem., Results: We conducted 16S rRNA meta-taxonomic analysis on honey and beebread samples that were collected from 15 apiaries in the southeast of England to quantify the bacteria associated with different beehives. We observed that honeybee products carry a significant variety of bacterial groups that comprise bee commensals, environmental bacteria and symbionts and pathogens of plants and animals. Remarkably, this bacterial diversity differs not only amongst apiaries, but also between the beehives of the same apiary. In particular, the levels of the bee commensals varied significantly, and their fluctuations correlated with the presence of different environmental bacteria and various apiculture practices., Conclusions: Our results show that every hive possesses their own distinct microbiome and that this very defined fingerprint is affected by multiple factors such as the nectar and pollen gathered from local plants, the management of the apiaries and the bacterial communities living around the beehives. Based on our findings, we suggest that the microbiome of beehives could be used as a valuable biosensor informing of the health of the honeybees and their surrounding environment., (© 2023. The Author(s).)

There is conflicting evidence as to whether Porifera (sponges) or Ctenophora (comb jellies) comprise the root of the animal phylogeny. Support for either a Porifera-sister or Ctenophore-sister tree has been extensively examined in the context of model selection, taxon sampling, and outgroup selection. The influence of dataset construction is comparatively understudied. We re-examine five animal phylogeny datasets that have supported either root hypothesis using an approach designed to enrich orthologous signal in phylogenomic datasets. We find that many component orthogroups in animal datasets fail to recover major lineages as monophyletic with the exception of Ctenophora, regardless of the supported root. Enriching these datasets to retain orthogroups recovering ≥3 major lineages reduces dataset size by up to 50% while retaining underlying phylogenetic information and taxon sampling. Site-heterogeneous phylogenomic analysis of these enriched datasets recovers both Porifera-sister and Ctenophora-sister positions, even with additional constraints on outgroup sampling. Two datasets which previously supported Ctenophora-sister support Porifera-sister upon enrichment. All enriched datasets display improved model fitness under posterior predictive analysis. While not conclusively rooting animals at either Porifera or Ctenophora, we do see an increase in signal for Porifera-sister and a decrease in signal for Ctenophore-sister when data are filtered for orthologous signal. Our results indicate that dataset size and construction as well as model fit influence animal root inference., (© The Author(s) 2023. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Conflicting studies place a group of bilaterian invertebrates containing xenoturbellids and acoelomorphs, the Xenacoelomorpha, as either the primary emerging bilaterian phylum 1 , 2 , 3 , 4 , 5 , 6 or within Deuterostomia, sister to Ambulacraria. 7 , 8 , 9 , 10 , 11 Although their placement as sister to the rest of Bilateria supports relatively simple morphology in the ancestral bilaterian, their alternative placement within Deuterostomia suggests a morphologically complex ancestral bilaterian along with extensive loss of major phenotypic traits in the Xenacoelomorpha. Recent studies have questioned whether Deuterostomia should be considered monophyletic at all. 10 , 12 , 13 Hidden paralogy and poor phylogenetic signal present a major challenge for reconstructing species phylogenies. 14 , 15 , 16 , 17 , 18 Here, we assess whether these issues have contributed to the conflict over the placement of Xenacoelomorpha. We reanalyzed published datasets, enriching for orthogroups whose gene trees support well-resolved clans elsewhere in the animal tree. 16 We find that most genes in previously published datasets violate incontestable clans, suggesting that hidden paralogy and low phylogenetic signal affect the ability to reconstruct branching patterns at deep nodes in the animal tree. We demonstrate that removing orthogroups that cannot recapitulate incontestable relationships alters the final topology that is inferred, while simultaneously improving the fit of the model to the data. We discover increased, but ultimately not conclusive, support for the existence of Xenambulacraria in our set of filtered orthogroups. At a time when we are progressing toward sequencing all life on the planet, we argue that long-standing contentious issues in the tree of life will be resolved using smaller amounts of better quality data that can be modeled adequately. 19 ., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 The Author(s). Published by Elsevier Inc. All rights reserved.)

Microbiomes are rife for biotechnological exploitation, particularly the rumen microbiome, due to their complexicity and diversity. In this study, antimicrobial peptides (AMPs) from the rumen microbiome (Lynronne 1, 2, 3 and P15s) were assessed for their therapeutic potential against seven clinical strains of Pseudomonas aeruginosa. All AMPs exhibited antimicrobial activity against all strains, with minimum inhibitory concentrations (MICs) ranging from 4-512 µg/mL. Time-kill kinetics of all AMPs at 3× MIC values against strains PAO1 and LES431 showed complete kill within 10 min to 4 h, although P15s was not bactericidal against PAO1. All AMPs significantly inhibited biofilm formation by strains PAO1 and LES431, and induction of resistance assays showed no decrease in activity against these strains. AMP cytotoxicity against human lung cells was also minimal. In terms of mechanism of action, the AMPs showed affinity towards PAO1 and LES431 bacterial membrane lipids, efficiently permeabilising the P. aeruginosa membrane. Transcriptome and metabolome analysis revealed increased catalytic activity at the cell membrane and promotion of β-oxidation of fatty acids. Finally, tests performed with the Galleria mellonella infection model showed that Lynronne 1 and 2 were efficacious in vivo, with a 100% survival rate following treatment at 32 mg/kg and 128 mg/kg, respectively. This study illustrates the therapeutic potential of microbiome-derived AMPs against P. aeruginosa infections., (© 2022. The Author(s).)

Here we report two antimicrobial peptides (AMPs), HG2 and HG4 identified from a rumen microbiome metagenomic dataset, with activity against multidrug-resistant (MDR) bacteria, especially methicillin-resistant Staphylococcus aureus (MRSA) strains, a major hospital and community-acquired pathogen. We employed the classifier model design to analyse, visualise, and interpret AMP activities. This approach allowed in silico discrimination of promising lead AMP candidates for experimental evaluation. The lead AMPs, HG2 and HG4, are fast-acting and show anti-biofilm and anti-inflammatory activities in vitro and demonstrated little toxicity to human primary cell lines. The peptides were effective in vivo within a Galleria mellonella model of MRSA USA300 infection. In terms of mechanism of action, HG2 and HG4 appear to interact with the cytoplasmic membrane of target cells and may inhibit other cellular processes, whilst preferentially binding to bacterial lipids over human cell lipids. Therefore, these AMPs may offer additional therapeutic templates for MDR bacterial infections., (© 2022. The Author(s).)

Antimicrobial resistance (AMR) is a serious threat to public health globally; it is estimated that AMR bacteria caused 1.27 million deaths in 2019, and this is set to rise to 10 million deaths annually. Agricultural and soil environments act as antimicrobial resistance gene (ARG) reservoirs, operating as a link between different ecosystems and enabling the mixing and dissemination of resistance genes. Due to the close interactions between humans and agricultural environments, these AMR gene reservoirs are a major risk to both human and animal health. In this study, we aimed to identify the resistance gene reservoirs present in four microbiomes: poultry, ruminant, swine gastrointestinal (GI) tracts coupled with those from soil. This large study brings together every poultry, swine, ruminant, and soil shotgun metagenomic sequence available on the NCBI sequence read archive for the first time. We use the ResFinder database to identify acquired antimicrobial resistance genes in over 5,800 metagenomes. ARGs were diverse and widespread within the metagenomes, with 235, 101, 167, and 182 different resistance genes identified in the poultry, ruminant, swine, and soil microbiomes, respectively. The tetracycline resistance genes were the most widespread in the livestock GI microbiomes, including tet (W)_1, tet (Q)_1, tet (O)_1, and tet (44)_1. The tet (W)_1 resistance gene was found in 99% of livestock GI tract microbiomes, while tet (Q)_1 was identified in 93%, tet (O)_1 in 82%, and finally tet (44)_1 in 69%. Metatranscriptomic analysis confirmed these genes were "real" and expressed in one or more of the livestock GI tract microbiomes, with tet (40)_1 and tet (O)_1 expressed in all three livestock microbiomes. In soil, the most abundant ARG was the oleandomycin resistance gene, ole (B)_1. A total of 55 resistance genes were shared by the four microbiomes, with 11 ARGs actively expressed in two or more microbiomes. By using all available metagenomes we were able to mine a large number of samples and describe resistomes in 37 countries. This study provides a global insight into the diverse and abundant antimicrobial resistance gene reservoirs present in both livestock and soil microbiomes., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Lawther, Santos, Oyama, Rubino, Morrison, Creevey, McGrath and Huws.)

With an increasing human population access to ruminant products is an important factor in global food supply. While ruminants contribute to climate change, climate change could also affect ruminant production. Here we investigated how the plant response to climate change affects forage quality and subsequent rumen fermentation. Models of near future climate change (2050) predict increases in temperature, CO 2 , precipitation and altered weather systems which will produce stress responses in field crops. We hypothesised that pre-exposure to altered climate conditions causes compositional changes and also primes plant cells such that their post-ingestion metabolic response to the rumen is altered. This "stress memory" effect was investigated by screening ten forage grass varieties in five differing climate scenarios, including current climate (2020), future climate (2050), or future climate plus flooding, drought or heat shock. While varietal differences in fermentation were detected in terms of gas production, there was little effect of elevated temperature or CO 2 compared with controls (2020). All varieties consistently showed decreased digestibility linked to decreased methane production as a result of drought or an acute flood treatment. These results indicate that efforts to breed future forage varieties should target tolerance of acute stress rather than long term climate., (© 2022. The Author(s).)

Antimicrobial peptides (AMPs) can efficiently control different microbial pathogens and show the potential to be applied in clinical practice and livestock production. In this work, the aim was to isolate AMP-producing ruminal streptococci and to characterize their genetic features through whole-genome sequencing. We cultured 463 bacterial isolates from the rumen of Nelore bulls, 81 of which were phenotypically classified as being Streptococcaceae . Five isolates with broad-range activity were genome sequenced and confirmed as being Streptococcus lutetiensis . The genetic features linked to their antimicrobial activity or adaptation to the rumen environment were characterized through comparative genomics. The genome of S. lutetiensis UFV80 harbored a putative CRISPR-Cas9 system (Type IIA). Computational tools were used to discover novel biosynthetic clusters linked to the production of bacteriocins. All bacterial genomes harbored genetic clusters related to the biosynthesis of class I and class II bacteriocins. SDS-PAGE confirmed the results obtained in silico and demonstrated that the class II bacteriocins predicted in the genomes of three S. lutetiensis strains had identical molecular mass (5197 Da). These results demonstrate that ruminal bacteria of the Streptococcus bovis / equinus complex represent a promising source of novel antimicrobial peptides.

Bacteriophages (phages) are viruses that target bacteria, with the ability to lyse and kill host bacterial cells. Due to this, they have been of some interest as a therapeutic since their discovery in the early 1900s, but with the recent increase in antibiotic resistance, phages have seen a resurgence in attention. Current methods of isolation and purification of phages can be long and tedious, with caesium chloride concentration gradients the gold standard for purifying a phage fraction. Isolation of novel phages requires centrifugation and ultrafiltration of mixed samples, such as water sources, effluent or faecal samples etc, to prepare phage filtrates for further testing. We propose countercurrent chromatography as a novel and alternative approach to use when studying phages, as a scalable and high-yield method for obtaining phage fractions. However, the full extent of the usefulness and resolution of separation with this technique has not been researched; it requires optimization and ample testing before this can be revealed. Here we present an initial study to determine survivability of two phages, T4 and ϕX174, using only water as a mobile phase in a Spectrum Series 20 HPCCC. Both phages were found to remain active once eluted from the column. Phages do not fully elute from the column and sodium hydroxide is necessary to flush the column between runs to deactivate remaining phages., Competing Interests: D.R. and C.B. are employed by Dynamic Extractions Ltd. J.C.A.F. has been an employee of Dynamic Extractions Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (© 2022 The Authors.)

Motivation: The biases in CoDing Sequence (CDS) prediction tools, which have been based on historic genomic annotations from model organisms, impact our understanding of novel genomes and metagenomes. This hinders the discovery of new genomic information as it results in predictions being biased towards existing knowledge. To date, users have lacked a systematic and replicable approach to identify the strengths and weaknesses of any CDS prediction tool and allow them to choose the right tool for their analysis., Results: We present an evaluation framework (ORForise) based on a comprehensive set of 12 primary and 60 secondary metrics that facilitate the assessment of the performance of CDS prediction tools. This makes it possible to identify which performs better for specific use-cases. We use this to assess 15 ab initio- and model-based tools representing those most widely used (historically and currently) to generate the knowledge in genomic databases. We find that the performance of any tool is dependent on the genome being analysed, and no individual tool ranked as the most accurate across all genomes or metrics analysed. Even the top-ranked tools produced conflicting gene collections, which could not be resolved by aggregation. The ORForise evaluation framework provides users with a replicable, data-led approach to make informed tool choices for novel genome annotations and for refining historical annotations., Availability and Implementation: Code and datasets for reproduction and customisation are available at https://github.com/NickJD/ORForise., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author(s) 2021. Published by Oxford University Press.)

Background: Schistosomiasis is a neglected tropical parasitic and chronic disease affecting hundreds of millions of people. Adult schistosomes reside in the blood stream of the definitive mammalian host. These helminth parasites possess two epithelial surfaces, the tegument and the gastrodermis, both of which interact with the host during immune evasion and in nutrient uptake., Methods: Female ARC Swiss mice (4-6 weeks old) were infected percutaneously with Schistosoma japonicum cercariae freshly shed from Oncomelania hupensis quadrasi snails (Philippines strain). Fluorescent in situ hybridisation (FISH) was performed by using fresh adult S. japonicum perfused from those infected mice. Adult S. japonicum worms were processed to isolate the tegument from the carcass containing the gastrodermis; blood and bile were collected individually from infected and uninfected mice. Total DNA extracted from all those samples were used for microbiome profiling., Results: FISH and microbiome profiling showed the presence of bacterial populations on two epithelial surfaces of adult worms, suggesting they were distinct not only from the host blood but also from each other. Whereas microbial diversity was reduced overall in the parasite epithelial tissues when compared with that of host blood, specific bacterial taxa, including Anoxybacillus and Escherichia, were elevated on the tegument. Minimal differences were evident in the microbiome of host blood during an active infection, compared with that of control uninfected blood. However, sampling of bile from infected animals identified some differences compared with controls, including elevated levels of Limnohabitans, Clostridium and Curvibacter., Conclusions: Using FISH and microbial profiling, we were able to demonstrate, for the first time, that bacteria are presented on the epithelial surfaces of adult schistosomes. These schistosome surface-associated bacteria, which are distinct from the host blood microenvironment, should be considered as a new and important component of the host-schistosome interaction. The importance of individual bacterial species in relation to schistosome parasitism needs further elucidation., Competing Interests: The authors have declared that no competing interests exist.

Persistent mating-induced endometritis is a major cause of poor fertility rates in the mare. Endometritis can be investigated using an ex vivo equine endometrial explant system which measures uterine inflammation using prostaglandin F 2α as a biomarker. However, this model has yet to undergo a wide-ranging assessment through transcriptomics. In this study, we assessed the transcriptomes of cultured endometrial explants and the optimal temporal window for their use. Endometrium harvested immediately post-mortem from native pony mares ( n = 8) were sampled (0 h) and tissue explants were cultured for 24, 48 and 72 h. Tissues were stored in RNALater, total RNA was extracted and sequenced. Differentially expressed genes (DEGs) were defined using DESeq2 (R/Bioconductor). Principal component analysis indicated that the greatest changes in expression occurred in the first 24 h of culture when compared to autologous biopsies at 0 h. Fewer DEGs were seen between 24 and 48 h of culture suggesting the system was more stable than during the first 24 h. No genes were differentially expressed between 48 and 72 h but the low number of background gene expression suggested that explant viability was compromised after 48 h. ESR1 , MMP9 , PTGS2 , PMAIP1 , TNF , GADD45B and SELE genes were used as biomarkers of endometrial function, cell death and inflammation across tissue culture timepoints. STRING assessments of gene ontology suggested that DEGs between 24 and 48 h were linked to inflammation, immune system, cellular processes, environmental information processing and signal transduction, with an upregulation of most biomarker genes at 24 h. Taken together our observations indicated that 24-48 h is the optimal temporal window when the explant model can be used, as explants restore microcirculation, perform wound healing and tackle inflammation during this period. This key observation will facilitate the appropriate use of this as a model for further research into the equine endometrium and potentially the progression of mating-induced endometritis to persistent inflammation between 24 and 48 h.

Background: Gut microbiomes, such as the rumen, greatly influence host nutrition due to their feed energy-harvesting capacity. We investigated temporal ecological interactions facilitating energy harvesting at the fresh perennial ryegrass (PRG)-biofilm interface in the rumen using an in sacco approach and prokaryotic metatranscriptomic profiling., Results: Network analysis identified two distinct sub-microbiomes primarily representing primary (≤ 4 h) and secondary (≥ 4 h) colonisation phases and the most transcriptionally active bacterial families (i.e Fibrobacteriaceae, Selemondaceae and Methanobacteriaceae) did not interact with either sub-microbiome, indicating non-cooperative behaviour. Conversely, Prevotellaceae had most transcriptional activity within the primary sub-microbiome (focussed on protein metabolism) and Lachnospiraceae within the secondary sub-microbiome (focussed on carbohydrate degradation). Putative keystone taxa, with low transcriptional activity, were identified within both sub-microbiomes, highlighting the important synergistic role of minor bacterial families; however, we hypothesise that they may be 'cheating' in order to capitalise on the energy-harvesting capacity of other microbes. In terms of chemical cues underlying transition from primary to secondary colonisation phases, we suggest that AI-2-based quorum sensing plays a role, based on LuxS gene expression data, coupled with changes in PRG chemistry., Conclusions: In summary, we show that fresh PRG-attached prokaryotes are resilient and adapt quickly to changing niches. This study provides the first major insight into the complex temporal ecological interactions occurring at the plant-biofilm interface within the rumen. The study also provides valuable insights into potential plant breeding strategies for development of the utopian plant, allowing optimal sustainable production of ruminants. Video Abstract.

Motivation: Population-level genetic variation enables competitiveness and niche specialization in microbial communities. Despite the difficulty in culturing many microbes from an environment, we can still study these communities by isolating and sequencing DNA directly from an environment (metagenomics). Recovering the genomic sequences of all isoforms of a given gene across all organisms in a metagenomic sample would aid evolutionary and ecological insights into microbial ecosystems with potential benefits for medicine and biotechnology. A significant obstacle to this goal arises from the lack of a computationally tractable solution that can recover these sequences from sequenced read fragments. This poses a problem analogous to reconstructing the two sequences that make up the genome of a diploid organism (i.e. haplotypes) but for an unknown number of individuals and haplotypes., Results: The problem of single individual haplotyping was first formalized by Lancia et al. in 2001. Now, nearly two decades later, we discuss the complexity of 'haplotyping' metagenomic samples, with a new formalization of Lancia et al.'s data structure that allows us to effectively extend the single individual haplotype problem to microbial communities. This work describes and formalizes the problem of recovering genes (and other genomic subsequences) from all individuals within a complex community sample, which we term the metagenomic individual haplotyping problem. We also provide software implementations for a pairwise single nucleotide variant (SNV) co-occurrence matrix and greedy graph traversal algorithm., Availability and Implementation: Our reference implementation of the described pairwise SNV matrix (Hansel) and greedy haplotype path traversal algorithm (Gretel) is open source, MIT licensed and freely available online at github.com/samstudio8/hansel and github.com/samstudio8/gretel, respectively., (© The Author(s) 2021. Published by Oxford University Press.)

Dietary lipids increase energy density in dairy cow diets and in some cases can increase beneficial fatty acids (FA) in milk and dairy products. However, the degree of FA saturation may affect the rumen microbiome. The objective of this study was to determine the long-term effects of feeding saturated (hydrogenated vegetable oil; HVO) or unsaturated (olive oil; OO) fatty acid (FA) sources on the rumen microbiome of dairy cows. For 63 days, 15 mid-lactating cows were fed with either a basal diet (no fat supplement), or the basal diet supplemented with 3% dry matter (DM), either HVO or OO. Rumen contents were collected on days 21, 42 and 63 for 16S rRNA gene sequencing using the Illumina MiSeq platform. The results reveal dominance of the phyla Firmicutes (71.5%) and Bacteroidetes (26.2%), and their respective prevalent genera Succiniclasticum (19.4%) and Prevotella (16.6%). Succiniclasticum increased with both treatments at all time points. Prevotella was reduced on day 42 in both diets. Bacterial diversity alpha or beta were not affected by diets. Predicted bacterial functions by CowPI showed changes in energy and protein metabolism. Overall, 3% DM of lipid supplementation over 63 days can be used in dairy cow diets without major impacts on global bacterial community structure.

Background: The rise of microbial antibiotic resistance is a leading threat to the health of the human population. As such, finding new approaches to tackle these microbes, including development of novel antibiotics is vital., Results: In this study, we mined a rumen eukaryotic metatranscriptomic library for novel Antimicrobial peptides (AMPs) using computational approaches and thereafter characterised the therapeutic potential of the AMPs. We identified a total of 208 potentially novel AMPs from the ruminal eukaryotome, and characterised one of those, namely Lubelisin. Lubelisin (GIVAWFWRLAR) is an α-helical peptide, 11 amino acid long with theoretical molecular weight of 1373.76 D. In the presence of Lubelisin, strains of methicillin-resistant Staphylococcus aureus (MRSA) USA300 and EMRSA-15 were killed within 30 min of exposure with ≥10 3 and 10 4  CFU/mL reduction in viable cells respectively. Cytotoxicity of Lubelisin against both human and sheep erythrocytes was low resulting in a therapeutic index of 0.43. Membrane permeabilisation assays using propidium iodide alongside transmission electron microscopy revealed that cytoplasmic membrane damage may contribute to the antimicrobial activities of Lubelisin., Conclusions: We demonstrate that the rumen eukaryotome is a viable source for the discovery of antimicrobial molecules for the treatment of bacterial infections and further development of these may provide part of the potential solution to the ongoing problem of antimicrobial resistance. The role of these AMPs in the ecological warfare within the rumen is also currently unknown.

Preimplantation factor (PIF) is an embryo derived peptide which exerts an immune modulatory effect on human endometrium, promoting immune tolerance to the embryo whilst maintaining the immune response to invading pathogens. While bovine embryos secrete PIF, the effect on the bovine endometrium is unknown. Maternal recognition of pregnancy is driven by an embryo-maternal cross talk, however the process differs between humans and cattle. As many embryos are lost during the early part of pregnancy in cattle, a greater knowledge of factors affecting the embryo-maternal crosstalk, such as PIF, is needed to improve fertility. Therefore, for the first time, we demonstrate the effect of synthetic PIF (sPIF) on the bovine transcriptome in an ex vivo bovine endometrial tissue culture model. Explants were cultured for 30h with sPIF (100nM) or in control media. Total RNA was analysed via RNA-sequencing. As a result of sPIF treatment, 102 genes were differentially expressed compared to the control (Padj<0.1), although none by more than 2-fold. The majority of genes (78) were downregulated. Pathway analysis revealed targeting of several immune based pathways. Genes for the TNF, NF-κB, IL-17, MAPK and TLR signalling pathways were down-regulated by sPIF. However, some immune genes were demonstrated to be upregulated following sPIF treatment, including C3. Steroid biosynthesis was the only over-represented pathway with all genes upregulated. We demonstrate that sPIF can modulate the bovine endometrial transcriptome in an immune modulatory manner, like that in the human endometrium, however, the regulation of genes was much weaker than in previous human work., Competing Interests: No authors have competing interests.

Background: Gene expression profiles can provide insights into the molecular machinery behind tissue functions and, in turn, can further our understanding of environmental responses, and developmental and evolutionary processes. During vertebrate evolution, the skin has played a crucial role, displaying a wide diversity of essential functions. To unravel the molecular basis of skin specialisations and adaptations, we compared gene expression in the skin with eight other tissues in a phylogenetically and ecologically diverse species sample of one of the most neglected vertebrate groups, the caecilian amphibians (order Gymnophiona)., Results: The skin of the five studied caecilian species showed a distinct gene expression profile reflecting its developmental origin and showing similarities to other epithelial tissues. We identified 59 sequences with conserved enhanced expression in the skin that might be associated with caecilian dermal specialisations. Some of the up-regulated genes shared expression patterns with human skin and potentially are involved in skin functions across vertebrates. Variation trends in gene expression were detected between mid and posterior body skin suggesting different functions between body regions. Several candidate biologically active peptides were also annotated., Conclusions: Our study provides the first atlas of differentially expressed sequences in caecilian tissues and a baseline to explore the molecular basis of the skin functions in caecilian amphibians, and more broadly in vertebrates.

Although the prokaryotic communities of the rumen microbiome are being uncovered through genome sequencing, little is known about the resident viral populations. Whilst temperate phages can be predicted as integrated prophages when analyzing bacterial and archaeal genomes, the genetics underpinning lytic phages remain poorly characterized. To the five genomes of bacteriophages isolated from rumen-associated samples sequenced and analyzed previously, this study adds a further five novel genomes and predictions gleaned from them to further the understanding of the rumen phage population. Lytic bacteriophages isolated from fresh ovine and bovine fecal and rumen fluid samples were active against the predominant fibrolytic ruminal bacterium Butyrivibrio fibrisolvens . The double stranded DNA genomes were sequenced and reconstructed into single circular complete contigs. Based on sequence similarity and genome distances, the five phages represent four species from three separate genera, consisting of: (1) Butyrivibrio phages Arian and Bo-Finn; (2) Butyrivibrio phages Idris and Arawn; and (3) Butyrivibrio phage Ceridwen. They were predicted to all belong to the Siphoviridae family, based on evidence in the genomes such as size, the presence of the tail morphogenesis module, genes that share similarity to those in other siphovirus isolates and phylogenetic analysis using phage proteomes. Yet, phylogenomic analysis and sequence similarity of the entire phage genomes revealed that these five phages are unique and novel. These phages have only been observed undergoing the lytic lifecycle, but there is evidence in the genomes of phages Arawn and Idris for the potential to be temperate. However, there is no evidence in the genome of the bacterial host Butyrivibrio fibrisolvens of prophage genes or genes that share similarity with the phage genomes., (Copyright © 2020 Friedersdorff, Kingston-Smith, Pachebat, Cookson, Rooke and Creevey.)

The rumen protozoa, alongside fungi, comprise the eukaryotic portion of the rumen microbiome. Rumen protozoa may account for up to 50% of biomass, yet their role in this ecosystem remains unclear. Early experiments inferred a role in carbohydrate and protein metabolism, but due to their close association with bacteria, definitively attributing these functions to the protozoa was challenging. The advent of 'omic technologies has created opportunities to broaden our understanding of the rumen protozoa. This study aimed to utilize these methods to further our understanding of the role that protozoa play in the rumen in terms of their metabolic capacities, and in doing so, contribute valuable sequence data to reduce the chance of mis or under-representation of the rumen protozoa in meta'omic datasets. Rumen protozoa were isolated and purified using glucose-based sedimentation and differential centrifugation, extracted RNA was Poly(A) fraction enriched and DNase treated before use in a phage-based, cDNA metatranscriptomic library. Biochemical activity testing of the phage library showed 6 putatively positive plaques in response to carboxymethyl cellulose agar (indicative of cellulose activity), and no positive results for tributyrin (indicative of esterase/lipase activity) or egg yolk agar (indicative of proteolysis). Direct sequencing of the cDNA was also conducted using the Illumina HiSeq 2500. The metatranscriptome identified a wealth of carbohydrate-active enzymes which accounted for 8% of total reads. The most highly expressed carbohydrate-active enzymes were glycosyl hydrolases 5 and 11, polysaccharide lyases and deacetylases, xylanases and enzymes active against pectin, mannan and chitin; the latter likely used to digest rumen fungi which contain a chitin-rich cell membrane. Codon usage analysis of expressed genes also showed evidence of horizontal gene transfer, suggesting that many of these enzymes were acquired from the rumen bacteria in an evolutionary response to the carbohydrate-rich environment of the rumen. This study provides evidence of the significant contribution that the protozoa make to carbohydrate breakdown in the rumen, potentially using horizontally acquired genes, and highlights their predatory capacity., (Copyright © 2020 Williams, Thomas, McEwan, Rees Stevens, Creevey and Huws.)

The rumen contains a multi-kingdom, commensal microbiome, including protozoa, bacteria, archaea, fungi and viruses, which enables ruminant herbivores to ferment and utilize plant feedstuffs that would be otherwise indigestible. Within the rumen, virus populations are diverse and highly abundant, often out-numbering the microbial populations that they both predate on and co-exist with. To date the research effort devoted to understanding rumen-associated viral populations has been considerably less than that given to the other microbial populations, yet their contribution to maintaining microbial population balance, intra-ruminal microbial lysis, fiber breakdown, nutrient cycling and genetic transfer may be highly significant. This review follows the technological advances which have contributed to our current understanding of rumen viruses and drawing on knowledge from other environmental and animal-associated microbiomes, describes the known and potential roles and impacts viruses have on rumen function and speculates on the future directions of rumen viral research., (Copyright © 2020 Gilbert, Townsend, Crew, Hitch, Friedersdorff, Creevey, Pope, Ouwerkerk and Jameson.)

Background: The rumen contains a myriad of microbes whose primary role is to degrade and ferment dietary nutrients, which then provide the host with energy and nutrients. Rumen microbes commonly attach to ingested plant materials and form biofilms for effective plant degradation. Quorum sensing (QS) is a well-recognised form of bacterial communication in most biofilm communities, with homoserine lactone (AHL)-based QS commonly being used by Gram-negative bacteria alone and AI-2 Lux-based QS communication being used to communicate across Gram-negative and Gram-positive bacteria. However, bacterial cell to cell communication in the rumen is poorly understood. In this study, rumen bacterial genomes from the Hungate collection and Genbank were prospected for QS-related genes. To check that the discovered QS genes are actually expressed in the rumen, we investigated expression levels in rumen metatranscriptome datasets., Results: A total of 448 rumen bacterial genomes from the Hungate collection and Genbank, comprised of 311 Gram-positive, 136 Gram-negative and 1 Gram stain variable bacterium, were analysed. Abundance and distribution of AHL and AI-2 signalling genes showed that only one species (Citrobacter sp. NLAE-zl-C269) of a Gram-negative bacteria appeared to possess an AHL synthase gene, while the Lux-based genes (AI-2 QS) were identified in both Gram-positive and Gram-positive bacteria (191 genomes representing 38.2% of total genomes). Of these 192 genomes, 139 are from Gram-positive bactreetteria and 53 from Gram-negative bacteria. We also found that the genera Butyrivibrio, Prevotella, Ruminococcus and Pseudobutyrivibrio, which are well known as the most abundant bacterial genera in the rumen, possessed the most lux-based AI-2 QS genes. Gene expression levels within the metatranscriptome dataset showed that Prevotella, in particular, expressed high levels of LuxS synthase suggesting that this genus plays an important role in QS within the rumen., Conclusion: This is the most comprehensive study of QS in the rumen microbiome to date. This study shows that AI-2-based QS is rife in the rumen. These results allow a greater understanding on plant-microbe interactions in the rumen.

Genomic and transcriptomic analyses were performed to investigate nonribosomal peptide synthetases (NRPS) and polyketide synthases (PKS) in 310 genomes of ruminal/fecal microorganisms. A total of 119 biosynthetic genes potentially encoding distinct nonribosomal peptides (NRPs) and polyketides (PKs) were predicted in the ruminal microbial genomes and functional annotation separated these genes into 19 functional categories. The phylogenetic reconstruction of the 16S rRNA sequences coupled to the distribution of the three 'backbone' genes involved in NRPS and PKS biosyntheses suggested that these genes were not acquired through horizontal gene transfer. Metatranscriptomic analyses revealed that the predominant genes involved in the synthesis of NRPs and PKs were more abundant in sheep rumen datasets. Reads mapping to the NRPS and PKS biosynthetic genes were represented in the active ruminal microbial community, with transcripts being highly expressed in the bacterial community attached to perennial ryegrass, and following the main changes occurring between primary and secondary colonization of the forage incubated with ruminal fluid. This study is the first comprehensive characterization demonstrating the rich genetic capacity for NRPS and PKS biosyntheses within rumen bacterial genomes, which highlights the potential functional roles of secondary metabolites in the rumen ecosystem., (© FEMS 2019.)

Most animal species reproduce sexually and fully parthenogenetic lineages are usually short lived in evolution. Still, parthenogenesis may be advantageous as it avoids the cost of sex and permits colonization by single individuals. Panagrolaimid nematodes have colonized environments ranging from arid deserts to Arctic and Antarctic biomes. Many are obligatory meiotic parthenogens, and most have cryptobiotic abilities, being able to survive repeated cycles of complete desiccation and freezing. To identify systems that may contribute to these striking abilities, we sequenced and compared the genomes and transcriptomes of parthenogenetic and outcrossing panagrolaimid species, including cryptobionts and non-cryptobionts. The parthenogens are triploids, most likely originating through hybridization. Adaptation to cryptobiosis shaped the genomes of panagrolaimid nematodes and is associated with the expansion of gene families and signatures of selection on genes involved in cryptobiosis. All panagrolaimids have acquired genes through horizontal gene transfer, some of which are likely to contribute to cryptobiosis., (Copyright © 2019 The Authors. Published by Elsevier Inc. All rights reserved.)

As previous studies have demonstrated a link between the porcine intestinal microbiome and feed efficiency (FE), microbiota manipulation may offer a means of improving FE in pigs. A fecal microbiota transplantation procedure (FMTp), using fecal extracts from highly feed-efficient pigs, was performed in pregnant sows ( n  = 11), with a control group ( n  = 11) receiving no FMTp. At weaning, offspring were allocated, within sow treatment, to (i) control ( n  = 67; no dietary supplement) or (ii) inulin ( n  = 65; 6-week dietary inulin supplementation) treatments. The sow FMTp, alone or in combination with inulin supplementation in offspring, reduced offspring body weight by 8.1 to 10.6 kg at ∼140 days of age, but there was no effect on feed intake. It resulted in better FE, greater bacterial diversity, and higher relative abundances of potentially beneficial bacterial taxa ( Fibrobacter and Prevotella ) in offspring. Due to the FMTp and/or inulin supplementation, relative abundances of potential pathogens ( Chlamydia and Treponema ) in the ileum and cecal concentrations of butyric acid were significantly lower. The maternal FMTp led to a greater number of jejunal goblet cells in offspring. Inulin supplementation alone did not affect growth or FE but upregulated duodenal genes linked to glucose and volatile fatty acid homeostasis and increased the mean platelet volume but reduced ileal propionic acid concentrations, granulocyte counts, and serum urea concentrations. Overall, the FMTp in pregnant sows, with or without dietary inulin supplementation in offspring, beneficially modulated offspring intestinal microbiota (albeit mostly low-relative-abundance taxa) and associated physiological parameters. Although FE was improved, the detrimental effect on growth limits the application of this FMTp-inulin strategy in commercial pig production. IMPORTANCE As previous research suggests a link between microbiota and FE, modulation of the intestinal microbiome may be effective in improving FE in pigs. The FMTp in gestating sows, alone or in combination with postweaning dietary inulin supplementation in offspring, achieved improvements in FE and resulted in a higher relative abundance of intestinal bacteria associated with fiber degradation and a lower relative abundance of potential pathogens. However, there was a detrimental effect on growth, although this may not be wholly attributable to microbiota transplantation, as antibiotic and other interventions were also part of the FMT regimen. Therefore, further work with additional control groups is needed to disentangle the effects of each component of the FMTp in order to develop a regimen with practical applications in pig production. Additional research based on findings from this study may also identify specific dietary supplements for the promotion/maintenance of the microbiota transferred via the maternal FMTp, thereby optimizing pig growth and FE., (Copyright © 2019 American Society for Microbiology.)

Gene fusion occurs when two or more individual genes with independent open reading frames becoming juxtaposed under the same open reading frame creating a new fused gene. A small number of gene fusions described in detail have been associated with novel functions, for example, the hominid-specific PIPSL gene, TNFSF12, and the TWE-PRIL gene family. We use Sequence Similarity Networks and species level comparisons of great ape genomes to identify 45 new genes that have emerged by transcriptional readthrough, that is, transcription-derived gene fusion. For 35 of these putative gene fusions, we have been able to assess available RNAseq data to determine whether there are reads that map to each breakpoint. A total of 29 of the putative gene fusions had annotated transcripts (9/29 of which are human-specific). We carried out RT-qPCR in a range of human tissues (placenta, lung, liver, brain, and testes) and found that 23 of the putative gene fusion events were expressed in at least one tissue. Examining the available ribosome foot-printing data, we find evidence for translation of three of the fused genes in human. Finally, we find enrichment for transcription-derived gene fusions in regions of known segmental duplication in human. Together, our results implicate chromosomal structural variation brought about by segmental duplication with the emergence of novel transcripts and translated protein products., (© The Author(s) 2019. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Increasingly, large phylogenomic data sets include transcriptomic data from nonmodel organisms. This not only has allowed controversial and unexplored evolutionary relationships in the tree of life to be addressed but also increases the risk of inadvertent inclusion of paralogs in the analysis. Although this may be expected to result in decreased phylogenetic support, it is not clear if it could also drive highly supported artifactual relationships. Many groups, including the hyperdiverse Lissamphibia, are especially susceptible to these issues due to ancient gene duplication events and small numbers of sequenced genomes and because transcriptomes are increasingly applied to resolve historically conflicting taxonomic hypotheses. We tested the potential impact of paralog inclusion on the topologies and timetree estimates of the Lissamphibia using published and de novo sequencing data including 18 amphibian species, from which 2,656 single-copy gene families were identified. A novel paralog filtering approach resulted in four differently curated data sets, which were used for phylogenetic reconstructions using Bayesian inference, maximum likelihood, and quartet-based supertrees. We found that paralogs drive strongly supported conflicting hypotheses within the Lissamphibia (Batrachia and Procera) and older divergence time estimates even within groups where no variation in topology was observed. All investigated methods, except Bayesian inference with the CAT-GTR model, were found to be sensitive to paralogs, but with filtering convergence to the same answer (Batrachia) was observed. This is the first large-scale study to address the impact of orthology selection using transcriptomic data and emphasizes the importance of quality over quantity particularly for understanding relationships of poorly sampled taxa., (© The Author(s) 2019. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Background: Evolution leaves an imprint in species through genetic change. At the molecular level, evolutionary changes can be explored by studying ratios of nucleotide substitutions. The interplay among molecular evolution, derived phenotypes, and ecological ranges can provide insights into adaptive radiations. Caecilians (order Gymnophiona), probably the least known of the major lineages of vertebrates, are limbless tropical amphibians, with adults of most species burrowing in soils (fossoriality). This enigmatic order of amphibians are very distinct phenotypically from other extant amphibians and likely from the ancestor of Lissamphibia, but little to nothing is known about the molecular changes underpinning their radiation. We hypothesised that colonization of various depths of tropical soils and of freshwater habitats presented new ecological opportunities to caecilians., Results: A total of 8540 candidate groups of orthologous genes from transcriptomic data of five species of caecilian amphibians and the genome of the frog Xenopus tropicalis were analysed in order to investigate the genetic machinery behind caecilian diversification. We found a total of 168 protein-coding genes with signatures of positive selection at different evolutionary times during the radiation of caecilians. The majority of these genes were related to functional elements of the cell membrane and extracellular matrix with expression in several different tissues. The first colonization of the tropical soils was connected to the largest number of protein-coding genes under positive selection in our analysis. From the results of our study, we highlighted molecular changes in genes involved in perception, reduction-oxidation processes, and aging that likely were involved in the adaptation to different soil strata., Conclusions: The genes inferred to have been under positive selection provide valuable insights into caecilian evolution, potentially underpin adaptations of caecilians to their extreme environments, and contribute to a better understanding of fossorial adaptations and molecular evolution in vertebrates.

RNA sequencing (RNA-seq) has become one of the most powerful tools to unravel the genomic basis of biological adaptation and diversity. Although challenging, RNA-seq is particularly promising for research on non-model, secretive species that cannot be observed in nature easily and therefore remain comparatively understudied. Among such animals, the caecilians (order Gymnophiona) likely constitute the least known group of vertebrates, despite being an old and remarkably distinct lineage of amphibians. Here, we characterize multi-tissue transcriptomes for five species of caecilians that represent a broad level of diversity across the order. We identified vertebrate homologous elements of caecilian functional genes of varying tissue specificity that reveal a great number of unclassified gene families, especially for the skin. We annotated several protein domains for those unknown candidate gene families to investigate their function. We also conducted supertree analyses of a phylogenomic dataset of 1,955 candidate orthologous genes among five caecilian species and other major lineages of vertebrates, with the inferred tree being in agreement with current views of vertebrate evolution and systematics. Our study provides insights into the evolution of vertebrate protein-coding genes, and a basis for future research on the molecular elements underlying the particular biology and adaptations of caecilian amphibians., (© The Author(s) 2018. Published by Oxford University Press on behalf of Kazusa DNA Research Institute.)

Background: Microorganisms are responsible for fermentation within the rumen and have been reported to contribute to the variation in feed efficiency of cattle. However, to what extent the breed affects the rumen microbiome and its association with host feed efficiency is unknown. Here, rumen microbiomes of beef cattle (n = 48) from three breeds (Angus, Charolais, Kinsella composite hybrid) with high and low feed efficiency were explored using metagenomics and metatranscriptomics, aiming to identify differences between functional potentials and activities of same rumen microbiomes and to evaluate the effects of host breed and feed efficiency on the rumen microbiome., Results: Rumen metagenomes were more closely clustered together and thus more conserved among individuals than metatranscriptomes, suggesting that inter-individual functional variations at the RNA level were higher than those at the DNA level. However, while mRNA enrichment significantly increased the sequencing depth of mRNA and generated similar functional profiles to total RNA-based metatranscriptomics, it led to biased abundance estimation of several transcripts. We observed divergent rumen microbial composition (metatranscriptomic level) and functional potentials (metagenomic level) among three breeds, but differences in functional activity (metatranscriptomic level) were less apparent. Differential rumen microbial features (e.g., taxa, diversity indices, functional categories, and genes) were detected between cattle with high and low feed efficiency, and most of them were breed-specific., Conclusions: Metatranscriptomes represent real-time functional activities of microbiomes and have the potential to better associate rumen microorganisms with host performances compared to metagenomics. As total RNA-based metatranscriptomics seem to avoid potential biases caused by mRNA enrichment and allow simultaneous use of rRNA for generation of compositional profiles, we suggest their use for linking the rumen microbiome with host phenotypes in future studies. However, if exploration of specific lowly expressed genes is desired, mRNA enrichment is recommended as it will enhance the resolution of mRNA. Finally, the differential microbial features observed between efficient and inefficient steers tended to be specific to breeds, suggesting that interactions between host breed genotype and the rumen microbiome contribute to the variations in feed efficiency observed. These breed-associated differences represent an opportunity to engineer specific rumen microbiomes through selective breeding of the hosts.

This work aimed to gain insight into the transition from milk to solid feeding at weaning combining genomics and metabolomics on rumen contents from goat kids treated with a methanogenic inhibitor (bromochloromethane, BCM). Sixteen goats giving birth to two kids were used. Eight does were treated (D+) with BCM after giving birth and over 2 months. One kid per doe in both groups was treated with BCM (k+) for 3 months while the other was untreated (k-). Rumen samples were collected from kids at weaning (W) and 1 (W + 1) and 4 (W + 4) months after and from does at weaning and subjected to 16S pyrosequencing and metabolomics analyses combining GC/LC-MS. Results from pyrosequencing showed a clear effect of age of kids, with more diverse bacterial community as solid feed becomes more important after weaning. A number of specific OTUs were significantly different as a result of BCM treatment of the kid at W while at W + 1 and W + 4 less OTUs were significantly changed. At W + 1, Prevotella was increased and Butyrivibrio decreased in BCM treated kids. At W + 4 only the effect of treating mothers resulted in significant changes in the abundance of some OTUs: Ruminococcus, Butyrivibrio and Prevotella . The analysis of the OTUs shared by different treatments revealed that kids at weaning had the largest number of unique OTUs compared with kids at W + 1 (137), W + 4 (238), and does (D) (23). D + k+ kids consistently shared more OTUs with mothers than the other three groups at the three sampling times. The metalobomic study identified 473 different metabolites. In does, lipid super pathway included the highest number of metabolites that were modified by BCM, while in kids all super-pathways were evenly affected. The metabolomic profile of samples from kids at W was different in composition as compared to W + 1 and W + 4, which may be directly ascribed to the process of rumen maturation and changes in the solid diet. This study shows the complexity of the bacterial community and metabolome in the rumen before weaning, which clearly differ from that after weaning and highlight the importance of the dam in transmitting the primary bacterial community after birth.

The rumen is a complex ecosystem composed of anaerobic bacteria, protozoa, fungi, methanogenic archaea and phages. These microbes interact closely to breakdown plant material that cannot be digested by humans, whilst providing metabolic energy to the host and, in the case of archaea, producing methane. Consequently, ruminants produce meat and milk, which are rich in high-quality protein, vitamins and minerals, and therefore contribute to food security. As the world population is predicted to reach approximately 9.7 billion by 2050, an increase in ruminant production to satisfy global protein demand is necessary, despite limited land availability, and whilst ensuring environmental impact is minimized. Although challenging, these goals can be met, but depend on our understanding of the rumen microbiome. Attempts to manipulate the rumen microbiome to benefit global agricultural challenges have been ongoing for decades with limited success, mostly due to the lack of a detailed understanding of this microbiome and our limited ability to culture most of these microbes outside the rumen. The potential to manipulate the rumen microbiome and meet global livestock challenges through animal breeding and introduction of dietary interventions during early life have recently emerged as promising new technologies. Our inability to phenotype ruminants in a high-throughput manner has also hampered progress, although the recent increase in "omic" data may allow further development of mathematical models and rumen microbial gene biomarkers as proxies. Advances in computational tools, high-throughput sequencing technologies and cultivation-independent "omics" approaches continue to revolutionize our understanding of the rumen microbiome. This will ultimately provide the knowledge framework needed to solve current and future ruminant livestock challenges.

The rumen is a complex ecosystem. It is the primary site for microbial fermentation of ingested feed allowing conversion of a low nutritional feed source into high quality meat and milk products. However, digestive inefficiencies lead to production of high amounts of environmental pollutants; methane and nitrogenous waste. These inefficiencies could be overcome by development of forages which better match the requirements of the rumen microbial population. Although challenging, the application of meta-proteomics has potential for a more complete understanding of the rumen ecosystem than sequencing approaches alone. Here, we have implemented a meta-proteomic approach to determine the association between taxonomies of microbial sources of the most abundant proteins in the rumens of forage-fed dairy cows, with taxonomic abundances typical of those previously described by metagenomics. Reproducible proteome profiles were generated from rumen samples. The most highly abundant taxonomic phyla in the proteome were Bacteriodetes, Firmicutes and Proteobacteria, which corresponded with the most abundant taxonomic phyla determined from 16S rRNA studies. Meta-proteome data indicated differentiation between metabolic pathways of the most abundant phyla, which is in agreement with the concept of diversified niches within the rumen microbiota.

Metataxonomic 16S rDNA based studies are a commonplace and useful tool in the research of the microbiome, but they do not provide the full investigative power of metagenomics and metatranscriptomics for revealing the functional potential of microbial communities. However, the use of metagenomic and metatranscriptomic technologies is hindered by high costs and skills barrier necessary to generate and interpret the data. To address this, a tool for Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt) was developed for inferring the functional potential of an observed microbiome profile, based on 16S data. This allows functional inferences to be made from metataxonomic 16S rDNA studies with little extra work or cost, but its accuracy relies on the availability of completely sequenced genomes of representative organisms from the community being investigated. The rumen microbiome is an example of a community traditionally underrepresented in genome and sequence databases, but recent efforts by projects such as the Global Rumen Census and Hungate 1000 have resulted in a wide sampling of 16S rDNA profiles and almost 500 fully sequenced microbial genomes from this environment. Using this information, we have developed "CowPI," a focused version of the PICRUSt tool provided for use by the wider scientific community in the study of the rumen microbiome. We evaluated the accuracy of CowPI and PICRUSt using two 16S datasets from the rumen microbiome: one generated from rDNA and the other from rRNA where corresponding metagenomic and metatranscriptomic data was also available. We show that the functional profiles predicted by CowPI better match estimates for both the meta-genomic and transcriptomic datasets than PICRUSt, and capture the higher degree of genetic variation and larger pangenomes of rumen organisms. Nonetheless, whilst being closer in terms of predictive power for the rumen microbiome, there were differences when compared to both the metagenomic and metatranscriptome data and so we recommend, where possible, functional inferences from 16S data should not replace metagenomic and metatranscriptomic approaches. The tool can be accessed at http://www.cowpi.org and is provided to the wider scientific community for use in the study of the rumen microbiome.

Antibiotic resistance is an increasingly important environmental pollutant with direct consequences for human health. Identification of environmental sources of antibiotic resistance genes (ARGs) makes it possible to follow their evolution and prevent their entry into the clinical setting. ARGs have been found in environmental sources exogenous to the original source and previous studies have shown that these genes are capable of being transferred from livestock to humans. Due to the nature of farming and the slaughter of ruminants for food, humans interact with these animals in close proximity, and for this reason it is important to consider the risks to human health. In this study, we characterised the ARG populations in the ovine rumen, termed the resistome. This was done using the Comprehensive Antibiotic Resistance Database (CARD) to identify the presence of genes conferring resistance to antibiotics within the rumen. Genes were successfully mapped to those that confer resistance to a total of 30 different antibiotics. Daptomycin was identified as the most common antibiotic for which resistance is present, suggesting that ruminants may be a source of daptomycin ARGs. Colistin resistance, conferred by the gene pmrE, was also found to be present within all samples, with an average abundance of 800 counts. Due to the high abundance of some ARGs (against daptomycin) and the presence of rare ARGs (against colistin), we suggest further study and monitoring of the rumen resistome as a possible source of clinically relevant ARGs., (Copyright © 2018 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Productivity of ruminant livestock depends on the rumen microbiota, which ferment indigestible plant polysaccharides into nutrients used for growth. Understanding the functions carried out by the rumen microbiota is important for reducing greenhouse gas production by ruminants and for developing biofuels from lignocellulose. We present 410 cultured bacteria and archaea, together with their reference genomes, representing every cultivated rumen-associated archaeal and bacterial family. We evaluate polysaccharide degradation, short-chain fatty acid production and methanogenesis pathways, and assign specific taxa to functions. A total of 336 organisms were present in available rumen metagenomic data sets, and 134 were present in human gut microbiome data sets. Comparison with the human microbiome revealed rumen-specific enrichment for genes encoding de novo synthesis of vitamin B 12 , ongoing evolution by gene loss and potential vertical inheritance of the rumen microbiome based on underrepresentation of markers of environmental stress. We estimate that our Hungate genome resource represents ∼75% of the genus-level bacterial and archaeal taxa present in the rumen.

Previous studies suggest a link between intestinal microbiota and porcine feed efficiency (FE). Therefore, we investigated whether fecal microbiota transplantation (FMT) in sows and/or neonatal offspring, using inocula derived from highly feed-efficient pigs, could improve offspring FE. Pregnant sows were assigned to control or FMT treatments and the subsequent offspring to control treatment, FMT once (at birth), or FMT four times (between birth and weaning). FMT altered sow fecal and colostrum microbiota compositions and resulted in lighter offspring body weight at 70 and 155 days of age when administered to sows and/or offspring. This was accompanied by FMT-associated changes within the offspring's intestinal microbiota, mostly in the ileum. These included transiently higher fecal bacterial diversity and load and numerous compositional differences at the phylum and genus levels (e.g., Spirochaetes and Bacteroidetes at high relative abundances and mostly members of Clostridia , respectively), as well as differences in the abundances of predicted bacterial pathways. In addition, intestinal morphology was negatively impacted, duodenal gene expression altered, and serum protein and cholesterol concentrations reduced due to FMT in sows and/or offspring. Taken together, the results suggest poorer absorptive capacity and intestinal health, most likely explaining the reduced body weight. An additive effect of FMT in sows and offspring also occurred for some parameters. Although these findings have negative implications for the practical use of the FMT regime used here for improving FE in pigs, they nonetheless demonstrate the enormous impact of early-life intestinal microbiota on the host phenotype. IMPORTANCE Here, for the first time, we investigate FMT as a novel strategy to modulate the porcine intestinal microbiota in an attempt to improve FE in pigs. However, reprogramming the maternal and/or offspring microbiome by using fecal transplants derived from highly feed-efficient pigs did not recapitulate the highly efficient phenotype in the offspring and, in fact, had detrimental effects on lifetime growth. Although these findings may not be wholly attributable to microbiota transplantation, as antibiotic and purgative were also part of the regime in sows, similar effects were also seen in offspring, in which these interventions were not used. Nonetheless, additional work is needed to unravel the effects of each component of the FMT regime and to provide additional mechanistic insights. This may lead to the development of an FMT procedure with practical applications for the improvement of FE in pigs, which could in turn improve the profitability of pig production.

Background: The consensus emerging from the study of microbiomes is that they are far more complex than previously thought, requiring better assemblies and increasingly deeper sequencing. However, current metagenomic assembly techniques regularly fail to incorporate all, or even the majority in some cases, of the sequence information generated for many microbiomes, negating this effort. This can especially bias the information gathered and the perceived importance of the minor taxa in a microbiome., Results: We propose a simple but effective approach, implemented in Python, to address this problem. Based on an iterative methodology, our workflow (called Spherical) carries out successive rounds of assemblies with the sequencing reads not yet utilised. This approach also allows the user to reduce the resources required for very large datasets, by assembling random subsets of the whole in a "divide and conquer" manner., Conclusions: We demonstrate the accuracy of Spherical using simulated data based on completely sequenced genomes and the effectiveness of the workflow at retrieving lost information for taxa in three published metagenomics studies of varying sizes. Our results show that Spherical increased the amount of reads utilized in the assembly by up to 109% compared to the base assembly. The additional contigs assembled by the Spherical workflow resulted in a significant (P < 0.05) changes in the predicted taxonomic profile of all datasets analysed. Spherical is implemented in Python 2.7 and freely available for use under the MIT license. Source code and documentation is hosted publically at: https://github.com/thh32/Spherical .

Antimicrobial peptides (AMPs) are promising drug candidates to target multi-drug resistant bacteria. The rumen microbiome presents an underexplored resource for the discovery of novel microbial enzymes and metabolites, including AMPs. Using functional screening and computational approaches, we identified 181 potentially novel AMPs from a rumen bacterial metagenome. Here, we show that three of the selected AMPs (Lynronne-1, Lynronne-2 and Lynronne-3) were effective against numerous bacterial pathogens, including methicillin-resistant Staphylococcus aureus (MRSA). No decrease in MRSA susceptibility was observed after 25 days of sub-lethal exposure to these AMPs. The AMPs bound preferentially to bacterial membrane lipids and induced membrane permeability leading to cytoplasmic leakage. Topical administration of Lynronne-1 (10% w/v) to a mouse model of MRSA wound infection elicited a significant reduction in bacterial counts, which was comparable to treatment with 2% mupirocin ointment. Our findings indicate that the rumen microbiome may provide viable alternative antimicrobials for future therapeutic application.

Background: Negative energy balance (NEB) is an imbalance between energy intake and energy requirements for lactation and body maintenance affecting high-yielding dairy cows and is of considerable economic importance due to its negative impact on fertility and health in dairy herds. It is anticipated that the cow hypothalamus experiences extensive biochemical changes during the early post partum period in an effort to re-establish metabolic homeostasis. However, there is variation in the tolerance to NEB between individual cows. In order to understand the genomic regulation of ovulation in hypothalamic tissue during NEB, mRNA transcriptional patterns between tolerant and sensitive animals were examined. A short term dietary restriction heifer model was developed which induced abrupt onset of anoestrus in some animals (Restricted Anovulatory; RA) while others maintained oestrous cyclicity (Restricted Ovulatory; RO). A third control group (C) received a higher level of normal feeding., Results: A total of 15,295 genes were expressed in hypothalamic tissue. Between RA and C groups 137 genes were differentially expressed, whereas between RO and C, 32 genes were differentially expressed. Differentially expressed genes were involved in the immune response and cellular motility in RA and RO groups, respectively, compared to C group. The largest difference between groups was observed in the comparison between RA and RO heifers, with 1094 genes shown to be significantly differentially expressed (SDE). Pathway analysis showed that these SDE genes were associated with 6 canonical pathways (P < 0.01), of which neuroactive ligand-receptor interaction was the most significant. Within the comparisons the main over-represented pathway functions were immune response including neuroprotection (CXCL10, Q1KLR3, IFIH1, IL1 and IL8; RA v C and RA v RO); energy homeostasis (AgRP and NPY; RA v RO); cell motility (CADH1, DSP and TSP4; RO v C) and prevention of GnRH release (NTSR1 IL1α, IL1β, NPY and PACA; RA v RO)., Conclusions: This information will assist in understanding the genomic factors regulating the influence of diet restriction on fertility and may assist in optimising nutritional and management systems for the improvement in reproductive performance.

Twelve mature (aged 5-16 years) horses and ponies of mixed breed and type were fed restricted (1.25% BM Dry matter) quantities of one of two fiber based diets formulated to be iso-caloric. Diet 1 comprised of 0.8% body mass (BM) of chaff based complete feed plus 0.45% BM low energy grass hay (the same hay used for both diets). Diet 2 comprised 0.1% BM of a nutrient balancer plus 1.15% BM grass hay. Fecal samples were collected at week 10 and week 16. DNA was extracted and the V1-V2 regions of 16SrDNA were 454-pyrosequenced to investigate the bacterial microbiome of the horse. The two most abundant phyla found in both diets and sampling periods were the Firmicutes and Bacteroidetes. There was a clear reduction in Bacteroidetes with a concordant increase in Firmicutes over time. There was a limited degree of stability within the bacterial community of the hindgut of horses, with 65% of bacteria retained, over a 6 week period whilst on a uniform diet. The presence of a core community defined by being present in all samples (each animal/diet combination) included in the study and being present at 0.1% relative abundance (or greater) was identified. In total 65 operational taxonomic units (OTUs) were identified that fit the definition of core making up 21-28% of the total sequences recovered. As with total population the most abundant phyla were the Bacteroidetes followed by the Firmicutes , however there was no obvious shift in phyla due to period. Indeed, when the relative abundance of OTUs was examined across diets and periods there was no significant effect of diet or period alone or in combination on the relative abundance of the core OTUs.