Search

Abstract Isomalto/malto-polysaccharides (IMMPs) are a novel type of soluble dietary fibres with a prebiotic potential promoting growth of beneficial microbes in the gut. However, the mode of action of IMMPs remains unknown. Previous studies on IMMPs showed an increase in total bacteria, especially lactobacilli, and higher production of short chain fatty acids (SCFA) when IMMPs were fed to rats or used during in vitro fermentation. Here we used metatranscriptomics to investigate how IMMPs with different amounts of α − (1 → 6) glycosidic linkages affected microbial function during incubation with human fecal inoculum. We showed that active microbial community dynamics during fermentation varied depending on the type of IMMP used and that the observed changes were reflected in the community gene expression profiles. Based on metatranscriptome analysis, members of Bacteroides, Lactobacillus and Bifidobacterium were the predominant degraders of IMMPs, and the increased gene expression in these bacteria correlated with high amounts of α − (1 → 6) glycosidic linkages. We also noted an increase in relative abundance of these bacteria and an activation of pathways involved in SCFA synthesis. Our findings could provide a baseline for more targeted approaches in designing prebiotics for specific bacteria and to achieve more controlled modulation of microbial activity towards desired health outcomes.

Abstract The rhizosphere, the region of soil surrounding roots of plants, is colonized by a unique population of Plant Growth Promoting Rhizobacteria (PGPR). Many important PGPR as well as plant pathogens belong to the genus Pseudomonas. There is, however, uncertainty on the divide between beneficial and pathogenic strains as previously thought to be signifying genomic features have limited power to separate these strains. Here we used the Genome properties (GP) common biological pathways annotation system and Machine Learning (ML) to establish the relationship between the genome wide GP composition and the plant-associated lifestyle of 91 Pseudomonas strains isolated from the rhizosphere and the phyllosphere representing both plant-associated phenotypes. GP enrichment analysis, Random Forest model fitting and feature selection revealed 28 discriminating features. A test set of 75 new strains confirmed the importance of the selected features for classification. The results suggest that GP annotations provide a promising computational tool to better classify the plant-associated lifestyle.

Abstract Background Microbial production of propionate from diluted streams of ethanol (e.g., deriving from syngas fermentation) is a sustainable alternative to the petrochemical production route. Yet, few ethanol-fermenting propionigenic bacteria are known, and understanding of their metabolism is limited. Anaerotignum neopropionicum is a propionate-producing bacterium that uses the acrylate pathway to ferment ethanol and CO2 to propionate and acetate. In this work, we used computational and experimental methods to study the metabolism of A. neopropionicum and, in particular, the pathway for conversion of ethanol into propionate. Results Our work describes iANEO_SB607, the first genome-scale metabolic model (GEM) of A. neopropionicum. The model was built combining the use of automatic tools with an extensive manual curation process, and it was validated with experimental data from this and published studies. The model predicted growth of A. neopropionicum on ethanol, lactate, sugars and amino acids, matching observed phenotypes. In addition, the model was used to implement a dynamic flux balance analysis (dFBA) approach that accurately predicted the fermentation profile of A. neopropionicum during batch growth on ethanol. A systematic analysis of the metabolism of A. neopropionicum combined with model simulations shed light into the mechanism of ethanol fermentation via the acrylate pathway, and revealed the presence of the electron-transferring complexes NADH-dependent reduced ferredoxin:NADP+ oxidoreductase (Nfn) and acryloyl-CoA reductase-EtfAB, identified for the first time in this bacterium. Conclusions The realisation of the GEM iANEO_SB607 is a stepping stone towards the understanding of the metabolism of the propionate-producer A. neopropionicum. With it, we have gained insight into the functioning of the acrylate pathway and energetic aspects of the cell, with focus on the fermentation of ethanol. Overall, this study provides a basis to further exploit the potential of propionigenic bacteria as microbial cell factories.

One-carbon (C1) compounds are promising feedstocks for the sustainable production of commodity chemicals. CO2 is a particularly advantageous C1-feedstock since it is an unwanted industrial off-gas that can be converted into valuable products while reducing its atmospheric levels. Acetogens are microorganisms that can grow on CO2/H2 gas mixtures and syngas converting these substrates into ethanol and acetate. Co-cultivation of acetogens with other microbial species that can further process such products, can expand the variety of products to, for example, medium chain fatty acids (MCFA) and longer chain alcohols. Solventogens are microorganisms known to produce MCFA and alcohols via the acetone-butanol-ethanol (ABE) fermentation in which acetate is a key metabolite. Thus, co-cultivation of an acetogen and a solventogen in a consortium provides a potential platform to produce valuable chemicals from CO2. In this study, metabolic modeling was implemented to design a new co-culture of an acetogen and a solventogen to produce butyrate from CO2/H2 mixtures. The model-driven approach suggested the ability of the studied solventogenic species to grow on lactate/glycerol with acetate as co-substrate. This ability was confirmed experimentally by cultivation of Clostridium beijerinckii on these substrates in batch serum bottles and subsequently in pH-controlled bioreactors. Community modeling also suggested that a novel microbial consortium consisting of the acetogen Clostridium autoethanogenum, and the solventogen C. beijerinckii would be feasible and stable. On the basis of this prediction, a co-culture was experimentally established. C. autoethanogenum grew on CO2/H2 producing acetate and traces of ethanol. Acetate was in turn, consumed by C. beijerinckii together with lactate, producing butyrate. These results show that community modeling of metabolism is a valuable tool to guide the design of microbial consortia for the tailored production of chemicals from renewable resources.

