BridGE: a pathway-based analysis tool for detecting genetic interactions from GWAS

Genetic interactions have the potential to modulate phenotypes, including human disease. In principle, genome-wide association studies (GWAS) provide a platform for detecting genetic interactions; however, traditional methods for identifying them, which tend to focus on testing individual variant pairs, lack statistical power. In this protocol, we describe a novel computational approach, called Bridging Gene sets with Epistasis (BridGE), for discovering genetic interactions between biological pathways from GWAS data. We present a Python-based implementation of BridGE along with instructions for its application to a typical human GWAS cohort. The major stages include initial data processing and quality control, construction of a variant-level genetic interaction network, measurement of pathway-level genetic interactions, evaluation of statistical significance using sample permutations and generation of results in a standardized output format. The BridGE software pipeline includes options for running the analysis on multiple cores and multiple nodes for users who have access to computing clusters or a cloud computing environment. In a cluster computing environment with 10 nodes and 100 GB of memory per node, the method can be run in less than 24 h for typical human GWAS cohorts. Using BridGE requires knowledge of running Python programs and basic shell script programming experience.

The PRC2.1 Subcomplex Opposes G1 Progression through Regulation of CCND1 and CCND2

Progression through the G1 phase of the cell cycle is the most highly regulated step in cellular division. We employed a chemogenomics approach to discover novel cellular networks that regulate cell cycle progression. This approach uncovered functional clusters of genes that altered sensitivity of cells to inhibitors of the G1/S transition. Mutation of components of the Polycomb Repressor Complex 2 rescued growth inhibition caused by the CDK4/6 inhibitor palbociclib, but not to inhibitors of S phase or mitosis. In addition to its core catalytic subunits, mutation of the PRC2.1 accessory protein MTF2, but not the PRC2.2 protein JARID2, rendered cells resistant to palbociclib treatment. We found that PRC2.1 (MTF2), but not PRC2.2 (JARID2), was critical for promoting H3K27me3 deposition at CpG islands genome-wide and in promoters. This included the CpG islands in the promoter of the CDK4/6 cyclins CCND1 and CCND2, and loss of MTF2 lead to upregulation of both CCND1 and CCND2. Our results demonstrate a role for PRC2.1, but not PRC2.2, in promoting G1 progression.

Decomposition of dynamic transcriptomic responses during effector-triggered immunity reveals conserved responses in two distinct plant cell populations

Rapid plant immune responses in the appropriate cells are needed for effective defense against pathogens. Although transcriptome analysis is often used to describe overall immune responses, collection of transcriptome data with sufficient resolution in both space and time is challenging. We reanalyzed public Arabidopsis time-course transcriptome data obtained after low-dose inoculation with a Pseudomonas syringae strain expressing the effector AvrRpt2, which induces effector-triggered immunity in Arabidopsis. Double-peak time-course patterns are prevalent among thousands of upregulated genes. We implemented a multi-compartment modeling approach to decompose the double-peak pattern into two single-peak patterns for each gene. The decomposed peaks reveal an "echoing" pattern: the peak times of the first and second peaks correlate well across most upregulated genes. We demonstrated that the two peaks likely represent responses of two distinct cell populations that respond either cell autonomously or indirectly to AvrRpt2. Thus, the peak decomposition has extracted spatial information from the time-course data. The echoing pattern also indicates a conserved transcriptome response with different initiation times between the two cell populations despite different elicitor types. A gene set highly overlapping with the conserved gene set is also upregulated with similar kinetics during pattern-triggered immunity. Activation of a WRKY network via different entry-point WRKYs can explain the similar but not identical transcriptome responses elicited by different elicitor types. We discuss potential benefits of the properties of the WRKY activation network as an immune signaling network in light of pressure from rapidly evolving pathogens.

RNF4 prevents genomic instability caused by chronic DNA under-replication

Eukaryotic genome stability is maintained by a complex and diverse set of molecular processes. One class of enzymes that promotes proper DNA repair, replication and cell cycle progression comprises small ubiquitin-like modifier (SUMO)-targeted E3 ligases, or STUbLs. Previously, we reported a role for the budding yeast STUbL synthetically lethal with sgs1 (Slx) 5/8 in preventing G2/M-phase arrest in a minichromosome maintenance protein 10 (Mcm10)-deficient model of replication stress. Here, we extend these studies to human cells, examining the requirement for the human STUbL RING finger protein 4 (RNF4) in MCM10 mutant cancer cells. We find that MCM10 and RNF4 independently promote origin firing but regulate DNA synthesis epistatically and, unlike in yeast, the negative genetic interaction between RNF4 and MCM10 causes cells to accumulate in G1-phase. When MCM10 is deficient, RNF4 prevents excessive DNA under-replication at hard-to-replicate regions that results in large DNA copy number alterations and severely reduced viability. Overall, our findings highlight that STUbLs participate in species-specific mechanisms to maintain genome stability, and that human RNF4 is required for origin activation in the presence of chronic replication stress.

A scalable platform for efficient CRISPR-Cas9 chemical-genetic screens of DNA damage-inducing compounds

Current approaches to define chemical-genetic interactions (CGIs) in human cell lines are resource-intensive. We designed a scalable chemical-genetic screening platform by generating a DNA damage response (DDR)-focused custom sgRNA library targeting 1011 genes with 3033 sgRNAs. We performed five proof-of-principle compound screens and found that the compounds' known modes-of-action (MoA) were enriched among the compounds' CGIs. These scalable screens recapitulated expected CGIs at a comparable signal-to-noise ratio (SNR) relative to genome-wide screens. Furthermore, time-resolved CGIs, captured by sequencing screens at various time points, suggested an unexpected, late interstrand-crosslinking (ICL) repair pathway response to camptothecin-induced DNA damage. Our approach can facilitate screening compounds at scale with 20-fold fewer resources than commonly used genome-wide libraries and produce biologically informative CGI profiles.

Dimensionality reduction methods for extracting functional networks from large-scale CRISPR screens

CRISPR-Cas9 screens facilitate the discovery of gene functional relationships and phenotype-specific dependencies. The Cancer Dependency Map (DepMap) is the largest compendium of whole-genome CRISPR screens aimed at identifying cancer-specific genetic dependencies across human cell lines. A mitochondria-associated bias has been previously reported to mask signals for genes involved in other functions, and thus, methods for normalizing this dominant signal to improve co-essentiality networks are of interest. In this study, we explore three unsupervised dimensionality reduction methods-autoencoders, robust, and classical principal component analyses (PCA)-for normalizing the DepMap to improve functional networks extracted from these data. We propose a novel "onion" normalization technique to combine several normalized data layers into a single network. Benchmarking analyses reveal that robust PCA combined with onion normalization outperforms existing methods for normalizing the DepMap. Our work demonstrates the value of removing low-dimensional signals from the DepMap before constructing functional gene networks and provides generalizable dimensionality reduction-based normalization tools.

N-terminal acetylation shields proteins from degradation and promotes age-dependent motility and longevity

