Drug-specific T-cell responses in patients with liver injury following treatment with the BACE inhibitor atabecestat.

BACKGROUND: Atabecestat is an orally administered BACE inhibitor developed to treat Alzheimer's disease. Elevations in hepatic enzymes were detected in a number of in trial patients, which resulted in termination of the drug development programme. Immunohistochemical characterization of liver tissue from an index case of atabecestat-mediated liver injury revealed an infiltration of T-lymphocytes in areas of hepatocellular damage. This coupled with the fact that liver injury had a delayed onset suggests that the adaptive immune system may be involved in the pathogenesis. The aim of this study was to generate and characterize atabecestat(metabolite)-responsive T-cell clones from patients with liver injury. METHODS: Peripheral blood mononuclear cells were cultured with atabecestat and its metabolites (diaminothiazine [DIAT], N-acetyl DIAT & epoxide) and cloning was attempted in a number of patients. Atabecestat(metabolite)-responsive clones were analysed in terms of T-cell phenotype, function, pathways of T-cell activation and cross-reactivity with structurally related compounds. RESULTS: CD4+ T-cell clones activated with the DIAT metabolite were detected in 5 out of 8 patients (up to 4.5% cloning efficiency). Lower numbers of CD4+ and CD8+ clones displayed reactivity against atabecestat. Clones proliferated and secreted IFN-γ, IL-13 and cytolytic molecules following atabecestat or DIAT stimulation. Certain atabecestat and DIAT-responsive clones cross-reacted with N-acetyl DIAT; however, no cross-reactivity was observed between atabecestat and DIAT. CD4+ clones were activated through a direct, reversible compound-HLA class II interaction with no requirement for protein processing. CONCLUSION: The detection of atabecestat metabolite-responsive T-cell clones activated via a pharmacological interactions pathway in patients with liver injury is indicative of an immune-based mechanism for the observed hepatic enzyme elevations.

A scalable SCENIC workflow for single-cell gene regulatory network analysis.

This protocol explains how to perform a fast SCENIC analysis alongside standard best practices steps on single-cell RNA-sequencing data using software containers and Nextflow pipelines. SCENIC reconstructs regulons (i.e., transcription factors and their target genes) assesses the activity of these discovered regulons in individual cells and uses these cellular activity patterns to find meaningful clusters of cells. Here we present an improved version of SCENIC with several advances. SCENIC has been refactored and reimplemented in Python (pySCENIC), resulting in a tenfold increase in speed, and has been packaged into containers for ease of use. It is now also possible to use epigenomic track databases, as well as motifs, to refine regulons. In this protocol, we explain the different steps of SCENIC: the workflow starts from the count matrix depicting the gene abundances for all cells and consists of three stages. First, coexpression modules are inferred using a regression per-target approach (GRNBoost2). Next, the indirect targets are pruned from these modules using cis-regulatory motif discovery (cisTarget). Lastly, the activity of these regulons is quantified via an enrichment score for the regulon's target genes (AUCell). Nonlinear projection methods can be used to display visual groupings of cells based on the cellular activity patterns of these regulons. The results can be exported as a loom file and visualized in the SCope web application. This protocol is illustrated on two use cases: a peripheral blood mononuclear cell data set and a panel of single-cell RNA-sequencing cancer experiments. For a data set of 10,000 genes and 50,000 cells, the pipeline runs in <2 h.

IAMBEE: a web-service for the identification of adaptive pathways from parallel evolved clonal populations.

IAMBEE is a web server designed for the Identification of Adaptive Mutations in Bacterial Evolution Experiments (IAMBEE). Input data consist of genotype information obtained from independently evolved clonal populations or strains that show the same adapted behavior (phenotype). To distinguish adaptive from passenger mutations, IAMBEE searches for neighborhoods in an organism-specific interaction network that are recurrently mutated in the adapted populations. This search for recurrently mutated network neighborhoods, as proxies for pathways is driven by additional information on the functional impact of the observed genetic changes and their dynamics during adaptive evolution. In addition, the search explicitly accounts for the differences in mutation rate between the independently evolved populations. Using this approach, IAMBEE allows exploiting parallel evolution to identify adaptive pathways. The web-server is freely available at http://bioinformatics.intec.ugent.be/iambee/ with no login requirement.

The Crabtree Effect Shapes the Saccharomyces cerevisiae Lag Phase during the Switch between Different Carbon Sources.

When faced with environmental changes, microbes often enter a temporary growth arrest during which they reprogram the expression of specific genes to adapt to the new conditions. A prime example of such a lag phase occurs when microbes need to switch from glucose to other, less-preferred carbon sources. Despite its industrial relevance, the genetic network that determines the duration of the lag phase has not been studied in much detail. Here, we performed a genome-wide Bar-Seq screen to identify genetic determinants of the Saccharomyces cerevisiae glucose-to-galactose lag phase. The results show that genes involved in respiration, and specifically those encoding complexes III and IV of the electron transport chain, are needed for efficient growth resumption after the lag phase. Anaerobic growth experiments confirmed the importance of respiratory energy conversion in determining the lag phase duration. Moreover, overexpression of the central regulator of respiration, HAP4, leads to significantly shorter lag phases. Together, these results suggest that the glucose-induced repression of respiration, known as the Crabtree effect, is a major determinant of microbial fitness in fluctuating carbon environments.IMPORTANCE The lag phase is arguably one of the prime characteristics of microbial growth. Longer lag phases result in lower competitive fitness in variable environments, and the duration of the lag phase is also important in many industrial processes where long lag phases lead to sluggish, less efficient fermentations. Despite the immense importance of the lag phase, surprisingly little is known about the exact molecular processes that determine its duration. Our study uses the molecular toolbox of S. cerevisiae combined with detailed growth experiments to reveal how the transition from fermentative to respirative metabolism is a key bottleneck for cells to overcome the lag phase. Together, our findings not only yield insight into the key molecular processes and genes that influence lag duration but also open routes to increase the efficiency of industrial fermentations and offer an experimental framework to study other types of lag behavior.

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.

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

MOTIVATION: Subtyping 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-clustering RESULTS: We 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. AVAILABILITY AND IMPLEMENTATION: The implementation is available at: https://github.com/rankmatrixfactorisation/SRF CONTACT: kathleen.marchal@intec.ugent.be, siegfried.nijssen@cs.kuleuven.be SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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.

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.

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. PheNetic's 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/.

Dynamics of the Saccharomyces cerevisiae transcriptome during bread dough fermentation.

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 .

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.