Abstract Background The genus Xanthomonas has long been considered to consist predominantly of plant pathogens, but over the last decade there has been an increasing number of reports on non-pathogenic and endophytic members. As Xanthomonas species are prevalent pathogens on a wide variety of important crops around the world, there is a need to distinguish between these plant-associated phenotypes. To date a large number of Xanthomonas genomes have been sequenced, which enables the application of machine learning (ML) approaches on the genome content to predict this phenotype. Until now such approaches to the pathogenomics of Xanthomonas strains have been hampered by the fragmentation of information regarding pathogenicity of individual strains over many studies. Unification of this information into a single resource was therefore considered to be an essential step. Results Mining of 39 papers considering both plant-associated phenotypes, allowed for a phenotypic classification of 578 Xanthomonas strains. For 65 plant-pathogenic and 53 non-pathogenic strains the corresponding genomes were available and de novo annotated for the presence of Pfam protein domains used as features to train and compare three ML classification algorithms; CART, Lasso and Random Forest. Conclusion The literature resource in combination with recursive feature extraction used in the ML classification algorithms provided further insights into the virulence enabling factors, but also highlighted domains linked to traits not present in pathogenic strains.

Abstract Background Akkermansia muciniphila is a member of the human gut microbiota where it resides in the mucus layer and uses mucin as the sole carbon, nitrogen and energy source. A. muciniphila is the only representative of the Verrucomicrobia phylum in the human gut. However, A. muciniphila 16S rRNA gene sequences have also been found in the intestines of many vertebrates. Results We detected A. muciniphila-like bacteria in the intestines of animals belonging to 15 out of 16 mammalian orders. In addition, other species belonging to the Verrucomicrobia phylum were detected in fecal samples. We isolated 10 new A. muciniphila strains from the feces of chimpanzee, siamang, mouse, pig, reindeer, horse and elephant. The physiology and genome of these strains were highly similar in comparison to the type strain A. muciniphila MucT. Overall, the genomes of the new strains showed high average nucleotide identity (93.9 to 99.7%). In these genomes, we detected considerable conservation of at least 75 of the 78 mucin degradation genes that were previously detected in the genome of the type strain MucT. Conclusions The low genomic divergence observed in the new strains may indicate that A. muciniphila favors mucosal colonization independent of the differences in hosts. In addition, the conserved mucus degradation capability points towards a similar beneficial role of the new strains in regulating host metabolic health.

Rational lifestyle engineering has improved performance of microbial cell factories in Industrial Biotechnology. However, few innovations make it through the Valley of Death to market implementation. To gain insights into the views of industry and academia on key bottlenecks and opportunities to reach market implementation, interviews were conducted with industrial and academic participants, helping us gathering the characteristics that any cell factory and bioprocess must have as well as commonly recognised opportunities.We found that academics are limited by technical factors, whereas industry is restricted by technical, sector-dependent and societal ones, leading to misalignments of interest which often result in miscommunication and missed cooperation opportunities. Although both consider that academia must perform curiosity-driven research, there is pressure for short-term industrial applications, which widens the Valley of Death in Biotechnology. In this context, start-ups could be the answer to traverse this valley more effectively, particularly when embedded in adequate innovation ecosystems.

Aspergillus niger is a robust microbial cell factory for organic acid production. However, the regulation of many industrially important pathways is still poorly understood. The regulation of the glucose oxidase (Gox) expression system, involved in the biosynthesis of gluconic acid, has recently been uncovered. The results of that study show hydrogen peroxide, a by-product of the extracellular conversion of glucose to gluconate, has a pivotal role as a signaling molecule in the induction of this system. In this study, the facilitated diffusion of hydrogen peroxide via aquaporin water channels (AQPs) was studied. AQPs are transmembrane proteins of the major intrinsic proteins (MIPs) superfamily. In addition to water and glycerol, they may also transport small solutes such as hydrogen peroxide. The genome sequence of A. niger N402 was screened for putative AQPs. Seven AQPs were found and could be classified into three main groups. One protein (AQPA) belonged to orthodox AQP, three (AQPB, AQPD, and AQPE) were grouped in aquaglyceroporins (AQGP), two (AQPC and AQPF) were in X-intrinsic proteins (XIPs), and the other (AQPG) could not be classified. Their ability to facilitate diffusion of hydrogen peroxide was identified using yeast phenotypic growth assays and by studying AQP gene knock-outs in A. niger. The X-intrinsic protein AQPF appears to play roles in facilitating hydrogen peroxide transport across the cellular membrane in both Saccharomyces cerevisiae and A. niger experiments.

Abstract Background Pseudomonas putida KT2440 is a metabolically versatile, HV1-certified, genetically accessible, and thus interesting microbial chassis for biotechnological applications. However, its obligate aerobic nature hampers production of oxygen sensitive products and drives up costs in large scale fermentation. The inability to perform anaerobic fermentation has been attributed to insufficient ATP production and an inability to produce pyrimidines under these conditions. Addressing these bottlenecks enabled growth under micro-oxic conditions but does not lead to growth or survival under anoxic conditions. Results Here, a data-driven approach was used to develop a rational design for a P. putida KT2440 derivative strain capable of anaerobic respiration. To come to the design, data derived from a genome comparison of 1628 Pseudomonas strains was combined with genome-scale metabolic modelling simulations and a transcriptome dataset of 47 samples representing 14 environmental conditions from the facultative anaerobe Pseudomonas aeruginosa. Conclusions The results indicate that the implementation of anaerobic respiration in P. putida KT2440 would require at least 49 additional genes of known function, at least 8 genes encoding proteins of unknown function, and 3 externally added vitamins.

