Ethanol-lactate transition of Lachancea thermotolerans is linked to nitrogen metabolism.

Climate change increases sugar content in grapes, resulting in unwanted increase in ethanol content of wine. Lachancea thermotolerans ferments glucose and fructose into both ethanol and lactate, decreasing final ethanol content and positively affecting wine acidity. Reported Lachancea thermotolerans strains show big variation in lactate production during fermentation. However, a mechanistic understanding of this lactate producing phenotype is currently lacking. Through a combination of metabolomics, transcriptomics, genomics and computational methods we show that the lactate production is induced by amino acid limitation in a high lactate producing strain. We found in fermentations in synthetic grape juice media that lactate production starts in the last stages of growth, marked by decreased growth rate and increased expression levels of stress related genes. This onset of lactate production is specific for the high lactate producing strain and independent of oxygen availability. The onset of lactate production was changed by increased amino acid content of the media, and it is shown by both computational methods and amino acid measurements that at the onset of lactate production amino acids become limiting for growth. This study shows that lactate production of Lachancea thermotolerans is directly linked to nitrogen availability in the media, an insight that can further aid in the improvement of wine quality.

Proteome constraints reveal targets for improving microbial fitness in nutrient-rich environments.

Cells adapt to different conditions via gene expression that tunes metabolism for maximal fitness. Constraints on cellular proteome may limit such expression strategies and introduce trade-offs. Resource allocation under proteome constraints has explained regulatory strategies in bacteria. It is unclear, however, to what extent these constraints can predict evolutionary changes, especially for microorganisms that evolved under nutrient-rich conditions, i.e., multiple available nitrogen sources, such as Lactococcus lactis. Here, we present a proteome-constrained genome-scale metabolic model of L. lactis (pcLactis) to interpret growth on multiple nutrients. Through integration of proteomics and flux data, in glucose-limited chemostats, the model predicted glucose and arginine uptake as dominant constraints at low growth rates. Indeed, glucose and arginine catabolism were found upregulated in evolved mutants. At high growth rates, pcLactis correctly predicted the observed shutdown of arginine catabolism because limited proteome availability favored lactate for ATP production. Thus, our model-based analysis is able to identify and explain the proteome constraints that limit growth rate in nutrient-rich environments and thus form targets of fitness improvement.

Genetic Elements Orchestrating Lactobacillus crispatus Glycogen Metabolism in the Vagina.

Glycogen in the female lower reproductive tract is a major carbon source for colonization and acidification by common vaginal Lactobacillus species, such as Lactobacillus crispatus. Previously, we identified the amylopullulanase encoding gene pulA of Lactobacillus crispatus to correlate with the ability to autonomously utilize glycogen for growth. Here, we further characterize genetic variation and differential regulation of pulA affecting the presence of its gene product on the outer surface layer. We show that alpha-glucan degrading activity dissipates when Lactobacillus crispatus is grown on glucose, maltose and maltotriose, in agreement with carbon catabolite repression elements flanking the pulA gene. Proteome analysis of the S-layer confirmed that the amylopullulanase protein is highly abundant in an S-layer enriched fraction, but not in a strain with a defective amylopullulanase variant or in an amylopullulanase-sufficient strain grown on glucose. In addition, we provide evidence that Lactobacillus crispatus pulA mutants are relevant in vivo, as they are commonly observed in metagenome datasets of human vaginal microbial communities. Analysis of the largest publicly available dataset of 1507 human vaginal metagenomes indicates that among the 270 samples that contain a Lactobacillus crispatuspulA gene, 62 samples (23%) had a defective variant of this gene. Taken together, these results demonstrate that both environmental, as well as genetic factors explain the variation of Lactobacillus crispatus alpha-glucosidases in the vaginal environment.

Adaption to glucose limitation is modulated by the pleotropic regulator CcpA, independent of selection pressure strength.

