Hydrogen peroxide sensitivity connects the activity of COX5A and NPR3 to the regulation of YAP1 expression

Reactive oxygen species (ROS) are among the most severe types of cellular stressors with the ability to damage essential cellular biomolecules. Excess levels of ROS are correlated with multiple pathophysiological conditions including neurodegeneration, diabetes, atherosclerosis, and cancer. Failure to regulate the severely imbalanced levels of ROS can ultimately lead to cell death, highlighting the importance of investigating the molecular mechanisms involved in the detoxification procedures that counteract the effects of these compounds in living organisms. One of the most abundant forms of ROS is H2 O2 , mainly produced by the electron transport chain in the mitochondria. Numerous genes have been identified as essential to the process of cellular detoxification. Yeast YAP1, which is homologous to mammalian AP-1 type transcriptional factors, has a key role in oxidative detoxification by upregulating the expression of antioxidant genes in yeast. The current study reveals novel functions for COX5A and NPR3 in H2 O2 -induced stress by demonstrating that their deletions result in a sensitive phenotype. Our follow-up investigations indicate that COX5A and NPR3 regulate the expression of YAP1 through an alternative mode of translation initiation. These novel gene functions expand our understanding of the regulation of gene expression and defense mechanism of yeast against oxidative stress.

The Involvement of YNR069C in Protein Synthesis in the Baker's Yeast, Saccharomyces cerevisiae

Maintaining translation fidelity is a critical step within the process of gene expression. It requires the involvement of numerous regulatory elements to ensure the synthesis of functional proteins. The efficient termination of protein synthesis can play a crucial role in preserving this fidelity. Here, we report on investigating a protein of unknown function, YNR069C (also known as BSC5), for its activity in the process of translation. We observed a significant increase in the bypass of premature stop codons upon the deletion of YNR069C. Interestingly, the genomic arrangement of this ORF suggests a compatible mode of expression reliant on translational readthrough, incorporating the neighboring open reading frame. We also showed that the deletion of YNR069C results in an increase in the rate of translation. Based on our results, we propose that YNR069C may play a role in translation fidelity, impacting the overall quantity and quality of translation. Our genetic interaction analysis supports our hypothesis, associating the role of YNR069C to the regulation of protein synthesis.

Differential gene expression provides leads to environmentally regulated soybean seed protein content

Soybean is an important global source of plant-based protein. A persistent trend has been observed over the past two decades that soybeans grown in western Canada have lower seed protein content than soybeans grown in eastern Canada. In this study, 10 soybean genotypes ranging in average seed protein content were grown in an eastern location (control) and three western locations (experimental) in Canada. Seed protein and oil contents were measured for all lines in each location. RNA-sequencing and differential gene expression analysis were used to identify differentially expressed genes that may account for relatively low protein content in western-grown soybeans. Differentially expressed genes were enriched for ontologies and pathways that included amino acid biosynthesis, circadian rhythm, starch metabolism, and lipid biosynthesis. Gene ontology, pathway mapping, and quantitative trait locus (QTL) mapping collectively provide a close inspection of mechanisms influencing nitrogen assimilation and amino acid biosynthesis between soybeans grown in the East and West. It was found that western-grown soybeans had persistent upregulation of asparaginase (an asparagine hydrolase) and persistent downregulation of asparagine synthetase across 30 individual differential expression datasets. This specific difference in asparagine metabolism between growing environments is almost certainly related to the observed differences in seed protein content because of the positive correlation between seed protein content at maturity and free asparagine in the developing seed. These results provided pointed information on seed protein-related genes influenced by environment. This information is valuable for breeding programs and genetic engineering of geographically optimized soybeans.

Genome-wide analysis of cold imbibition stress in soybean, Glycine max

In Canada, the length of the frost-free season necessitates planting crops as early as possible to ensure that the plants have enough time to reach full maturity before they are harvested. Early planting carries inherent risks of cold water imbibition (specifically less than 4°C) affecting seed germination. A marker dataset developed for a previously identified Canadian soybean GWAS panel was leveraged to investigate the effect of cold water imbibition on germination. Seed from a panel of 137 soybean elite cultivars, grown in the field at Ottawa, ON, over three years, were placed on filter paper in petri dishes and allowed to imbibe water for 16 hours at either 4°C or 20°C prior to being transferred to a constant 20°C. Observations on seed germination, defined as the presence of a 1 cm radicle, were done from day two to seven. A three-parameter exponential rise to a maximum equation (3PERM) was fitted to estimate germination, time to the one-half maximum germination, and germination uniformity for each cultivar. Genotype-by-sequencing was used to identify SNPs in 137 soybean lines, and using genome-wide association studies (GWAS - rMVP R package, with GLM, MLM, and FarmCPU as methods), haplotype block analysis, and assumed linkage blocks of ±100 kbp, a threshold for significance was established using the qvalue package in R, and five significant SNPs were identified on chromosomes 1, 3, 4, 6, and 13 for maximum germination after cold water imbibition. Percent of phenotypic variance explained (PVE) and allele substitution effect (ASE) eliminated two of the five candidate SNPs, leaving three QTL regions on chromosomes 3, 6, and 13 (Chr3-3419152, Chr6-5098454, and Chr13-29649544). Based on the gene ontology (GO) enrichment analysis, 14 candidate genes whose function is predicted to include germination and cold tolerance related pathways were identified as candidate genes. The identified QTLs can be used to select future soybean cultivars tolerant to cold water imbibition and mitigate risks associated with early soybean planting.

Large-scale data mining pipeline for identifying novel soybean genes involved in resistance against the soybean cyst nematode

The soybean cyst nematode (SCN) [Heterodera glycines Ichinohe] is a devastating pathogen of soybean [Glycine max (L.) Merr.] that is rapidly becoming a global economic issue. Two loci conferring SCN resistance have been identified in soybean, Rhg1 and Rhg4; however, they offer declining protection. Therefore, it is imperative that we identify additional mechanisms for SCN resistance. In this paper, we develop a bioinformatics pipeline to identify protein-protein interactions related to SCN resistance by data mining massive-scale datasets. The pipeline combines two leading sequence-based protein-protein interaction predictors, the Protein-protein Interaction Prediction Engine (PIPE), PIPE4, and Scoring PRotein INTeractions (SPRINT) to predict high-confidence interactomes. First, we predicted the top soy interacting protein partners of the Rhg1 and Rhg4 proteins. Both PIPE4 and SPRINT overlap in their predictions with 58 soybean interacting partners, 19 of which had GO terms related to defense. Beginning with the top predicted interactors of Rhg1 and Rhg4, we implement a "guilt by association" in silico proteome-wide approach to identify novel soybean genes that may be involved in SCN resistance. This pipeline identified 1,082 candidate genes whose local interactomes overlap significantly with the Rhg1 and Rhg4 interactomes. Using GO enrichment tools, we highlighted many important genes including five genes with GO terms related to response to the nematode (GO:0009624), namely, Glyma.18G029000, Glyma.11G228300, Glyma.08G120500, Glyma.17G152300, and Glyma.08G265700. This study is the first of its kind to predict interacting partners of known resistance proteins Rhg1 and Rhg4, forming an analysis pipeline that enables researchers to focus their search on high-confidence targets to identify novel SCN resistance genes in soybean.

