Search

Abstract This study aimed to highlight the molecular and biochemical changes induced by methylglyoxal (MGO) exposure in SH-SY5Y human neuroblastoma cells, and to explore how these changes contribute to its neurotoxicity, utilizing an integrated proteomics and metabolomics approach. Using label-free quantitative nanoLC-MS/MS proteomics and targeted LC-TQ-MS/MS-based metabolomics, the results revealed that MGO exposure, particularly at cytotoxic levels, significantly altered the proteome and metabolome of SH-SY5Y cells. Analysis of proteomics data showed significant alterations in cellular functions including protein synthesis, cellular structural integrity, mitochondrial function, and oxidative stress responses. Analysis of metabolomics and integration of metabolomics and proteomics data highlighted significant changes in key metabolic pathways including arginine biosynthesis, glutathione metabolism, cysteine and methionine metabolism, and the tricarboxylic acid cycle. These results suggest that MGO exposure induced both toxic effects and adaptive responses in cells. MGO exposure led to increased endoplasmic reticulum stress, disruptions in cellular adhesion and extracellular matrix integrity, mitochondrial dysfunction, and amino acid metabolism disruption, contributing to cellular toxicity. Conversely, cells exhibited adaptive responses by upregulating protein synthesis, activating the Nrf2 pathway, and reprogramming metabolism to counteract dicarbonyl stress and maintain energy levels. Furthermore, a set of key proteins and metabolites associated with these changes were shown to exhibit a significant concentration-dependent decrease or increase in their expression levels with increasing MGO concentrations, suggesting their potential as biomarkers for MGO exposure. Taken together, these findings provide insight into the molecular mechanisms underlying MGO-induced neurotoxicity and potential targets for therapeutic intervention.

Abstract Proteins and peptides found in human milk have bioactive potential to benefit the newborn and support healthy development. Research has been carried out on the health benefits of proteins and peptides, but many questions still need to be answered about the nature of these components, how they are formed, and how they end up in the milk. This study explored and elucidated the complexity of the human milk proteome and peptidome. Proteins and peptides were analyzed with non-targeted nanoLC-Orbitrap-MS/MS in a selection of 297 milk samples from the CHILD Cohort Study. Protein and peptide abundances were determined, and a network was inferred using Gaussian graphical modeling (GGM), allowing an investigation of direct associations. This study showed that signatures of (1) specific mechanisms of transport of different groups of proteins, (2) proteolytic degradation by proteases and aminopeptidases, and (3) coagulation and complement activation are present in human milk. These results show the value of an integrated approach in evaluating large-scale omics data sets and provide valuable information for studies that aim to associate protein or peptide profiles from biofluids such as milk with specific physiological characteristics.

During early lactation, dairy cows have a negative energy balance since their energy demands exceed their energy intake: in this study, we aimed to investigate the association between diet and plasma metabolomics profiles and how these relate to energy unbalance of course in the early-lactation stage. Holstein-Friesian cows were randomly assigned to a glucogenic (n = 15) or lipogenic (n = 15) diet in early lactation. Blood was collected in week 2 and week 4 after calving. Plasma metabolite profiles were detected using liquid chromatography–mass spectrometry (LC-MS), and a total of 39 metabolites were identified. Two plasma metabolomic profiles were available every week for each cow. Metabolite abundance and metabolite ratios were used for the analysis using the XGboost algorithm to discriminate between diet treatment and lactation week. Using metabolite ratios resulted in better discrimination performance compared with the metabolite abundances in assigning cows to a lipogenic diet or a glucogenic diet. The quality of the discrimination of performance of lipogenic diet and glucogenic diet effects improved from 0.606 to 0.753 and from 0.696 to 0.842 in week 2 and week 4 (as measured by area under the curve, AUC), when the metabolite abundance ratios were used instead of abundances. The top discriminating ratios for diet were the ratio of arginine to tyrosine and the ratio of aspartic acid to valine in week 2 and week 4, respectively. For cows fed the lipogenic diet, choline and the ratio of creatinine to tryptophan were top features to discriminate cows in week 2 vs. week 4. For cows fed the glucogenic diet, methionine and the ratio of 4-hydroxyproline to choline were top features to discriminate dietary effects in week 2 or week 4. This study shows the added value of using metabolite abundance ratios to discriminate between lipogenic and glucogenic diet and lactation weeks in early-lactation cows when using metabolomics data. The application of this research will help to accurately regulate the nutrition of lactating dairy cows and promote sustainable agricultural development.