BACKGROUND: A central theme in (micro)biology is understanding the molecular basis of fitness i.e. which strategies are successful under which conditions; how do organisms implement such strategies at the molecular level; and which constraints shape the trade-offs between alternative strategies. Highly standardized microbial laboratory evolution experiments are ideally suited to approach these questions. For example, prolonged chemostats provide a constant environment in which the growth rate can be set, and the adaptive process of the organism to such environment can be subsequently characterized. RESULTS: We performed parallel laboratory evolution of Lactococcus lactis in chemostats varying the quantitative value of the selective pressure by imposing two different growth rates. A mutation in one specific amino acid residue of the global transcriptional regulator of carbon metabolism, CcpA, was selected in all of the evolution experiments performed. We subsequently showed that this mutation confers predictable fitness improvements at other glucose-limited growth rates as well. In silico protein structural analysis of wild type and evolved CcpA, as well as biochemical and phenotypic assays, provided the underpinning molecular mechanisms that resulted in the specific reprogramming favored in constant environments. CONCLUSION: This study provides a comprehensive understanding of a case of microbial evolution and hints at the wide dynamic range that a single fitness-enhancing mutation may display. It demonstrates how the modulation of a pleiotropic regulator can be used by cells to improve one trait while simultaneously work around other limiting constraints, by fine-tuning the expression of a wide range of cellular processes.

Training for translation between disciplines: a philosophy for life and data sciences curricula.

Motivation: Our society has become data-rich to the extent that research in many areas has become impossible without computational approaches. Educational programmes seem to be lagging behind this development. At the same time, there is a growing need not only for strong data science skills, but foremost for the ability to both translate between tools and methods on the one hand, and application and problems on the other. Results: Here we present our experiences with shaping and running a masters' programme in bioinformatics and systems biology in Amsterdam. From this, we have developed a comprehensive philosophy on how translation in training may be achieved in a dynamic and multidisciplinary research area, which is described here. We furthermore describe two requirements that enable translation, which we have found to be crucial: sufficient depth and focus on multidisciplinary topic areas, coupled with a balanced breadth from adjacent disciplines. Finally, we present concrete suggestions on how this may be implemented in practice, which may be relevant for the effectiveness of life science and data science curricula in general, and of particular interest to those who are in the process of setting up such curricula. Supplementary information: Supplementary data are available at Bioinformatics online.

Protein costs do not explain evolution of metabolic strategies and regulation of ribosomal content: does protein investment explain an anaerobic bacterial Crabtree effect?

Protein investment costs are considered a major driver for the choice of alternative metabolic strategies. We tested this premise in Lactococcus lactis, a bacterium that exhibits a distinct, anaerobic version of the bacterial Crabtree/Warburg effect; with increasing growth rates it shifts from a high yield metabolic mode [mixed-acid fermentation; 3 adenosine triphosphate (ATP) per glucose] to a low yield metabolic mode (homolactic fermentation; 2 ATP per glucose). We studied growth rate-dependent relative transcription and protein ratios, enzyme activities, and fluxes of L. lactis in glucose-limited chemostats, providing a high-quality and comprehensive data set. A three- to fourfold higher growth rate rerouted metabolism from acetate to lactate as the main fermentation product. However, we observed hardly any changes in transcription, protein levels and enzyme activities. Even levels of ribosomal proteins, constituting a major investment in cellular machinery, changed only slightly. Thus, contrary to the original hypothesis, central metabolism in this organism appears to be hardly regulated at the level of gene expression, but rather at the metabolic level. We conclude that L. lactis is either poorly adapted to growth at low and constant glucose concentrations, or that protein costs play a less important role in fitness than hitherto assumed.

Molecular assessment of bacterial vaginosis by Lactobacillus abundance and species diversity.

BACKGROUND: To date, women are most often diagnosed with bacterial vaginosis (BV) using microscopy based Nugent scoring or Amsel criteria. However, the accuracy is less than optimal. The aim of the present study was to confirm the identity of known BV-associated composition profiles and evaluate indicators for BV using three molecular methods. METHODS: Evaluation of indicators for BV was carried out by 16S rRNA amplicon sequencing of the V5-V7 region, a tailor-made 16S rRNA oligonucleotide-based microarray, and a PCR-based profiling technique termed IS-profiling, which is based on fragment variability of the 16S-23S rRNA intergenic spacer region. An inventory of vaginal bacterial species was obtained from 40 females attending a Dutch sexually transmitted infection outpatient clinic, of which 20 diagnosed with BV (Nugent score 7-10), and 20 BV negative (Nugent score 0-3). RESULTS: Analysis of the bacterial communities by 16S rRNA amplicon sequencing revealed two clusters in the BV negative women, dominated by either Lactobacillus iners or Lactobacillus crispatus and three distinct clusters in the BV positive women. In the former, there was a virtually complete, negative correlation between L. crispatus and L. iners. BV positive subjects showed cluster profiles that were relatively high in bacterial species diversity and dominated by anaerobic species, including Gardnerella vaginalis, and those belonging to the Families of Lachnospiraceae and Leptotrichiaceae. Accordingly, the Gini-Simpson index of species diversity, and the relative abundance Lactobacillus species appeared consistent indicators for BV. Under the conditions used, only the 16S rRNA amplicon sequencing method was suitable to assess species diversity, while all three molecular composition profiling methods were able to indicate Lactobacillus abundance in the vaginal microbiota. CONCLUSION: An affordable and simple molecular test showing a depletion of the genus Lactobacillus in combination with an increased species diversity of vaginal microbiota could serve as an alternative and practical diagnostic method for the assessment of BV.

