Dries De Maeyer

Adaptation to High Ethanol Reveals Complex Evolutionary Pathways

Tolerance to high levels of ethanol is an ecologically and industrially relevant phenotype of microbes, but the molecular mechanisms underlying this complex trait remain largely unknown. Here, we use long-term experimental evolution of isogenic yeast populations of different initial ploidy to study adaptation to increasing levels of ethanol. Whole-genome sequencing of more than 30 evolved populations and over 100 adapted clones isolated throughout this two-year evolution experiment revealed how a complex interplay of de novo single nucleotide mutations, copy number variation, ploidy changes, mutator phenotypes, and clonal interference led to a significant increase in ethanol tolerance. Although the specific mutations differ between different evolved lineages, application of a novel computational pipeline, PheNetic, revealed that many mutations target functional modules involved in stress response, cell cycle regulation, DNA repair and respiration. Measuring the fitness effects of selected mutations introduced in non-evolved ethanol-sensitive cells revealed several adaptive mutations that had previously not been implicated in ethanol tolerance, including mutations in PRT1, VPS70 and MEX67. Interestingly, variation in VPS70 was recently identified as a QTL for ethanol tolerance in an industrial bio-ethanol strain. Taken together, our results show how, in contrast to adaptation to some other stresses, adaptation to a continuous complex and severe stress involves interplay of different evolutionary mechanisms. In addition, our study reveals functional modules involved in ethanol resistance and identifies several mutations that could help to improve the ethanol tolerance of industrial yeasts.

Identification of a complex genetic network underlying Saccharomyces cerevisiae colony morphology

When grown on solid substrates, different microorganisms often form colonies with very specific morphologies. Whereas the pioneers of microbiology often used colony morphology to discriminate between species and strains, the phenomenon has not received much attention recently. In this study, we use a genome‐wide assay in the model yeast Saccharomyces cerevisiae to identify all genes that affect colony morphology. We show that several major signalling cascades, including the MAPK, TORC, SNF1 and RIM101 pathways play a role, indicating that morphological changes are a reaction to changing environments. Other genes that affect colony morphology are involved in protein sorting and epigenetic regulation. Interestingly, the screen reveals only few genes that are likely to play a direct role in establishing colony morphology, with one notable example being FLO11, a gene encoding a cell‐surface adhesin that has already been implicated in colony morphology, biofilm formation, and invasive and pseudohyphal growth. Using a series of modified promoters for fine‐tuning FLO11 expression, we confirm the central role of Flo11 and show that differences in FLO11 expression result in distinct colony morphologies. Together, our results provide a first comprehensive look at the complex genetic network that underlies the diversity in the morphologies of yeast colonies.

Dynamics of the Saccharomyces cerevisiae transcriptome during bread dough fermentation.

ABSTRACT The behavior of yeast cells during industrial processes such as the production of beer, wine, and bioethanol has been extensively studied. In contrast, our knowledge about yeast physiology during solid-state processes, such as bread dough, cheese, or cocoa fermentation, remains limited. We investigated changes in the transcriptomes of three genetically distinct Saccharomyces cerevisiae strains during bread dough fermentation. Our results show that regardless of the genetic background, all three strains exhibit similar changes in expression patterns. At the onset of fermentation, expression of glucose-regulated genes changes dramatically, and the osmotic stress response is activated. The middle fermentation phase is characterized by the induction of genes involved in amino acid metabolism. Finally, at the latest time point, cells suffer from nutrient depletion and activate pathways associated with starvation and stress responses. Further analysis shows that genes regulated by the high-osmolarity glycerol (HOG) pathway, the major pathway involved in the response to osmotic stress and glycerol homeostasis, are among the most differentially expressed genes at the onset of fermentation. More importantly, deletion of HOG1 and other genes of this pathway significantly reduces the fermentation capacity. Together, our results demonstrate that cells embedded in a solid matrix such as bread dough suffer severe osmotic stress and that a proper induction of the HOG pathway is critical for optimal fermentation.

PheNetic: Network-based interpretation of unstructured gene lists in E. coli