Peptides of a Feather: How Computation Is Taking Peptide Therapeutics under Its Wing

Leveraging computation in the development of peptide therapeutics has garnered increasing recognition as a valuable tool to generate novel therapeutics for disease-related targets. To this end, computation has transformed the field of peptide design through identifying novel therapeutics that exhibit enhanced pharmacokinetic properties and reduced toxicity. The process of in-silico peptide design involves the application of molecular docking, molecular dynamics simulations, and machine learning algorithms. Three primary approaches for peptide therapeutic design including structural-based, protein mimicry, and short motif design have been predominantly adopted. Despite the ongoing progress made in this field, there are still significant challenges pertaining to peptide design including: enhancing the accuracy of computational methods; improving the success rate of preclinical and clinical trials; and developing better strategies to predict pharmacokinetics and toxicity. In this review, we discuss past and present research pertaining to the design and development of in-silico peptide therapeutics in addition to highlighting the potential of computation and artificial intelligence in the future of disease therapeutics.

Synthetic virology approaches to improve the safety and efficacy of oncolytic virus therapies

The large coding potential of vaccinia virus (VV) vectors is a defining feature. However, limited regulatory switches are available to control viral replication as well as timing and dosing of transgene expression in order to facilitate safe and efficacious payload delivery. Herein, we adapt drug-controlled gene switches to enable control of virally encoded transgene expression, including systems controlled by the FDA-approved rapamycin and doxycycline. Using ribosome profiling to characterize viral promoter strength, we rationally design fusions of the operator element of different drug-inducible systems with VV promoters to produce synthetic promoters yielding robust inducible expression with undetectable baseline levels. We also generate chimeric synthetic promoters facilitating additional regulatory layers for VV-encoded synthetic transgene networks. The switches are applied to enable inducible expression of fusogenic proteins, dose-controlled delivery of toxic cytokines, and chemical regulation of VV replication. This toolbox enables the precise modulation of transgene circuitry in VV-vectored oncolytic virus design.

A Correlation between 3'-UTR of OXA1 Gene and Yeast Mitochondrial Translation

Mitochondria possess their own DNA (mtDNA) and are capable of carrying out their transcription and translation. Although protein synthesis can take place in mitochondria, the majority of the proteins in mitochondria have nuclear origin. 3' and 5' untranslated regions of mRNAs (3'-UTR and 5'-UTR, respectively) are thought to play key roles in directing and regulating the activity of mitochondria mRNAs. Here we investigate the association between the presence of 3'-UTR from OXA1 gene on a prokaryotic reporter mRNA and mitochondrial translation in yeast. OXA1 is a nuclear gene that codes for mitochondrial inner membrane insertion protein and its 3'-UTR is shown to direct its mRNA toward mitochondria. It is not clear, however, if this mRNA may also be translated by mitochondria. In the current study, using a β-galactosidase reporter gene, we provide genetic evidence for a correlation between the presence of 3'-UTR of OXA1 on an mRNA and mitochondrial translation in yeast.

IRES-mediated translation in bacteria

Despite the similarity in fundamental goals of translation initiation between different domains of life, it is one of the most phylogenetically diverse steps of the central dogma of molecular biology. In a classical view, the translation signals for prokaryotes and eukaryotes are distinct from each other. This idea was challenged by the finding that the Internal Ribosome Entry Site (IRES) belonging to Plautia stali intestine virus (PSIV) could bypass the domain-specific boundaries and effectively initiate translation in E. coli. This finding led us to investigate whether the ability of PSIV IRES to initiate translation in E. coli is specific to this IRES and also to study features that allow this viral IRES to mediate prokaryotic translation initiation. We observed that certain IRESs may also possess the ability to initiate E. coli translation. Our results also indicated that the structural integrity of the PSIV IRES in translation in prokaryotes does not appear to be as critical as it is in eukaryotes. We also demonstrated that two regions of the PSIV IRES with complementarity to 16S ribosomal RNA are important for the ability of this IRES to initiate translation in E. coli.

DBP7 and YRF1-6 Are Involved in Cell Sensitivity to LiCl by Regulating the Translation of PGM2 mRNA

Lithium chloride (LiCl) has been widely researched and utilized as a therapeutic option for bipolar disorder (BD). Several pathways, including cell signaling and signal transduction pathways in mammalian cells, are shown to be regulated by LiCl. LiCl can negatively control the expression and activity of PGM2, a phosphoglucomutase that influences sugar metabolism in yeast. In the presence of galactose, when yeast cells are challenged by LiCl, the phosphoglucomutase activity of PGM2p is decreased, causing an increase in the concentration of toxic galactose metabolism intermediates that result in cell sensitivity. Here, we report that the null yeast mutant strains DBP7∆ and YRF1-6∆ exhibit increased LiCl sensitivity on galactose-containing media. Additionally, we demonstrate that DBP7 and YRF1-6 modulate the translational level of PGM2 mRNA, and the observed alteration in translation seems to be associated with the 5'-untranslated region (UTR) of PGM2 mRNA. Furthermore, we observe that DBP7 and YRF1-6 influence, to varying degrees, the translation of other mRNAs that carry different 5'-UTR secondary structures.

P-TarPmiR accurately predicts plant-specific miRNA targets

microRNAs (miRNAs) are small non-coding ribonucleic acids that post-transcriptionally regulate gene expression through the targeting of messenger RNA (mRNAs). Most miRNA target predictors have focused on animal species and prediction performance drops substantially when applied to plant species. Several rule-based miRNA target predictors have been developed in plant species, but they often fail to discover new miRNA targets with non-canonical miRNA-mRNA binding. Here, the recently published TarDB database of plant miRNA-mRNA data is leveraged to retrain the TarPmiR miRNA target predictor for application on plant species. Rigorous experiment design across four plant test species demonstrates that animal-trained predictors fail to sustain performance on plant species, and that the use of plant-specific training data improves accuracy depending on the quantity of plant training data used. Surprisingly, our results indicate that the complete exclusion of animal training data leads to the most accurate plant-specific miRNA target predictor indicating that animal-based data may detract from miRNA target prediction in plants. Our final plant-specific miRNA prediction method, dubbed P-TarPmiR, is freely available for use at http://ptarpmir.cu-bic.ca . The final P-TarPmiR method is used to predict targets for all miRNA within the soybean genome. Those ranked predictions, together with GO term enrichment, are shared with the research community.