Availability of public goods shapes the evolution of competing metabolic strategies.

Tradeoffs provide a rationale for the outcome of natural selection. A prominent example is the negative correlation between the growth rate and the biomass yield in unicellular organisms. This tradeoff leads to a dilemma, where the optimization of growth rate is advantageous for an individual, whereas the optimization of the biomass yield would be advantageous for a population. High-rate strategies are observed in a broad variety of organisms such as Escherichia coli, yeast, and cancer cells. Growth in suspension cultures favors fast-growing organisms, whereas spatial structure is of importance for the evolution of high-yield strategies. Despite this realization, experimental methods to directly select for increased yield are lacking. We here show that the serial propagation of a microbial population in a water-in-oil emulsion allows selection of strains with increased biomass yield. The propagation in emulsion creates a spatially structured environment where the growth-limiting substrate is privatized for populations founded by individual cells. Experimental evolution of several isogenic Lactococcus lactis strains demonstrated the existence of a tradeoff between growth rate and biomass yield as an apparent Pareto front. The underlying mutations altered glucose transport and led to major shifts between homofermentative and heterofermentative metabolism, accounting for the changes in metabolic efficiency. The results demonstrated the impact of privatizing a public good on the evolutionary outcome between competing metabolic strategies. The presented approach allows the investigation of fundamental questions in biology such as the evolution of cooperation, cell-cell interactions, and the relationships between environmental and metabolic constraints.

Genome-wide transposon mutagenesis indicates that Mycobacterium marinum customizes its virulence mechanisms for survival and replication in different hosts.

The interaction of environmental bacteria with unicellular eukaryotes is generally considered a major driving force for the evolution of intracellular pathogens, allowing them to survive and replicate in phagocytic cells of vertebrate hosts. To test this hypothesis on a genome-wide level, we determined for the intracellular pathogen Mycobacterium marinum whether it uses conserved strategies to exploit host cells from both protozoan and vertebrate origin. Using transposon-directed insertion site sequencing (TraDIS), we determined differences in genetic requirements for survival and replication in phagocytic cells of organisms from different kingdoms. In line with the general hypothesis, we identified a number of general virulence mechanisms, including the type VII protein secretion system ESX-1, biosynthesis of polyketide lipids, and utilization of sterols. However, we were also able to show that M. marinum contains an even larger set of host-specific virulence determinants, including proteins involved in the modification of surface glycolipids and, surprisingly, the auxiliary proteins of the ESX-1 system. Several of these factors were in fact counterproductive in other hosts. Therefore, M. marinum contains different sets of virulence factors that are tailored for specific hosts. Our data imply that although amoebae could function as a training ground for intracellular pathogens, they do not fully prepare pathogens for crossing species barriers.

Microbial domestication signatures of Lactococcus lactis can be reproduced by experimental evolution.