We present the OrganelX e-Science Web Server that provides a user-friendly implementation of the In-Pero and In-Mito classifiers for sub-peroxisomal and sub-mitochondrial localization of peroxisomal and mitochondrial proteins and the Is-PTS1 algorithm for detecting and validating potential peroxisomal proteins carrying a PTS1 signal sequence. The OrganelX e-Science Web Server is available at https://organelx.hpc.rug.nl/fasta/.

Abstract The ability to detect and characterize bacteria within a biological sample is crucial for the monitoring of infections and epidemics, as well as for the study of human health and its relationship with commensal microorganisms. To this aim, a commonly used technique is the 16S rRNA gene targeted sequencing. PCR-amplified 16S sequences derived from the sample of interest are usually clustered into the so-called Operational Taxonomic Units (OTUs) based on pairwise similarities. Then, representative OTU sequences are compared with reference (human-made) databases to derive their phylogeny and taxonomic classification. Here, we propose a new reference-free approach to define the phylogenetic distance between bacteria based on protein domains, which are the evolving units of proteins. We extract the protein domain profiles of 3368 bacterial genomes and we use an ecological approach to model their Relative Species Abundance distribution. Based on the model parameters, we then derive a new measurement of phylogenetic distance. Finally, we show that such model-based distance is capable of detecting differences between bacteria in cases in which the 16S rRNA-based method fails, providing a possibly complementary approach , which is particularly promising for the analysis of bacterial populations measured by shotgun sequencing.

Abstract Background Necrotising soft tissue infections (NSTIs) are rapidly progressing bacterial infections usually caused by either several pathogens in unison (polymicrobial infections) or Streptococcus pyogenes (mono-microbial infection). These infections are rare and are associated with high mortality rates. However, the underlying pathogenic mechanisms in this heterogeneous group remain elusive. Methods In this study, we built interactomes at both the population and individual levels consisting of host-pathogen interactions inferred from dual RNA-Seq gene transcriptomic profiles of the biopsies from NSTI patients. Results NSTI type-specific responses in the host were uncovered. The S. pyogenes mono-microbial subnetwork was enriched with host genes annotated with involved in cytokine production and regulation of response to stress. The polymicrobial network consisted of several significant associations between different species (S. pyogenes, Porphyromonas asaccharolytica and Escherichia coli) and host genes. The host genes associated with S. pyogenes in this subnetwork were characterised by cellular response to cytokines. We further found several virulence factors including hyaluronan synthase, Sic1, Isp, SagF, SagG, ScfAB-operon, Fba and genes upstream and downstream of EndoS along with bacterial housekeeping genes interacting with the human stress and immune response in various subnetworks between host and pathogen. Conclusions At the population level, we found aetiology-dependent responses showing the potential modes of entry and immune evasion strategies employed by S. pyogenes, congruent with general cellular processes such as differentiation and proliferation. After stratifying the patients based on the subject-specific networks to study the patient-specific response, we observed different patient groups with different collagens, cytoskeleton and actin monomers in association with virulence factors, immunogenic proteins and housekeeping genes which we utilised to postulate differing modes of entry and immune evasion for different bacteria in relationship to the patients’ phenotype.

In the scientific literature data analysis results are often presented when samples from different experiments or different conditions, technical replicates or times series are merged to increase the sample size before calculating the correlation coefficient. This way of proceeding violates two basic assumptions underlying the use of the correlation coefficient: sampling from one population and independence of the observations (independence of errors). Since correlations are used to measure and infer associations between biological entities, this has tremendous implications on the reliability of scientific results, as the violation of these assumption leads to wrong and biased results. In this technical note, I review some basic properties of the Pearson’s correlation coefficient and illustrate some exemplary problems with simulated and experimental data, taking a didactic approach with the use of supporting graphical examples.

