A systems biology approach reveals the role of a novel methyltransferase in response to chemical stress and lipid homeostasis

Transcriptional response of Saccharomyces cerevisiae to lactic acid enantiomers

The model yeast, Saccharomyces cerevisiae, is a popular object for both fundamental and applied research, including the development of biosensors and industrial production of pharmaceutical compounds. However, despite multiple studies exploring S. cerevisiae transcriptional response to various substances, this response is unknown for some substances produced in yeast, such as D-lactic acid (DLA). Here, we explore the transcriptional response of the BY4742 strain to a wide range of DLA concentrations (from 0.05 to 45 mM), and compare it to the response to 45 mM L-lactic acid (LLA). We recorded a response to 5 and 45 mM DLA (125 and 113 differentially expressed genes (DEGs), respectively; > 50% shared) and a less pronounced response to 45 mM LLA (63 DEGs; > 30% shared with at least one DLA treatment). Our data did not reveal natural yeast promoters quantitatively sensing DLA but provide the first description of the transcriptome-wide response to DLA and enrich our understanding of the LLA response. Some DLA-activated genes were indeed related to lactate metabolism, as well as iron uptake and cell wall structure. Additional analyses showed that at least some of these genes were activated only by acidic form of DLA but not its salt, revealing the role of pH. The list of LLA-responsive genes was similar to those published previously and also included iron uptake and cell wall genes, as well as genes responding to other weak acids. These data might be instrumental for optimization of lactate production in yeast and yeast co-cultivation with lactic acid bacteria. KEY POINTS: • We present the first dataset on yeast transcriptional response to DLA. • Differential gene expression was correlated with yeast growth inhibition. • The transcriptome response to DLA was richer in comparison to LLA.

Profiling The Compendium Of Changes In Saccharomyces cerevisiae Due To Mutations That Alter Availability Of The Main Methyl Donor S-Adenosylmethionine

The SAM1 and SAM2 genes encode for S-AdenosylMethionine (AdoMet) synthetase enzymes, with AdoMet serving as the main methyl donor. We have previously shown that independent deletion of these genes alters chromosome stability and AdoMet concentrations in opposite ways in S. cerevisiae. To characterize other changes occurring in these mutants, we grew wildtype, sam1∆/sam1∆, and sam2∆/sam2∆ strains in 15 different Phenotypic Microarray plates with different components, equal to 1440 wells, and measured for growth variations. RNA-Sequencing was also carried out on these strains and differential gene expression determined for each mutant. In this study, we explore how the phenotypic growth differences are linked to the altered gene expression, and thereby predict the mechanisms by which loss of the SAM genes and subsequent AdoMet level changes, impact S. cerevisiae pathways and processes. We present six stories, discussing changes in sensitivity or resistance to azoles, cisplatin, oxidative stress, arginine biosynthesis perturbations, DNA synthesis inhibitors, and tamoxifen, to demonstrate the power of this novel methodology to broadly profile changes due to gene mutations. The large number of conditions that result in altered growth, as well as the large number of differentially expressed genes with wide-ranging functionality, speaks to the broad array of impacts that altering methyl donor abundance can impart, even when the conditions tested were not specifically selected as targeting known methyl involving pathways. Our findings demonstrate that some cellular changes are directly related to AdoMet-dependent methyltransferases and AdoMet availability, some are directly linked to the methyl cycle and its role is production of several important cellular components, and others reveal impacts of SAM gene mutations on previously unconnected pathways.

Cantharidin downregulates PSD1 expression and inhibits autophagic flux in yeast cells

Cantharidin is a terpenoid compound of insect origin, naturally produced by male blister beetles as an anti-predatory mechanism. Cantharidin has anticancer properties, which are attributed to its ability to induce cell cycle arrest, DNA damage, MAPK signalling pathway and apoptosis. Cantharidin has been reported to induce apoptosis in triple-negative breast cancer cells by suppressing autophagy via downregulation of Beclin 1 expression and autophagosome formation. However, it remains unclear which stage of the autophagic pathway is targeted by cantharidin. Herein, we report that yeast cells are sensitive to cantharidin, and external supplementation of ethanolamine (ETA) ameliorates the cytotoxicity. In addition, cantharidin downregulates phosphatidylserine decarboxylase1 (PSD1) expression. We also report that cantharidin inhibits autophagic flux, and external administration of ETA could rescue this inhibition. Additionally, co-treatment with chloroquine sensitized the autophagy inhibitory effects of cantharidin. We conclude that yeast cells are sensitive to cantharidin due to inhibition of autophagic flux.