Abstract Background Cutaneotrichosporon oleaginosus ATCC 20509 is a fast-growing oleaginous basidiomycete yeast that is able to grow in a wide range of low-cost carbon sources including crude glycerol, a byproduct of biodiesel production. When glycerol is used as a carbon source, this yeast can accumulate more than 50% lipids (w/w) with high concentrations of mono-unsaturated fatty acids. Results To increase our understanding of this yeast and to provide a knowledge base for further industrial use, a FAIR re-annotated genome was used to build a genome-scale, constraint-based metabolic model containing 1553 reactions involving 1373 metabolites in 11 compartments. A new description of the biomass synthesis reaction was introduced to account for massive lipid accumulation in conditions with high carbon-to-nitrogen (C/N) ratio in the media. This condition-specific biomass objective function is shown to better predict conditions with high lipid accumulation using glucose, fructose, sucrose, xylose, and glycerol as sole carbon source. Conclusion Contributing to the economic viability of biodiesel as renewable fuel, C. oleaginosus ATCC 20509 can effectively convert crude glycerol waste streams in lipids as a potential bioenergy source. Performance simulations are essential to identify optimal production conditions and to develop and fine tune a cost-effective production process. Our model suggests ATP-citrate lyase as a possible target to further improve lipid production.

Abstract Background Pandemics, even more than other medical problems, require swift integration of knowledge. When caused by a new virus, understanding the underlying biology may help finding solutions. In a setting where there are a large number of loosely related projects and initiatives, we need common ground, also known as a “commons.” Wikidata, a public knowledge graph aligned with Wikipedia, is such a commons and uses unique identifiers to link knowledge in other knowledge bases. However, Wikidata may not always have the right schema for the urgent questions. In this paper, we address this problem by showing how a data schema required for the integration can be modeled with entity schemas represented by Shape Expressions. Results As a telling example, we describe the process of aligning resources on the genomes and proteomes of the SARS-CoV-2 virus and related viruses as well as how Shape Expressions can be defined for Wikidata to model the knowledge, helping others studying the SARS-CoV-2 pandemic. How this model can be used to make data between various resources interoperable is demonstrated by integrating data from NCBI (National Center for Biotechnology Information) Taxonomy, NCBI Genes, UniProt, and WikiPathways. Based on that model, a set of automated applications or bots were written for regular updates of these sources in Wikidata and added to a platform for automatically running these updates. Conclusions Although this workflow is developed and applied in the context of the COVID-19 pandemic, to demonstrate its broader applicability it was also applied to other human coronaviruses (MERS, SARS, human coronavirus NL63, human coronavirus 229E, human coronavirus HKU1, human coronavirus OC4).

Abstract Graphene nanomaterials have attracted a great interest during the last years for different applications, but their possible impact on different biological systems remains unclear. Here, an assessment to understand the toxicity of commercial polycarboxylate functionalized graphene nanoplatelets (GN) on the unicellular fungal model Saccharomyces cerevisiae was performed. While cell proliferation was not negatively affected even in the presence of 800 mg L−1 of the nanomaterial for 24 hours, oxidative stress was induced at a lower concentration (160 mg L−1), after short exposure periods (2 and 4 hours). No DNA damage was observed under a comet assay analysis under the studied conditions. In addition, to pinpoint the molecular mechanisms behind the early oxidative damage induced by GN and to identify possible toxicity pathways, the transcriptome of S. cerevisiae exposed to 160 and 800 mg L−1 of GN was studied. Both GN concentrations induced expression changes in a common group of genes (337), many of them related to the fungal response to reduce the nanoparticles toxicity and to maintain cell homeostasis. Also, a high number of genes were only differentially expressed in the GN800 condition (3254), indicating that high GN concentrations can induce severe changes in the physiological state of the yeast.

Microbial fermentation of synthesis gas (syngas) is becoming more attractive for sustainable production of commodity chemicals. To date, syngas fermentation focuses mainly on the use of Clostridium species for the production of small organic molecules such as ethanol and acetate. The co-cultivation of syngas-fermenting microorganisms with chain-elongating bacteria can expand the range of possible products, allowing, for instance, the production of medium-chain fatty acids (MCFA) and alcohols from syngas. To explore these possibilities, we report herein a genome-scale, constraint-based metabolic model to describe growth of a co-culture of Clostridium autoethanogenum and Clostridium kluyveri on syngas for the production of valuable compounds. Community flux balance analysis was used to gain insight into the metabolism of the two strains and their interactions, and to reveal potential strategies enabling production of butyrate and hexanoate.The model suggests that one strategy to optimize the production of medium-chain fatty-acids from syngas would be the addition of succinate. According to the prediction, addition of succinate would increase the pool of crotonyl-CoA and the ethanol/acetate uptake ratio in C. kluyveri, resulting in a flux of up to 60% of electrons into hexanoate. Another potential way to further optimize butyrate and hexanoate production would be an increase of C. autoethanogenum ethanol production. Blocking either acetaldehyde dehydrogenase or formate dehydrogenase (ferredoxin) activity or formate transport, in the C. autoethanogenum metabolic model could potentially lead to an up to 150% increase in ethanol production.