BackgroundThe human milk proteome comprises a vast number of proteins with immunomodulatory functions, but it is not clear how this relates to allergy of the mother or allergy development in the breastfed infant. This study aimed to explore the relation between the human milk proteome and allergy of both mother and child.MethodsProteins were analyzed in milk samples from a subset of 300 mother-child dyads from the Canadian CHILD Cohort Study, selected based on maternal and child allergy phenotypes. For this selection, the definition of “allergy” included food allergy, eczema, allergic rhinitis, and asthma. Proteins were analyzed with non-targeted shotgun proteomics using filter-aided sample preparation (FASP) and nanoLC-Orbitrap-MS/MS. Protein abundances, based on label-free quantification, were compared using multiple statistical approaches, including univariate, multivariate, and network analyses.ResultsUsing univariate analysis, we observed a trend that milk for infants who develop an allergy by 3 years of age contains higher abundances of immunoglobulin chains, irrespective of the allergy status of the mother. This observation suggests a difference in the milk’s immunological potential, which might be related to the development of the infant’s immune system. Furthermore, network analysis showed overall increased connectivity of proteins in the milk of allergic mothers and milk for infants who ultimately develop an allergy. This difference in connectivity was especially noted for proteins involved in the protein translation machinery and may be due to the physiological status of the mother, which is reflected in the interconnectedness of proteins in her milk. In addition, it was shown that network analysis complements the other methods for data analysis by revealing complex associations between the milk proteome and mother-child allergy status.ConclusionTogether, these findings give new insights into how the human milk proteome, through differences in the abundance of individual proteins and protein-protein associations, relates to the allergy status of mother and child. In addition, these results inspire new research directions into the complex interplay of the mother-milk-infant triad and allergy.

Systems and Synthetic Biology are complementary fields emerging side-by-side into mainstream scientific research. Whilst systems biologists focus on understanding natural systems, synthetic biologists wish to modify, adapt and re-purpose biological systems towards certain desired goals, for example enhancing efficiency and robustness of desired biological traits. In both fields, data analysis, predictive mathematical modelling, experimental design, and controlled experimentation are crucial to obtain reproducible results and understand how applications can be scaled to larger systems and processes. As such, students from Life Sciences, Engineering, and Mathematics backgrounds must be taught fundamentals in biological systems, experimental techniques, mathematics, and data analysis/statistics. In addition, students must be trained for future multidisciplinary careers, where the interaction and communication between experimental and modelling researchers is fundamental. With the acceleration of technological developments (both computational and experimental) continuing unabated, educators need to bridge the increasing gap between fundamentally-required knowledge and skills that students need to pursue future academic or industrial research projects. In this paper, we will discuss how we have re-designed an introductory course in Systems and Synthetic Biology at Wageningen University and Research (Netherlands) that is targeted simultaneously to mathematical/computational students with an interest in biology and experimental methods, and to Life Science students interested in learning how biological systems can be mathematically analysed and modelled. The course highlights the links between fundamental methodologies and recently developed technologies within the Systems and Synthetic Biology fields. The course was re-designed for the 2021/22 academic year, we report that students from biology and biotechnology programmes graded their satisfaction of the course as 4.4 out of 5. We discuss how the course can act as a gateway to advanced courses in Systems Biology-oriented curricula (comprising: data infrastructure, modelling, and experimental synthetic biology), and towards future research projects.

Abstract Background Staphylococcus and Streptococcus species can cause many different diseases, ranging from mild skin infections to life-threatening necrotizing fasciitis. Both genera consist of commensal species that colonize the skin and nose of humans and animals, and of which some can display a pathogenic phenotype. Results We compared 235 Staphylococcus and 315 Streptococcus genomes based on their protein domain content. We show the relationships between protein persistence and essentiality by integrating essentiality predictions from two metabolic models and essentiality measurements from six large-scale transposon mutagenesis experiments. We identified clusters of strains within species based on proteins associated to similar biological processes. We built Random Forest classifiers that predicted the zoonotic potential. Furthermore, we identified shared attributes between of Staphylococcus aureus and Streptococcus pyogenes that allow them to cause necrotizing fasciitis. Conclusions Differences observed in clustering of strains based on functional groups of proteins correlate with phenotypes such as host tropism, capability to infect multiple hosts and drug resistance. Our method provides a solid basis towards large-scale prediction of phenotypes based on genomic information.