Application of digital holographic tomography in antitumor effect of cantharides complex on 4T1 breast cancer cells

The study focuses on a methodology providing noninvasive monitoring and evaluation of the antitumor effect of traditional Chinese medicine, cantharides complex (canth), on 4T1 breast tumor cells. Digital holographic tomography (DHT) and developed data post-processing algorithms were used for quantitative estimation of changes in optical and morphological parameters of cells. We calculated and compared data on the refractive index, thickness, and projected area of 4T1 breast tumor cells in control untreated specimens and those treated with doxorubicin hydrochloride (DOX), canth, and their combinations. Post-treatment changes in cellular morphology recorded by DHT demonstrated that the two drugs led to noticeably different morphological changes in cells that can be presumably associated with different pathways of their death, apoptosis, or necrosis. The effect of combined treatment with these two drugs strongly depended on their relative concentrations and could lead to changes characteristic either for DOX or for canth; however, being more profound than those obtained when using each drug solely. The results obtained by DHT are in a good correspondence with commonly used cell viability analysis and immunofluorescent analysis of changes in cellular cytoskeleton.

ABC transporter Pdr5 is required for cantharidin resistance in Saccharomyces cerevisiae

Cantharidin is a potent anti-cancer drug and is known to exert its cytotoxic effects in several cancer cell lines. Although we have ample knowledge about its mode of action, we still know a little about cantharidin associated drug resistance mechanisms which dictates the efficacy and cytotoxic potential of this drug. In this direction, in the present study we employed Sacharomyces cerevisiae as a model organism and screened mutants of pleiotropic drug resistance network of genes for their susceptibility to cantharidin. We show that growth of pdr1Δ and pdr1Δpdr3Δ was severely reduced in presence of cantharidin whereas that of pdr3Δ remain unaffected when compared to wildtype. Loss of one of the PDR1 target genes PDR5, encoding an ABC membrane efflux pump, rendered the cells hypersensitive whereas overexpression of it conferred resistance. Additionally, cantharidin induced the upregulation of both PDR1 and PDR5 genes. Interestingly, pdr1Δpdr5Δ double deletion mutants were hypersensitive to cantharidin showing a synergistic effect in its cellular detoxification. Furthermore, transcriptional activation of PDR5 post cantharidin treatment was majorly dependent on the presence of Pdr1 and less significantly of Pdr3 transcription factors. Altogether our findings suggest that Pdr1 acts to increase cantharidin resistance by elevating the level of Pdr5 which serves as a major detoxification safeguard under CAN stress.

The Crg1 N-Terminus Is Essential for Methyltransferase Activity and Cantharidin Resistance in Saccharomyces cerevisiae

Crg1 is an S-adenosylmethionine (SAM)-dependent methyltransferase required for cantharidin resistance in yeast. Crg1 has a well-characterized methyltransferase domain that inactivates cantharidin by methylation. However, the remaining part of the Crg1 protein is yet to be functionally characterized. In this study, we identified an essential role of the N-terminus of Crg1 in methyltransferase activity and cantharidin resistance. Yeast cells lacking 41 residues of the N-terminus of Crg1 ( crg1ΔN) showed hypersensitivity to cantharidin as same as the null mutant, crg1. The mass spectrometry-based biochemical enzyme assay revealed a loss of methyltransferase activity in Crg1ΔN, which justifies the loss of cantharidin resistance, as well. The subcellular distribution of Crg1ΔN-daGFP showed cytoplasmic aggregates, whereas wild-type Crg1-daGFP was distributed normally in the cytoplasm. Interestingly, the Crg1-methyltransferase domain point mutants (D44A, D67A, and E105A/D108A) also showed the same cytoplasmic aggregates as Crg1ΔN-daGFP. In silico prediction of the tertiary structures of these mutants indicated an altered protein conformation. Altogether, these observations suggest that the N-terminal truncation, as well as the point mutations in the methyltransferase domain, alters the native folding of Crg1 methyltransferase, resulting in a loss of enzyme activity. Furthermore, the crg1ΔN mutant showed the same phenotypes as the crg1 null mutant in the presence of cantharidin, i.e., lethal cell growth, PE auxotrophy, temperature sensitivity, endoplasmic reticulum stress, GPI anchor missorting, and cell wall damage. Overall, this study identifies an essential role of the N-terminus of Crg1 in methyltransferase activity and cantharidin resistance.

The natural anticancer agent cantharidin alters GPI-anchored protein sorting by targeting Cdc1-mediated remodeling in endoplasmic reticulum