Most eukaryotic proteins are N-terminally acetylated, but the functional impact on a global scale has remained obscure. Using genome-wide CRISPR knockout screens in human cells, we reveal a strong genetic dependency between a major N-terminal acetyltransferase and specific ubiquitin ligases. Biochemical analyses uncover that both the ubiquitin ligase complex UBR4-KCMF1 and the acetyltransferase NatC recognize proteins bearing an unacetylated N-terminal methionine followed by a hydrophobic residue. NatC KO-induced protein degradation and phenotypes are reversed by UBR knockdown, demonstrating the central cellular role of this interplay. We reveal that loss of Drosophila NatC is associated with male sterility, reduced longevity, and age-dependent loss of motility due to developmental muscle defects. Remarkably, muscle-specific overexpression of UbcE2M, one of the proteins targeted for NatC KO-mediated degradation, suppresses defects of NatC deletion. In conclusion, NatC-mediated N-terminal acetylation acts as a protective mechanism against protein degradation, which is relevant for increased longevity and motility.

RNF4 and USP7 cooperate in ubiquitin-regulated steps of DNA replication

DNA replication requires precise regulation achieved through post-translational modifications, including ubiquitination and SUMOylation. These modifications are linked by the SUMO-targeted E3 ubiquitin ligases (STUbLs). Ring finger protein 4 (RNF4), one of only two mammalian STUbLs, participates in double-strand break repair and resolving DNA-protein cross-links. However, its role in DNA replication has been poorly understood. Using CRISPR/Cas9 genetic screens, we discovered an unexpected dependency of RNF4 mutants on ubiquitin specific peptidase 7 (USP7) for survival in TP53-null retinal pigment epithelial cells. TP53-/-/RNF4-/-/USP7-/- triple knockout (TKO) cells displayed defects in DNA replication that cause genomic instability. These defects were exacerbated by the proteasome inhibitor bortezomib, which limited the nuclear ubiquitin pool. A shortage of free ubiquitin suppressed the ataxia telangiectasia and Rad3-related (ATR)-mediated checkpoint response, leading to increased cell death. In conclusion, RNF4 and USP7 work cooperatively to sustain a functional level of nuclear ubiquitin to maintain the integrity of the genome.

CD69 marks a subpopulation of acute myeloid leukemia with enhanced colony forming capacity and a unique signaling activation state

In acute myeloid leukemia (AML), leukemia stem cells (LSCs) have self-renewal potential and are responsible for relapse. We previously showed that, in Mll-AF9/NRASG12V murine AML, CD69 expression marks an LSC-enriched subpopulation with enhanced in vivo self-renewal capacity. Here, we used CyTOF to define activated signaling pathways in LSC subpopulations in Mll-AF9/NRASG12V AML. Furthermore, we compared the signaling activation states of CD69High and CD36High subsets of primary human AML. The human CD69High subset expresses low levels of Ki67 and high levels of NFκB and pMAPKAPKII. Additionally, the human CD69High AML subset also has enhanced colony-forming capacity. We applied Bayesian network modeling to compare the global signaling network within the human AML subsets. We find that distinct signaling states, distinguished by NFκB and pMAPKAPKII levels, correlate with divergent functional subsets, defined by CD69 and CD36 expression, in human AML. Targeting NFκB with proteasome inhibition diminished colony formation.

Reproducibility metrics for context-specific CRISPR screens

CRISPR screens are used extensively to systematically interrogate the phenotype-to-genotype problem. In contrast to early CRISPR screens, which defined core cell fitness genes, most current efforts now aim to identify context-specific phenotypes that differentiate a cell line, genetic background, or condition of interest, such as a drug treatment. While CRISPR-related technologies have shown great promise and a fast pace of innovation, a better understanding of standards and methods for quality assessment of CRISPR screen results is crucial to guide technology development and application. Specifically, many commonly used metrics for quantifying screen quality do not accurately measure the reproducibility of context-specific hits. We highlight the importance of reporting reproducibility statistics that directly relate to the purpose of the screen and suggest the use of metrics that are sensitive to context-specific signal. A record of this paper's transparent peer review process is included in the supplemental information.

Fanconi anemia-associated chromosomal radial formation is dependent on POLθ-mediated alternative end joining

Activation of the Fanconi anemia (FA) pathway after treatment with mitomycin C (MMC) is essential for preventing chromosome translocations termed "radials." When replication forks stall at MMC-induced interstrand crosslinks (ICLs), the FA pathway is activated to orchestrate ICL unhooking and repair of the DNA break intermediates. However, in FA-deficient cells, how ICL-associated breaks are resolved in a manner that leads to radials is unclear. Here, we demonstrate that MMC-induced radials are dependent on DNA polymerase theta (POLθ)-mediated alternative end joining (A-EJ). Specifically, we show that radials observed in FANCD2^-/- cells are dependent on POLθ and DNA ligase III and occur independently of classical non-homologous end joining. Furthermore, treatment of FANCD2^-/- cells with POLθ inhibitors abolishes radials and leads to the accumulation of breaks co-localizing with common fragile sites. Uniformly, these observations implicate A-EJ in radial formation and provide mechanistic insights into the treatment of FA pathway-deficient cancers with POLθ inhibitors.

Systems-level transcriptional regulation of Caenorhabditis elegans metabolism

Metabolism is controlled to ensure organismal development and homeostasis. Several mechanisms regulate metabolism, including allosteric control and transcriptional regulation of metabolic enzymes and transporters. So far, metabolism regulation has mostly been described for individual genes and pathways, and the extent of transcriptional regulation of the entire metabolic network remains largely unknown. Here, we find that three-quarters of all metabolic genes are transcriptionally regulated in the nematode Caenorhabditis elegans. We find that many annotated metabolic pathways are coexpressed, and we use gene expression data and the iCEL1314 metabolic network model to define coregulated subpathways in an unbiased manner. Using a large gene expression compendium, we determine the conditions where subpathways exhibit strong coexpression. Finally, we develop "WormClust," a web application that enables a gene-by-gene query of genes to view their association with metabolic (sub)-pathways. Overall, this study sheds light on the ubiquity of transcriptional regulation of metabolism and provides a blueprint for similar studies in other organisms, including humans.

Dimensionality reduction methods for extracting functional networks from large-scale CRISPR screens

CRISPR-Cas9 screens facilitate the discovery of gene functional relationships and phenotype-specific dependencies. The Cancer Dependency Map (DepMap) is the largest compendium of whole-genome CRISPR screens aimed at identifying cancer-specific genetic dependencies across human cell lines. A mitochondria-associated bias has been previously reported to mask signals for genes involved in other functions, and thus, methods for normalizing this dominant signal to improve co-essentiality networks are of interest. In this study, we explore three unsupervised dimensionality reduction methods - autoencoders, robust, and classical principal component analyses (PCA) - for normalizing the DepMap to improve functional networks extracted from these data. We propose a novel "onion" normalization technique to combine several normalized data layers into a single network. Benchmarking analyses reveal that robust PCA combined with onion normalization outperforms existing methods for normalizing the DepMap. Our work demonstrates the value of removing low-dimensional signals from the DepMap before constructing functional gene networks and provides generalizable dimensionality reduction-based normalization tools.