BackgroundSince sequencing techniques have become less expensive, larger sample sizes are applicable for microbiota studies. The aim of this study is to show how, and to what extent, different diversity metrics and different compositions of the microbiota influence the needed sample size to observe dissimilar groups. Empirical 16S rRNA amplicon sequence data obtained from animal experiments, observational human data, and simulated data were used to perform retrospective power calculations. A wide variation of alpha diversity and beta diversity metrics were used to compare the different microbiota datasets and the effect on the sample size.ResultsOur data showed that beta diversity metrics are the most sensitive to observe differences as compared with alpha diversity metrics. The structure of the data influenced which alpha metrics are the most sensitive. Regarding beta diversity, the Bray–Curtis metric is in general the most sensitive to observe differences between groups, resulting in lower sample size and potential publication bias.ConclusionWe recommend performing power calculations and to use multiple diversity metrics as an outcome measure. To improve microbiota studies, awareness needs to be raised on the sensitivity and bias for microbiota research outcomes created by the used metrics rather than biological differences. We have seen that different alpha and beta diversity metrics lead to different study power: because of this, one could be naturally tempted to try all possible metrics until one or more are found that give a statistically significant test result, i.e., p-value < α. This way of proceeding is one of the many forms of the so-called p-value hacking. To this end, in our opinion, the only way to protect ourselves from (the temptation of) p-hacking would be to publish a statistical plan before experiments are initiated, describing the outcomes of interest and the corresponding statistical analyses to be performed.

Colorectal cancer (CRC), one of the most prevalent and deadly cancers worldwide, generally evolves from adenomatous polyps. The understanding of the molecular mechanisms underlying this pathological evolution is crucial for diagnostic and prognostic purposes. Integrative systems biology approaches offer an optimal point of view to analyze CRC and patients with polyposis. The present study analyzed the association networks constructed from a publicly available array of 113 serum metabolites measured on a cohort of 234 subjects from three groups (66 CRC patients, 76 patients with polyposis, and 92 healthy controls), which concentrations were obtained via targeted liquid chromatography-tandem mass spectrometry. In terms of architecture, topology, and connectivity, the metabolite-metabolite association network of CRC patients appears to be completely different with respect to patients with polyposis and healthy controls. The most relevant nodes in the CRC network are those related to energy metabolism. Interestingly, phenylalanine, tyrosine, and tryptophan metabolism are found to be involved in both CRC and polyposis. Our results demonstrate that the characterization of metabolite–metabolite association networks is a promising and powerful tool to investigate molecular aspects of CRC.

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.

Abstract In early lactation, dairy cows typically have a negative energy balance which has been related to metabolic disorders, compromised health and fertility, and reduced productive lifespan. Assessment of the energy balance, however, is not easy on the farm. Our aims were to investigate the milk metabolic profiles of dairy cows in early lactation, and to obtain models to estimate energy balance from milk metabolomics data and milk production traits. Milk samples were collected in week 2 and 7 after calving from 31 dairy cows. For each cow, the energy balance was calculated from energy intake, milk production traits and body weight. A total of 52 milk metabolites were detected using LC-QQQ-MS. Data from different lactation weeks was analysed by partial least squares analysis, the top 15 most relevant variables from the metabolomics data related to energy balance were used to develop reduced linear models to estimate energy balance by forward selection regression. Milk fat yield, glycine, choline and carnitine were important variables to estimate energy balance (adjusted R 2: 0.53 to 0.87, depending on the model). The relationship of these milk metabolites with energy balance is proposed to be related to their roles in cell renewal.