Cantharidin (CTD) is a potent anticancer small molecule produced by several species of blister beetle. It has been a traditional medicine for the management of warts and tumors for many decades. CTD suppresses tumor growth by inducing apoptosis, cell cycle arrest, and DNA damage and inhibits protein phosphatase 2 phosphatase activator (PP2A) and protein phosphatase 1 (PP1). CTD also alters lipid homeostasis, cell wall integrity, endocytosis, adhesion, and invasion in yeast cells. In this study, we identified additional molecular targets of CTD using a Saccharomyces cerevisiae strain that expresses a cantharidin resistance gene (CRG1), encoding a SAM-dependent methyltransferase that methylates and inactivates CTD. We found that CTD specifically affects phosphatidylethanolamine (PE)-associated functions that can be rescued by supplementing the growth media with ethanolamine (ETA). CTD also perturbed endoplasmic reticulum (ER) homeostasis and cell wall integrity by altering the sorting of glycosylphosphatidylinositol (GPI)-anchored proteins. A CTD-dependent genetic interaction profile of CRG1 revealed that the activity of the lipid phosphatase cell division control protein 1 (Cdc1) in GPI-anchor remodeling is the key target of CTD, independently of PP2A and PP1 activities. Moreover, experiments with human cells further suggested that CTD functions through a conserved mechanism in higher eukaryotes. Altogether, we conclude that CTD induces cytotoxicity by targeting Cdc1 activity in GPI-anchor remodeling in the ER.

Deciphering the Genic Basis of Yeast Fitness Variation by Simultaneous Forward and Reverse Genetics

The budding yeast Saccharomyces cerevisiae is the best studied eukaryote in molecular and cell biology, but its utility for understanding the genetic basis of phenotypic variation in natural populations is limited by inefficient association mapping due to strong and complex population structure. To overcome this challenge, we generated genome sequences for 85 strains and performed a comprehensive population genomic survey of a total of 190 diverse strains. We identified considerable variation in population structure among chromosomes and identified 181 genes that are absent from the reference genome. Many of these nonreference genes are expressed and we functionally confirmed that two of these genes confer increased resistance to antifungals. Next, we simultaneously measured the growth rates of over 4,500 laboratory strains, each of which lacks a nonessential gene, and 81 natural strains across multiple environments using unique DNA barcode present in each strain. By combining the genome-wide reverse genetic information gained from the gene deletion strains with a genome-wide association analysis from the natural strains, we identified genomic regions associated with fitness variation in natural populations. To experimentally validate a subset of these associations, we used reciprocal hemizygosity tests, finding that while the combined forward and reverse genetic approaches can identify a single causal gene, the phenotypic consequences of natural genetic variation often follow a complicated pattern. The resources and approach provided outline an efficient and reliable route to association mapping in yeast and significantly enhance its value as a model for understanding the genetic mechanisms underlying phenotypic variation and evolution in natural populations.

Transcriptional Responses of the Bacterium Burkholderia terrae BS001 to the Fungal Host Lyophyllum sp. Strain Karsten under Soil-Mimicking Conditions

In this study, the mycosphere isolate Burkholderia terrae BS001 was confronted with the soil fungus Lyophyllum sp. strain Karsten on soil extract agar plates in order to examine its transcriptional responses over time. At the initial stages of the experiment (T1-day 3; T2-day 5), contact between both partner organisms was absent, whereas in the final stage (T3-day 8), the two populations made intimate physical contact. Overall, a strong modulation of the strain BS001 gene expression patterns was found. First, the stationary-phase sigma factor RpoS, and numerous genes under its control, were strongly expressed as a response to the soil extract agar, and this extended over the whole temporal regime. In the system, B. terrae BS001 apparently perceived the presence of the fungal hyphae already at the early experimental stages (T1, T2), by strongly upregulating a suite of chemotaxis and flagellar motility genes. With respect to specific metabolism and energy generation, a picture of differential involvement in different metabolic routes was obtained. Initial (T1, T2) up- or downregulation of ethanolamine and mandelate uptake and utilization pathways was substituted by a strong investment, in the presence of the fungus, in the expression of putative metabolic gene clusters (T3). Specifically at T3, five clustered genes that are potentially involved in energy generation coupled to an oxidative stress response, and two genes encoding short-chain dehydrogenases/oxidoreductases (SDR), were highly upregulated. In contrast, the dnaE2 gene (related to general stress response; encoding error-prone DNA polymerase) was transcriptionally downregulated at this stage. This study revealed that B. terrae BS001, from a stress-induced state, resulting from the soil extract agar milieu, responds positively to fungal hyphae that encroach upon it, in a temporally dynamic manner. The response is characterized by phases in which the modulation of (1) chemotaxis, (2) metabolic activity, and (3) oxidative stress responses are key mechanisms.