At the present time, omics experiments are commonly used in wet lab practice to identify leads involved in interesting phenotypes. These omics experiments often result in unstructured gene lists, the interpretation of which in terms of pathways or the mode of action is challenging. To aid in the interpretation of such gene lists, we developed PheNetic, a decision theoretic method that exploits publicly available information, captured in a comprehensive interaction network to obtain a mechanistic view of the listed genes. PheNetic selects from an interaction network the sub-networks highlighted by these gene lists. We applied PheNetic to an Escherichia coli interaction network to reanalyse a previously published KO compendium, assessing gene expression of 27 E. coli knock-out mutants under mild acidic conditions. Being able to unveil previously described mechanisms involved in acid resistance demonstrated both the performance of our method and the added value of our integrated E. coli network. PheNetic is available at .

Elucidation of the Mode of Action of a New Antibacterial Compound Active against Staphylococcus aureus and Pseudomonas aeruginosa

Nosocomial and community-acquired infections caused by multidrug resistant bacteria represent a major human health problem. Thus, there is an urgent need for the development of antibiotics with new modes of action. In this study, we investigated the antibacterial characteristics and mode of action of a new antimicrobial compound, SPI031 (N-alkylated 3, 6-dihalogenocarbazol 1-(sec-butylamino)-3-(3,6-dichloro-9H-carbazol-9-yl)propan-2-ol), which was previously identified in our group. This compound exhibits broad-spectrum antibacterial activity, including activity against the human pathogens Staphylococcus aureus and Pseudomonas aeruginosa. We found that SPI031 has rapid bactericidal activity (7-log reduction within 30 min at 4x MIC) and that the frequency of resistance development against SPI031 is low. To elucidate the mode of action of SPI031, we performed a macromolecular synthesis assay, which showed that SPI031 causes non-specific inhibition of macromolecular biosynthesis pathways. Liposome leakage and membrane permeability studies revealed that SPI031 rapidly exerts membrane damage, which is likely the primary cause of its antibacterial activity. These findings were supported by a mutational analysis of SPI031-resistant mutants, a transcriptome analysis and the identification of transposon mutants with altered sensitivity to the compound. In conclusion, our results show that SPI031 exerts its antimicrobial activity by causing membrane damage, making it an interesting starting point for the development of new antibacterial therapies.

Network-Based Analysis of eQTL Data to Prioritize Driver Mutations

In clonal systems, interpreting driver genes in terms of molecular networks helps understanding how these drivers elicit an adaptive phenotype. Obtaining such a network-based understanding depends on the correct identification of driver genes. In clonal systems, independent evolved lines can acquire a similar adaptive phenotype by affecting the same molecular pathways, a phenomenon referred to as parallelism at the molecular pathway level. This implies that successful driver identification depends on interpreting mutated genes in terms of molecular networks. Driver identification and obtaining a network-based understanding of the adaptive phenotype are thus confounded problems that ideally should be solved simultaneously. In this study, a network-based eQTL method is presented that solves both the driver identification and the network-based interpretation problem. As input the method uses coupled genotype-expression phenotype data (eQTL data) of independently evolved lines with similar adaptive phenotypes and an organism-specific genome-wide interaction network. The search for mutational consistency at pathway level is defined as a subnetwork inference problem, which consists of inferring a subnetwork from the genome-wide interaction network that best connects the genes containing mutations to differentially expressed genes. Based on their connectivity with the differentially expressed genes, mutated genes are prioritized as driver genes. Based on semisynthetic data and two publicly available data sets, we illustrate the potential of the network-based eQTL method to prioritize driver genes and to gain insights in the molecular mechanisms underlying an adaptive phenotype. The method is available at http://bioinformatics.intec.ugent.be/phenetic_eqtl/index.html

PheNetic: network-based interpretation of molecular profiling data