Abstract Intestinal epithelial cells, like Caco-2, are commonly used to study the interaction between food, other luminal factors and the host, often supported by microarray analysis to study the changes in gene expression as a result of the exposure. However, no compiled dataset for Caco-2 has ever been initiated and Caco-2-dedicated gene expression networks are barely available. Here, 341 Caco-2-specific microarray samples were collected from public databases and from in-house experiments pertaining to Caco-2 cells exposed to pathogens, probiotics and several food compounds. Using these datasets, a gene functional association network specific for Caco-2 was generated containing 8937 nodes 129711 edges. Two in silico methods, a modified version of biclustering and the new Differential Expression Correlation Analysis, were developed to identify Caco-2-specific gene targets within a pathway of interest. These methods were subsequently applied to the AhR and Nrf2 signalling pathways and altered expression of the predicted target genes was validated by qPCR in Caco-2 cells exposed to coffee extracts, known to activate both AhR and Nrf2 pathways. The datasets and in silico method(s) to identify and predict responsive target genes can be used to more efficiently design experiments to study Caco-2/intestinal epithelial-relevant biological processes.

In this study, we investigated blood lipoprotein and lipid fraction profiles, quantified using nuclear magnetic resonance, in a cohort of 844 healthy blood donors, integrating standard univariate and multivariate analysis with predictive modeling and network analysis. We observed a strong association of lipoprotein and lipid main fraction profiles with sex and age. Our results suggest an age-dependent remodulation of lipase lipoprotein activity in men and a change in the mechanisms controlling the ratio between esterified and non-esterified cholesterol in both men and women.

Principal Component Analysis (PCA) is a powerful statistical technique for reducing the complexity of data and making patterns and relationships within the data more easily understandable. By using PCA, students can learn to identify the most important features of a data set, visualize relationships between variables, and make informed decisions based on the data. As such, PCA can be an effective tool to increase students data literacy by providing a visual and intuitive way to understand and work with data. This article outlines a teaching strategy to introduce and explain PCA using basic mathematics and statistics together with visual demonstrations.

Metabolite differential connectivity analysis has been successful in investigating potential molecular mechanisms underlying different conditions in biological systems. Correlation and Mutual Information (MI) are two of the most common measures to quantify association and for building metabolite—metabolite association networks and to calculate differential connectivity. In this study, we investigated the performance of correlation and MI to identify significantly differentially connected metabolites. These association measures were compared on (i) 23 publicly available metabolomic data sets and 7 data sets from other fields, (ii) simulated data with known correlation structures, and (iii) data generated using a dynamic metabolic model to simulate real-life observed metabolite concentration profiles. In all cases, we found more differentially connected metabolites when using correlation indices as a measure for association than MI. We also observed that different MI estimation algorithms resulted in difference in performance when applied to data generated using a dynamic model. We concluded that there is no significant benefit in using MI as a replacement for standard Pearson’s or Spearman’s correlation when the application is to quantify and detect differentially connected metabolites.

DNA based methods have been widely used to study the complexity of the rumen microbiota, and it is well known that the method of DNA extraction is a critical step in enabling accurate assessment of this complexity. Rumen fluid (RF) and fibrous content (FC) fractions differ substantially in terms of their physical nature and associated microorganisms. The aim of this study was therefore to assess the effect of four DNA extraction methods (RBB, PBB, FDSS, PQIAmini) differing in cell lysis and/or DNA recovery methods on the observed microbial diversity in RF and FC fractions using samples from four rumen cannulated dairy cows fed 100% grass silage (GS100), 67% GS and 33% maize silage (GS67MS33), 33% GS and 67% MS (GS33MS67), or 100% MS (MS100). An ANOVA statistical test was applied on DNA quality and yield measurements, and it was found that the DNA yield was significantly affected by extraction method (p < 0.001) and fraction (p < 0.001). The 260/280 ratio was not affected by extraction (p = 0.08) but was affected by fraction (p = 0.03). On the other hand, the 260/230 ratio was affected by extraction method (p < 0.001) but not affected by fraction (p = 0.8). However, all four extraction procedures yielded DNA suitable for further analysis of bacterial, archaeal and anaerobic fungal communities using quantitative PCR and pyrosequencing of relevant taxonomic markers. Redundancy analysis (RDA) of bacterial 16S rRNA gene sequence data at the family level showed that there was a significant effect of rumen fraction (p = 0.012), and that PBB (p = 0.012) and FDSS (p = 0.024) also significantly contributed to explaining the observed variation in bacterial community composition. Whilst the DNA extraction method affected the apparent bacterial community composition, no single extraction method could be concluded to be ineffective. No obvious effect of DNA extraction method on the anaerobic fungi or archaea was observed, although fraction effects were evident for both. In summary, the comprehensive assessment of observed communities of bacteria, archaea and anaerobic fungi described here provides insight into a rational basis for selecting an optimal methodology to obtain a representative picture of the rumen microbiota.

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.

