Intergenomic evolution and metabolic cross-talk between rumen and thermophilic autotrophic methanogenic archaea

Discovery of novel virulence mechanisms in Clostridium botulinum type A3 using genome-wide analysis

OBJECTIVE

Clostridium botulinum type A is a neurotoxin-producing, spore-forming anaerobic bacterium that causes botulism in humans. The evolutionary genomic context of this organism is not yet known to understand its molecular virulence mechanisms in the human intestinal tract. Hence, this study aimed to investigate the mechanisms underlying virulence and pathogenesis by comparing the genomic contexts across species, serotypes, and subtypes.

METHODS

A comparative genomic approach was used to analyze evolutionary genomic relationships, intergenomic distances, syntenic blocks, replication origins, and gene abundance with phylogenomic neighbors.

RESULTS

Type A strains have shown genomic proximity to group I strains with distinct accessory genes and vary even within subtypes. Phylogenomic data showed that type C and D strains were distantly related to a group I and group II strains. Synthetic plots indicated that orthologous genes might have evolved from Clostridial ancestry to subtype A3 strains, whereas syntonic out-paralogs might have emerged between subtypes A3 and A1 through α-events. Gene abundance analysis revealed the key roles of genes involved in biofilm formation, cell-cell communication, human diseases, and drug resistance compared to the pathogenic Clostridia. Moreover, we identified 43 unique genes in the type A3 genome, of which 29 were involved in the pathophysiological processes and other genes contributed to amino acid metabolism. The C. botulinum type A3 genome contains 14 new virulence proteins that can provide the ability to confer antibiotic resistance, virulence exertion and adherence to host cells, the host immune system, and mobility of extrachromosomal genetic elements.

CONCLUSION

The results of our study provide insight into the understanding of new virulence mechanisms to discover new therapeutics for the treatment of human diseases caused by type A3 strains.

Metagenomic and Metabolomic Insights Into the Mechanism Underlying the Disparity in Milk Yield of Holstein Cows

This study was conducted to investigate the metabolic mechanism underlying the disparity in the milk yield of Holstein cows. Eighteen lactating Holstein cows in their second parity and 56 (±14.81 SD) days in milking (DIM) were selected from 94 cows. Based on the milk yield of the cows, they were divided into two groups of nine cows each, the high milk yield group (HP) (44.57 ± 2.11 kg/day) and the low milk yield group (LP) (26.71 ± 0.70 kg/day). The experimental cows were fed the same diet and kept under the same management system for more than 60 days. Rumen metagenomics revealed that two Archaea genera, one Bacteria genus, eight Eukaryota genera, and two Virus genera differ between the HP and LP groups. The analysis of metabolites in the rumen fluid, milk, and serum showed that several metabolites differed between the HP and LP groups. Correlation analysis between the predominant microbiota and milk yield-associated metabolites (MP-metabolites) revealed that four Bacteria and two Eukaryota genera have a positive relationship with MP-metabolites. Pathway enrichment analysis of the differential metabolites revealed that five pathways were enriched in all the samples (two pathways in the milk, two pathways in the serum, and one pathway in the rumen fluid). Further investigation revealed that the low milk yield observed in the LP group might be due to an upregulation in dopamine levels in the rumen fluid and milk, which could inhibit the release of prolactin or suppress the action of oxytocin in the udder resulting in reduced milk yield. On the other hand, the high milk yield in the HP group is attributed to an upregulation in citrulline, and N-acetylornithine, which could be used as substrates for energy metabolism in the citric acid cycle and ultimately gluconeogenesis.

Comparative genomic analysis of hyper-ammonia producing Acetoanaerobium sticklandii DSM 519 with purinolytic Gottschalkia acidurici 9a and pathogenic Peptoclostridium difficile 630

Acetoanaerobium sticklandii DSM519 (CST) is a hype-ammonia producing non-pathogenic anaerobe that can use amino acids as important carbon and energy sources through the Stickland reactions. Biochemical aspects of this organism have been extensively studied, but systematic studies addressing its metabolic discrepancy remain scant. In this perspective, we have intensively analyzed its genomic and metabolic characteristics to comprehend the evolutionary conservation of amino acid catabolism by a comparative genomic approach. The whole-genome data indicated that CST has shown a phylogenomic similarity with hyper-ammonia producing, purinolytic, and proteolytic pathogenic Clostridia. CST has shown to common genomic context sharing across the purinolytic Gottschalkia acidurici 9a and pathogenic Peptoclostridium difficile 630. Genome syntenic analysis described that syntenic orthologs might be originated from the recent ancestor at a slow evolution rate and syntenic-out paralogs evolved from either CDF or CAC via α-event and β-event. Collinearity of either gene orders or gene families was adjusted with syntenic out-paralogs across these genomes. The genome-wide metabolic analysis predicted 11 unique putative metabolic subsystems from the CST genome for amino acid catabolism and hydrogen production. The in silico analysis of our study revealed that a characteristic system for amino acid catabolism-directed biofuel synthesis might have slowly evolved and established as a core genomic content of CST.