A Multi-Year, Multi-Cultivar Approach to Differential Expression Analysis of High- and Low-Protein Soybean (Glycine max)

Soybean (Glycine max (L.) Merr.) is among the most valuable crops based on its nutritious seed protein and oil. Protein quality, evaluated as the ratio of glycinin (11S) to β-conglycinin (7S), can play a role in food and feed quality. To help uncover the underlying differences between high and low protein soybean varieties, we performed differential expression analysis on high and low total protein soybean varieties and high and low 11S soybean varieties grown in four locations across Eastern and Western Canada over three years (2018-2020). Simultaneously, ten individual differential expression datasets for high vs. low total protein soybeans and ten individual differential expression datasets for high vs. low 11S soybeans were assessed, for a total of 20 datasets. The top 15 most upregulated and the 15 most downregulated genes were extracted from each differential expression dataset and cross-examination was conducted to create shortlists of the most consistently differentially expressed genes. Shortlisted genes were assessed for gene ontology to gain a global appreciation of the commonly differentially expressed genes. Genes with roles in the lipid metabolic pathway and carbohydrate metabolic pathway were differentially expressed in high total protein and high 11S soybeans in comparison to their low total protein and low 11S counterparts. Expression differences were consistent between East and West locations with the exception of one, Glyma.03G054100. These data are important for uncovering the genes and biological pathways responsible for the difference in seed protein between high and low total protein or 11S cultivars.

GmSWEET29 and Paralog GmSWEET34 Are Differentially Expressed between Soybeans Grown in Eastern and Western Canada

Over the past two decades soybeans grown in western Canada have persistently had lower seed protein than those grown in eastern Canada. To understand the discrepancy in seed protein content between eastern- and western-grown soybeans, RNA-seq and differential expression analysis have been investigated. Ten soybean genotypes, ranging from low to high in seed protein content, were grown in four locations across eastern (Ottawa) and western (Morden, Brandon, and Saskatoon) Canada. Differential expression analysis revealed 34 differentially expressed genes encoding Glycine max Sugars Will Eventually be Exported Transporters (GmSWEETs), including paralogs GmSWEET29 and GmSWEET34 (AtSWEET2 homologs) that were consistently upregulated across all ten genotypes in each of the western locations over three years. GmSWEET29 and GmSWEET34 are likely candidates underlying the lower seed protein content of western soybeans. GmSWEET20 (AtSWEET12 homolog) was downregulated in the western locations and may also play a role in lower seed protein content. These findings are valuable for improving soybean agriculture in western growing regions, establishing more strategic and efficient agricultural practices.

Discovery and identification of genes involved in DNA damage repair in yeast

DNA repair defects are common in tumour cells and can lead to misrepair of double-strand breaks (DSBs), posing a significant challenge to cellular integrity. The overall mechanisms of DSB have been known for decades. However, the list of the genes that affect the efficiency of DSB repair continues to grow. Additional factors that play a role in DSB repair pathways have yet to be identified. In this study, we present a computational approach to identify novel gene functions that are involved in DNA damage repair in Saccharomyces cerevisiae. Among the primary candidates, GAL7, YMR130W, and YHI9 were selected for further analysis since they had not previously been identified as being active in DNA repair pathways. Originally, GAL7 was linked to galactose metabolism. YHI9 and YMR130W encode proteins of unknown functions. Laboratory testing of deletion strains gal7Δ, ymr130wΔ, and yhi9Δ implicated all 3 genes in Homologous Recombination (HR) and/or Non-Homologous End Joining (NHEJ) repair pathways, and enhanced sensitivity to DNA damage-inducing drugs suggested involvement in the broader DNA damage repair machinery. A subsequent genetic interaction analysis revealed interconnections of these three genes, most strikingly through SIR2, SIR3 and SIR4 that are involved in chromatin regulation and DNA damage repair network.

A computational approach to rapidly design peptides that detect SARS-CoV-2 surface protein S

The coronavirus disease 19 (COVID-19) caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) prompted the development of diagnostic and therapeutic frameworks for timely containment of this pandemic. Here, we utilized our non-conventional computational algorithm, InSiPS, to rapidly design and experimentally validate peptides that bind to SARS-CoV-2 spike (S) surface protein. We previously showed that this method can be used to develop peptides against yeast proteins, however, the applicability of this method to design peptides against other proteins has not been investigated. In the current study, we demonstrate that two sets of peptides developed using InSiPS method can detect purified SARS-CoV-2 S protein via ELISA and Surface Plasmon Resonance (SPR) approaches, suggesting the utility of our strategy in real time COVID-19 diagnostics. Mass spectrometry-based salivary peptidomics shortlist top SARS-CoV-2 peptides detected in COVID-19 patients’ saliva, rendering them attractive SARS-CoV-2 diagnostic targets that, when subjected to our computational platform, can streamline the development of potent peptide diagnostics of SARS-CoV-2 variants of concern. Our approach can be rapidly implicated in diagnosing other communicable diseases of immediate threat.

Novel QTL for Low Seed Cadmium Accumulation in Soybean

Soybean is a valuable crop, used in animal feed and for human consumption. Selecting soybean cultivars with low seed cadmium (Cd) concentration is important for the purpose of minimizing the transfer of Cd into the human body. To ensure international trade, farmers need to produce soybean that meets the European Union (EU) Cd limit of 0.2 mg kg-1. In this study, we evaluated two populations of recombinant inbred lines (RILs), X5154 and X4050, for seed Cd accumulation. Linkage maps were constructed with 325 and 280 polymorphic simple sequence repeat (SSR) markers, respectively, and used to identify a novel minor quantitative trait locus (QTL) on chromosome 13 in the X4050 population between SSR markers Satt522 and Satt218. Based on a gene ontology search within the QTL region, seven genes were identified as candidates responsible for low seed Cd accumulation, including Glyma.13G308700 and Glyma.13G309100. In addition, we confirmed the known major gene, Cda1, in the X5154 population and developed KASP and CAPS/dCAPS allele-specific markers for efficient marker-assisted breeding for Cda1.

A Broad Review of Soybean Research on the Ongoing Race to Overcome Soybean Cyst Nematode