Over the past two decades, nuclear magnetic resonance (NMR) has emerged as one of the three principal analytical techniques used in metabolomics (the other two being gas chromatography coupled to mass spectrometry (GC-MS) and liquid chromatography coupled with single-stage mass spectrometry (LC-MS)). The relative ease of sample preparation, the ability to quantify metabolite levels, the high level of experimental reproducibility, and the inherently nondestructive nature of NMR spectroscopy have made it the preferred platform for long-term or large-scale clinical metabolomic studies. These advantages, however, are often outweighed by the fact that most other analytical techniques, including both LC-MS and GC-MS, are inherently more sensitive than NMR, with lower limits of detection typically being 10 to 100 times better. This review is intended to introduce readers to the field of NMR-based metabolomics and to highlight both the advantages and disadvantages of NMR spectroscopy for metabolomic studies. It will also explore some of the unique strengths of NMR-based metabolomics, particularly with regard to isotope selection/detection, mixture deconvolution via 2D spectroscopy, automation, and the ability to noninvasively analyze native tissue specimens. Finally, this review will highlight a number of emerging NMR techniques and technologies that are being used to strengthen its utility and overcome its inherent limitations in metabolomic applications.

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.

The ability to restore homeostasis upon environmental challenges has been proposed as a measure for health. Metabolic profiling of plasma samples during the challenge response phase should offer a profound view on the flexibility of a phenotype to cope with daily stressors. Current data modeling approaches, however, struggle to extract biological descriptors from time-resolved metabolite profiles that are able to discriminate between different phenotypes. Thus, for the case of oxylipin responses in plasma upon an oral glucose tolerance test we developed a modeling approach that incorporates a priori biological pathway knowledge. The degradation pathways of arachidonic and linoleic acids were modeled using a regression model based on a pseudo-steady-state approximated model, resulting in a parameter A that summarizes the relative enzymatic activity in these pathways. Analysis of the phenotypic parameters As suggests that different phenotypes can be discriminated according to preferred relative activity of the arachidonic and linoleic pathway. Correlation analysis shows that there is little or no competition between the arachidonic and linoleic acid pathways, although they share the same enzymes.

While in situ chemical oxidation is often used to remediate tetrachloroethene (PCE) contaminated locations, very little is known about its influence on microbial composition and organohalide respiration (OHR) activity. Here, we investigate the impact of oxidation with permanganate on OHR rates, the abundance of organohalide respiring bacteria (OHRB) and reductive dehalogenase (rdh) genes using quantitative PCR, and microbial community composition through sequencing of 16S rRNA genes. A PCE degrading enrichment was repeatedly treated with low (25 μmol), medium (50 μmol), or high (100 μmol) permanganate doses, or no oxidant treatment (biotic control). Low and medium treatments led to higher OHR rates and enrichment of several OHRB and rdh genes, as compared to the biotic control. Improved degradation rates can be attributed to enrichment of (1) OHRB able to also utilize Mn oxides as a terminal electron acceptor and (2) non-dechlorinating community members of the Clostridiales and Deltaproteobacteria possibly supporting OHRB by providing essential co-factors. In contrast, high permanganate treatment disrupted dechlorination beyond cis-dichloroethene and caused at least a 2-4 orders of magnitude reduction in the abundance of all measured OHRB and rdh genes, as compared to the biotic control. High permanganate treatments resulted in a notably divergent microbial community, with increased abundances of organisms affiliated with Campylobacterales and Oceanospirillales capable of dissimilatory Mn reduction, and decreased abundance of presumed supporters of OHRB. Although OTUs classified within the OHR-supportive order Clostridiales and OHRB increased in abundance over the course of 213 days following the final 100 μmol permanganate treatment, only limited regeneration of PCE dechlorination was observed in one of three microcosms, suggesting strong chemical oxidation treatments can irreversibly disrupt OHR. Overall, this detailed investigation into dose-dependent changes of microbial composition and activity due to permanganate treatment provides insight into the mechanisms of OHR stimulation or disruption upon chemical oxidation.