Genetic architecture of protein and oil content in soybean seed and meal

Soybean is grown primarily for the protein and oil extracted from its seed and its value is influenced by these components. The objective of this study was to map marker-trait associations (MTAs) for the concentration of seed protein, oil, and meal protein using the soybean nested association mapping (SoyNAM) population. The composition traits were evaluated on seed harvested from over 5000 inbred lines of the SoyNAM population grown in 10 field locations across 3 years. Estimated heritabilities were at least 0.85 for all three traits. The genotyping of lines with single nucleotide polymorphism markers resulted in the identification of 107 MTAs for the three traits. When MTAs for the three traits that mapped within 5 cM intervals were binned together, the MTAs were mapped to 64 intervals on 19 of the 20 soybean chromosomes. The majority of the MTA effects were small and of the 107 MTAs, 37 were for protein content, 39 for meal protein, and 31 for oil content. For cases where a protein and oil MTAs mapped to the same interval, most (94%) significant effects were opposite for the two traits, consistent with the negative correlation between these traits. A coexpression analysis identified candidate genes linked to MTAs and 18 candidate genes were identified. The large number of small effect MTAs for the composition traits suggest that genomic prediction would be more effective in improving these traits than marker-assisted selection.

Context-dependent regulation of ferroptosis sensitivity

Ferroptosis is an important mediator of pathophysiological cell death and an emerging target for cancer therapy. Whether ferroptosis sensitivity is governed by a single regulatory mechanism is unclear. Here, based on the integration of 24 published chemical genetic screens combined with targeted follow-up experimentation, we find that the genetic regulation of ferroptosis sensitivity is highly variable and context-dependent. For example, the lipid metabolic gene acyl-coenzyme A (CoA) synthetase long chain family member 4 (ACSL4) appears far more essential for ferroptosis triggered by direct inhibition of the lipid hydroperoxidase glutathione peroxidase 4 (GPX4) than by cystine deprivation. Despite this, distinct pro-ferroptotic stimuli converge upon a common lethal effector mechanism: accumulation of lipid peroxides at the plasma membrane. These results indicate that distinct genetic mechanisms regulate ferroptosis sensitivity, with implications for the initiation and analysis of this process in vivo.

The Parkinson's disease protein alpha-synuclein is a modulator of processing bodies and mRNA stability

Alpha-synuclein (αS) is a conformationally plastic protein that reversibly binds to cellular membranes. It aggregates and is genetically linked to Parkinson's disease (PD). Here, we show that αS directly modulates processing bodies (P-bodies), membraneless organelles that function in mRNA turnover and storage. The N terminus of αS, but not other synucleins, dictates mutually exclusive binding either to cellular membranes or to P-bodies in the cytosol. αS associates with multiple decapping proteins in close proximity on the Edc4 scaffold. As αS pathologically accumulates, aberrant interaction with Edc4 occurs at the expense of physiologic decapping-module interactions. mRNA decay kinetics within PD-relevant pathways are correspondingly disrupted in PD patient neurons and brain. Genetic modulation of P-body components alters αS toxicity, and human genetic analysis lends support to the disease-relevance of these interactions. Beyond revealing an unexpected aspect of αS function and pathology, our data highlight the versatility of conformationally plastic proteins with high intrinsic disorder.

Acute aerobic exercise reveals FAHFAs distinguish the metabolomes of overweight and normal weight runners

Background

Responses of the metabolome to acute aerobic exercise may predict maximum oxygen consumption (VO₂max) and longer-term outcomes, including the development of diabetes and its complications.

Methods

Serum samples were collected from overweight/obese trained (OWT) and normal-weight trained (NWT) runners prior to and immediately after a supervised 90-minute treadmill run at 60% VO₂max (NWT = 14, OWT = 11) in a cross-sectional study. We applied a liquid chromatography high-resolution–mass spectrometry–based untargeted metabolomics platform to evaluate the effect of acute aerobic exercise on the serum metabolome.

Results

NWT and OWT metabolic profiles shared increased circulating acylcarnitines and free fatty acids (FFAs) with exercise, while intermediates of adenine metabolism, inosine, and hypoxanthine were strongly correlated with body fat percentage and VO₂max. Untargeted metabolomics-guided follow-up quantitative lipidomic analysis revealed that baseline levels of fatty acid esters of hydroxy fatty acids (FAHFAs) were generally diminished in the OWT group. FAHFAs negatively correlated with visceral fat mass and HOMA-IR. Strikingly, a 4-fold decrease in FAHFAs was provoked by acute aerobic running in NWT participants, an effect that negatively correlated with circulating IL-6; these effects were not observed in the OWT group. Machine learning models based on a preexercise metabolite profile that included FAHFAs, FFAs, and adenine intermediates predicted VO₂max.

Conclusion

These findings in overweight human participants and healthy controls indicate that exercise-provoked changes in FAHFAs distinguish normal-weight from overweight participants and could predict VO₂max. These results support the notion that FAHFAs could modulate the inflammatory response, fuel utilization, and insulin resistance.

Trial registration

ClinicalTrials.gov, NCT02150889.

Funding

NIH DK091538, AG069781, DK098203, TR000114, UL1TR002494.

Probing ion channel functional architecture and domain recombination compatibility by massively parallel domain insertion profiling

Protein domains are the basic units of protein structure and function. Comparative analysis of genomes and proteomes showed that domain recombination is a main driver of multidomain protein functional diversification and some of the constraining genomic mechanisms are known. Much less is known about biophysical mechanisms that determine whether protein domains can be combined into viable protein folds. Here, we use massively parallel insertional mutagenesis to determine compatibility of over 300,000 domain recombination variants of the Inward Rectifier K⁺ channel Kir2.1 with channel surface expression. Our data suggest that genomic and biophysical mechanisms acted in concert to favor gain of large, structured domain at protein termini during ion channel evolution. We use machine learning to build a quantitative biophysical model of domain compatibility in Kir2.1 that allows us to derive rudimentary rules for designing domain insertion variants that fold and traffic to the cell surface. Positional Kir2.1 responses to motif insertion clusters into distinct groups that correspond to contiguous structural regions of the channel with distinct biophysical properties tuned towards providing either folding stability or gating transitions. This suggests that insertional profiling is a high-throughput method to annotate function of ion channel structural regions.

Leveraging machine learning essentiality predictions and chemogenomic interactions to identify antifungal targets

Fungal pathogens pose a global threat to human health, with Candida albicans among the leading killers. Systematic analysis of essential genes provides a powerful strategy to discover potential antifungal targets. Here, we build a machine learning model to generate genome-wide gene essentiality predictions for C. albicans and expand the largest functional genomics resource in this pathogen (the GRACE collection) by 866 genes. Using this model and chemogenomic analyses, we define the function of three uncharacterized essential genes with roles in kinetochore function, mitochondrial integrity, and translation, and identify the glutaminyl-tRNA synthetase Gln4 as the target of N-pyrimidinyl-β-thiophenylacrylamide (NP-BTA), an antifungal compound.