Plant pathogens greatly impact food security of the ever-growing human population. Breeding resistant crops is one of the most sustainable strategies to overcome the negative effects of these biotic stressors. In order to efficiently breed for resistant plants, the specific plant-pathogen interactions should be understood. Soybean is a short-day legume that is a staple in human food and animal feed due to its high nutritional content. Soybean cyst nematode (SCN) is a major soybean stressor infecting soybean worldwide including in China, Brazil, Argentina, USA and Canada. There are many Quantitative Trait Loci (QTLs) conferring resistance to SCN that have been identified; however, only two are widely used: rhg1 and Rhg4. Overuse of cultivars containing these QTLs/genes can lead to SCN resistance breakdown, necessitating the use of additional strategies. In this manuscript, a literature review is conducted on research related to soybean resistance to SCN. The main goal is to provide a current understanding of the mechanisms of SCN resistance and list the areas of research that could be further explored.

Manganese-Induced Neurotoxicity through Impairment of Cross-Talk Pathways in Human Neuroblastoma Cell Line SH-SY5Y Differentiated with Retinoic Acid

Manganese (Mn) is an important element; yet acute and/or chronic exposure to this metal has been linked to neurotoxicity and neurodegenerative illnesses such as Parkinson’s disease and others via an unknown mechanism. To better understand it, we exposed a human neuroblastoma cell model (SH-SY5Y) to two Mn chemical species, MnCl₂ and Citrate of Mn(II) (0–2000 µM), followed by a cell viability assay, transcriptomics, and bioinformatics. Even though these cells have been chemically and genetically modified, which may limit the significance of our findings, we discovered that by using RA-differentiated cells instead of undifferentiated SH-SY5Y cell line, both chemical species induce a similar toxicity, potentially governed by disruption of protein metabolism, with some differences. The MnCl₂ altered amino acid metabolism, which affects RNA metabolism and protein synthesis. Citrate of Mn(II), however, inhibited the E3 ubiquitin ligases–target protein degradation pathway, which can lead to the buildup of damaged/unfolded proteins, consistent with histone modification. Finally, we discovered that Mn(II)-induced cytotoxicity in RA-SH-SY5Y cells shared 84 percent of the pathways involved in neurodegenerative diseases.

HIF-1α Regulation of Cytokine Production following TLR3 Engagement in Murine Bone Marrow-Derived Macrophages Is Dependent on Viral Nucleic Acid Length and Glucose Availability

Hypoxia-inducible factor-1α (HIF-1α) is an important regulator of glucose metabolism and inflammatory cytokine production in innate immune responses. Viruses modulate HIF-1α to support viral replication and the survival of infected cells, but it is unclear if this transcription factor also plays an important role in regulating antiviral immune responses. In this study, we found that short and long dsRNA differentially engage TLR3, inducing distinct levels of proinflammatory cytokine production (TNF-α and IL-6) in bone marrow-derived macrophages from C57BL/6 mice. These responses are associated with differential accumulation of HIF-1α, which augments NF-κB activation. Unlike TLR4 responses, increased HIF-1α following TLR3 engagement is not associated with significant alterations in glycolytic activity and was more pronounced in low glucose conditions. We also show that the mechanisms supporting HIF-1α stabilization may differ following stimulation with short versus long dsRNA and that pyruvate kinase M2 and mitochondrial reactive oxygen species play a central role in these processes. Collectively, this work suggests that HIF-1α may fine-tune proinflammatory cytokine production during early antiviral immune responses, particularly when there is limited glucose availability or under other conditions of stress. Our findings also suggest we may be able to regulate the magnitude of proinflammatory cytokine production during antiviral responses by targeting proteins or molecules that contribute to HIF-1α stabilization.

Human-Soybean Allergies: Elucidation of the Seed Proteome and Comprehensive Protein-Protein Interaction Prediction

The soybean crop, Glycine max (L.) Merr., is consumed by humans, Homo sapiens, worldwide. While the respective bodies of literature and -omics data for each of these organisms are extensive, comparatively few studies investigate the molecular biological processes occurring between the two. We are interested in elucidating the network of protein-protein interactions (PPIs) involved in human-soybean allergies. To this end, we leverage state-of-the-art sequence-based PPI predictors amenable to predicting the enormous comprehensive interactome between human and soybean. A network-based analytical approach is proposed, leveraging similar interaction profiles to identify candidate allergens and proteins involved in the allergy response. Interestingly, the predicted interactome can be explored from two complementary perspectives: which soybean proteins are predicted to interact with specific human proteins and which human proteins are predicted to interact with specific soybean proteins. A total of eight proteins (six specific to the human proteome and two to the soy proteome) have been identified and supported by the literature to be involved in human health, specifically related to immunological and neurological pathways. This study, beyond generating the most comprehensive human-soybean interactome to date, elucidated a soybean seed interactome and identified several proteins putatively consequential to human health.

Propionic acid disrupts endocytosis, cell cycle, and cellular respiration in yeast

Objective

We previously identified propionic acid as a microbially-produced volatile organic compound with fungicidal activity against several pathogenic fungi. The purpose of this work is to better understand how propionic acid affects fungi by examining some of the effects of this compound on the yeast cell.

Results

We show that propionic acid causes a dramatic increase in the uptake of lucifer yellow in yeast cells, which is consistent with enhanced endocytosis. Additionally, using a propidium iodide assay, we show that propionic acid treatment causes a significant increase in the proportion of yeast cells in G₁ and a significant decrease in the proportion of cells in G₂, suggesting that propionic acid causes a cell cycle arrest in yeast. Finally, we show that the reduction of MTT is attenuated in yeast cells treated with propionic acid, indicating that propionic acid disrupts cellular respiration. Understanding the effects of propionic acid on the yeast cell may aid in assessing the broader utility of this compound.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13104-021-05752-z.

The conserved Tpk1 regulates non-homologous end joining double-strand break repair by phosphorylation of Nej1, a homolog of the human XLF

The yeast cyclic AMP-dependent protein kinase A (PKA) is a ubiquitous serine-threonine kinase, encompassing three catalytic (Tpk1-3) and one regulatory (Bcy1) subunits. Evidence suggests PKA involvement in DNA damage checkpoint response, but how DNA repair pathways are regulated by PKA subunits remains inconclusive. Here, we report that deleting the tpk1 catalytic subunit reduces non-homologous end joining (NHEJ) efficiency, whereas tpk2-3 and bcy1 deletion does not. Epistatic analyses revealed that tpk1, as well as the DNA damage checkpoint kinase (dun1) and NHEJ factor (nej1), co-function in the same pathway, and parallel to the NHEJ factor yku80. Chromatin immunoprecipitation and resection data suggest that tpk1 deletion influences repair protein recruitments and DNA resection. Further, we show that Tpk1 phosphorylation of Nej1 at S298 (a Dun1 phosphosite) is indispensable for NHEJ repair and nuclear targeting of Nej1 and its binding partner Lif1. In mammalian cells, loss of PRKACB (human homolog of Tpk1) also reduced NHEJ efficiency, and similarly, PRKACB was found to phosphorylate XLF (a Nej1 human homolog) at S263, a corresponding residue of the yeast Nej1 S298. Together, our results uncover a new and conserved mechanism for Tpk1 and PRKACB in phosphorylating Nej1 (or XLF), which is critically required for NHEJ repair.