Abstract Background The genus Trichococcus currently contains nine species: T. flocculiformis, T. pasteurii, T. palustris, T. collinsii, T. patagoniensis, T. ilyis, T. paludicola, T. alkaliphilus, and T. shcherbakoviae. In general, Trichococcus species can degrade a wide range of carbohydrates. However, only T. pasteurii and a non-characterized strain of Trichococcus, strain ES5, have the capacity of converting glycerol to mainly 1,3-propanediol. Comparative genomic analysis of Trichococcus species provides the opportunity to further explore the physiological potential and uncover novel properties of this genus. Results In this study, a genotype-phenotype comparative analysis of Trichococcus strains was performed. The genome of Trichococcus strain ES5 was sequenced and included in the comparison with the other nine type strains. Genes encoding functions related to e.g. the utilization of different carbon sources (glycerol, arabinan and alginate), antibiotic resistance, tolerance to low temperature and osmoregulation could be identified in all the sequences analysed. T. pasteurii and Trichococcus strain ES5 contain a operon with genes encoding necessary enzymes for 1,3-PDO production from glycerol. All the analysed genomes comprise genes encoding for cold shock domains, but only five of the Trichococcus species can grow at 0 °C. Protein domains associated to osmoregulation mechanisms are encoded in the genomes of all Trichococcus species, except in T. palustris, which had a lower resistance to salinity than the other nine studied Trichococcus strains. Conclusions Genome analysis and comparison of ten Trichococcus strains allowed the identification of physiological traits related to substrate utilization and environmental stress resistance (e.g. to cold and salinity). Some substrates were used by single species, e.g. alginate by T. collinsii and arabinan by T. alkaliphilus. Strain ES5 may represent a subspecies of Trichococcus flocculiformis and contrary to the type strain (DSM 2094T), is able to grow on glycerol with the production of 1,3-propanediol.

Abstract Background Pseudomonas putida is a metabolically versatile, genetically accessible, and stress-robust species with outstanding potential to be used as a workhorse for industrial applications. While industry recognises the importance of robustness under micro-oxic conditions for a stable production process, the obligate aerobic nature of P. putida, attributed to its inability to produce sufficient ATP and maintain its redox balance without molecular oxygen, severely limits its use for biotechnology applications. Results Here, a combination of genome-scale metabolic modelling and comparative genomics is used to pinpoint essential $$\text {O}_{2}$$ O2 -dependent processes. These explain the inability of the strain to grow under anoxic conditions: a deficient ATP generation and an inability to synthesize essential metabolites. Based on this, several P. putida recombinant strains were constructed harbouring acetate kinase from Escherichia coli for ATP production, and a class I dihydroorotate dehydrogenase and a class III anaerobic ribonucleotide triphosphate reductase from Lactobacillus lactis for the synthesis of essential metabolites. Initial computational designs were fine-tuned by means of adaptive laboratory evolution. Conclusions We demonstrated the value of combining in silico approaches, experimental validation and adaptive laboratory evolution for microbial design by making the strictly aerobic Pseudomonas putida able to grow under micro-oxic conditions.

Lactic acid (LA) is a promising bio-based chemical that has broad applications in food, nutraceutical, and bioplastic industries. However, production of the D-form of LA (D-LA) from fermentative organisms is lacking. In this study, Saccharomyces cerevisiae harboring the D-lactate dehydrogenase (DLDH) gene from Leuconostoc mesenteroides was constructed (CEN.PK2_DLDH). To increase D-LA production, the CRISPR/Cas12a system was used for the deletion of gpd1, gpd2, and adh1 to minimize glycerol and ethanol production. Although an improved D-LA titer was observed for both CEN.PK2_DLDHΔgpd and CEN.PK2_DLDHΔgpdΔadh1, growth impairment was observed. To enhance the D-LA productivity, CEN.PK2_DLDHΔgpd was crossed with the weak acid-tolerant S. cerevisiae BCC39850. The isolated hybrid2 showed a maximum D-LA concentration of 23.41 ± 1.65 g/L, equivalent to the improvement in productivity and yield by 2.2 and 1.5 folds, respectively. The simultaneous saccharification and fermentation using alkaline pretreated sugarcane bagasse by the hybrid2 led to an improved D-LA conversion yield on both the washed solid and whole slurry (0.33 and 0.24 g/g glucan). Our findings show the exploitation of natural yeast diversity and the potential strategy of gene editing combined with conventional breeding on improving the performance of S. cerevisiae for the production of industrially potent products.

Ca. Neoehrlichia mikurensis is widely prevalent in I. ricinus across Europe and has been associated with human disease. However, diagnostic modalities are limited, and much is still unknown about its biology. Here, we present the first complete Ca. Neoehrlichia mikurensis genomes directly derived from wildlife reservoir host tissues, using both long- and short-read sequencing technologies. This pragmatic approach provides an alternative to obtaining sufficient material from clinical cases, a difficult task for emerging infectious diseases, and to expensive and challenging bacterial isolation and culture methods. Both genomes exhibit a larger chromosome than the currently available Ca. Neoehrlichia mikurensis genomes and expand the ability to find new targets for the development of supportive laboratory diagnostics in the future. Moreover, this method could be utilized for other tick-borne pathogens that are difficult to culture.