Analysis of combinatorial CRISPR screens with the Orthrus scoring pipeline

The continued improvement of combinatorial CRISPR screening platforms necessitates the development of new computational pipelines for scoring combinatorial screening data. Unlike for single-guide RNA (sgRNA) pooled screening platforms, combinatorial scoring for multiplexed systems is confounded by guide design parameters such as the number of gRNAs per construct, the position of gRNAs along constructs, and additional features that may impact gRNA expression, processing or capture. In this protocol we describe Orthrus, an R package for processing, scoring and analyzing combinatorial CRISPR screening data that addresses these challenges. This protocol walks through the application of Orthrus to previously published combinatorial screening data from the CHyMErA experimental system, a platform we recently developed that pairs Cas9 with Cas12a gRNAs and enables programmed targeting of multiple genomic sites. We demonstrate Orthrus' features for screen quality assessment and two distinct scoring modes for dual guide RNAs (dgRNAs) that target the same gene twice or dgRNAs that target two different genes. Running Orthrus requires basic R programming experience, ~5-10 min of computational time and 15-60 min total.

Trans-acting genetic variation affects the expression of adjacent genes

Gene expression differences among individuals are shaped by trans-acting expression quantitative trait loci (eQTLs). Most trans-eQTLs map to hotspot locations that influence many genes. The molecular mechanisms perturbed by hotspots are often assumed to involve "vertical" cascades of effects in pathways that can ultimately affect the expression of thousands of genes. Here, we report that trans-eQTLs can affect the expression of adjacent genes via "horizontal" mechanisms that extend along a chromosome. Genes affected by trans-eQTL hotspots in the yeast Saccharomyces cerevisiae were more likely to be located next to each other than expected by chance. These paired hotspot effects tended to occur at adjacent genes that also show coexpression in response to genetic and environmental perturbations, suggesting shared mechanisms. Physical proximity and shared chromatin state, in addition to regulation of adjacent genes by similar transcription factors, were independently associated with paired hotspot effects among adjacent genes. Paired effects of trans-eQTLs can occur at neighboring genes even when these genes do not share a common function. This phenomenon could result in unexpected connections between regulatory genetic variation and phenotypes.

Systematic analysis of bypass suppression of essential genes

Essential genes tend to be highly conserved across eukaryotes, but, in some cases, their critical roles can be bypassed through genetic rewiring. From a systematic analysis of 728 different essential yeast genes, we discovered that 124 (17%) were dispensable essential genes. Through whole-genome sequencing and detailed genetic analysis, we investigated the genetic interactions and genome alterations underlying bypass suppression. Dispensable essential genes often had paralogs, were enriched for genes encoding membrane-associated proteins, and were depleted for members of protein complexes. Functionally related genes frequently drove the bypass suppression interactions. These gene properties were predictive of essential gene dispensability and of specific suppressors among hundreds of genes on aneuploid chromosomes. Our findings identify yeast's core essential gene set and reveal that the properties of dispensable essential genes are conserved from yeast to human cells, correlating with human genes that display cell line-specific essentiality in the Cancer Dependency Map (DepMap) project.

Improving Measures of Chemical Structural Similarity Using Machine Learning on Chemical-Genetic Interactions

A common strategy for identifying molecules likely to possess a desired biological activity is to search large databases of compounds for high structural similarity to a query molecule that demonstrates this activity, under the assumption that structural similarity is predictive of similar biological activity. However, efforts to systematically benchmark the diverse array of available molecular fingerprints and similarity coefficients have been limited by a lack of large-scale datasets that reflect biological similarities of compounds. To elucidate the relative performance of these alternatives, we systematically benchmarked 11 different molecular fingerprint encodings, each combined with 13 different similarity coefficients, using a large set of chemical–genetic interaction data from the yeast Saccharomyces cerevisiae as a systematic proxy for biological activity. We found that the performance of different molecular fingerprints and similarity coefficients varied substantially and that the all-shortest path fingerprints paired with the Braun-Blanquet similarity coefficient provided superior performance that was robust across several compound collections. We further proposed a machine learning pipeline based on support vector machines that offered a fivefold improvement relative to the best unsupervised approach. Our results generally suggest that using high-dimensional chemical–genetic data as a basis for refining molecular fingerprints can be a powerful approach for improving prediction of biological functions from chemical structures.

Comparing Host Module Activation Patterns and Temporal Dynamics in Infection by Influenza H1N1 Viruses

Influenza is a serious global health threat that shows varying pathogenicity among different virus strains. Understanding similarities and differences among activated functional pathways in the host responses can help elucidate therapeutic targets responsible for pathogenesis. To compare the types and timing of functional modules activated in host cells by four influenza viruses of varying pathogenicity, we developed a new DYNAmic MOdule (DYNAMO) method that addresses the need to compare functional module utilization over time. This integrative approach overlays whole genome time series expression data onto an immune-specific functional network, and extracts conserved modules exhibiting either different temporal patterns or overall transcriptional activity. We identified a common core response to influenza virus infection that is temporally shifted for different viruses. We also identified differentially regulated functional modules that reveal unique elements of responses to different virus strains. Our work highlights the usefulness of combining time series gene expression data with a functional interaction map to capture temporal dynamics of the same cellular pathways under different conditions. Our results help elucidate conservation of the immune response both globally and at a granular level, and provide mechanistic insight into the differences in the host response to infection by influenza strains of varying pathogenicity.

Environmental robustness of the global yeast genetic interaction network

Phenotypes associated with genetic variants can be altered by interactions with other genetic variants (GxG), with the environment (GxE), or both (GxGxE). Yeast genetic interactions have been mapped on a global scale, but the environmental influence on the plasticity of genetic networks has not been examined systematically. To assess environmental rewiring of genetic networks, we examined 14 diverse conditions and scored 30,000 functionally representative yeast gene pairs for dynamic, differential interactions. Different conditions revealed novel differential interactions, which often uncovered functional connections between distantly related gene pairs. However, the majority of observed genetic interactions remained unchanged in different conditions, suggesting that the global yeast genetic interaction network is robust to environmental perturbation and captures the fundamental functional architecture of a eukaryotic cell.

A method for benchmarking genetic screens reveals a predominant mitochondrial bias

We present FLEX (Functional evaluation of experimental perturbations), a pipeline that leverages several functional annotation resources to establish reference standards for benchmarking human genome-wide CRISPR screen data and methods for analyzing them. FLEX provides a quantitative measurement of the functional information captured by a given gene-pair dataset and a means to explore the diversity of functions captured by the input dataset. We apply FLEX to analyze data from the diverse cell line screens generated by the DepMap project. We identify a predominant mitochondria-associated signal within co-essentiality networks derived from these data and explore the basis of this signal. Our analysis and time-resolved CRISPR screens in a single cell line suggest that the variable phenotypes associated with mitochondria genes across cells may reflect screen dynamics and protein stability effects rather than genetic dependencies. We characterize this functional bias and demonstrate its relevance for interpreting differential hits in any CRISPR screening context. More generally, we demonstrate the utility of the FLEX pipeline for performing robust comparative evaluations of CRISPR screens or methods for processing them.