Identifying new variation at the J locus, previously identified as e6, in long juvenile 'Paranagoiana' soybean

KEY MESSAGE

A previously identified soybean maturity locus, E6, is discovered to be J, with the long juvenile allele in Paranagoiana now deemed j-x. Soybean grown at latitudes of ~20° or lower can produce lower grain yields due to the short days. This limitation can be overcome by using the long juvenile trait (LJ) which delays flowering under short day conditions. Two LJ loci have been mapped to the same location on Gm04, J and E6. The objective of this research was to investigate the e6 allele in 'Paranagoiana' and determine if E6 and J are the same locus or linked loci. KASP markers showed that e6 lines did not have the j-1 allele of LJ PI 159925. A population fixed for E1 but segregating for E6, with e6 introgressed from Paranagoiana, showed single gene control for flowering and maturity under short days. Sequencing Glyma.04G050200, the J gene, with long amplification Taq found that the e6 line 'Paranagoiana' contains a Ty1-copia retrotransposon of ~10,000 bp, inserted within exon 4. PCR amplification of the cDNA of Glyma.04G050200 also showed differences between the mRNA sequences (presence of insertion in j-x). Hence, we conclude that the loci E6 and J are one locus and deem this new variation found in Paranagoiana as j-x.

Manganese-induced neurotoxicity in cerebellar granule neurons due to perturbation of cell network pathways with potential implications for neurodegenerative disorders

Manganese (Mn) is essential for living organisms, playing an important role in nervous system function. Nevertheless, chronic and/or acute exposure to this metal, especially during early life stages, can lead to neurotoxicity and dementia by unclear mechanisms. Thus, based on previous works of our group with yeast and zebrafish, we hypothesized that the mechanisms mediating manganese-induced neurotoxicity can be associated with the alteration of protein metabolism. These mechanisms may also depend on the chemical speciation of manganese. Therefore, the current study aimed at investigating the mechanisms mediating the toxic effects of manganese in primary cultures of cerebellar granule neurons (CGNs). By exposing cultured CGNs to different chemical species of manganese ([[2-[(dithiocarboxy)amino]ethyl]carbamodithioato]](2-)-kS,kS']manganese, named maneb (MB), and [[1,2-ethanediylbis[carbamodithioato]](2-)]manganese mixture with [[1,2-ethanediylbis[carbamodithioato]](2-)]zinc, named mancozeb (MZ), and manganese chloride (MnCl₂)), and using the MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide) assay, we observed that both MB and MZ induced similar cytotoxicity (LC₅₀∼ 7-9 μM), which was higher than that of MnCl₂ (LC₅₀∼ 27 μM). Subsequently, we applied systems biology approaches, including metallomics, proteomics, gene expression and bioinformatics, and revealed that independent of chemical speciation, for non-cytotoxic concentrations (0.3-3 μM), Mn-induced neurotoxicity in CGNs is associated with metal dyshomeostasis and impaired protein metabolism. In this way, we verified that MB induced more post-translational alterations than MnCl₂, which can be a plausible explanation for cytotoxic differences between both chemical species. The metabolism of proteins is one of the most energy consuming cellular processes and its impairment appears to be a key event of some cellular stress processes reported separately in other studies such as cell cycle arrest, energy impairment, cell signaling, excitotoxicity, immune response, potential protein accumulation and apoptosis. Interestingly, we verified that Mn-induced neurotoxicity shares pathways associated with the development of Alzheimer's disease, Amyotrophic Lateral Sclerosis, Huntington's disease, and Parkinson's disease. This has been observed in baker's yeast and zebrafish suggesting that the mode of action of Mn may be evolutionarily conserved.

Lithium Chloride Sensitivity in Yeast and Regulation of Translation

For decades, lithium chloride (LiCl) has been used as a treatment option for those living with bipolar disorder (BD). As a result, many studies have been conducted to examine its mode of action, toxicity, and downstream cellular responses. We know that LiCl is able to affect cell signaling and signaling transduction pathways through protein kinase C and glycogen synthase kinase-3, which are considered to be important in regulating gene expression at the translational level. However, additional downstream effects require further investigation, especially in translation pathway. In yeast, LiCl treatment affects the expression, and thus the activity, of PGM2, a phosphoglucomutase involved in sugar metabolism. Inhibition of PGM2 leads to the accumulation of intermediate metabolites of galactose metabolism causing cell toxicity. However, it is not fully understood how LiCl affects gene expression in this matter. In this study, we identified three genes, NAM7, PUS2, and RPL27B, which increase yeast LiCl sensitivity when deleted. We further demonstrate that NAM7, PUS2, and RPL27B influence translation and exert their activity through the 5′-Untranslated region (5′-UTR) of PGM2 mRNA in yeast.

Sensitivity of yeast to lithium chloride connects the activity of YTA6 and YPR096C to translation of structured mRNAs

Lithium Chloride (LiCl) toxicity, mode of action and cellular responses have been the subject of active investigations over the past decades. In yeast, LiCl treatment is reported to reduce the activity and alters the expression of PGM2, a gene that encodes a phosphoglucomutase involved in sugar metabolism. Reduced activity of phosphoglucomutase in the presence of galactose causes an accumulation of intermediate metabolites of galactose metabolism leading to a number of phenotypes including growth defect. In the current study, we identify two understudied yeast genes, YTA6 and YPR096C that when deleted, cell sensitivity to LiCl is increased when galactose is used as a carbon source. The 5’-UTR of PGM2 mRNA is structured. Using this region, we show that YTA6 and YPR096C influence the translation of PGM2 mRNA.

PIPE4: Fast PPI Predictor for Comprehensive Inter- and Cross-Species Interactomes