Experimental evolution is a powerful approach to unravel how selective forces shape microbial genotypes and phenotypes. To this date, the available examples focus on the adaptation to conditions specific to the laboratory. The lactic acid bacterium Lactococcus lactis naturally occurs on plants and in dairy environments, and it is proposed that dairy strains originate from the plant niche. Here we investigate the adaptation of a L. lactis strain isolated from a plant to a dairy niche by propagating it for 1000 generations in milk. Two out of three independently evolved strains displayed significantly increased acidification rates and biomass yields in milk. Genome resequencing, revealed six, seven, and 28 mutations in the three strains, including point mutations in loci related to amino acid biosynthesis and transport and in the gene encoding MutL, which is involved in DNA mismatch repair. Two strains lost a conjugative transposon containing genes important in the plant niche but dispensable in milk. A plasmid carrying an extracellular protease was introduced by transformation. Although improving growth rate and growth yield significantly, the plasmid was rapidly lost. Comparative transcriptome and phenotypic analyses confirmed that major physiological changes associated with improved growth in milk relate to nitrogen metabolism and the loss or down-regulation of several pathways involved in the utilization of complex plant polymers. Reproducing the transition from the plant to the dairy niche through experimental evolution revealed several genome, transcriptome, and phenotype signatures that resemble those seen in strains isolated from either niche.

Bioinformatics and systems biology: bridging the gap between heterogeneous student backgrounds.

Teaching students with very diverse backgrounds can be extremely challenging. This article uses the Bioinformatics and Systems Biology MSc in Amsterdam as a case study to describe how the knowledge gap for students with heterogeneous backgrounds can be bridged. We show that a mix in backgrounds can be turned into an advantage by creating a stimulating learning environment for the students. In the MSc Programme, conversion classes help to bridge differences between students, by mending initial knowledge and skill gaps. Mixing students from different backgrounds in a group to solve a complex task creates an opportunity for the students to reflect on their own abilities. We explain how a truly interdisciplinary approach to teaching helps students of all backgrounds to achieve the MSc end terms. Moreover, transferable skills obtained by the students in such a mixed study environment are invaluable for their later careers.

Genome instability in Lactobacillus rhamnosus GG.

We describe here a comparative genome analysis of three dairy product isolates of Lactobacillus rhamnosus GG (LGG) and the ATCC 53103 reference strain to the published genome sequence of L. rhamnosus GG. The analysis showed that in two of three isolates, major DNA segments were missing from the genomic islands LGGISL1,2. The deleted DNA segments consist of 34 genes in one isolate and 84 genes in the other and are flanked by identical insertion elements. Among the missing genes are the spaCBA genes, which encode pilin subunits involved in adhesion to mucus and persistence of the strains in the human intestinal tract. Subsequent quantitative PCR analyses of six commercial probiotic products confirmed that two more products contain a heterogeneous population of L. rhamnosus GG variants, including genotypes with or without spaC. These results underline the relevance for quality assurance and control measures targeting genome stability in probiotic strains and justify research assessing the effect of genetic rearrangements in probiotics on the outcome of in vitro and in vivo efficacy studies.

Vesicle trafficking via the Spitzenkörper during hyphal tip growth in Rhizoctonia solani.

Growing hyphae of Rhizoctonia solani were stained with the endocytic marker dye FM4-64 and imaged by confocal microscopy. Staining of the plasma membrane was followed by labeling of organelles in the cytoplasm (after ~1 min) and of the Spitzenkörper (Spk; after ~2 min). Fluorescence recovery after photobleaching (FRAP) of the stained Spk demonstrated the vectorial flow of secretory vesicles from the apical cytoplasm to the Spk. This flux was modelled in a two-compartment model. The turnover time of the vesicles of the Spk was estimated to be 1.3-2.5 min. These results are roughly consistent with the expected flux of vesicles through the Spk based on the number of secretory vesicles within the Spk and the number of secretory vesicles that would be necessary to fuse with the apical plasma membrane to maintain hyphal extension rates. These results suggest that membrane retrieval via endocytosis is not as significant as previously suggested.

Identifying antimicrobials and their metabolites in wastewater and surface water with effect-directed analysis.

This study aimed to identify antimicrobial contaminants in the aquatic environment with effect-directed analysis. Wastewater influent, effluent, and surface water (up- and downstream of the discharge location) were sampled at two study sites. The samples were enriched, subjected to high-resolution fractionation, and the resulting 80 fractions were tested in an antibiotics bioassay. The resulting bioactive fractions guided the suspect and nontargeted identification strategy in the high-resolution mass spectrometry data that was recorded in parallel. Chemical features were annotated with reference databases, assessed on annotation quality, and assigned identification confidence levels. To identify antibiotic metabolites, Phase I metabolites were predicted in silico for over 500 antibiotics and included as a suspect list. Predicted retention times and fragmentation patterns reduced the number of annotations to consider for confirmation testing. Overall, the bioactivity of three fractions could be explained by the identified antibiotics (clarithromycin and azithromycin) and an antibiotic metabolite (14-OH(R) clarithromycin), explaining 78% of the bioactivity measured at one study site. The applied identification strategy successfully identified antibiotic metabolites in the aquatic environment, emphasizing the need to include the toxic effects of bioactive metabolites in environmental risk assessments.