Identification of Candidate Susceptibility Genes to Puccinia graminis f. sp. tritici in Wheat

Wheat stem rust disease caused by Puccinia graminis f. sp. tritici (Pgt) is a global threat to wheat production. Fast evolving populations of Pgt limit the efficacy of plant genetic resistance and constrain disease management strategies. Understanding molecular mechanisms that lead to rust infection and disease susceptibility could deliver novel strategies to deploy crop resistance through genetic loss of disease susceptibility. We used comparative transcriptome-based and orthology-guided approaches to characterize gene expression changes associated with Pgt infection in susceptible and resistant Triticum aestivum genotypes as well as the non-host Brachypodium distachyon. We targeted our analysis to genes with differential expression in T. aestivum and genes suppressed or not affected in B. distachyon and report several processes potentially linked to susceptibility to Pgt, such as cell death suppression and impairment of photosynthesis. We complemented our approach with a gene co-expression network analysis to identify wheat targets to deliver resistance to Pgt through removal or modification of putative susceptibility genes.

Revealing Functional Crosstalk between Distinct Bioprocesses through Reciprocal Functional Tests of Genetically Interacting Genes

To systematically explore the genes mediating functional crosstalk between metazoan biological processes, we apply comparative genetic interaction (GI) mapping in Saccharomyces cerevisiae and Caenorhabditis elegans to generate an inter-bioprocess network consisting of 178 C. elegans GIs. The GI network spans six annotated biological processes including aging, intracellular transport, microtubule-based processes, cytokinesis, lipid metabolic processes, and anatomical structure development. By proposing a strategy called "reciprocal functional test" for interacting gene pairs, we discover a group of genes that mediate crosstalk between distinct biological processes. In particular, we identify the ribosomal S6 Kinase/RSKS-1, previously characterized as an mTOR (mechanistic target of rapamycin) effector, as a regulator of DAF-2 endosomal recycling transport, which traces a functional correlation between endocytic recycling and aging processes. Together, our results provide an alternative and effective strategy for identifying genes and pathways that mediate crosstalk between bioprocesses with little prior knowledge.

EXO1 resection at G-quadruplex structures facilitates resolution and replication

G-quadruplexes represent unique roadblocks to DNA replication, which tends to stall at these secondary structures. Although G-quadruplexes can be found throughout the genome, telomeres, due to their G-richness, are particularly predisposed to forming these structures and thus represent difficult-to-replicate regions. Here, we demonstrate that exonuclease 1 (EXO1) plays a key role in the resolution of, and replication through, telomeric G-quadruplexes. When replication forks encounter G-quadruplexes, EXO1 resects the nascent DNA proximal to these structures to facilitate fork progression and faithful replication. In the absence of EXO1, forks accumulate at stabilized G-quadruplexes and ultimately collapse. These collapsed forks are preferentially repaired via error-prone end joining as depletion of EXO1 diverts repair away from error-free homology-dependent repair. Such aberrant repair leads to increased genomic instability, which is exacerbated at chromosome termini in the form of dysfunction and telomere loss.

Exploring whole-genome duplicate gene retention with complex genetic interaction analysis

Whole-genome duplication has played a central role in the genome evolution of many organisms, including the human genome. Most duplicated genes are eliminated, and factors that influence the retention of persisting duplicates remain poorly understood. We describe a systematic complex genetic interaction analysis with yeast paralogs derived from the whole-genome duplication event. Mapping of digenic interactions for a deletion mutant of each paralog, and of trigenic interactions for the double mutant, provides insight into their roles and a quantitative measure of their functional redundancy. Trigenic interaction analysis distinguishes two classes of paralogs: a more functionally divergent subset and another that retained more functional overlap. Gene feature analysis and modeling suggest that evolutionary trajectories of duplicated genes are dictated by combined functional and structural entanglement factors.

Systematic mapping of genetic interactions for de novo fatty acid synthesis identifies C12orf49 as a regulator of lipid metabolism

The de novo synthesis of fatty acids has emerged as a therapeutic target for various diseases, including cancer. Because cancer cells are intrinsically buffered to combat metabolic stress, it is important to understand how cells may adapt to the loss of de novo fatty acid biosynthesis. Here, we use pooled genome-wide CRISPR screens to systematically map genetic interactions (GIs) in human HAP1 cells carrying a loss-of-function mutation in fatty acid synthase (FASN), whose product catalyses the formation of long-chain fatty acids. FASN-mutant cells show a strong dependence on lipid uptake that is reflected in negative GIs with genes involved in the LDL receptor pathway, vesicle trafficking and protein glycosylation. Further support for these functional relationships is derived from additional GI screens in query cell lines deficient in other genes involved in lipid metabolism, including LDLR, SREBF1, SREBF2 and ACACA. Our GI profiles also identify a potential role for the previously uncharacterized gene C12orf49 (which we call LUR1) in regulation of exogenous lipid uptake through modulation of SREBF2 signalling in response to lipid starvation. Overall, our data highlight the genetic determinants underlying the cellular adaptation associated with loss of de novo fatty acid synthesis and demonstrate the power of systematic GI mapping for uncovering metabolic buffering mechanisms in human cells.

Identification of nodulation-related genes in Medicago truncatula using genome-wide association studies and co-expression networks

Genome-wide association studies (GWAS) have proven to be a valuable approach for identifying genetic intervals associated with phenotypic variation in Medicago truncatula. These intervals can vary in size, depending on the historical local recombination. Typically, significant intervals span numerous gene models, limiting the ability to resolve high-confidence candidate genes underlying the trait of interest. Additional genomic data, including gene co-expression networks, can be combined with the genetic mapping information to successfully identify candidate genes. Co-expression network analysis provides information about the functional relationships of each gene through its similarity of expression patterns to other well-defined clusters of genes. In this study, we integrated data from GWAS and co-expression networks to pinpoint candidate genes that may be associated with nodule-related phenotypes in M. truncatula. We further investigated a subset of these genes and confirmed that several had existing evidence linking them nodulation, including MEDTR2G101090 (PEN3-like), a previously validated gene associated with nodule number.

Genetic interaction mapping and exon-resolution functional genomics with a hybrid Cas9-Cas12a platform

Systematic mapping of genetic interactions (GIs) and interrogation of the functions of sizable genomic segments in mammalian cells represent important goals of biomedical research. To advance these goals, we present a CRISPR (clustered regularly interspaced short palindromic repeats)-based screening system for combinatorial genetic manipulation that employs coexpression of CRISPR-associated nucleases 9 and 12a (Cas9 and Cas12a) and machine-learning-optimized libraries of hybrid Cas9-Cas12a guide RNAs. This system, named Cas Hybrid for Multiplexed Editing and screening Applications (CHyMErA), outperforms genetic screens using Cas9 or Cas12a editing alone. Application of CHyMErA to the ablation of mammalian paralog gene pairs reveals extensive GIs and uncovers phenotypes normally masked by functional redundancy. Application of CHyMErA in a chemogenetic interaction screen identifies genes that impact cell growth in response to mTOR pathway inhibition. Moreover, by systematically targeting thousands of alternative splicing events, CHyMErA identifies exons underlying human cell line fitness. CHyMErA thus represents an effective screening approach for GI mapping and the functional analysis of sizable genomic regions, such as alternative exons.