Mycoplasma hyopneumoniae (M. hyopneumoniae) causes enzootic pneumonia in pigs but it is still largely unknown which host-pathogen interactions enable persistent infection and cause disease. In this study, we analyzed the host and bacterial transcriptomes during infection using RNA sequencing. Comparison of the transcriptome of lung lesion tissue from infected pigs with lung tissue from non-infected animals, identified 424 differentially expressed genes (FDR < 0.01 and fold change > 1.5LOG2). These genes were part of the following major pathways of the immune system: interleukin signaling (type 4, 10, 13, and 18), regulation of Toll-like receptors by endogenous ligand and activation of C3 and C5 in the complement system. Besides analyzing the lung transcriptome, a sampling protocol was developed to obtain enough bacterial mRNA from infected lung tissue for RNA sequencing. This was done by flushing infected lobes in the lung, and subsequently enriching for bacterial RNA. On average, 2.2 million bacterial reads were obtained per biological replicate to analyze the bacterial in vivo transcriptome. We compared the in vivo bacterial transcriptome with the transcriptome of bacteria grown in vitro and identified 22 up-regulated and 30 down-regulated genes (FDR < 0.01 and fold change > 2LOG2). Six out of seven genes in the operon encoding the mycoplasma specific F1-like ATPase (MHP_RS02445-MHP_RS02475) and all genes in the operon MHP_RS01965-MHP_RS01990 with functions related to nucleotide metabolism, spermidine transport and glycerol-3-phoshate transport were up-regulated in vivo. Down-regulated in vivo were genes related to glycerol uptake, cilium adhesion (P102), cell division and myo-inositol metabolism. In addition to providing a novel method to isolate bacterial mRNA from infected lung, this study provided insights into changes in gene expression during infection, which could help development of novel treatment strategies against enzootic pneumonia caused by M. hyopneumoniae.

NG-Tax 2.0 is a semantic framework for FAIR high-throughput analysis and classification of marker gene amplicon sequences including bacterial and archaeal 16S ribosomal RNA (rRNA), eukaryotic 18S rRNA and ribosomal intergenic transcribed spacer sequences. It can directly use single or merged reads, paired-end reads and unmerged paired-end reads from long range fragments as input to generate de novo amplicon sequence variants (ASV). Using the RDF data model, ASV’s can be automatically stored in a graph database as objects that link ASV sequences with the full data-wise and element-wise provenance, thereby achieving the level of interoperability required to utilize such data to its full potential. The graph database can be directly queried, allowing for comparative analyses of over thousands of samples and is connected with an interactive Rshiny toolbox for analysis and visualization of (meta) data. Additionally, NG-Tax 2.0 exports an extended BIOM 1.0 (JSON) file as starting point for further analyses by other means. The extended BIOM file contains new attribute types to include information about the command arguments used, the sequences of the ASVs formed, classification confidence scores and is backwards compatible. The performance of NG-Tax 2.0 was compared with DADA2, using the plugin in the QIIME 2 analysis pipeline. Fourteen 16S rRNA gene amplicon mock community samples were obtained from the literature and evaluated. Precision of NG-Tax 2.0 was significantly higher with an average of 0.95 vs 0.58 for QIIME2-DADA2 while recall was comparable with an average of 0.85 and 0.77, respectively. NG-Tax 2.0 is written in Java. The code, the ontology, a Galaxy platform implementation, the analysis toolbox, tutorials and example SPARQL queries are freely available at http://wurssb.gitlab.io/ngtax under the MIT License.

Abstract Humans are not autonomous entities. We are all living in a complex environment, interacting not only with our peers, but as true holobionts; we are also very much in interaction with our coexisting microbial ecosystems living on and especially within us, in the intestine. Intestinal microorganisms, often collectively referred to as intestinal microbiota, contribute significantly to our daily energy uptake by breaking down complex carbohydrates into simple sugars, which are fermented to short-chain fatty acids and subsequently absorbed by human cells. They also have an impact on our immune system, by suppressing or enhancing the growth of malevolent and beneficial microbes. Our lifestyle can have a large influence on this ecosystem. What and how much we consume can tip the ecological balance in the intestine. A “western diet” containing mainly processed food will have a different effect on our health than a balanced diet fortified with pre- and probiotics. In recent years, new technologies have emerged, which made a more detailed understanding of microbial communities and ecosystems feasible. This includes progress in the sequencing of PCR-amplified phylogenetic marker genes as well as the collective microbial metagenome and metatranscriptome, allowing us to determine with an increasing level of detail, which microbial species are in the microbiota, understand what these microorganisms do and how they respond to changes in lifestyle and diet. These new technologies also include the use of synthetic and in vitro systems, which allow us to study the impact of substrates and addition of specific microbes to microbial communities at a high level of detail, and enable us to gather quantitative data for modelling purposes. Here, we will review the current state of microbiome research, summarizing the computational methodologies in this area and highlighting possible outcomes for personalized nutrition and medicine.