Vertebrate host phylogeny influences gut archaeal diversity

Commonly used 16S rRNA gene primers do not detect the full range of archaeal diversity present in the vertebrate gut. As a result, several questions regarding the archaeal component of the gut microbiota remain, including which Archaea are host-associated, the specificities of such associations and the major factors influencing archaeal diversity. Using 16S rRNA gene amplicon sequencing with primers that specifically target Archaea, we obtained sufficient sequence data from 185 gastrointestinal samples collected from 110 vertebrate species that span five taxonomic classes (Mammalia, Aves, Reptilia, Amphibia and Actinopterygii), of which the majority were wild. We provide evidence for previously undescribed Archaea-host associations, including Bathyarchaeia and Methanothermobacter, the latter of which was prevalent among Aves and relatively abundant in species with higher body temperatures, although this association could not be decoupled from host phylogeny. Host phylogeny explained archaeal diversity more strongly than diet, while specific taxa were associated with both factors, and cophylogeny was significant and strongest for mammalian herbivores. Methanobacteria was the only class predicted to be present in the last common ancestors of mammals and all host species. Further analysis indicated that Archaea-Bacteria interactions have a limited effect on archaeal diversity. These findings expand our current understanding of Archaea-vertebrate associations.

Comparative genomic analysis reveals starvation survival systems in Methanothermobacter thermautotrophicus ΔH

Methanothermobacter thermautotrophicus ΔH (MTH) is a thermophilic hydrogenotrophic methanogenic archaeon capable of reducing CO₂ with H₂ to produce methane gas. It is the potential candidate in the biomethanation of CO₂ and CO in anaerobic reactors and biogas upgrading process. However, systematic studies addressing its genome conservation and function remain scant in this genome. In this study, we have evaluated its evolutionary resemblance and metabolic discrepancy, particularly in starvation survival systems by comparing the genomic contexts with Methanothermobacter marburgensis str. Marburg (MMG) and Methanobacterium formicicum DSM 1535 (MFO). The phylogenomic analysis of this study indicated that there was a strong phylogenomic signal among MTH, MMG, and MFO in the whole-genome tree. DNA replication machinery was conserved in the MTH genome and might have evolved at different evolution rates. Genome synteny analysis observed collinearity of either gene orders or gene families has to be maintained with syntenic blocks located in the syntenic out-paralogs. A genome-wide metabolic analysis identified some unique putative metabolic subsystems in MTH, which are proposed to determine its growth characteristics in diverse environments. MTH genome comprised of 93 unique genes-coding for starvation survival and stress-response proteins. These proteins confer its adaptation to nutritional deprivation and other abiotic stresses. MTH has a typical system to withstand its growth and cell viability during stable operation and recovery after prolonged starvation. Thus, the present work will provide an insight to improve the genome refinement and metabolic reconstruction in parallel to other closely related species.

Comparative analysis of differential proteome-wide protein-protein interaction network of Methanobrevibacter ruminantium M1

A proteome-wide protein-protein interaction (PPI) network of Methanobrevibacter ruminantium M1 (MRU), a predominant rumen methanogen, was constructed from its metabolic genes using a gene neighborhood algorithm and then compared with closely related rumen methanogens Using proteome-wide PPI approach, we constructed network encompassed 2194 edges and 637 nodes interacting with 634 genes. Network quality and robustness of functional modules were assessed with gene ontology terms. A structure-function-metabolism mapping for each protein has been carried out with efforts to extract experimental PPI concomitant information from the literature. The results of our study revealed that some topological properties of its network were robust for sharing homologous protein interactions across heterotrophic and hydrogenotrophic methanogens. MRU proteome has shown to establish many PPI sub-networks for associated metabolic subsystems required to survive in the rumen environment. MRU genome found to share interacting proteins from its PPI network involved in specific metabolic subsystems distinct to heterotrophic and hydrogenotrophic methanogens. Across these proteomes, the interacting proteins from differential PPI networks were shared in common for the biosynthesis of amino acids, nucleosides, and nucleotides and energy metabolism in which more fractions of protein pairs shared with Methanosarcina acetivorans. Our comparative study expedites our knowledge to understand a complex proteome network associated with typical metabolic subsystems of MRU and to improve its genome-scale reconstruction in the future.

Performance and microbial community profiles in pilot-scale biofilter for the simultaneous removal of ammonia, iron and manganese at different manganese concentrations

To accelerate extensive application of biological manganese removal technology, a pilot-scale biofilter for ammonia, iron and manganese removal was constructed to investigate the removal performance and microbial community profiles at different manganese concentrations. When manganese in influent increased from 1 to 10 mg/L, the pollutants were completely removed. Ammonia and iron was slightly changed along the filter depth, while manganese obviously increased. In 0 m of the filter depth, the abundance of Gallionella (iron oxidizing bacteria, IOB) increased, while Crenothrix (IOB) decreased. The abundance of Gallionella (manganese oxidizing bacteria, MnOB) in 0.4 and 0.8 m increased to 16.82% and 12.37%, respectively; and Crenothrix (MnOB) in 0.8 m increased to 19.95%, but decreased to 25.08% in 0.4 m. The abundance of ammonia oxidizing bacteria (AOB, Nitrosococcus) decreased in 0.4 and 0.8 m. The biofilter presented a high ability to remove manganese, and had a broad application prospect.