Molecular profiling experiments have become standard in current wet-lab practices. Classically, enrichment analysis has been used to identify biological functions related to these experimental results. Combining molecular profiling results with the wealth of currently available interactomics data, however, offers the opportunity to identify the molecular mechanism behind an observed molecular phenotype. In this paper, we therefore introduce ‘PheNetic’, a user-friendly web server for inferring a sub-network based on probabilistic logical querying. PheNetic extracts from an interactome, the sub-network that best explains genes prioritized through a molecular profiling experiment. Depending on its run mode, PheNetic searches either for a regulatory mechanism that gave explains to the observed molecular phenotype or for the pathways (in)activated in the molecular phenotype. The web server provides access to a large number of interactomes, making sub-network inference readily applicable to a wide variety of organisms. The inferred sub-networks can be interactively visualized in the browser. PheNetics method and use are illustrated using an example analysis of differential expression results of ampicillin treated Escherichia coli cells. The PheNetic web service is available at http://bioinformatics.intec.ugent.be/phenetic/.

Simultaneous discovery of cancer subtypes and subtype features by molecular data integration.

MOTIVATIONnSubtyping cancer is key to an improved and more personalized prognosis/treatment. The increasing availability of tumor related molecular data provides the opportunity to identify molecular subtypes in a data-driven way. Molecular subtypes are defined as groups of samples that have a similar molecular mechanism at the origin of the carcinogenesis. The molecular mechanisms are reflected by subtype-specific mutational and expression features. Data-driven subtyping is a complex problem as subtyping and identifying the molecular mechanisms that drive carcinogenesis are confounded problems. Many current integrative subtyping methods use global mutational and/or expression tumor profiles to group tumor samples in subtypes but do not explicitly extract the subtype-specific features. We therefore present a method that solves both tasks of subtyping and identification of subtype-specific features simultaneously. Hereto our method integrates` mutational and expression data while taking into account the clonal properties of carcinogenesis. Key to our method is a formalization of the problem as a rank matrix factorization of ranked data that approaches the subtyping problem as multi-view bi-clusteringnnnRESULTSnWe introduce a novel integrative framework to identify subtypes by combining mutational and expression features. The incomparable measurement data is integrated by transformation into ranked data and subtypes are defined as multi-view bi-clusters We formalize the model using rank matrix factorization, resulting in the SRF algorithm. Experiments on simulated data and the TCGA breast cancer data demonstrate that SRF is able to capture subtle differences that existing methods may miss.nnnAVAILABILITY AND IMPLEMENTATIONnThe implementation is available at: https://github.com/rankmatrixfactorisation/SRF CONTACT: kathleen.marchal@intec.ugent.be, siegfried.nijssen@cs.kuleuven.bennnSUPPLEMENTARY INFORMATIONnSupplementary data are available at Bioinformatics online.

SSA-ME Detection of cancer driver genes using mutual exclusivity by small subnetwork analysis

Because of its clonal evolution a tumor rarely contains multiple genomic alterations in the same pathway as disrupting the pathway by one gene often is sufficient to confer the complete fitness advantage. As a result, many cancer driver genes display mutual exclusivity across tumors. However, searching for mutually exclusive gene sets requires analyzing all possible combinations of genes, leading to a problem which is typically too computationally complex to be solved without a stringent a priori filtering, restricting the mutations included in the analysis. To overcome this problem, we present SSA-ME, a network-based method to detect cancer driver genes based on independently scoring small subnetworks for mutual exclusivity using a reinforced learning approach. Because of the algorithmic efficiency, no stringent upfront filtering is required. Analysis of TCGA cancer datasets illustrates the added value of SSA-ME: well-known recurrently mutated but also rarely mutated drivers are prioritized. We show that using mutual exclusivity to detect cancer driver genes is complementary to state-of-the-art approaches. This framework, in which a large number of small subnetworks are being analyzed in order to solve a computationally complex problem (SSA), can be generically applied to any problem in which local neighborhoods in a network hold useful information.

Path Finding in Biological Networks

Understanding the cellular behavior from a systems perspective requires the identification of functional and physical interactions among diverse molecular entities in a cell (i.e., DNA/RNA, proteins, and metabolites). The most straightforward way to represent such datasets is by means of molecular networks of which nodes correspond to molecular entities and edges to the interactions amongst those entities. Nowadays with large amounts of interaction data being generated, genome-wide networks can be created for an increasing number of organisms. These networks can be exploited to study a molecular entity like a protein in a wider context than just in isolation and provide a way of representing our knowledge of the system as a whole. On the other hand, viewing a single entity or an experimental dataset in the light of an interaction network can reveal previous unknown insights in biological processes.