The need for larger-scale and increasingly complex protein-protein interaction (PPI) prediction tasks demands that state-of-the-art predictors be highly efficient and adapted to inter- and cross-species predictions. Furthermore, the ability to generate comprehensive interactomes has enabled the appraisal of each PPI in the context of all predictions leading to further improvements in classification performance in the face of extreme class imbalance using the Reciprocal Perspective (RP) framework. We here describe the PIPE4 algorithm. Adaptation of the PIPE3/MP-PIPE sequence preprocessing step led to upwards of 50x speedup and the new Similarity Weighted Score appropriately normalizes for window frequency when applied to any inter- and cross-species prediction schemas. Comprehensive interactomes for three prediction schemas are generated: (1) cross-species predictions, where Arabidopsis thaliana is used as a proxy to predict the comprehensive Glycine max interactome, (2) inter-species predictions between Homo sapiens-HIV1, and (3) a combined schema involving both cross- and inter-species predictions, where both Arabidopsis thaliana and Caenorhabditis elegans are used as proxy species to predict the interactome between Glycine max (the soybean legume) and Heterodera glycines (the soybean cyst nematode). Comparing PIPE4 with the state-of-the-art resulted in improved performance, indicative that it should be the method of choice for complex PPI prediction schemas.

Differential remodeling of the electron transport chain is required to support TLR3 and TLR4 signaling and cytokine production in macrophages

Increasing evidence suggests that mitochondria play a critical role in driving innate immune responses against bacteria and viruses. However, it is unclear if differential reprogramming of mitochondrial function contributes to the fine tuning of pathogen specific immune responses. Here, we found that TLR3 and TLR4 engagement on murine bone marrow derived macrophages was associated with differential remodeling of electron transport chain complex expression. This remodeling was associated with differential accumulation of mitochondrial and cytosolic ROS, which were required to support ligand specific inflammatory and antiviral cytokine production. We also found that the magnitude of TLR3, but not TLR4, responses were modulated by glucose availability. Under conditions of low glucose, TLR3 engagement was associated with increased ETC complex III expression, increased mitochondrial and cytosolic ROS and increased inflammatory and antiviral cytokine production. This amplification was selectively reversed by targeting superoxide production from the outer Q-binding site of the ETC complex III. These results suggest that ligand specific modulation of the ETC may act as a rheostat that fine tunes innate immune responses via mitochondrial ROS production. Modulation of these processes may represent a novel mechanism to modulate the nature as well as the magnitude of antiviral vs. inflammatory immune responses.

Mode of action of nisin on Escherichia coli

Nisin is a class I polycyclic bacteriocin produced by the bacterium Lactococcus lactis, which is used extensively as a food additive to inhibit the growth of foodborne Gram-positive bacteria. Nisin also inhibits growth of Gram-negative bacteria when combined with membrane-disrupting chelators such as citric acid. To gain insight into nisin's mode of action, we analyzed chemical-genetic interactions and identified nisin-sensitive Escherichia coli strains in the Keio library of knockout mutants. The most sensitive mutants fell into two main groups. The first group accords with the previously proposed mode of action based on studies with Gram-positive bacteria, whereby nisin interacts with factors involved in cell wall, membrane, envelope biogenesis. We identified an additional, novel mode of action for nisin based on the second group of sensitive mutants that involves cell cycle and DNA replication, recombination, and repair. Further analyses supported these two distinct modes of action.

In Silico Engineering of Synthetic Binding Proteins from Random Amino Acid Sequences

Synthetic proteins with high affinity and selectivity for a protein target can be used as research tools, biomarkers, and pharmacological agents, but few methods exist to design such proteins de novo. To this end, the In-Silico Protein Synthesizer (InSiPS) was developed to design synthetic binding proteins (SBPs) that bind pre-determined targets while minimizing off-target interactions. InSiPS is a genetic algorithm that refines a pool of random sequences over hundreds of generations of mutation and selection to produce SBPs with pre-specified binding characteristics. As a proof of concept, we design SBPs against three yeast proteins and demonstrate binding and functional inhibition of two of three targets in vivo. Peptide SPOT arrays confirm binding sites, and a permutation array demonstrates target specificity. Our foundational approach will support the field of de novo design of small binding polypeptide motifs and has robust applicability while offering potential advantages over the limited number of techniques currently available.

•

InSiPS engineers synthetic binding proteins (SBPs) using primary protein sequence

•

SBPs are designed to a bind a target protein and avoid “off-target” interactions

•

Binding and functional inhibition of two of three target proteins in yeast is demonstrated

•

Our new approach offers advantages over alternative tools that rely on 3D models

Zinc oxide and silver nanoparticles toxicity in the baker's yeast, Saccharomyces cerevisiae

Engineered nanomaterials (ENMs) are increasingly incorporated into a variety of commercial applications and consumer products; however, ENMs may possess cytotoxic properties due to their small size. This study assessed the effects of two commonly used ENMs, zinc oxide nanoparticles (ZnONPs) and silver nanoparticles (AgNPs), in the model eukaryote Saccharomyces cerevisiae. A collection of ≈4600 S. cerevisiae deletion mutant strains was used to deduce the genes, whose absence makes S. cerevisiae more prone to the cytotoxic effects of ZnONPs or AgNPs. We demonstrate that S. cerevisiae strains that lack genes involved in transmembrane and membrane transport, cellular ion homeostasis, and cell wall organization or biogenesis exhibited the highest sensitivity to ZnONPs. In contrast, strains that lack genes involved in transcription and RNA processing, cellular respiration, and endocytosis and vesicular transport exhibited the highest sensitivity to AgNPs. Secondary assays confirmed that ZnONPs affected cell wall function and integrity, whereas AgNPs exposure decreased transcription, reduced endocytosis, and led to a dysfunctional electron transport system. This study supports the use of S. cerevisiae Gene Deletion Array as an effective high-throughput technique to determine cellular targets of ENM toxicity.

Systematic Genetic Screens Reveal the Dynamic Global Functional Organization of the Bacterial Translation Machinery

Bacterial protein synthesis is an essential, conserved, and environmentally responsive process. Yet, many of its components and dependencies remain unidentified. To address this gap, we used quantitative synthetic genetic arrays to map functional relationships among >48,000 gene pairs in Escherichia coli under four culture conditions differing in temperature and nutrient availability. The resulting data provide global functional insights into the roles and associations of genes, pathways, and processes important for efficient translation, growth, and environmental adaptation. We predict and independently verify the requirement of unannotated genes for normal translation, including a previously unappreciated role of YhbY in 30S biogenesis. Dynamic changes in the patterns of genetic dependencies across the four growth conditions and data projections onto other species reveal overarching functional and evolutionary pressures impacting the translation system and bacterial fitness, underscoring the utility of systematic screens for investigating protein synthesis, adaptation, and evolution.

The sensitivity of the yeast, Saccharomyces cerevisiae, to acetic acid is influenced by DOM34 and RPL36A