Functional annotation of operome from Methanothermobacter thermautotrophicus ΔH: An insight to metabolic gap filling

Methanothermobacter thermautotrophicus ΔH (MTH) is a potential methanogen known to reduce CO₂ with H₂ for producing methane biofuel in thermophilic digesters. The genome of this organism contains ~50.5% conserved hypothetical proteins (HPs; operome) whose function is still not determined precisely. Here, we employed a combined bioinformatics approach to annotate a precise function to HPs and categorize them as enzymes, binding proteins, and transport proteins. Results of our study show that 315 (35.6%) HPs have exhibited well-defined functions contributing imperative roles in diverse cellular metabolism. Some of them are responsible for stress-response mechanisms and cell cycle, membrane transport, and regulatory processes. The genome-neighborhood analysis found five important gene clusters (dsr, ehb, kaiC, cmr, and gas) involving in the energetic metabolism and defense systems. MTH operome contains 223 enzymes with 15 metabolic subsystems, 15 cell cycle proteins, 17 transcriptional regulators and 33 binding proteins. Functional annotation of its operome is thus more fundamental to a profound understanding of the molecular and cellular machinery at systems-level.

Methanobacterium formicicum as a target rumen methanogen for the development of new methane mitigation interventions: A review

Methanobacterium formicicum (Methanobacteriaceae family) is an endosymbiotic methanogenic Archaean found in the digestive tracts of ruminants and elsewhere. It has been significantly implicated in global CH₄ emission during enteric fermentation processes. In this review, we discuss current genomic and metabolic aspects of this microorganism for the purpose of the discovery of novel veterinary therapeutics. This microorganism encompasses a typical H₂ scavenging system, which facilitates a metabolic symbiosis across the H₂ producing cellulolytic bacteria and fumarate reducing bacteria. To date, five genome-scale metabolic models (iAF692, iMG746, iMB745, iVS941 and iMM518) have been developed. These metabolic reconstructions revealed the cellular and metabolic behaviors of methanogenic archaea. The characteristics of its symbiotic behavior and metabolic crosstalk with competitive rumen anaerobes support understanding of the physiological function and metabolic fate of shared metabolites in the rumen ecosystem. Thus, systems biological characterization of this microorganism may provide a new insight to realize its metabolic significance for the development of a healthy microbiota in ruminants. An in-depth knowledge of this microorganism may allow us to ensure a long term sustainability of ruminant-based agriculture.

Recent development of Ori-Finder system and DoriC database for microbial replication origins

DNA replication begins at replication origins in all three domains of life. Identification and characterization of replication origins are important not only in providing insights into the structure and function of the replication origins but also in understanding the regulatory mechanisms of the initiation step in DNA replication. The Z-curve method has been used in the identification of replication origins in archaeal genomes successfully since 2002. Furthermore, the Web servers of Ori-Finder and Ori-Finder 2 have been developed to predict replication origins in both bacterial and archaeal genomes based on the Z-curve method, and the replication origins with manual curation have been collected into an online database, DoriC. Ori-Finder system and DoriC database are currently used in the research field of DNA replication origins in prokaryotes, including: (i) identification of oriC regions in bacterial and archaeal genomes; (ii) discovery and analysis of the conserved sequences within oriC regions; and (iii) strand-biased analysis of bacterial genomes.

Up to now, more and more predicted results by Ori-Finder system were supported by subsequent experiments, and Ori-Finder system has been used to identify the replication origins in > 100 newly sequenced prokaryotes in their genome reports. In addition, the data in DoriC database have been widely used in the large-scale analyses of replication origins and strand bias in prokaryotic genomes. Here, we review the development of Ori-Finder system and DoriC database as well as their applications. Some future directions and aspects for extending the application of Ori-Finder and DoriC are also presented.

Distribution and genetic diversity of microbial populations in the pilot-scale biofilter for simultaneous removal of ammonia, iron and manganese from real groundwater

A pilot-scale biofilter treating real groundwater was developed in this study, which showed that ammonia, iron and manganese were mainly removed at 0.4, 0.4 and 0.8 m of the filter bed, respectively, and the corresponding removal efficiencies were 90.82%, 95.48% and 95.90% in steady phase, respectively. The variation of microbial populations in the biofilter during start-up process was also investigated using high-throughput pyrosequencing (HTP). Results indicated that the main functional microbes for ammonia, iron and manganese removal were Nitrosomonas, Crenothrix and Crenothrix, respectively, which was mainly distributed at 0.8, 0, and 0.8 m of the filter bed with a corresponding abundance of 8.7%, 28.12% and 11.33% in steady phase, respectively. Kinds of other bacteria which may be related to methane, hydrogen sulfide and organic matter removal, were also found. In addition, small part of archaea was also detected, such as Candidatus Nitrososphaera, which plays a role in nitritation.