Little is known about the metabolic state of Mycobacterium tuberculosis (Mtb) inside the phagosome, a compartment inside phagocytes for killing pathogens and other foreign substances. We have developed a combined model of Mtb and human metabolism, sMtb-RECON and used this model to predict the metabolic state of Mtb during infection of the host. Amino acids are predicted to be used for energy production as well as biomass formation. Subsequently we assessed the effect of increasing dosages of drugs targeting metabolism on the metabolic state of the pathogen and predict resulting metabolic adaptations and flux rerouting through various pathways. In particular, the TCA cycle becomes more important upon drug application, as well as alanine, aspartate, glutamate, proline, arginine and porphyrin metabolism, while glycine, serine, and threonine metabolism become less important. We modeled the effect of 11 metabolically active drugs. Notably, the effect of eight could be recreated and two major profiles of the metabolic state were predicted. The profiles of the metabolic states of Mtb affected by the drugs BTZ043, cycloserine and its derivative terizidone, ethambutol, ethionamide, propionamide, and isoniazid were very similar, while TMC207 is predicted to have quite a different effect on metabolism as it inhibits ATP synthase and therefore indirectly interferes with a multitude of metabolic pathways.

Aspergillus niger is the major industrial citrate producer worldwide. Export as well as uptake of citric acid are believed to occur by active, proton-dependent, symport systems. Both are major bottlenecks for industrial citrate production. Therefore, we assessed the consequences of deleting the citT gene encoding the A. niger citrate exporter, effectively blocking active citrate export. We followed the consumption of glucose and citrate as carbon sources, monitored the secretion of organic acids and carried out a thorough transcriptome pathway enrichment analysis. Under controlled cultivation conditions that normally promote citrate secretion, the knock-out strain secreted negligible amounts of citrate. Blocking active citrate export in this way led to a reduced glucose uptake and a reduced expression of high-affinity glucose transporter genes, mstG and mstH. The glyoxylate shunt was strongly activated and an increased expression of the OAH gene was observed, resulting in a more than two-fold higher concentration of oxalate in the medium. Deletion of citT did not affect citrate uptake suggesting that citrate export and citrate uptake are uncoupled from the system.

Background: Massive high-throughput sequencing of short, hypervariable segments of the 16S ribosomal RNA (rRNA) gene has transformed the methodological landscape describing microbial diversity within and across complex biomes. However, several studies have shown that the methodology rather than the biological variation is responsible for the observed sample composition and distribution. This compromises meta-analyses, although this fact is often disregarded. Results: To facilitate true meta-analysis of microbiome studies, we developed NG-Tax, a pipeline for 16S rRNA gene amplicon sequence analysis that was validated with different mock communities and benchmarked against QIIME as a frequently used pipeline. The microbial composition of 49 independently amplified mock samples was characterized by sequencing two variable 16S rRNA gene regions, V4 and V5-V6, in three separate sequencing runs on Illumina’s HiSeq2000 platform. This allowed for the evaluation of important causes of technical bias in taxonomic classification: 1) run-to-run sequencing variation, 2) PCR–error, and 3) region/primer specific amplification bias. Despite the short read length (~140 nt) and all technical biases, the average specificity of the taxonomic assignment for the phylotypes included in the mock communities was 97.78%. On average 99.95% and 88.43% of the reads could be assigned to at least family or genus level, respectively, while assignment to ‘spurious genera’ represented on average only 0.21% of the reads per sample. Analysis of α- and β-diversity confirmed conclusions guided by biology rather than the aforementioned methodological aspects, which was not achieved with QIIME. Conclusions: Different biological outcomes are commonly observed due to 16S rRNA region-specific performance. NG-Tax demonstrated high robustness against choice of region and other technical biases associated with 16S rRNA gene amplicon sequencing studies, diminishing their impact and providing accurate qualitative and quantitative representation of the true sample composition. This will improve comparability between studies and facilitate efforts towards standardization.

Aspergillus niger is an industrially important source for gluconic acid and glucose oxidase (GOx), a secreted commercially important flavoprotein which catalyses the oxidation of β-D-glucose by molecular oxygen to D-glucolactone and hydrogen peroxide. Expression of goxC, the GOx encoding gene and the concomitant two step conversion of glucose to gluconic acid requires oxygen and the presence of significant amounts of glucose in the medium and is optimally induced at pH 5.5. The molecular mechanisms underlying regulation of goxC expression are, however, still enigmatic. Genetic studies aimed at understanding GOx induction have indicated the involvement of at least seven complementation groups, for none of which the molecular basis has been resolved. In this study, a mapping-by-sequencing forward genetics approach was used to uncover the molecular role of the goxB locus in goxC expression. Using the Illumina and PacBio sequencing platforms a hybrid high quality draft genome assembly of laboratory strain N402 was obtained and used as a reference for mapping of genomic reads obtained from the derivative NW103:goxB mutant strain. The goxB locus encodes a thioredoxin reductase. A deletion of the encoding gene in the N402 parent strain led to a high constitutive expression level of the GOx and the lactonase encoding genes required for the two-step conversion of glucose in gluconic acid and of the catR gene encoding catalase R. This high constitutive level of expression was observed to be irrespective of the carbon source and oxidative stress applied. A model clarifying the role of GoxB in the regulation of the expression of goxC involving hydrogen peroxide as second messenger is presented.