Single-Cell Gene Expression Analyses Reveal Distinct Self-Renewing and Proliferating Subsets in the Leukemia Stem Cell Compartment in Acute Myeloid Leukemia

Standard chemotherapy for acute myeloid leukemia (AML) targets proliferative cells and efficiently induces complete remission; however, many patients relapse and die of their disease. Relapse is caused by leukemia stem cells (LSC), the cells with self-renewal capacity. Self-renewal and proliferation are separate functions in normal hematopoietic stem cells (HSC) in steady-state conditions. If these functions are also separate functions in LSCs, then antiproliferative therapies may fail to target self-renewal, allowing for relapse. We investigated whether proliferation and self-renewal are separate functions in LSCs as they often are in HSCs. Distinct transcriptional profiles within LSCs of Mll-AF9/NRAS^G12V murine AML were identified using single-cell RNA sequencing. Single-cell qPCR revealed that these genes were also differentially expressed in primary human LSCs and normal human HSPCs. A smaller subset of these genes was upregulated in LSCs relative to HSPCs; this subset of genes constitutes "LSC-specific" genes in human AML. To assess the differences between these profiles, we identified cell surface markers, CD69 and CD36, whose genes were differentially expressed between these profiles. In vivo mouse reconstitution assays resealed that only CD69^High LSCs were capable of self-renewal and were poorly proliferative. In contrast, CD36^High LSCs were unable to transplant leukemia but were highly proliferative. These data demonstrate that the transcriptional foundations of self-renewal and proliferation are distinct in LSCs as they often are in normal stem cells and suggest that therapeutic strategies that target self-renewal, in addition to proliferation, are critical to prevent relapse and improve survival in AML. SIGNIFICANCE: These findings define and functionally validate a self-renewal gene profile of leukemia stem cells at the single-cell level and demonstrate that self-renewal and proliferation are distinct in AML. GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/80/3/458/F1.large.jpg.

Integrating yeast chemical genomics and mammalian cell pathway analysis

Chemical genomics has been applied extensively to evaluate small molecules that modulate biological processes in Saccharomyces cerevisiae. Here, we use yeast as a surrogate system for studying compounds that are active against metazoan targets. Large-scale chemical-genetic profiling of thousands of synthetic and natural compounds from the Chinese National Compound Library identified those with high-confidence bioprocess target predictions. To discover compounds that have the potential to function like therapeutic agents with known targets, we also analyzed a reference library of approved drugs. Previously uncharacterized compounds with chemical-genetic profiles resembling existing drugs that modulate autophagy and Wnt/β-catenin signal transduction were further examined in mammalian cells, and new modulators with specific modes of action were validated. This analysis exploits yeast as a general platform for predicting compound bioactivity in mammalian cells.

Using BEAN-counter to quantify genetic interactions from multiplexed barcode sequencing experiments

The construction of genome-wide mutant collections has enabled high-throughput, high-dimensional quantitative characterization of gene and chemical function, particularly via genetic and chemical-genetic interaction experiments. As the throughput of such experiments increases with improvements in sequencing technology and sample multiplexing, appropriate tools must be developed to handle the large volume of data produced. Here, we describe how to apply our approach to high-throughput, fitness-based profiling of pooled mutant yeast collections using the BEAN-counter software pipeline (Barcoded Experiment Analysis for Next-generation sequencing) for analysis. The software has also successfully processed data from Schizosaccharomyces pombe, Escherichia coli, and Zymomonas mobilis mutant collections. We provide general recommendations for the design of large-scale, multiplexed barcode sequencing experiments. The procedure outlined here was used to score interactions for ~4 million chemical-by-mutant combinations in our recently published chemical-genetic interaction screen of nearly 14,000 chemical compounds across seven diverse compound collections. Here we selected a representative subset of these data on which to demonstrate our analysis pipeline. BEAN-counter is open source, written in Python, and freely available for academic use. Users should be proficient at the command line; advanced users who wish to analyze larger datasets with hundreds or more conditions should also be familiar with concepts in analysis of high-throughput biological data. BEAN-counter encapsulates the knowledge we have accumulated from, and successfully applied to, our multiplexed, pooled barcode sequencing experiments. This protocol will be useful to those interested in generating their own high-dimensional, quantitative characterizations of gene or chemical function in a high-throughput manner.

Domain insertion permissibility-guided engineering of allostery in ion channels

Allostery is a fundamental principle of protein regulation that remains hard to engineer, particularly in membrane proteins such as ion channels. Here we use human Inward Rectifier K⁺ Channel Kir2.1 to map site-specific permissibility to the insertion of domains with different biophysical properties. We find that permissibility is best explained by dynamic protein properties, such as conformational flexibility. Several regions in Kir2.1 that are equivalent to those regulated in homologs, such as G-protein-gated inward rectifier K⁺ channels (GIRK), have differential permissibility; that is, for these sites permissibility depends on the structural properties of the inserted domain. Our data and the well-established link between protein dynamics and allostery led us to propose that differential permissibility is a metric of latent allosteric capacity in Kir2.1. In support of this notion, inserting light-switchable domains into sites with predicted latent allosteric capacity renders Kir2.1 activity sensitive to light.

Allostery is a fundamental principle of protein regulation that remains challenging to engineer. Here authors screen human Inward Rectifier K + Channel Kir2.1 for permissibility to domain insertions and propose that differential permissibility is a metric of latent allosteric capacity in Kir2.1.

MOSAIC: a chemical-genetic interaction data repository and web resource for exploring chemical modes of action

Summary

Chemical-genomic approaches that map interactions between small molecules and genetic perturbations offer a promising strategy for functional annotation of uncharacterized bioactive compounds. We recently developed a new high-throughput platform for mapping chemical-genetic (CG) interactions in yeast that can be scaled to screen large compound collections, and we applied this system to generate CG interaction profiles for more than 13 000 compounds. When integrated with the existing global yeast genetic interaction network, CG interaction profiles can enable mode-of-action prediction for previously uncharacterized compounds as well as discover unexpected secondary effects for known drugs. To facilitate future analysis of these valuable data, we developed a public database and web interface named MOSAIC. The website provides a convenient interface for querying compounds, bioprocesses (Gene Ontology terms) and genes for CG information including direct CG interactions, bioprocesses and gene-level target predictions. MOSAIC also provides access to chemical structure information of screened molecules, chemical-genomic profiles and the ability to search for compounds sharing structural and functional similarity. This resource will be of interest to chemical biologists for discovering new small molecule probes with specific modes-of-action as well as computational biologists interested in analysing CG interaction networks.

Availability and implementation

MOSAIC is available at http://mosaic.cs.umn.edu.

Contact