The presence of acetic acid during industrial alcohol fermentation reduces the yield of fermentation by imposing additional stress on the yeast cells. The biology of cellular responses to stress has been a subject of vigorous investigations. Although much has been learned, details of some of these responses remain poorly understood. Members of heat shock chaperone HSP proteins have been linked to acetic acid and heat shock stress responses in yeast. Both acetic acid and heat shock have been identified to trigger different cellular responses including reduction of global protein synthesis and induction of programmed cell death. Yeast HSC82 and HSP82 code for two important heat shock proteins that together account for 1–2% of total cellular proteins. Both proteins have been linked to responses to acetic acid and heat shock. In contrast to the overall rate of protein synthesis which is reduced, the expression of HSC82 and HSP82 is induced in response to acetic acid stress. In the current study we identified two yeast genes DOM34 and RPL36A that are linked to acetic acid and heat shock sensitivity. We investigated the influence of these genes on the expression of HSP proteins. Our observations suggest that Dom34 and RPL36A influence translation in a CAP-independent manner.

Evolution of protein-protein interaction networks in yeast

Interest in the evolution of protein-protein and genetic interaction networks has been rising in recent years, but the lack of large-scale high quality comparative datasets has acted as a barrier. Here, we carried out a comparative analysis of computationally predicted protein-protein interaction (PPI) networks from five closely related yeast species. We used the Protein-protein Interaction Prediction Engine (PIPE), which uses a database of known interactions to make sequence-based PPI predictions, to generate high quality predicted interactomes. Simulated proteomes and corresponding PPI networks were used to provide null expectations for the extent and nature of PPI network evolution. We found strong evidence for conservation of PPIs, with lower than expected levels of change in PPIs for about a quarter of the proteome. Furthermore, we found that changes in predicted PPI networks are poorly predicted by sequence divergence. Our analyses identified a number of functional classes experiencing fewer PPI changes than expected, suggestive of purifying selection on PPIs. Our results demonstrate the added benefit of considering predicted PPI networks when studying the evolution of closely related organisms.

Mapping and identification of a potential candidate gene for a novel maturity locus, E10, in soybean

KEY MESSAGE

E10 is a new maturity locus in soybean and FT4 is the predicted/potential functional gene underlying the locus. Flowering and maturity time traits play crucial roles in economic soybean production. Early maturity is critical for north and west expansion of soybean in Canada. To date, 11 genes/loci have been identified which control time to flowering and maturity; however, the molecular bases of almost half of them are not yet clear. We have identified a new maturity locus called "E10" located at the end of chromosome Gm08. The gene symbol E10e10 has been approved by the Soybean Genetics Committee. The e10e10 genotype results in 5-10 days earlier maturity than E10E10. A set of presumed E10E10 and e10e10 genotypes was used to identify contrasting SSR and SNP haplotypes. These haplotypes, and their association with maturity, were maintained through five backcross generations. A functional genomics approach using a predicted protein-protein interaction (PPI) approach (Protein-protein Interaction Prediction Engine, PIPE) was used to investigate approximately 75 genes located in the genomic region that SSR and SNP analyses identified as the location of the E10 locus. The PPI analysis identified FT4 as the most likely candidate gene underlying the E10 locus. Sequence analysis of the two FT4 alleles identified three SNPs, in the 5'UTR, 3'UTR and fourth exon in the coding region, which result in differential mRNA structures. Allele-specific markers were developed for this locus and are available for soybean breeders to efficiently develop earlier maturing cultivars using molecular marker assisted breeding.

Elevated levels of ribosomal proteins eL36 and eL42 control expression of Hsp90 in rhabdomyosarcoma

Mammalian 90 kDa heat shock protein (Hsp90) is a ubiquitous molecular chaperone whose expression is selectively upregulated during stress, although the precise control mechanism of this increase is yet to be fully elucidated. We used polysome profiling to show that Hsp90α mRNA is selectively translated, while global translation is inhibited during heat stress. Furthermore, we have identified 2 ribosomal proteins, eL36 and eL42 that modulate Hsp90α expression under both normal and heat shock conditions. Importantly, we noted that expression of eL36 and eL42 is elevated in a panel of human rhabdomyosarcomas where it drives high expression of Hsp90 and modulates sensitivity of these cells to an Hsp90 inhibitor 17-AAG.

The genome of an Encephalitozoon cuniculi type III strain reveals insights into the genetic diversity and mode of reproduction of a ubiquitous vertebrate pathogen

Encephalitozoon cuniculi is a model microsporidian species with a mononucleate nucleus and a genome that has been extensively studied. To date, analyses of genome diversity have revealed the existence of four genotypes in E. cuniculi (EcI, II, III and IV). Genome sequences are available for EcI, II and III, and are all very divergent, possibly diploid and genetically homogeneous. The mechanisms that cause low genetic diversity in E. cuniculi (for example, selfing, inbreeding or a combination of both), as well as the degree of genetic variation in their natural populations, have been hard to assess because genome data have been so far gathered from laboratory-propagated strains. In this study, we aim to tackle this issue by analyzing the complete genome sequence of a natural strain of E. cuniculi isolated in 2013 from a steppe lemming. The strain belongs to the EcIII genotype and has been designated EcIII-L. The EcIII-L genome sequence harbors genomic features intermediate to known genomes of II and III lab strains, and we provide primers that differentiate the three E. cuniculi genotypes using a single PCR. Surprisingly, the EcIII-L genome is also highly homogeneous, harbors signatures of heterozygosity and also one strain-specific single-nucleotide polymorphism (SNP) that introduces a stop codon in a key meiosis gene, Spo11. Functional analyses using a heterologous system demonstrate that this SNP leads to a deficient meiosis in a model fungus. This indicates that EcIII-L meiotic machinery may be presently broken. Overall, our findings reveal previously unsuspected genome diversity in E. cuniculi, some of which appears to affect genes of primary importance for the biology of this pathogen.

Conditional Epistatic Interaction Maps Reveal Global Functional Rewiring of Genome Integrity Pathways in Escherichia coli

As antibiotic resistance is increasingly becoming a public health concern, an improved understanding of the bacterial DNA damage response (DDR), which is commonly targeted by antibiotics, could be of tremendous therapeutic value. Although the genetic components of the bacterial DDR have been studied extensively in isolation, how the underlying biological pathways interact functionally remains unclear. Here, we address this by performing systematic, unbiased, quantitative synthetic genetic interaction (GI) screens and uncover widespread changes in the GI network of the entire genomic integrity apparatus of Escherichia coli under standard and DNA-damaging growth conditions. The GI patterns of untreated cultures implicated two previously uncharacterized proteins (YhbQ and YqgF) as nucleases, whereas reorganization of the GI network after DNA damage revealed DDR roles for both annotated and uncharacterized genes. Analyses of pan-bacterial conservation patterns suggest that DDR mechanisms and functional relationships are near universal, highlighting a modular and highly adaptive genomic stress response.