Genome-scale metabolic models of Mycobacterium tuberculosis (Mtb), the causative agent of tuberculosis, have been envisioned as a platform for drug discovery. By systematically probing the networks that underpin such models, the reactions that are essential for Mtb are identified. A majority of these reactions are catalyzed by enzymes and thus represent candidate drug targets to fight an Mtb infection. Nevertheless, this is complicated by the limited knowledge on the environment that Mtb encounters during infection. Modeling the behavior of the bacteria during infection requires knowledge of the so-called biomass reaction that represents bacterial biomass composition. This composition varies in different environments or bacterial growth phases. Accurate modeling of the metabolic state requires a precise biomass reaction for the described condition. In recent years, additional insights in the in-host environment occupied by Mtb have been gained as transcript abundance data of interacting host and pathogen have become available. Therefore, we used transcript abundance data and developed a straightforward and systematic method to obtain a condition-specific biomass reaction for Mtb during in vitro growth and during infection of its host. The method described herein is virtually free of any pre-set assumptions on uptake rates of nutrients, making it suitable for exploring environments with limited accessibility. The condition-specific biomass reaction represents the “metabolic objective” of Mtb in a given environment (in-host growth and growth on defined medium) at a specific time point, and as such allows modeling the bacterial metabolic state in these environments. Five different biomass reactions were used to predict nutrient uptake rates and gene essentiality. Predictions were subsequently compared to available experimental data. Our results show that nutrient uptake can accurately be predicted. Gene essentiality can also be predicted but accurate predictions remain difficult to obtain. In conclusion, a viable strategy to model Mtb metabolism in hard-to-access environments that is virtually free of pre-set assumptions is provided.

Background The microbiota in the small intestine relies on their capacity to rapidly import and ferment available carbohydrates to survive in a complex and highly competitive ecosystem. Understanding how these communities function requires elucidating the role of its key players, the interactions among them and with their environment/host. Methods The genome of the gut bacterium Romboutsia ilealis CRIBT was sequenced with multiple technologies (Illumina paired-end, mate-pair and PacBio). The transcriptome was sequenced (Illumina HiSeq) after growth on three different carbohydrate sources, and short chain fatty acids were measured via HPLC. Results We present the complete genome of Romboutsia ilealis CRIBT, a natural inhabitant and key player of the small intestine of rats. R. ilealis CRIBT possesses a circular chromosome of 2,581,778 bp and a plasmid of 6,145 bp, carrying 2,351 and eight predicted protein coding sequences, respectively. Analysis of the genome revealed limited capacity to synthesize amino acids and vitamins, whereas multiple and partially redundant pathways for the utilization of different relatively simple carbohydrates are present. Transcriptome analysis allowed identification of the key components in the degradation of glucose, L-fucose and fructo-oligosaccharides. Discussion This revealed that R. ilealis CRIBT is adapted to a nutrient-rich environment where carbohydrates, amino acids and vitamins are abundantly available.

Aspergillus niger has an innate ability to secrete various organic acids, including citrate. The conditions required for A. niger citrate overproduction are well described, but the physiological reasons underlying extracellular citrate accumulation are not yet fully understood. One of the less understood culture conditions is the requirement of growth-limiting iron concentrations. While this has been attributed to iron-dependent citrate metabolizing enzymes, this straightforward relationship does not always hold true. Here, we show that an increase in citrate secretion under iron limited conditions is a physiological response consistent with a role of citrate as A. niger iron siderophore. We found that A. niger citrate secretion increases with decreasing amounts of iron added to the culture medium and, in contrast to previous findings, this response is independent of the nitrogen source. Differential transcriptomics analyses of the two A. niger mutants NW305 (gluconate non-producer) and NW186 (gluconate and oxalate non-producer) revealed up-regulation of the citrate biosynthesis gene citA under iron limited conditions compared to iron replete conditions. In addition, we show that A. niger can utilize Fe(III) citrate as iron source. Finally, we discuss our findings in the general context of the pH-dependency of A. niger organic acid production, offering an explanation, besides competition, for why A. niger organic acid production is a sequential process influenced by the external pH of the culture medium.

A functional comparative genome analysis is essential to understand the mechanisms underlying bacterial evolution and adaptation. Detection of functional orthologs using standard global sequence similarity methods faces several problems; the need for defining arbitrary acceptance thresholds for similarity and alignment length, lateral gene acquisition and the high computational cost for finding bi-directional best matches at a large scale. We investigated the use of protein domain architectures for large scale functional comparative analysis as an alternative method. The performance of both approaches was assessed through functional comparison of 446 bacterial genomes sampled at different taxonomic levels. We show that protein domain architectures provide a fast and efficient alternative to methods based on sequence similarity to identify groups of functionally equivalent proteins within and across taxonomic boundaries, and it is suitable for large scale comparative analysis. Running both methods in parallel pinpoints potential functional adaptations that may add to bacterial fitness.

The filamentous fungus Rhizopus delemar naturally accumulates relatively high amounts of fumarate. Although the culture conditions that increase fumarate yields are well established, the network underlying the accumulation of fumarate is not yet fully understood. We set out to increase the knowledge about fumarate accumulation in R. delemar. To this end, we combined a transcriptomics and proteomics approach to identify key metabolic pathways involved in fumarate production in R. delemar, and propose that a substantial part of the fumarate accumulated in R. delemar during nitrogen starvation results from the urea cycle due to amino acid catabolism.

The cupin-type phosphoglucose isomerase (PfPGI) from the hyperthermophilic archaeon Pyrococcus furiosus catalyzes the reversible isomerization of glucose-6-phosphate to fructose-6-phosphate. We investigated PfPGI using protein-engineering bioinformatics tools to select functionally-important residues based on correlated mutation analyses. A pair of amino acids in the periphery of PfPGI was found to be the dominant co-evolving mutation. The position of these selected residues was found to be non-obvious to conventional protein engineering methods. We designed a small smart library of variants by substituting the co-evolved pair and screened their biochemical activity, which revealed their functional relevance. Four mutants were further selected from the library for purification, measurement of their specific activity, crystal structure determination, and metal cofactor coordination analysis. Though the mutant structures and metal cofactor coordination were strikingly similar, variations in their activity correlated with their fine-tuned dynamics and solvent access regulation. Alternative, small smart libraries for enzyme optimization are suggested by our approach, which is able to identify non-obvious yet beneficial mutations.