Understanding the diversity and robustness of the metabolism of bacteria is fundamental for understanding how bacteria evolve and adapt to different environments. In this study, we characterised 121 Streptococcus strains and studied metabolic diversity from a protein domain perspective. Metabolic pathways were described in terms of the promiscuity of domains participating in metabolic pathways that were inferred to be functional. Promiscuity was defined by adapting existing measures based on domain abundance and versatility. The approach proved to be successful in capturing bacterial metabolic flexibility and species diversity, indicating that it can be described in terms of reuse and sharing functional domains in different proteins involved in metabolic activity. Additionally, we showed striking differences among metabolic organisation of the pathogenic serotype 2 Streptococcus suis and other strains.

Metabolomics is an “omic” science that is now emerging with the purpose of elaborating a comprehensive analysis of the metabolome, which is the complete set of metabolites (i.e., small molecules intermediates) in an organism, tissue, cell, or biofluid. In the past decade, metabolomics has already proved to be useful for the characterization of several pathological conditions and offers promises as a clinical tool. A metabolomics investigation of coeliac disease (CD) revealed that a metabolic fingerprint for CD can be defined, which accounts for three different but complementary components: malabsorption, energy metabolism, and alterations in gut microflora and/or intestinal permeability. In this review, we will discuss the major advancements in metabolomics of CD, in particular with respect to the role of gut microbiome and energy metabolism.

BACKGROUND: Metabolomics has attracted the interest of the medical community for its potential in predicting early derangements from a healthy to a diseased metabolic phenotype. One key issue is the diversity observed in metabolic profiles of different healthy individuals, commonly attributed to the variation of intrinsic (such as (epi)genetic variation, gut microbiota, etc.) and extrinsic factors (such as dietary habits, life-style and environmental conditions). Understanding the relative contributions of these factors is essential to establish the robustness of the healthy individual metabolic phenotype. METHODS: To assess the relative contribution of intrinsic and extrinsic factors we compared multilevel analysis results obtained from subjects of Homo sapiens and Macaca mulatta, the latter kept in a controlled environment with a standardized diet by making use of previously published data and results. RESULTS: We observed similarities for the two species and found the diversity of urinary metabolic phenotypes as identified by nuclear magnetic resonance (NMR) spectroscopy could be ascribed to the complex interplay of intrinsic factors and, to a lesser extent, of extrinsic factors in particular minimizing the role played by diet in shaping the metabolic phenotype. Moreover, we show that despite the standardization of diet as the most relevant extrinsic factor, a clear individual and discriminative metabolic fingerprint also exists for monkeys. We investigate the metabolic phenotype both at the static (i.e., at the level of the average metabolite concentration) and at the dynamic level (i.e., concerning their variation over time), and we show that these two components sum up to the overall phenotype with different relative contributions of about 1/4 and 3/4, respectively, for both species. Finally, we show that the great degree diversity observed in the urinary metabolic phenotype of both species can be attributed to differences in both the static and dynamic part of their phenotype.

One of the first steps in analyzing high-dimensional functional genomics data is an exploratory analysis of such data. Cluster Analysis and Principal Component Analysis are then usually the method of choice. Despite their versatility they also have a severe drawback: they do not always generate simple and interpretable solutions. On the basis of the observation that functional genomics data often contain both informative and non-informative variation, we propose a method that finds sets of variables containing informative variation. This informative variation is subsequently expressed in easily interpretable simplivariate components.We present a new implementation of the recently introduced simplivariate models. In this implementation, the informative variation is described by multiplicative models that can adequately represent the relations between functional genomics data. Both a simulated and two real-life metabolomics data sets show good performance of the method.