Regional vulnerability of brain white matter in vanishing white matter.

Vanishing white matter (VWM) is a leukodystrophy that primarily manifests in young children. In this disease, the brain white matter is differentially affected in a predictable pattern with telencephalic brain areas being most severely affected, while others remain allegedly completely spared. Using high-resolution mass spectrometry-based proteomics, we investigated the proteome patterns of the white matter in the severely affected frontal lobe and normal appearing pons in VWM and control cases to identify molecular bases underlying regional vulnerability. By comparing VWM patients to controls, we identified disease-specific proteome patterns. We showed substantial changes in both the VWM frontal and pons white matter at the protein level. Side-by-side comparison of brain region-specific proteome patterns further revealed regional differences. We found that different cell types were affected in the VWM frontal white matter than in the pons. Gene ontology and pathway analyses identified involvement of region specific biological processes, of which pathways involved in cellular respiratory metabolism were overarching features. In the VWM frontal white matter, proteins involved in glycolysis/gluconeogenesis and metabolism of various amino acids were decreased compared to controls. By contrast, in the VWM pons white matter, we found a decrease in proteins involved in oxidative phosphorylation. Taken together, our data show that brain regions are affected in parallel in VWM, but to different degrees. We found region-specific involvement of different cell types and discovered that cellular respiratory metabolism is likely to be differentially affected across white matter regions in VWM. These region-specific changes help explain regional vulnerability to pathology in VWM.

PhenoLink--a web-tool for linking phenotype to ~omics data for bacteria: application to gene-trait matching for Lactobacillus plantarum strains.

BACKGROUND: Linking phenotypes to high-throughput molecular biology information generated by ~omics technologies allows revealing cellular mechanisms underlying an organism's phenotype. ~Omics datasets are often very large and noisy with many features (e.g., genes, metabolite abundances). Thus, associating phenotypes to ~omics data requires an approach that is robust to noise and can handle large and diverse data sets. RESULTS: We developed a web-tool PhenoLink (http://bamics2.cmbi.ru.nl/websoftware/phenolink/) that links phenotype to ~omics data sets using well-established as well new techniques. PhenoLink imputes missing values and preprocesses input data (i) to decrease inherent noise in the data and (ii) to counterbalance pitfalls of the Random Forest algorithm, on which feature (e.g., gene) selection is based. Preprocessed data is used in feature (e.g., gene) selection to identify relations to phenotypes. We applied PhenoLink to identify gene-phenotype relations based on the presence/absence of 2847 genes in 42 Lactobacillus plantarum strains and phenotypic measurements of these strains in several experimental conditions, including growth on sugars and nitrogen-dioxide production. Genes were ranked based on their importance (predictive value) to correctly predict the phenotype of a given strain. In addition to known gene to phenotype relations we also found novel relations. CONCLUSIONS: PhenoLink is an easily accessible web-tool to facilitate identifying relations from large and often noisy phenotype and ~omics datasets. Visualization of links to phenotypes offered in PhenoLink allows prioritizing links, finding relations between features, finding relations between phenotypes, and identifying outliers in phenotype data. PhenoLink can be used to uncover phenotype links to a multitude of ~omics data, e.g., gene presence/absence (determined by e.g.: CGH or next-generation sequencing), gene expression (determined by e.g.: microarrays or RNA-seq), or metabolite abundance (determined by e.g.: GC-MS).

Standardized assay medium to measure Lactococcus lactis enzyme activities while mimicking intracellular conditions.

