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.

Multiomics analysis of tolerant interaction of potato with potato virus Y.

Potato virus Y (PVY) is the most economically important viral pathogen of potato worldwide. Different potato cultivars react to the pathogen differently, resulting in resistant, tolerant or disease outcome of the interaction. Here we focus on tolerant interaction between potato cv. Désirée and PVYNTN. To capture the response in its full complexity, we analyzed the dynamic changes on multiple molecular levels, including transcriptomics, sRNAomics, degradomics, proteomics and hormonomics. The analysis was complemented by the measurements of viral accumulation, photosynthetic activity and phenotypisation of the symptoms. Besides cv. Désirée we also studied its transgenic counterpart depleted for the accumulation of salicylic acid (NahG-Désirée). This multiomics analysis provides better insights into the mechanisms leading to tolerant response of potato to viral infection and can be used as a base in further studies of plant immunity regulation.

DiNAR: revealing hidden patterns of plant signalling dynamics using Differential Network Analysis in R.

BACKGROUND: Progress in high-throughput molecular methods accompanied by more complex experimental designs demands novel data visualisation solutions. To specifically answer the question which parts of the specifical biological system are responding in particular perturbation, integrative approach in which experimental data are superimposed on a prior knowledge network is shown to be advantageous. RESULTS: We have developed DiNAR, Differential Network Analysis in R, a user-friendly application with dynamic visualisation that integrates multiple condition high-throughput data and extensive biological prior knowledge. Implemented differential network approach and embedded network analysis allow users to analyse condition-specific responses in the context of topology of interest (e.g. immune signalling network) and extract knowledge concerning patterns of signalling dynamics (i.e. rewiring in network structure between two or more biological conditions). We validated the usability of software on the Arabidopsis thaliana and Solanum tuberosum datasets, but it is set to handle any biological instances. CONCLUSIONS: DiNAR facilitates detection of network-rewiring events, gene prioritisation for future experimental design and allows capturing dynamics of complex biological system. The fully cross-platform Shiny App is hosted and freely available at https://nib-si.shinyapps.io/DiNAR. The most recent version of the source code is available at https://github.com/NIB-SI/DiNAR/ with a DOI 10.5281/zenodo.1230523 of the archived version in Zenodo.

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.

Comparison between Proteome and Transcriptome Response in Potato (Solanum tuberosum L.) Leaves Following Potato Virus Y (PVY) Infection.

Plant diseases caused by viral infection are affecting all major crops. Being an obligate intracellular organisms, chemical control of these pathogens is so far not applied in the field except to control the insect vectors of the viruses. Understanding of molecular responses of plant immunity is therefore economically important, guiding the enforcement of crop resistance. To disentangle complex regulatory mechanisms of the plant immune responses, understanding system as a whole is a must. However, integrating data from different molecular analysis (transcriptomics, proteomics, metabolomics, smallRNA regulation etc.) is not straightforward. We evaluated the response of potato (Solanum tuberosum L.) following the infection with potato virus Y (PVY). The response has been analyzed on two molecular levels, with microarray transcriptome analysis and mass spectroscopy-based proteomics. Within this report, we performed detailed analysis of the results on both levels and compared two different approaches for analysis of proteomic data (spectral count versus MaxQuant). To link the data on different molecular levels, each protein was mapped to the corresponding potato transcript according to StNIB paralogue grouping. Only 33% of the proteins mapped to microarray probes in a one-to-one relation and additionally many showed discordance in detected levels of proteins with corresponding transcripts. We discussed functional importance of true biological differences between both levels and showed that the reason for the discordance between transcript and protein abundance lies partly in complexity and structure of biological regulation of proteome and transcriptome and partly in technical issues contributing to it.

Solanum venturii, a suitable model system for virus-induced gene silencing studies in potato reveals StMKK6 as an important player in plant immunity.

BACKGROUND: Virus-induced gene silencing (VIGS) is an optimal tool for functional analysis of genes in plants, as the viral vector spreads throughout the plant and causes reduced expression of selected gene over the whole plant. Potato (Solanum tuberosum) is one of the most important food crops, therefore studies performing functional analysis of its genes are very important. However, the majority of potato cultivars used in laboratory experimental setups are not well amenable to available VIGS systems, thus other model plants from Solanaceae family are used (usually Nicotiana benthamiana). Wild potato relatives can be a better choice for potato model, but their potential in this field was yet not fully explored. This manuscript presents the set-up of VIGS, based on Tobacco rattle virus (TRV) in wild potato relatives for functional studies in potato-virus interactions. RESULTS: Five different potato cultivars, usually used in our lab, did not respond to silencing of phytoene desaturase (PDS) gene with TRV-based vector. Thus screening of a large set of wild potato relatives (different Solanum species and their clones) for their susceptibility to VIGS was performed by silencing PDS gene. We identified several responsive species and further tested susceptibility of these genotypes to potato virus Y (PVY) strain NTN and N. In some species we observed that the presence of empty TRV vector restricted the movement of PVY. Fluorescently tagged PVY(N)-GFP spread systemically in only five of tested wild potato relatives. Based on the results, Solanum venturii (VNT366-2) was selected as the most suitable system for functional analysis of genes involved in potato-PVY interaction. The system was tested by silencing two different plant immune signalling-related kinases, StWIPK and StMKK6. Silencing of StMKK6 enabled faster spreading of the virus throughout the plant, while silencing of WIPK had no effect on spreading of the virus. CONCLUSIONS: The system employing S. venturii (VNT366-2) and PVY(N)-GFP is a suitable method for fast and simple functional analysis of genes involved in potato-PVY interactions. Additionally, a set of identified VIGS responsive species of wild potato relatives could serve as a tool for general studies of potato gene function.

Bimodal dynamics of primary metabolism-related responses in tolerant potato-Potato virus Y interaction.

