Stress Knowledge Map: A knowledge graph resource for systems biology analysis of plant stress responses.

Stress Knowledge Map (SKM; https://skm.nib.si) is a publicly available resource containing two complementary knowledge graphs that describe the current knowledge of biochemical, signaling, and regulatory molecular interactions in plants: a highly curated model of plant stress signaling (PSS; 543 reactions) and a large comprehensive knowledge network (488 390 interactions). Both were constructed by domain experts through systematic curation of diverse literature and database resources. SKM provides a single entry point for investigations of plant stress response and related growth trade-offs, as well as interactive explorations of current knowledge. PSS is also formulated as a qualitative and quantitative model for systems biology and thus represents a starting point for a plant digital twin. Here, we describe the features of SKM and show, through two case studies, how it can be used for complex analyses, including systematic hypothesis generation and design of validation experiments, or to gain new insights into experimental observations in plant biology.

Oxytetracycline hyper-production through targeted genome reduction of Streptomyces rimosus.

Most biosynthetic gene clusters (BGC) encoding the synthesis of important microbial secondary metabolites, such as antibiotics, are either silent or poorly expressed; therefore, to ensure a strong pipeline of novel antibiotics, there is a need to develop rapid and efficient strain development approaches. This study uses comparative genome analysis to instruct rational strain improvement, using Streptomyces rimosus, the producer of the important antibiotic oxytetracycline (OTC) as a model system. Sequencing of the genomes of two industrial strains M4018 and R6-500, developed independently from a common ancestor, identified large DNA rearrangements located at the chromosome end. We evaluated the effect of these genome deletions on the parental S. rimosus Type Strain (ATCC 10970) genome where introduction of a 145 kb deletion close to the OTC BGC in the Type Strain resulted in massive OTC overproduction, achieving titers that were equivalent to M4018 and R6-500. Transcriptome data supported the hypothesis that the reason for such an increase in OTC biosynthesis was due to enhanced transcription of the OTC BGC and not due to enhanced substrate supply. We also observed changes in the expression of other cryptic BGCs; some metabolites, undetectable in ATCC 10970, were now produced at high titers. This study demonstrated for the first time that the main force behind BGC overexpression is genome rearrangement. This new approach demonstrates great potential to activate cryptic gene clusters of yet unexplored natural products of medical and industrial value.IMPORTANCEThere is a critical need to develop novel antibiotics to combat antimicrobial resistance. Streptomyces species are very rich source of antibiotics, typically encoding 20-60 biosynthetic gene clusters (BGCs). However, under laboratory conditions, most are either silent or poorly expressed so that their products are only detectable at nanogram quantities, which hampers drug development efforts. To address this subject, we used comparative genome analysis of industrial Streptomyces rimosus strains producing high titers of a broad spectrum antibiotic oxytetracycline (OTC), developed during decades of industrial strain improvement. Interestingly, large-scale chromosomal deletions were observed. Based on this information, we carried out targeted genome deletions in the native strain S. rimosus ATCC 10970, and we show that a targeted deletion in the vicinity of the OTC BGC significantly induced expression of the OTC BGC, as well as some other silent BGCs, thus suggesting that this approach may be a useful way to identify new natural products.

Transcriptome-informed identification and characterization of Planococcus citri cis- and trans-isoprenyl diphosphate synthase genes.

Insect physiology and reproduction depend on several terpenoid compounds, whose biosynthesis is mainly unknown. One enigmatic group of insect monoterpenoids are mealybug sex pheromones, presumably resulting from the irregular coupling activity of unidentified isoprenyl diphosphate synthases (IDSs). Here, we performed a comprehensive search for IDS coding sequences of the pest mealybug Planococcus citri. We queried the available genomic and newly generated short- and long-read P. citri transcriptomic data and identified 18 putative IDS genes, whose phylogenetic analysis indicates several gene family expansion events. In vitro testing confirmed regular short-chain coupling activity with five gene products. With the candidate with highest IDS activity, we also detected low amounts of irregular coupling products, and determined amino acid residues important for chain-length preference and irregular coupling activity. This work therefore provides an important foundation for deciphering terpenoid biosynthesis in mealybugs, including the sex pheromone biosynthesis in P. citri.

Cold plasma within a stable supercavitation bubble - A breakthrough technology for efficient inactivation of viruses in water.

Water scarcity, one of the most pressing challenges we face today, has developed for many reasons, including the increasing number of waterborne pollutants that affect the safety of the water environment. Waterborne human, animal and plant viruses represent huge health, environmental, and financial burden and thus it is important to efficiently inactivate them. Therefore, the main objective of this study was to construct a unique device combining plasma with supercavitation and to evaluate its efficiency for water decontamination with the emphasis on inactivation of viruses. High inactivation (>5 log10 PFU/mL) of bacteriophage MS2, a human enteric virus surrogate, was achieved after treatment of 0.43 L of recirculating water for up to 4 min. The key factors in the inactivation were short-lived reactive plasma species that damaged viral RNA. Water treated with plasma for a short time required for successful virus inactivation did not cause cytotoxic effects in the in vitro HepG2 cell model system or adverse effects on potato plant physiology. Therefore, the combined plasma-supercavitation device represents an environmentally-friendly technology that could provide contamination-free and safe water.

Chloroplast redox state changes mark cell-to-cell signaling in the hypersensitive response.

Hypersensitive response (HR)-conferred resistance is associated with induction of programmed cell death and pathogen spread restriction in its proximity. The exact role of chloroplastic reactive oxygen species and its link with salicylic acid (SA) signaling in HR remain unexplained. To unravel this, we performed a detailed spatiotemporal analysis of chloroplast redox response in palisade mesophyll and upper epidermis to potato virus Y (PVY) infection in a resistant potato genotype and its transgenic counterpart with impaired SA accumulation and compromised resistance. Besides the cells close to the cell death zone, we detected individual cells with oxidized chloroplasts further from the cell death zone. These are rare in SA-deficient plants, suggesting their role in signaling for resistance. We confirmed that chloroplast redox changes play important roles in signaling for resistance, as blocking chloroplast redox changes affected spatial responses at the transcriptional level. Through spatiotemporal study of stromule induction after PVY infection, we show that stromules are induced by cell death and also as a response to PVY multiplication at the front of infection. Overall induction of stromules is attenuated in SA-deficient plants.

pISA-tree - a data management framework for life science research projects using a standardised directory tree.

We developed pISA-tree, a straightforward and flexible data management solution for organisation of life science project-associated research data and metadata. pISA-tree was initiated by end-user requirements thus its strong points are practicality and low maintenance cost. It enables on-the-fly creation of enriched directory tree structure (project/Investigation/Study/Assay) based on the ISA model, in a standardised manner via consecutive batch files. Templates-based metadata is generated in parallel at each level enabling guided submission of experiment metadata. pISA-tree is complemented by two R packages, pisar and seekr. pisar facilitates integration of pISA-tree datasets into bioinformatic pipelines and generation of ISA-Tab exports. seekr enables synchronisation with the FAIRDOMHub repository. Applicability of pISA-tree was demonstrated in several national and international multi-partner projects. The system thus supports findable, accessible, interoperable and reusable (FAIR) research and is in accordance with the Open Science initiative. Source code and documentation of pISA-tree are available at https://github.com/NIB-SI/pISA-tree .

Transcriptional deregulation of stress-growth balance in Nicotiana benthamiana biofactories producing insect sex pheromones.

Plant biofactories are a promising platform for sustainable production of high-value compounds, among which are insect sex pheromones, a green alternative to conventional insecticides in agriculture. Recently, we have constructed transgenic Nicotiana benthamiana plants ("Sexy Plants", SxP) that successfully produce a blend of moth (Lepidoptera) sex pheromone compounds (Z)-11-hexadecen-1-ol and (Z)-11-hexadecenyl acetate. However, efficient biosynthesis of sex pheromones resulted in growth and developmental penalty, diminishing the potential for commercial use of SxP in biomanufacturing. To gain insight into the underlying molecular responses, we analysed the whole-genome transcriptome and evaluated it in relation to growth and pheromone production in low- and high-producing transgenic plants of v1.0 and v1.2 SxP lines. In our study, high-producing SxPv1.2 plants accumulated the highest amounts of pheromones but still maintained better growth compared to v1.0 high producers. For an in-depth biological interpretation of the transcriptomic data, we have prepared a comprehensive functional N. benthamiana genome annotation as well as gene translations to Arabidopsis thaliana, enabling functional information transfer by using Arabidopsis knowledge networks. Differential gene expression analysis, contrasting pheromone producers to wild-type plants, revealed that while only a few genes were differentially regulated in low-producing plants, high-producing plants exhibited vast transcriptional reprogramming. They showed signs of stress-like response, manifested as downregulation of photosynthesis-related genes and significant differences in expression of hormonal signalling and secondary metabolism-related genes, the latter presumably leading to previously reported volatilome changes. Further network analyses confirmed stress-like response with activation of jasmonic acid and downregulation of gibberellic acid signalling, illuminating the possibility that the observed growth penalty was not solely a consequence of a higher metabolic burden imposed upon constitutive expression of a heterologous biosynthetic pathway, but rather the result of signalling pathway perturbation. Our work presents an example of comprehensive transcriptomic analyses of disadvantageous stress signalling in N. benthamiana biofactory that could be applied to other bioproduction systems.

Analysis of Virus Spread Around the Cell Death Zone at Spatiotemporal Resolution Using Confocal Microscopy.

The role of programmed cell death (PCD) in hypersensitive response (HR)-conferred resistance depends on the type of host-pathogen interaction and therefore has to be studied for each individual pathosystem. Here we present and explain the protocol for studying the role of PCD in HR-conferred resistance in potato plants in the interaction with the viral pathogen. As an experimental system, we use genotype Rywal, where the virus spread is restricted and HR PCD develops 3 days post potato virus Y (PVY) inoculation. As a control of virus multiplication and spread, we include its transgenic counterpart impaired in salicylic acid (SA) accumulation (NahG-Rywal), in which the HR-PCD occurs but the spread of the virus is not restricted. To follow the occurrence of virus-infected cells and/or virus multiplication outside the cell death zone, we use GFP-tagged PVY (PVY-N605(123)-GFP) which can be monitored by confocal microscopy. Any other plant-pathogen system which results in PCD development could be studied using a modified version of this protocol.

Insect pest management in the age of synthetic biology.

Arthropod crop pests are responsible for 20% of global annual crop losses, a figure predicted to increase in a changing climate where the ranges of numerous species are projected to expand. At the same time, many insect species are beneficial, acting as pollinators and predators of pest species. For thousands of years, humans have used increasingly sophisticated chemical formulations to control insect pests but, as the scale of agriculture expanded to meet the needs of the global population, concerns about the negative impacts of agricultural practices on biodiversity have grown. While biological solutions, such as biological control agents and pheromones, have previously had relatively minor roles in pest management, biotechnology has opened the door to numerous new approaches for controlling insect pests. In this review, we look at how advances in synthetic biology and biotechnology are providing new options for pest control. We discuss emerging technologies for engineering resistant crops and insect populations and examine advances in biomanufacturing that are enabling the production of new products for pest control.

RNA Sequencing Analyses for Deciphering Potato Molecular Responses.

Understanding the molecular mechanisms of potato development and responses to environmental stressors is of utmost importance for achieving stable crop yields. RNA sequencing (RNA-Seq) provides an insight into responses of all of the organism genes to the environmental and developmental cues and thus provides insights into underlying modes of action. In this chapter, we guide a researcher through some of the most important steps in the analysis of transcriptomics data. The initial topic of experimental design is followed by a more wet-lab-oriented section on RNA-Seq sample preparation. Next, we present intermediate steps of data retrieval, quality control, mapping, and differential expression of the dataset and a section on how to expose your data to the public (i.e., public repositories) and make it findable, accessible, interoperable, and reusable (FAIR). In the last four sections, we describe specific tools or Web applications, which ease the exploration of generated results in the context of their gene function and network-based visualizations, specifically GoMapMan, GSEA, DiNAR, and Biomine Explorer. All sections are accompanied by potato dataset examples and include general hints and tricks, as well as potato specificities that one should be aware of.

Potentially Toxic Planktic and Benthic Cyanobacteria in Slovenian Freshwater Bodies: Detection by Quantitative PCR.

Due to increased frequency of cyanobacterial blooms and emerging evidence of cyanotoxicity in biofilm, reliable methods for early cyanotoxin threat detection are of major importance for protection of human, animal and environmental health. To complement the current methods of risk assessment, this study aimed to evaluate selected qPCR assays for detection of potentially toxic cyanobacteria in environmental samples. In the course of one year, 25 plankton and 23 biofilm samples were collected from 15 water bodies in Slovenia. Three different analyses were performed and compared to each other; qPCR targeting mcyE, cyrJ and sxtA genes involved in cyanotoxin production, LC-MS/MS quantifying microcystin, cylindrospermopsin and saxitoxin concentration, and microscopic analyses identifying potentially toxic cyanobacterial taxa. qPCR analyses detected potentially toxic Microcystis in 10 lake plankton samples, and potentially toxic Planktothrix cells in 12 lake plankton and one lake biofilm sample. A positive correlation was observed between numbers of mcyE gene copies and microcystin concentrations. Potential cylindrospermopsin- and saxitoxin-producers were detected in three and seven lake biofilm samples, respectively. The study demonstrated a potential for cyanotoxin production that was left undetected by traditional methods in both plankton and biofilm samples. Thus, the qPCR method could be useful in regular monitoring of water bodies to improve risk assessment and enable timely measures.

Multiple copies of the oxytetracycline gene cluster in selected Streptomyces rimosus strains can provide significantly increased titers.

BACKGROUND: Natural products are a valuable source of biologically active compounds that have applications in medicine and agriculture. One disadvantage with natural products is the slow, time-consuming strain improvement regimes that are necessary to ensure sufficient quantities of target compounds for commercial production. Although great efforts have been invested in strain selection methods, many of these technologies have not been improved in decades, which might pose a serious threat to the economic and industrial viability of such important bioprocesses. RESULTS: In recent years, introduction of extra copies of an entire biosynthetic pathway that encodes a target product in a single microbial host has become a technically feasible approach. However, this often results in minor to moderate increases in target titers. Strain stability and process reproducibility are the other critical factors in the industrial setting. Industrial Streptomyces rimosus strains for production of oxytetracycline are one of the most economically efficient strains ever developed, and thus these represent a very good industrial case. To evaluate the applicability of amplification of an entire gene cluster in a single host strain, we developed and evaluated various gene tools to introduce multiple copies of the entire oxytetracycline gene cluster into three different Streptomyces rimosus strains: wild-type, and medium and high oxytetracycline-producing strains. We evaluated the production levels of these engineered S. rimosus strains with extra copies of the oxytetracycline gene cluster and their stability, and the oxytetracycline gene cluster expression profiles; we also identified the chromosomal integration sites. CONCLUSIONS: This study shows that stable and reproducible increases in target secondary metabolite titers can be achieved in wild-type and in high oxytetracycline-producing strains, which always reflects the metabolic background of each independent S. rimosus strain. Although this approach is technically very demanding and requires systematic effort, when combined with modern strain selection methods, it might constitute a very valuable approach in industrial process development.

Incorporating RNA-Seq transcriptomics into glycosylation-integrating metabolic network modelling kinetics: Multiomic Chinese hamster ovary (CHO) cell bioreactors.

In this work, the kinetic model based on the previously developed metabolic and glycan reaction networks of the ovarian cells of the Chinese hamster ovary (CHO) cell line was improved by the inclusion of transcriptomic data that took into account the values of the RPKM gene (Reads per Kilobase of Exon per Million Reads Mapped). The transcriptomic (RNASeq) data were obtained together with metabolic and glycan data from the literature, and the concentrations with RPKM values were collected at several points in time from two fed-batch processes. First, the fluxes were determined by regression analysis of the metabolic data, then these fluxes were corrected by using the fold change in gene expression as a measure of enzyme concentrations. Next, the corrected fluxes in the kinetic model were used to calculate the concentration profiles of the metabolites, and literature data were used to evaluate the predicted results of the model. Compared to other studies where the concentration profiles of CHO cell metabolites were described using a kinetic model without consideration of RNA-Seq data to correct the fluxes, this model is unique. The additional integration of transcriptomic data led to better predictions of metabolic concentrations in the fed-batch process, which is a significant improvement of the modelling technique used.

Deregulation of whole-transcriptome gene expression in zebrafish (Danio rerio) after chronic exposure to low doses of imatinib mesylate in a complete life cycle study.

Imatinib mesylate (IM) is an anticancer drug that belongs to tyrosine kinase inhibitors. We report the results of the first investigation of the chronic exposure of zebrafish (Danio rerio) to IM. The exposure to IM (0.01, 1 and 100 µg/L) was initiated in adult fish and continued through hatching and the offspring generation for seven months. In addition to standard toxicological endpoints, induction of genotoxic effects and whole-genome transcriptome of liver samples of offspring generation of zebrafish were analysed. Exposure to IM did not affect the survival and growth of zebrafish, did not cause any histopathological changes, but it induced a marginal increase in the chromosomal damage in blood cells. The whole-genome transcriptome analyses demonstrated dose-dependent increase in the number of differentially expressed genes with a significantly higher number of deregulated genes in female fish compared to male. Differentially expressed genes included genes involved in response to DNA damage, cell cycle control and regulation of circadian rhythm. Based on the low genotoxic activity and the pattern of the changes in DNA damage responsive genes we consider that at current environmental exposure levels, IM represents low risk for genotoxic effects in aquatic organisms. Exposure to IM also induced deregulation of the expression of genes associated with steroidogenesis and hormone metabolism and function, which indicates hormone-disrupting activity of IM that has not been studied so far. The study provide new information on the potential consequences of chronic exposure to the residues of tyrosine kinase inhibitors, which remain to be further explored.

Precision transcriptomics of viral foci reveals the spatial regulation of immune-signaling genes and identifies RBOHD as an important player in the incompatible interaction between potato virus Y and potato.

Whereas the activation of resistance (R) proteins has been intensively studied, the downstream signaling mechanisms leading to the restriction of the pathogen remain mostly unknown. We studied the immunity network response conditioned by the potato Ny-1 gene against potato virus Y. We analyzed the processes in the cell death zone and surrounding tissue on the biochemical and gene expression levels in order to reveal the spatiotemporal regulation of the immune response. We show that the transcriptional response in the cell death zone and surrounding tissue is dependent on salicylic acid (SA). For some genes the spatiotemporal regulation is completely lost in the SA-deficient line, whereas other genes show a different response, indicating multiple connections between hormonal signaling modules. The induction of NADPH oxidase RBOHD expression occurs specifically on the lesion border during the resistance response. In plants with silenced RBOHD, the functionality of the resistance response is perturbed and the spread of the virus is not arrested at the site of infection. RBOHD is required for the spatial accumulation of SA, and conversely RBOHD is under the transcriptional regulation of SA. Using spatially resolved RNA-seq, we also identified spatial regulation of an UDP-glucosyltransferase, another component in feedback activation of SA biosynthesis, thus deciphering a novel aspect of resistance signaling.

Roles of small RNAs in the establishment of tolerant interaction between plants and viruses.

In a tolerant plant-virus interaction, viral multiplication is sustained without substantial effects on plant growth or reproduction. Such interactions are, in natural environments, frequent and sometimes even beneficial for both interactors. Here we compiled evidence showing that small RNAs modulate plant immune responses and growth, hence adjusting its physiology to enable a tolerant interaction. Importantly, the role of small RNAs in tolerant interactions resembles that required for establishment of a mutualistic symbiosis. Tolerance can become a sustainable strategy for breeding for virus resistance as selection pressure for emergence of more aggressive strains is low. Understanding the processes underlying establishment of tolerance is, therefore, important for the development of future crops.

Plant Molecular Responses to Potato Virus Y: A Continuum of Outcomes from Sensitivity and Tolerance to Resistance.

Potato virus Y (PVY) is the most economically important virus affecting potato production. PVY manipulates the plant cell machinery in order to successfully complete the infecting cycle. On the other side, the plant activates a sophisticated multilayer immune defense response to combat viral infection. The balance between these mechanisms, depending on the plant genotype and environment, results in a specific outcome that can be resistance, sensitivity, or tolerance. In this review, we summarize and compare the current knowledge on molecular events, leading to different phenotypic outcomes in response to PVY and try to link them with the known molecular mechanisms.

The need for standardisation in life science research - an approach to excellence and trust.

Today, academic researchers benefit from the changes driven by digital technologies and the enormous growth of knowledge and data, on globalisation, enlargement of the scientific community, and the linkage between different scientific communities and the society. To fully benefit from this development, however, information needs to be shared openly and transparently. Digitalisation plays a major role here because it permeates all areas of business, science and society and is one of the key drivers for innovation and international cooperation. To address the resulting opportunities, the EU promotes the development and use of collaborative ways to produce and share knowledge and data as early as possible in the research process, but also to appropriately secure results with the European strategy for Open Science (OS). It is now widely recognised that making research results more accessible to all societal actors contributes to more effective and efficient science; it also serves as a boost for innovation in the public and private sectors. However  for research data to be findable, accessible, interoperable and reusable the use of standards is essential. At the metadata level, considerable efforts in standardisation have already been made (e.g. Data Management Plan and FAIR Principle etc.), whereas in context with the raw data these fundamental efforts are still fragmented and in some cases completely missing. The CHARME consortium, funded by the European Cooperation in Science and Technology (COST) Agency, has identified needs and gaps in the field of standardisation in the life sciences and also discussed potential hurdles for implementation of standards in current practice. Here, the authors suggest four measures in response to current challenges to ensure a high quality of life science research data and their re-usability for research and innovation.

Network Modeling Unravels Mechanisms of Crosstalk between Ethylene and Salicylate Signaling in Potato.

To develop novel crop breeding strategies, it is crucial to understand the mechanisms underlying the interaction between plants and their pathogens. Network modeling represents a powerful tool that can unravel properties of complex biological systems. In this study, we aimed to use network modeling to better understand immune signaling in potato (Solanum tuberosum). For this, we first built on a reliable Arabidopsis (Arabidopsis thaliana) immune signaling model, extending it with the information from diverse publicly available resources. Next, we translated the resulting prior knowledge network (20,012 nodes and 70,091 connections) to potato and superimposed it with an ensemble network inferred from time-resolved transcriptomics data for potato. We used different network modeling approaches to generate specific hypotheses of potato immune signaling mechanisms. An interesting finding was the identification of a string of molecular events illuminating the ethylene pathway modulation of the salicylic acid pathway through Nonexpressor of PR Genes1 gene expression. Functional validations confirmed this modulation, thus supporting the potential of our integrative network modeling approach for unraveling molecular mechanisms in complex systems. In addition, this approach can ultimately result in improved breeding strategies for potato and other sensitive crops.