Knowledge of how the activity of enzymes is affected under in vivo conditions is essential for analyzing their regulation and constructing models that yield an integrated understanding of cell behavior. Current kinetic parameters for Lactococcus lactis are scattered through different studies and performed under different assay conditions. Furthermore, assay conditions often diverge from conditions prevailing in the intracellular environment. To establish uniform assay conditions that resemble intracellular conditions, we analyzed the intracellular composition of anaerobic glucose-limited chemostat cultures of L. lactis subsp. cremoris MG 1363. Based on this, we designed a new assay medium for enzyme activity measurements of growing cells of L. lactis, mimicking as closely as practically possible its intracellular environment. Procedures were optimized to be carried out in 96-well plates, and the reproducibility and dynamic range were checked for all enzyme activity measurements. The effects of freezing and the carryover of ammonium sulfate from the addition of coupling enzymes were also established. Activities of all 10 glycolytic and 4 fermentative enzymes were measured. Remarkably, most in vivo-like activities were lower than previously published data. Yet, the ratios of V(max) over measured in vivo fluxes were above 1. With this work, we have developed and extensively validated standard protocols for enzyme activity measurements for L. lactis.

Metabolic shifts: a fitness perspective for microbial cell factories.

Performance of industrial microorganisms as cell factories is limited by the capacity to channel nutrients to desired products, of which optimal production usually requires careful manipulation of process conditions, or strain improvement. The focus in process improvement is often on understanding and manipulating the regulation of metabolism. Nonetheless, one encounters situations where organisms are remarkably resilient to further optimization or their properties become unstable. Therefore it is important to understand the origin of these apparent limitations to find whether and how they can be improved. We argue that by considering fitness effects of regulation, a more generic explanation for certain behaviour can be obtained. In this view, apparent process limitations arise from trade-offs that cells faced as they evolved to improve fitness. A deeper understanding of such trade-offs using a systems biology approach can ultimately enhance performance of cell factories.

Cortical Pathology in Vanishing White Matter.

Vanishing white matter (VWM) is classified as a leukodystrophy with astrocytes as primary drivers in its pathogenesis. Magnetic resonance imaging has documented the progressive thinning of cortices in long-surviving patients. Routine histopathological analyses, however, have not yet pointed to cortical involvement in VWM. Here, we provide a comprehensive analysis of the VWM cortex. We employed high-resolution-mass-spectrometry-based proteomics and immunohistochemistry to gain insight into possible molecular disease mechanisms in the cortices of VWM patients. The proteome analysis revealed 268 differentially expressed proteins in the VWM cortices compared to the controls. A majority of these proteins formed a major protein interaction network. A subsequent gene ontology analysis identified enrichment for terms such as cellular metabolism, particularly mitochondrial activity. Importantly, some of the proteins with the most prominent changes in expression were found in astrocytes, indicating cortical astrocytic involvement. Indeed, we confirmed that VWM cortical astrocytes exhibit morphological changes and are less complex in structure than control cells. Our findings also suggest that these astrocytes are immature and not reactive. Taken together, we provide insights into cortical involvement in VWM, which has to be taken into account when developing therapeutic strategies.

Volatile compound fingerprinting of mixed-culture fermentations.

With the advent of the -omics era, classical technology platforms, such as hyphenated mass spectrometry, are currently undergoing a transformation toward high-throughput application. These novel platforms yield highly detailed metabolite profiles in large numbers of samples. Such profiles can be used as fingerprints for the accurate identification and classification of samples as well as for the study of effects of experimental conditions on the concentrations of specific metabolites. Challenges for the application of these methods lie in the acquisition of high-quality data, data normalization, and data mining. Here, a high-throughput fingerprinting approach based on analysis of headspace volatiles using ultrafast gas chromatography coupled to time of flight mass spectrometry (ultrafast GC/TOF-MS) was developed and evaluated for classification and screening purposes in food fermentation. GC-MS mass spectra of headspace samples of milk fermented by different mixed cultures of lactic acid bacteria (LAB) were collected and preprocessed in MetAlign, a dedicated software package for the preprocessing and comparison of liquid chromatography (LC)-MS and GC-MS data. The Random Forest algorithm was used to detect mass peaks that discriminated combinations of species or strains used in fermentations. Many of these mass peaks originated from key flavor compounds, indicating that the presence or absence of individual strains or combinations of strains significantly influenced the concentrations of these components. We demonstrate that the approach can be used for purposes like the selection of strains from collections based on flavor characteristics and the screening of (mixed) cultures for the presence or absence of strains. In addition, we show that strain-specific flavor characteristics can be traced back to genetic markers when comparative genome hybridization (CGH) data are available.