Genome-scale metabolic models (GEMs) are manually curated repositories describing the metabolic capabilities of an organism. GEMs have been successfully used in different research areas, ranging from systems medicine to biotechnology. However, the different naming conventions (namespaces) of databases used to build GEMs limit model reusability and prevent the integration of existing models. This problem is known in the GEM community, but its extent has not been analyzed in depth. In this study, we investigate the name ambiguity and the multiplicity of non-systematic identifiers and we highlight the (in)consistency in their use in 11 biochemical databases of biochemical reactions and the problems that arise when mapping between different namespaces and databases. We found that such inconsistencies can be as high as 83.1%, thus emphasizing the need for strategies to deal with these issues. Currently, manual verification of the mappings appears to be the only solution to remove inconsistencies when combining models. Finally, we discuss several possible approaches to facilitate (future) unambiguous mapping.

A functional comparative genome analysis is essential to understand the mechanisms underlying bacterial evolution and adaptation. Detection of functional orthologs using standard global sequence similarity methods faces several problems; the need for defining arbitrary acceptance thresholds for similarity and alignment length, lateral gene acquisition and the high computational cost for finding bi-directional best matches at a large scale. We investigated the use of protein domain architectures for large scale functional comparative analysis as an alternative method. The performance of both approaches was assessed through functional comparison of 446 bacterial genomes sampled at different taxonomic levels. We show that protein domain architectures provide a fast and efficient alternative to methods based on sequence similarity to identify groups of functionally equivalent proteins within and across taxonomic boundaries. As the computational cost scales linearly, and not quadratically with the number of genomes, it is suitable for large scale comparative analysis. Running both methods in parallel pinpoints potential functional adaptations that may add to bacterial fitness.

A functional comparative genome analysis is essential to understand the mechanisms underlying bacterial evolution and adaptation. Detection of functional orthologs using standard global sequence similarity methods faces several problems; the need for defining arbitrary acceptance thresholds for similarity and alignment length, lateral gene acquisition and the high computational cost for finding bi-directional best matches at a large scale. We investigated the use of protein domain architectures for large scale functional comparative analysis as an alternative method. The performance of both approaches was assessed through functional comparison of 446 bacterial genomes sampled at different taxonomic levels. We show that protein domain architectures provide a fast and efficient alternative to methods based on sequence similarity to identify groups of functionally equivalent proteins within and across taxonomic bounderies. As the computational cost scales linearly, and not quadratically with the number of genomes, it is suitable for large scale comparative analysis. Running both methods in parallel pinpoints potential functional adaptations that may add to bacterial fitness.

Tuberculosis remains one of the deadliest diseases. Emergence of drug-resistant and multidrug-resistant M. tuberculosis strains makes treating tuberculosis increasingly challenging. In order to develop novel intervention strategies, detailed understanding of the molecular mechanisms behind the success of this pathogen is required. Here, we review recent literature to provide a systems level overview of the molecular and cellular components involved in divalent metal homeostasis and their role in regulating the three main virulence strategies of M. tuberculosis: immune modulation, dormancy and phagosomal rupture. We provide a visual and modular overview of these components and their regulation. Our analysis identified a single regulatory cascade for these three virulence strategies that respond to limited availability of divalent metals in the phagosome.

Background Massive high-throughput sequencing of short, hypervariable segments of the 16S ribosomal RNA (rRNA) gene has transformed the methodological landscape describing microbial diversity within and across complex biomes. However, several studies have shown that the methodology rather than the biological variation is responsible for the observed sample composition and distribution. This compromises true meta-analyses, although this fact is often disregarded. Results To facilitate true meta-analysis of microbiome studies, we developed NG-Tax, a pipeline for 16S rRNA gene amplicon sequence analysis that was validated with different mock communities and benchmarked against QIIME as the currently most frequently used pipeline. The microbial composition of 49 independently amplified mock samples was characterized by sequencing two variable 16S rRNA gene regions, V4 and V5-V6, in three separate sequencing runs on Illumina’s HiSeq2000 platform. This allowed evaluating important factors of technical bias in taxonomic classification: 1) run-to-run sequencing variation, 2) PCR–error, and 3) region/primer specific amplification bias. Despite the short read length (~140 nt) and all technical biases, the average specificity of the taxonomic assignment for the phylotypes included in the mock communities was 96%. On average 99.94% and 92.02% of the reads could be assigned to at least family or genus level, respectively, while assignment to ‘spurious genera’ represented on average only 0.02% of the reads per sample. Analysis of α- and β-diversity confirmed conclusions guided by biology rather than the aforementioned methodological aspects, which was not the case when samples were analysed using QIIME. Conclusions Different biological outcomes are commonly observed due to 16S rRNA region-specific performance. NG-Tax demonstrated high robustness against choice of region and other technical biases associated with 16S rRNA gene amplicon sequencing studies, diminishing their impact and providing accurate qualitative and quantitative representation of the true sample composition. This will improve comparability between studies and facilitate efforts towards standardization.