Recent advances in developing insect natural products as potential modern day medicines

Except for honey as food, and silk for clothing and pollination of plants, people give little thought to the benefits of insects in their lives. This overview briefly describes significant recent advances in developing insect natural products as potential new medicinal drugs. This is an exciting and rapidly expanding new field since insects are hugely variable and have utilised an enormous range of natural products to survive environmental perturbations for 100s of millions of years. There is thus a treasure chest of untapped resources waiting to be discovered. Insects products, such as silk and honey, have already been utilised for thousands of years, and extracts of insects have been produced for use in Folk Medicine around the world, but only with the development of modern molecular and biochemical techniques has it become feasible to manipulate and bioengineer insect natural products into modern medicines. Utilising knowledge gleaned from Insect Folk Medicines, this review describes modern research into bioengineering honey and venom from bees, silk, cantharidin, antimicrobial peptides, and maggot secretions and anticoagulants from blood-sucking insects into medicines. Problems and solutions encountered in these endeavours are described and indicate that the future is bright for new insect derived pharmaceuticals treatments and medicines.

Probabilistic approach to predicting substrate specificity of methyltransferases

We present a general probabilistic framework for predicting the substrate specificity of enzymes. We designed this approach to be easily applicable to different organisms and enzymes. Therefore, our predictive models do not rely on species-specific properties and use mostly sequence-derived data. Maximum Likelihood optimization is used to fine-tune model parameters and the Akaike Information Criterion is employed to overcome the issue of correlated variables. As a proof-of-principle, we apply our approach to predicting general substrate specificity of yeast methyltransferases (MTases). As input, we use several physico-chemical and biological properties of MTases: structural fold, isoelectric point, expression pattern and cellular localization. Our method accurately predicts whether a yeast MTase methylates a protein, RNA or another molecule. Among our experimentally tested predictions, 89% were confirmed, including the surprising prediction that YOR021C is the first known MTase with a SPOUT fold that methylates a substrate other than RNA (protein). Our approach not only allows for highly accurate prediction of functional specificity of MTases, but also provides insight into general rules governing MTase substrate specificity.

Our approach is easily applicable to different organisms, because it does not rely on species-specific properties and uses mostly sequence-derived and other readily available data (e.g. isoelectric point or predicted structural fold). Tests on yeast MTases indicate that the accuracy of our predictions is ∼90%. We show that knowledge of substrate binding sites or corresponding motifs is not crucial for highly accurate general substrate specificity predictions of enzymes, and provide new insights into how such specificities are achieved at the molecular level. We predict substrate specificities not yet observed for a given class of enzymes, and experimentally verify our predictions.

Lifespan extension conferred by endoplasmic reticulum secretory pathway deficiency requires induction of the unfolded protein response

Cells respond to accumulation of misfolded proteins in the endoplasmic reticulum (ER) by activating the unfolded protein response (UPR) signaling pathway. The UPR restores ER homeostasis by degrading misfolded proteins, inhibiting translation, and increasing expression of chaperones that enhance ER protein folding capacity. Although ER stress and protein aggregation have been implicated in aging, the role of UPR signaling in regulating lifespan remains unknown. Here we show that deletion of several UPR target genes significantly increases replicative lifespan in yeast. This extended lifespan depends on a functional ER stress sensor protein, Ire1p, and is associated with constitutive activation of upstream UPR signaling. We applied ribosome profiling coupled with next generation sequencing to quantitatively examine translational changes associated with increased UPR activity and identified a set of stress response factors up-regulated in the long-lived mutants. Besides known UPR targets, we uncovered up-regulation of components of the cell wall and genes involved in cell wall biogenesis that confer resistance to multiple stresses. These findings demonstrate that the UPR is an important determinant of lifespan that governs ER stress and identify a signaling network that couples stress resistance to longevity.

An evaluation of high-throughput approaches to QTL mapping in Saccharomyces cerevisiae