Impact of sea-level rise on cross-shore sediment transport on fetch-limited barrier reef island beaches under modal and cyclonic conditions

A one-dimensional wave model is combined with an analytical sediment transport model to investigate the likely influence of sea-level rise on net cross-shore sediment transport on fetch-limited barrier reef and lagoon island beaches. The modelling considers if changes in the nearshore wave height and wave period in the lagoon induced by different water levels over the reef flat are likely to lead to net offshore or onshore movement of sediment. The results indicate that the effects of SLR on net sediment movement are highly variable and controlled by the bathymetry of the reef and lagoon. A significant range of reef-lagoon bathymetry, and notably shallow and narrow reefs, appears to lead hydrodynamic conditions and beaches that are likely to be stable or even accrete under SLR. Loss of reef structural complexity, particularly on the reef flat, increases the chance of sediment transport away from beaches and offshore.

Advances and challenges in mucosal adjuvant technology

Adjuvants play attractive roles in enhancement of immune response during vaccination; however, due to several challenges, only a limited number of adjuvants are licensed by health authorities. The lack of an effective mucosal adjuvant is even more significant as none of the licensed adjuvants revealed a strong enhancement in immune system after mucosal administration. Over the past two decades, several mucosal adjuvants have been developed to deliver antigens to the target cells in the mucosal immune system and increase specific immune responses. However, the safety and efficacy of these adjuvants for testing in human trials is still an important issue, requiring further study. In this article, we briefly review the challenges associated with most common mucosal adjuvants and discuss potential strategies for targeting the mucosal immune system.

Efficient prediction of human protein-protein interactions at a global scale

BACKGROUND

Our knowledge of global protein-protein interaction (PPI) networks in complex organisms such as humans is hindered by technical limitations of current methods.

RESULTS

On the basis of short co-occurring polypeptide regions, we developed a tool called MP-PIPE capable of predicting a global human PPI network within 3 months. With a recall of 23% at a precision of 82.1%, we predicted 172,132 putative PPIs. We demonstrate the usefulness of these predictions through a range of experiments.

CONCLUSIONS

The speed and accuracy associated with MP-PIPE can make this a potential tool to study individual human PPI networks (from genomic sequences alone) for personalized medicine.

Novel function discovery with GeneMANIA: a new integrated resource for gene function prediction in Escherichia coli

MOTIVATION

The model bacterium Escherichia coli is among the best studied prokaryotes, yet nearly half of its proteins are still of unknown biological function. This is despite a wealth of available large-scale physical and genetic interaction data. To address this, we extended the GeneMANIA function prediction web application developed for model eukaryotes to support E.coli.

RESULTS

We integrated 48 distinct E.coli functional interaction datasets and used the GeneMANIA algorithm to produce thousands of novel functional predictions and prioritize genes for further functional assays. Our analysis achieved cross-validation performance comparable to that reported for eukaryotic model organisms, and revealed new functions for previously uncharacterized genes in specific bioprocesses, including components required for cell adhesion, iron-sulphur complex assembly and ribosome biogenesis. The GeneMANIA approach for network-based function prediction provides an innovative new tool for probing mechanisms underlying bacterial bioprocesses.

CONTACT

gary.bader@utoronto.ca; mohan.babu@uregina.ca

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Impact of sea-level rise and coral mortality on the wave dynamics and wave forces on barrier reefs

A one-dimensional wave model was used to investigate the reef top wave dynamics across a large suite of idealized reef-lagoon profiles, representing barrier coral reef systems under different sea-level rise (SLR) scenarios. The modeling shows that the impacts of SLR vary spatially and are strongly influenced by the bathymetry of the reef and coral type. A complex response occurs for the wave orbital velocity and forces on corals, such that the changes in the wave dynamics vary reef by reef. Different wave loading regimes on massive and branching corals also leads to contrasting impacts from SLR. For many reef bathymetries, wave orbital velocities increase with SLR and cyclonic wave forces are reduced for certain coral species. These changes may be beneficial to coral health and colony resilience and imply that predicting SLR impacts on coral reefs requires careful consideration of the reef bathymetry and the mix of coral species.

Quantitative genome-wide genetic interaction screens reveal global epistatic relationships of protein complexes in Escherichia coli

Large-scale proteomic analyses in Escherichia coli have documented the composition and physical relationships of multiprotein complexes, but not their functional organization into biological pathways and processes. Conversely, genetic interaction (GI) screens can provide insights into the biological role(s) of individual gene and higher order associations. Combining the information from both approaches should elucidate how complexes and pathways intersect functionally at a systems level. However, such integrative analysis has been hindered due to the lack of relevant GI data. Here we present a systematic, unbiased, and quantitative synthetic genetic array screen in E. coli describing the genetic dependencies and functional cross-talk among over 600,000 digenic mutant combinations. Combining this epistasis information with putative functional modules derived from previous proteomic data and genomic context-based methods revealed unexpected associations, including new components required for the biogenesis of iron-sulphur and ribosome integrity, and the interplay between molecular chaperones and proteases. We find that functionally-linked genes co-conserved among γ-proteobacteria are far more likely to have correlated GI profiles than genes with divergent patterns of evolution. Overall, examining bacterial GIs in the context of protein complexes provides avenues for a deeper mechanistic understanding of core microbial systems.

Genome-wide genetic interaction (GI) screens have been performed in yeast, but no analogous large-scale studies have yet been reported for bacteria. Here, we have used E. coli synthetic genetic array (eSGA) technology developed by our group to quantitatively map GIs to reveal epistatic dependencies and functional cross-talk among ∼600,000 digenic mutant combinations. By combining this epistasis information with functional modules derived by our group's earlier efforts from proteomic and genomic context (GC)-based methods, we identify several unexpected pathway-level dependencies, functional links between protein complexes, and biological roles of uncharacterized bacterial gene products. As part of the study, two of our pathway predictions from GI screens were validated experimentally, where we confirmed the role of these new components in iron-sulphur biogenesis and ribosome integrity. We also extrapolated the epistatic connectivity diagram of E. coli to 233 distantly related γ-proteobacterial species lacking GI information, and identified co-conserved genes and functional modules important for bacterial pathogenesis. Overall, this study describes the first genome-scale map of GIs in gram-negative bacterium, and through integrative analysis with previously derived protein-protein and GC-based interaction networks presents a number of novel insights into the architecture of bacterial pathways that could not have been discerned through either network alone.