hisyo@riken.jp, yoshidam@riken.jp, charlie.boone@utoronto.ca or chadm@umn.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

Integrating Coexpression Networks with GWAS to Prioritize Causal Genes in Maize

Genome-wide association studies (GWAS) have identified loci linked to hundreds of traits in many different species. Yet, because linkage equilibrium implicates a broad region surrounding each identified locus, the causal genes often remain unknown. This problem is especially pronounced in nonhuman, nonmodel species, where functional annotations are sparse and there is frequently little information available for prioritizing candidate genes. We developed a computational approach, Camoco, that integrates loci identified by GWAS with functional information derived from gene coexpression networks. Using Camoco, we prioritized candidate genes from a large-scale GWAS examining the accumulation of 17 different elements in maize (Zea mays) seeds. Strikingly, we observed a strong dependence in the performance of our approach based on the type of coexpression network used: expression variation across genetically diverse individuals in a relevant tissue context (in our case, roots that are the primary elemental uptake and delivery system) outperformed other alternative networks. Two candidate genes identified by our approach were validated using mutants. Our study demonstrates that coexpression networks provide a powerful basis for prioritizing candidate causal genes from GWAS loci but suggests that the success of such strategies can highly depend on the gene expression data context. Both the software and the lessons on integrating GWAS data with coexpression networks generalize to species beyond maize.

Predicting bioprocess targets of chemical compounds through integration of chemical-genetic and genetic interactions

Chemical-genetic interactions–observed when the treatment of mutant cells with chemical compounds reveals unexpected phenotypes–contain rich functional information linking compounds to their cellular modes of action. To systematically identify these interactions, an array of mutants is challenged with a compound and monitored for fitness defects, generating a chemical-genetic interaction profile that provides a quantitative, unbiased description of the cellular function(s) perturbed by the compound. Genetic interactions, obtained from genome-wide double-mutant screens, provide a key for interpreting the functional information contained in chemical-genetic interaction profiles. Despite the utility of this approach, integrative analyses of genetic and chemical-genetic interaction networks have not been systematically evaluated. We developed a method, called CG-TARGET (Chemical Genetic Translation via A Reference Genetic nETwork), that integrates large-scale chemical-genetic interaction screening data with a genetic interaction network to predict the biological processes perturbed by compounds. In a recent publication, we applied CG-TARGET to a screen of nearly 14,000 chemical compounds in Saccharomyces cerevisiae, integrating this dataset with the global S. cerevisiae genetic interaction network to prioritize over 1500 compounds with high-confidence biological process predictions for further study. We present here a formal description and rigorous benchmarking of the CG-TARGET method, showing that, compared to alternative enrichment-based approaches, it achieves similar or better accuracy while substantially improving the ability to control the false discovery rate of biological process predictions. Additional investigation of the compatibility of chemical-genetic and genetic interaction profiles revealed that one-third of observed chemical-genetic interactions contributed to the highest-confidence biological process predictions and that negative chemical-genetic interactions overwhelmingly formed the basis of these predictions. We also present experimental validations of CG-TARGET-predicted tubulin polymerization and cell cycle progression inhibitors. Our approach successfully demonstrates the use of genetic interaction networks in the high-throughput functional annotation of compounds to biological processes.

Understanding how chemical compounds affect biological systems is of paramount importance as pharmaceutical companies strive to develop life-saving medicines, governments seek to regulate the safety of consumer products and agrichemicals, and basic scientists continue to study the fundamental inner workings of biological organisms. One powerful approach to characterize the effects of chemical compounds in living cells is chemical-genetic interaction screening. Using this approach, a collection of cells–each with a different defined genetic perturbation–is tested for sensitivity or resistance to the presence of a compound, resulting in a quantitative profile describing the functional effects of that compound on the cells. The work presented here describes our efforts to integrate compounds’ chemical-genetic interaction profiles with reference genetic interaction profiles containing information on gene function to predict the cellular processes perturbed by the compounds. We focused on specifically developing a method that could scale to perform these functional predictions for large collections of thousands of screened compounds and robustly control the false discovery rate. With chemical-genetic and genetic interaction screens now underway in multiple species including human cells, the method described here can be generally applied to enable the characterization of compounds’ effects across the tree of life.

Integrating genetic and protein-protein interaction networks maps a functional wiring diagram of a cell

Systematic experimental approaches have led to construction of comprehensive genetic and protein-protein interaction networks for the budding yeast, Saccharomyces cerevisiae. Genetic interactions capture functional relationships between genes using phenotypic readouts, while protein-protein interactions identify physical connections between gene products. These complementary, and largely non-overlapping, networks provide a global view of the functional architecture of a cell, revealing general organizing principles, many of which appear to be evolutionarily conserved. Here, we focus on insights derived from the integration of large-scale genetic and protein-protein interaction networks, highlighting principles that apply to both unicellular and more complex systems, including human cells. Network integration reveals fundamental connections involving key functional modules of eukaryotic cells, defining a core network of cellular function, which could be elaborated to explore cell-type specificity in metazoans.

In Vitro Pluripotent Stem Cell Differentiation to Hepatocyte Ceases Further Maturation at an Equivalent Stage of E15 in Mouse Embryonic Liver Development

Hepatocyte-like cells (HLCs) can be derived from pluripotent stem cells (PSCs) by sequential treatment of chemical cues to mimic the microenvironment of embryonic liver development. However, these HLCs do not reach the full maturity level of primary hepatocytes. In this study, we carried out a meta-analysis of cross-species transcriptome data of in vitro differentiation of human PSCs to HLCs and in vivo mouse embryonic liver development to identify the developmental stage at which HLC maturation was blocked at. Systematic variations were found associated with the data source and removed by batch correction. Using principal component analysis, HLCs from different stages of differentiation were aligned with mouse embryonic liver development chronologically. A "unified developmental time" (DT) scale was developed after aligning in vitro HLC differentiation and in vivo embryonic liver development. HLCs were found to cease further maturation at an equivalent stage of mouse embryonic day (E)13-15. Genes with discordant time dynamics were identified by aligning in vivo and in vitro data set onto a common DT scale. These genes may be targets of genetic intervention for enhancing the maturity of PSC-derived HLCs.

Features of the Chaperone Cellular Network Revealed through Systematic Interaction Mapping

A comprehensive view of molecular chaperone function in the cell was obtained through a systematic global integrative network approach based on physical (protein-protein) and genetic (gene-gene or epistatic) interaction mapping. This allowed us to decipher interactions involving all core chaperones (67) and cochaperones (15) of Saccharomyces cerevisiae. Our analysis revealed the presence of a large chaperone functional supercomplex, which we named the naturally joined (NAJ) chaperone complex, encompassing Hsp40, Hsp70, Hsp90, AAA+, CCT, and small Hsps. We further found that many chaperones interact with proteins that form foci or condensates under stress conditions. Using an in vitro reconstitution approach, we demonstrate condensate formation for the highly conserved AAA+ ATPases Rvb1 and Rvb2, which are part of the R2TP complex that interacts with Hsp90. This expanded view of the chaperone network in the cell clearly demonstrates the distinction between chaperones having broad versus narrow substrate specificities in protein homeostasis.