Dissecting the molecular basis of quantitative traits is a significant challenge and is essential for understanding complex diseases. Even in model organisms, precisely determining causative genes and their interactions has remained elusive, due in part to difficulty in narrowing intervals to single genes and in detecting epistasis or linked quantitative trait loci. These difficulties are exacerbated by limitations in experimental design, such as low numbers of analyzed individuals or of polymorphisms between parental genomes. We address these challenges by applying three independent high-throughput approaches for QTL mapping to map the genetic variants underlying 11 phenotypes in two genetically distant Saccharomyces cerevisiae strains, namely (1) individual analysis of >700 meiotic segregants, (2) bulk segregant analysis, and (3) reciprocal hemizygosity scanning, a new genome-wide method that we developed. We reveal differences in the performance of each approach and, by combining them, identify eight polymorphic genes that affect eight different phenotypes: colony shape, flocculation, growth on two nonfermentable carbon sources, and resistance to two drugs, salt, and high temperature. Our results demonstrate the power of individual segregant analysis to dissect QTL and address the underestimated contribution of interactions between variants. We also reveal confounding factors like mutations and aneuploidy in pooled approaches, providing valuable lessons for future designs of complex trait mapping studies.

A novel small molecule methyltransferase is important for virulence in Candida albicans

Candida albicans is an opportunistic pathogen capable of causing life-threatening infections in immunocompromised individuals. Despite its significant health impact, our understanding of C. albicans pathogenicity is limited, particularly at the molecular level. One of the largely understudied enzyme families in C. albicans are small molecule AdoMet-dependent methyltransferases (smMTases), which are important for maintenance of cellular homeostasis by clearing toxic chemicals, generating novel cellular intermediates, and regulating intra- and interspecies interactions. In this study, we demonstrated that C. albicans Crg1 (CaCrg1) is a bona fide smMTase that interacts with the toxin in vitro and in vivo. We report that CaCrg1 is important for virulence-related processes such as adhesion, hyphal elongation, and membrane trafficking. Biochemical and genetic analyses showed that CaCrg1 plays a role in the complex sphingolipid pathway: it binds to exogenous short-chain ceramides in vitro and interacts genetically with genes of glucosylceramide pathway, and the deletion of CaCRG1 leads to significant changes in the abundance of phytoceramides. Finally we found that this novel lipid-related smMTase is required for virulence in the waxmoth Galleria mellonella, a model of infection.

Architectures, mechanisms and molecular evolution of natural product methyltransferases

The addition of a methyl moiety to a small chemical is a common transformation in the biosynthesis of natural products across all three domains of life. These methylation reactions are most often catalysed by S-adenosyl-L-methionine (SAM)-dependent methyltransferases (MTs). MTs are categorized based on the electron-rich, methyl accepting atom, usually O, N, C, or S. SAM-dependent natural product MTs (NPMTs) are responsible for the modification of a wide array of structurally distinct substrates, including signalling and host defense compounds, pigments, prosthetic groups, cofactors, cell membrane and cell wall components, and xenobiotics. Most notably, methylation modulates the bioavailability, bioactivity, and reactivity of acceptor molecules, and thus exerts a central role on the functional output of many metabolic pathways. Our current understanding of the structural enzymology of NPMTs groups these phylogenetically diverse enzymes into two MT-superfamily fold classes (class I and class III). Structural biology has also shed light on the catalytic mechanisms and molecular bases for substrate specificity for over fifty NPMTs. These biophysical-based approaches have contributed to our understanding of NPMT evolution, demonstrating how a widespread protein fold evolved to accommodate chemically diverse methyl acceptors and to catalyse disparate mechanisms suited to the physiochemical properties of the target substrates. This evolutionary diversity suggests that NPMTs may serve as starting points for generating new biocatalysts.

Evolutionary and Biochemical Aspects of Chemical Stress Resistance in Saccharomyces cerevisiae

Large-scale chemical genetics screens (chemogenomics) in yeast have been widely used to find drug targets, understand the mechanism-of-action of compounds, and unravel the biochemistry of drug resistance. Chemogenomics is based on the comparison of growth of gene deletants in the presence and absence of a chemical substance. Such studies showed that more than 90% of the yeast genes are required for growth in the presence of at least one chemical. Analysis of these data, using computational approaches, has revealed non-trivial features of the natural chemical tolerance systems. As a result two non-overlapping sets of genes are seen to respectively impart robustness and evolvability in the context of natural chemical resistance. The former is composed of multidrug-resistance genes, whereas the latter comprises genes sharing chemical genetic profiles with many others. Recent publications showing the potential applications chemogenomics in studying the pharmacological basis of various drugs are discussed, as well as the expansion of chemogenomics to other organisms. Finally, integration of chemogenomics with sensitive sequence analysis and ubiquitination/phosphorylation data led to the discovery of a new conserved domain and important post-translational modification pathways involved in stress resistance.