BACKGROUND: Potato virus Y (PVY) is a major pathogen that causes substantial economic losses in worldwide potato production. Different potato cultivars differ in resistance to PVY, from severe susceptibility, through tolerance, to complete resistance. The aim of this study was to better define the mechanisms underlying tolerant responses of potato to infection by the particularly aggressive PVY(NTN) strain. We focused on the dynamics of the primary metabolism-related processes during PVY(NTN) infection. RESULTS: A comprehensive analysis of the dynamic changes in primary metabolism was performed, which included whole transcriptome analysis, nontargeted proteomics, and photosynthetic activity measurements in potato cv. Désirée and its transgenic counterpart depleted for accumulation of salicylic acid (NahG-Désirée). Faster multiplication of virus occurred in the NahG-Désirée, with these plants developing strong disease symptoms. We show that while the dynamics of responses at the transcriptional level are extensive and bimodal, this is only partially translated to the protein level, and to the final functional outcome. Photosynthesis-related genes are transiently induced before viral multiplication is detected and it is down-regulated later on. This is reflected as a deficiency of the photosynthetic apparatus at the onset of viral multiplication only. Interestingly, specific and constant up-regulation of some RuBisCO transcripts was detected in Désirée plants, which might be important, as these proteins have been shown to interact with viral proteins. In SA-deficient and more sensitive NahG-Désirée plants, consistent down-regulation of photosynthesis-related genes was detected. A constant reduction in the photochemical efficiency from the onset of viral multiplication was identified; in nontransgenic plants this decrease was only transient. The transient reduction in net photosynthetic rate occurred in both genotypes with the same timing, and coincided with changes in stomatal conductivity. CONCLUSIONS: Down-regulation of photosynthesis-related gene expression and decreased photosynthetic activity is in line with other studies that have reported the effects of biotic stress on photosynthesis. Here, we additionally detected induction of light-reaction components in the early stages of PVY(NTN) infection of tolerant interaction. As some of these components have already been shown to interact with viral proteins, their overproduction might contribute to the absence of symptoms in cv. Désirée.

Plant defence model revisions through iterative minimisation of constraint violations.

Biologists have been investigating plant defence response to virus infections; however, a comprehensive mathematical model of this complex process has not been developed. One obstacle in developing a dynamic model, useful for simulation, is the lack of kinetic data from which the model parameters could be determined. We address this problem by proposing a methodology for iterative improvement of the model parameters until the simulation results come close to the expectation of biology experts. These expectations are formalised in the form of constraints to be satisfied by the model simulations. In three iterative steps the model converged to satisfy the biology experts. There are two results of our approach: individual simulations and optimised model parameters, which provide a deeper insight into the biological system. Our constraint-driven optimisation approach allows for an efficient exploration of the dynamic behaviour of biological models and, at the same time, increases their reliability.

Signalling network construction for modelling plant defence response.

Plant defence signalling response against various pathogens, including viruses, is a complex phenomenon. In resistant interaction a plant cell perceives the pathogen signal, transduces it within the cell and performs a reprogramming of the cell metabolism leading to the pathogen replication arrest. This work focuses on signalling pathways crucial for the plant defence response, i.e., the salicylic acid, jasmonic acid and ethylene signal transduction pathways, in the Arabidopsis thaliana model plant. The initial signalling network topology was constructed manually by defining the representation formalism, encoding the information from public databases and literature, and composing a pathway diagram. The manually constructed network structure consists of 175 components and 387 reactions. In order to complement the network topology with possibly missing relations, a new approach to automated information extraction from biological literature was developed. This approach, named Bio3graph, allows for automated extraction of biological relations from the literature, resulting in a set of (component1, reaction, component2) triplets and composing a graph structure which can be visualised, compared to the manually constructed topology and examined by the experts. Using a plant defence response vocabulary of components and reaction types, Bio3graph was applied to a set of 9,586 relevant full text articles, resulting in 137 newly detected reactions between the components. Finally, the manually constructed topology and the new reactions were merged to form a network structure consisting of 175 components and 524 reactions. The resulting pathway diagram of plant defence signalling represents a valuable source for further computational modelling and interpretation of omics data. The developed Bio3graph approach, implemented as an executable language processing and graph visualisation workflow, is publically available at http://ropot.ijs.si/bio3graph/and can be utilised for modelling other biological systems, given that an adequate vocabulary is provided.

Dynamics of responses in compatible potato-Potato virus Y interaction are modulated by salicylic acid.

To investigate the dynamics of the potato-Potato virus Y (PVY) compatible interaction in relation to salicylic acid-controlled pathways we performed experiments using non-transgenic potato cv. Désirée, transgenic NahG-Désirée, cv. Igor and PVY(NTN), the most aggressive strain of PVY. The importance of salicylic acid in viral multiplication and symptom development was confirmed by pronounced symptom development in NahG-Désirée, depleted in salicylic acid, and reversion of the effect after spraying with 2,6-dichloroisonicotinic acid (a salicylic acid-analogue). We have employed quantitative PCR for monitoring virus multiplication, as well as plant responses through expression of selected marker genes of photosynthetic activity, carbohydrate metabolism and the defence response. Viral multiplication was the slowest in inoculated potato of cv. Désirée, the only asymptomatic genotype in the study. The intensity of defence-related gene expression was much stronger in both sensitive genotypes (NahG-Désirée and cv. Igor) at the site of inoculation than in asymptomatic plants (cv. Désirée). Photosynthesis and carbohydrate metabolism gene expression differed between the symptomatic and asymptomatic phenotypes. The differential gene expression pattern of the two sensitive genotypes indicates that the outcome of the interaction does not rely simply on one regulatory component, but similar phenotypical features can result from distinct responses at the molecular level.