Systematic analysis of complex genetic interactions

To systematically explore complex genetic interactions, we constructed ~200,000 yeast triple mutants and scored negative trigenic interactions. We selected double-mutant query genes across a broad spectrum of biological processes, spanning a range of quantitative features of the global digenic interaction network and tested for a genetic interaction with a third mutation. Trigenic interactions often occurred among functionally related genes, and essential genes were hubs on the trigenic network. Despite their functional enrichment, trigenic interactions tended to link genes in distant bioprocesses and displayed a weaker magnitude than digenic interactions. We estimate that the global trigenic interaction network is ~100 times as large as the global digenic network, highlighting the potential for complex genetic interactions to affect the biology of inheritance, including the genotype-to-phenotype relationship.

Different Modes of Negative Regulation of Plant Immunity by Calmodulin-Related Genes

Plant immune responses activated through the perception of microbe-associated molecular patterns, leading to pattern-triggered immunity, are tightly regulated. This results in low immune responses in the absence of pathogens and a rapid return to the resting state following an activation event. Here, we show that two CALMODULIN-LIKE genes, CML46 and CML47, negatively regulate salicylic acid accumulation and immunity in Arabidopsis (Arabidopsis thaliana). The double mutant cml46 cml47 is highly resistant to the pathogen Pseudomonas syringae pv maculicola (Pma). The effects of cml46 cml47 on Pma growth are genetically additive to that of cbp60a, a known negative regulator in the CALMODULIN-BINDING PROTEIN60 (CBP60) family. Transcriptome profiling revealed the effects of cbp60a and cml46 cml47 on both common and separate sets of genes, with the majorities of these differentially expressed genes being Pma responsive. CBP60g, a positive regulator of immunity in the CBP60 family, was found to be transcriptionally regulated by CBP60a, CML46, and CML47 Analysis of the flg22-induced mRNA levels of CBP60g in cbp60a and cml46 cml47 revealed that cml46 cml47 plants have higher induced expression while cbp60a plants retain elevated levels longer than wild-type plants. Assays for the effect of flg22 treatment on Pma growth showed that the effect is stronger in cml46 cml47 plants and lasts longer in cbp60a plants. Thus, the expression pattern of CBP60g is reflected in flg22-induced resistance to Pma.

Translocon Declogger Ste24 Protects against IAPP Oligomer-Induced Proteotoxicity

Aggregates of human islet amyloid polypeptide (IAPP) in the pancreas of patients with type 2 diabetes (T2D) are thought to contribute to β cell dysfunction and death. To understand how IAPP harms cells and how this might be overcome, we created a yeast model of IAPP toxicity. Ste24, an evolutionarily conserved protease that was recently reported to degrade peptides stuck within the translocon between the cytoplasm and the endoplasmic reticulum, was the strongest suppressor of IAPP toxicity. By testing variants of the human homolog, ZMPSTE24, with varying activity levels, the rescue of IAPP toxicity proved to be directly proportional to the declogging efficiency. Clinically relevant ZMPSTE24 variants identified in the largest database of exomes sequences derived from T2D patients were characterized using the yeast model, revealing 14 partial loss-of-function variants, which were enriched among diabetes patients over 2-fold. Thus, clogging of the translocon by IAPP oligomers may contribute to β cell failure.

Chemical genomic guided engineering of gamma-valerolactone tolerant yeast

Background

Gamma valerolactone (GVL) treatment of lignocellulosic bomass is a promising technology for degradation of biomass for biofuel production; however, GVL is toxic to fermentative microbes. Using a combination of chemical genomics with the yeast (Saccharomyces cerevisiae) deletion collection to identify sensitive and resistant mutants, and chemical proteomics to monitor protein abundance in the presence of GVL, we sought to understand the mechanism toxicity and resistance to GVL with the goal of engineering a GVL-tolerant, xylose-fermenting yeast.

Results

Chemical genomic profiling of GVL predicted that this chemical affects membranes and membrane-bound processes. We show that GVL causes rapid, dose-dependent cell permeability, and is synergistic with ethanol. Chemical genomic profiling of GVL revealed that deletion of the functionally related enzymes Pad1p and Fdc1p, which act together to decarboxylate cinnamic acid and its derivatives to vinyl forms, increases yeast tolerance to GVL. Further, overexpression of Pad1p sensitizes cells to GVL toxicity. To improve GVL tolerance, we deleted PAD1 and FDC1 in a xylose-fermenting yeast strain. The modified strain exhibited increased anaerobic growth, sugar utilization, and ethanol production in synthetic hydrolysate with 1.5% GVL, and under other conditions. Chemical proteomic profiling of the engineered strain revealed that enzymes involved in ergosterol biosynthesis were more abundant in the presence of GVL compared to the background strain. The engineered GVL strain contained greater amounts of ergosterol than the background strain.

Conclusions

We found that GVL exerts toxicity to yeast by compromising cellular membranes, and that this toxicity is synergistic with ethanol. Deletion of PAD1 and FDC1 conferred GVL resistance to a xylose-fermenting yeast strain by increasing ergosterol accumulation in aerobically grown cells. The GVL-tolerant strain fermented sugars in the presence of GVL levels that were inhibitory to the unmodified strain. This strain represents a xylose fermenting yeast specifically tailored to GVL produced hydrolysates.

Electronic supplementary material

The online version of this article (10.1186/s12934-017-0848-9) contains supplementary material, which is available to authorized users.

Coculture of Marine Invertebrate-Associated Bacteria and Interdisciplinary Technologies Enable Biosynthesis and Discovery of a New Antibiotic, Keyicin

Advances in genomics and metabolomics have made clear in recent years that microbial biosynthetic capacities on Earth far exceed previous expectations. This is attributable, in part, to the realization that most microbial natural product (NP) producers harbor biosynthetic machineries not readily amenable to classical laboratory fermentation conditions. Such "cryptic" or dormant biosynthetic gene clusters (BGCs) encode for a vast assortment of potentially new antibiotics and, as such, have become extremely attractive targets for activation under controlled laboratory conditions. We report here that coculturing of a Rhodococcus sp. and a Micromonospora sp. affords keyicin, a new and otherwise unattainable bis-nitroglycosylated anthracycline whose mechanism of action (MOA) appears to deviate from those of other anthracyclines. The structure of keyicin was elucidated using high resolution MS and NMR technologies, as well as detailed molecular modeling studies. Sequencing of the keyicin BGC (within the Micromonospora genome) enabled both structural and genomic comparisons to other anthracycline-producing systems informing efforts to characterize keyicin. The new NP was found to be selectively active against Gram-positive bacteria including both Rhodococcus sp. and Mycobacterium sp. E. coli-based chemical genomics studies revealed that keyicin's MOA, in contrast to many other anthracyclines, does not invoke nucleic acid damage.

Errata: Functional annotation of chemical libraries across diverse biological processes

This corrects the article DOI: 10.1038/nchembio.2436.