Transcriptome Analysis of Developing Grains from Wheat Cultivars TAM 111 and TAM 112 Reveal Cultivar-Specific Regulatory Networks.

Wheat flour's end-use quality is tightly linked to the quantity and composition of storage proteins in the endosperm. TAM 111 and TAM 112 are two popular cultivars grown in the Southern US Great Plains with significantly different protein content. To investigate regulatory differences, transcriptome data were analyzed from developing grains at early- and mid-filling stages. At the mid-filling stage, TAM 111 preferentially upregulated starch metabolism-related pathways compared to TAM 112, whereas amino acid metabolism and transporter-related pathways were over-represented in TAM 112. Elemental analyses also indicated a higher N percentage in TAM 112 at the mid-filling stage. To explore the regulatory variation, weighted correlation gene network was constructed from publicly available RNAseq datasets to identify the modules differentially regulated in TAM 111 and TAM 112. Further, the potential transcription factors (TFs) regulating those modules were identified using graphical least absolute shrinkage and selection operator (GLASSO). Homologs of the OsNF-Y family members with known starch metabolism-related functions showed higher connectivities in TAM 111. Multiple TFs with high connectivity in TAM 112 had predicted functions associated with ABA response in grain. These results will provide novel targets for breeders to explore and further our understanding in mechanisms regulating grain development.

Network modeling and inference of peroxisome proliferator-activated receptor pathway in high fat diet-linked obesity.

Systems biology aims to understand how holistic systems theory can be used to explain the observable living system characteristics, and mathematical modeling tools have been successful in understanding the intricate relationships underlying cellular functions. Lately, researchers have been interested in understanding molecular mechanisms underlying obesity, which is a major health concern worldwide and has been linked to several diseases. Various mechanisms such as peroxisome proliferator-activated receptors (PPARs) are known to modulate obesity-induced inflammation and its consequences. In this study, we have modeled the PPAR pathway using a Bayesian model and inferred the sub-pathways that are potentially responsible for the activation of the output processes that are associated with high fat diet (HFD)-induced obesity. We examined a previously published dataset from a study that compared gene expression profiles of 40 mice maintained on HFD against 40 mice fed with chow diet (CD). Our simulations have highlighted that GPCR and FATCD36 sub-pathways were aberrantly active in HFD mice and are therefore favorable targets for anti-obesity strategies. We further cross-validated our observations with experimental results from the literature. We believe that mathematical models such as those presented in the present study can help in inferring other pathways and deducing significant biological relationships.

Bayesian modeling of plant drought resistance pathway.

BACKGROUND: Plants are sessile organisms and are unable to relocate to favorable locations under extreme environmental conditions. Hence they have no choice but to acclimate and eventually adapt to the severe conditions to ensure their survival. As traditional methods of bolstering plant defense against stressful conditions come to their biological limit, we require newer methods that can allow us to strengthen plants' internal defense mechanism. These factors motivated us to look into the genetic networks of plants. The WRKY transcription factors are well known for their role in plant defense against biotic stresses, but recent studies have shed light on their activities against abiotic stresses such as drought. We modeled this network of WRKY transcription factors using Bayesian networks and applied inference algorithm to find the best regulators of drought response. Biologically intervening (activating/inhibiting) these regulators can bolster the defense response of plants against droughts. RESULT: We used real world data from the NCBI GEO database and synthetic data generated from dependencies in the Bayesian network to learn the network parameters. These parameters were estimated using both a Bayesian and a frequentist approach. The two sets of parameters were used in a utility-based inference algorithm to determine the best regulator of plant drought response in the WRKY transcription factor network. CONCLUSION: Our analysis revealed that activating the transcription factor WRKY18 had the highest likelihood of inducing drought response among all the other elements of the WRKY transcription factor network. Our observation was also supported by biological literature, as WRKY18 is known to regulate drought responsive genes positively. We also found that activating the protein complex WRKY60-60 had the second highest likelihood of inducing drought defense response. Consistent with the existing biological literature, we also found the transcription factor WRKY40 and the protein complex WRKY40-40 to suppress drought response.

A Bayesian approach to determine the composition of heterogeneous cancer tissue.

BACKGROUND: Cancer Tissue Heterogeneity is an important consideration in cancer research as it can give insights into the causes and progression of cancer. It is known to play a significant role in cancer cell survival, growth and metastasis. Determining the compositional breakup of a heterogeneous cancer tissue can also help address the therapeutic challenges posed by heterogeneity. This necessitates a low cost, scalable algorithm to address the challenge of accurate estimation of the composition of a heterogeneous cancer tissue. METHODS: In this paper, we propose an algorithm to tackle this problem by utilizing the data of accurate, but high cost, single cell line cell-by-cell observation methods in low cost aggregate observation method for heterogeneous cancer cell mixtures to obtain their composition in a Bayesian framework. RESULTS: The algorithm is analyzed and validated using synthetic data and experimental data. The experimental data is obtained from mixtures of three separate human cancer cell lines, HCT116 (Colorectal carcinoma), A2058 (Melanoma) and SW480 (Colorectal carcinoma). CONCLUSION: The algorithm provides a low cost framework to determine the composition of heterogeneous cancer tissue which is a crucial aspect in cancer research.

Simulating variance heterogeneity in quantitative genome wide association studies.

BACKGROUND: Analyzing Variance heterogeneity in genome wide association studies (vGWAS) is an emerging approach for detecting genetic loci involved in gene-gene and gene-environment interactions. vGWAS analysis detects variability in phenotype values across genotypes, as opposed to typical GWAS analysis, which detects variations in the mean phenotype value. RESULTS: A handful of vGWAS analysis methods have been recently introduced in the literature. However, very little work has been done for evaluating these methods. To enable the development of better vGWAS analysis methods, this work presents the first quantitative vGWAS simulation procedure. To that end, we describe the mathematical framework and algorithm for generating quantitative vGWAS phenotype data from genotype profiles. Our simulation model accounts for both haploid and diploid genotypes under different modes of dominance. Our model is also able to simulate any number of genetic loci causing mean and variance heterogeneity. CONCLUSIONS: We demonstrate the utility of our simulation procedure through generating a variety of genetic loci types to evaluate common GWAS and vGWAS analysis methods. The results of this evaluation highlight the challenges current tools face in detecting GWAS and vGWAS loci.

An in-silico study examining the induction of apoptosis by Cryptotanshinone in metastatic melanoma cell lines.

BACKGROUND: Metastatic melanoma is an aggressive form of skin cancer that evades various anti-cancer treatments including surgery, radio-,immuno- and chemo-therapy. TRAIL-induced apoptosis is a desirable method to treat melanoma since, unlike other treatments, it does not harm non-cancerous cells. The pro-inflammatory response to melanoma by nF κB and STAT3 pathways makes the cancer cells resist TRAIL-induced apoptosis. We show that due to to its dual action on DR5, a death receptor for TRAIL and on STAT3, Cryptotanshinone can be used to increase sensitivity to TRAIL. METHODS: The development of chemoresistance and invasive properties in melanoma cells involves several biological pathways. The key components of these pathways are represented as a Boolean network with multiple inputs and multiple outputs. RESULTS: The possible mutations in genes that can lead to cancer are captured by faults in the combinatorial circuit and the model is used to theoretically predict the effectiveness of Cryptotanshinone for inducing apoptosis in melanoma cell lines. This prediction is experimentally validated by showing that Cryptotanshinone can cause enhanced cell death in A375 melanoma cells. CONCLUSION: The results presented in this paper facilitate a better understanding of melanoma drug resistance. Furthermore, this framework can be used to detect additional drug intervention points in the pathway that could amplify the action of Cryptotanshinone.

Towards targeted combinatorial therapy design for the treatment of castration-resistant prostate cancer.

BACKGROUND: Prostate cancer is one of the most prevalent cancers in males in the United States and amongst the leading causes of cancer related deaths. A particularly virulent form of this disease is castration-resistant prostate cancer (CRPC), where patients no longer respond to medical or surgical castration. CRPC is a complex, multifaceted and heterogeneous malady with limited standard treatment options. RESULTS: The growth and progression of prostate cancer is a complicated process that involves multiple pathways. The signaling network comprising the integral constituents of the signature pathways involved in the development and progression of prostate cancer is modeled as a combinatorial circuit. The failures in the gene regulatory network that lead to cancer are abstracted as faults in the equivalent circuit and the Boolean circuit model is then used to design therapies tailored to counteract the effect of each molecular abnormality and to propose potentially efficacious combinatorial therapy regimens. Furthermore, stochastic computational modeling is utilized to identify potentially vulnerable components in the network that may serve as viable candidates for drug development. CONCLUSION: The results presented herein can aid in the design of scientifically well-grounded targeted therapies that can be employed for the treatment of prostate cancer patients.

HiMMe: using genetic patterns as a proxy for genome assembly reliability assessment.

BACKGROUND: The information content of genomes plays a crucial role in the existence and proper development of living organisms. Thus, tremendous effort has been dedicated to developing DNA sequencing technologies that provide a better understanding of the underlying mechanisms of cellular processes. Advances in the development of sequencing technology have made it possible to sequence genomes in a relatively fast and inexpensive way. However, as with any measurement technology, there is noise involved and this needs to be addressed to reach conclusions based on the resulting data. In addition, there are multiple intermediate steps and degrees of freedom when constructing genome assemblies that lead to ambiguous and inconsistent results among assemblers. METHODS: Here we introduce HiMMe, an HMM-based tool that relies on genetic patterns to score genome assemblies. Through a Markov chain, the model is able to detect characteristic genetic patterns, while, by introducing emission probabilities, the noise involved in the process is taken into account. Prior knowledge can be used by training the model to fit a given organism or sequencing technology. RESULTS: Our results show that the method presented is able to recognize patterns even with relatively small k-mer size choices and limited computational resources. CONCLUSIONS: Our methodology provides an individual quality metric per contig in addition to an overall genome assembly score, with a time complexity well below that of an aligner. Ultimately, HiMMe provides meaningful statistical insights that can be leveraged by researchers to better select contigs and genome assemblies for downstream analysis.

Proliferation symmetry breaking in growing tissues.

Morphogenesis of developing tissues results from anisotropic growth, typically driven by polarized patterns of gene expression. Here we propose an alternative model of anisotropic growth driven by self-organized feed-back between cell polarity, mechanical pressure, and cell division rates. Specifically, cell polarity alignment can induce spontaneous symmetry breaking in proliferation, resulting from the anisotropic distribution of mechanical pressure in the tissue. We show that proliferation anisotropy can be controlled by cellular elasticity, motility and contact inhibition, thereby elucidating the design principles for anisotropic morphogenesis.

Boolean modeling and fault diagnosis in oxidative stress response.

BACKGROUND: Oxidative stress is a consequence of normal and abnormal cellular metabolism and is linked to the development of human diseases. The effective functioning of the pathway responding to oxidative stress protects the cellular DNA against oxidative damage; conversely the failure of the oxidative stress response mechanism can induce aberrant cellular behavior leading to diseases such as neurodegenerative disorders and cancer. Thus, understanding the normal signaling present in oxidative stress response pathways and determining possible signaling alterations leading to disease could provide us with useful pointers for therapeutic purposes. Using knowledge of oxidative stress response pathways from the literature, we developed a Boolean network model whose simulated behavior is consistent with earlier experimental observations from the literature. Concatenating the oxidative stress response pathways with the PI3-Kinase-Akt pathway, the oxidative stress is linked to the phenotype of apoptosis, once again through a Boolean network model. Furthermore, we present an approach for pinpointing possible fault locations by using temporal variations in the oxidative stress input and observing the resulting deviations in the apoptotic signature from the normally predicted pathway. Such an approach could potentially form the basis for designing more effective combination therapies against complex diseases such as cancer. RESULTS: In this paper, we have developed a Boolean network model for the oxidative stress response. This model was developed based on pathway information from the current literature pertaining to oxidative stress. Where applicable, the behaviour predicted by the model is in agreement with experimental observations from the published literature. We have also linked the oxidative stress response to the phenomenon of apoptosis via the PI3k/Akt pathway. CONCLUSIONS: It is our hope that some of the additional predictions here, such as those pertaining to the oscillatory behaviour of certain genes in the presence of oxidative stress, will be experimentally validated in the near future. Of course, it should be pointed out that the theoretical procedure presented here for pinpointing fault locations in a biological network with feedback will need to be further simplified before it can be even considered for practical biological validation.

Cancer therapy design based on pathway logic.

MOTIVATION: Cancer encompasses various diseases associated with loss of cell cycle control, leading to uncontrolled cell proliferation and/or reduced apoptosis. Cancer is usually caused by malfunction(s) in the cellular signaling pathways. Malfunctions occur in different ways and at different locations in a pathway. Consequently, therapy design should first identify the location and type of malfunction to arrive at a suitable drug combination. RESULTS: We consider the growth factor (GF) signaling pathways, widely studied in the context of cancer. Interactions between different pathway components are modeled using Boolean logic gates. All possible single malfunctions in the resulting circuit are enumerated and responses of the different malfunctioning circuits to a 'test' input are used to group the malfunctions into classes. Effects of different drugs, targeting different parts of the Boolean circuit, are taken into account in deciding drug efficacy, thereby mapping each malfunction to an appropriate set of drugs.

Tissue specific expression of UMAMIT amino acid transporters in wheat.

Wheat grain protein content and composition are important for its end-use quality. Protein synthesis during the grain filling phase is supported by the amino acids remobilized from the vegetative tissue, the process in which both amino acid importers and exporters are expected to be involved. Previous studies identified amino acid importers that might function in the amino acid remobilization in wheat. However, the amino acid exporters involved in this process have been unexplored so far. In this study, we have curated the Usually Multiple Amino acids Move In and out Transporter (UMAMIT) family of transporters in wheat. As expected, the majority of UMAMITs were found as triads in the A, B, and D genomes of wheat. Expression analysis using publicly available data sets identified groups of TaUMAMITs expressed in root, leaf, spike, stem and grain tissues, many of which were temporarily regulated. Strong expression of TaUMAMITs was detected in the late senescing leaves and transfer cells in grains, both of which are the expected site of apoplastic amino acid transport during grain filling. Biochemical characterization of selected TaUMAMITs revealed that TaUMAMIT17 shows a strong amino acid export activity and might play a role in amino acid transfer to the grains.

In Silico Modeling of the Induction of Apoptosis by Cryptotanshinone in Osteosarcoma Cell Lines.

Osteosarcoma (OS) is the most common primary malignant bone tumor of both children and pet canines. Its characteristic genomic instability and complexity coupled with the dearth of knowledge about its etiology has made improvement in the current treatment difficult. We use the existing literature about the biological pathways active in OS and combine it with the current research involving natural compounds to identify new targets and design more effective drug therapies. The key components of these pathways are modeled as a Boolean network with multiple inputs and multiple outputs. The combinatorial circuit is employed to theoretically predict the efficacies of various drugs in combination with Cryptotanshinone. We show that the action of the herbal drug, Cryptotanshinone on OS cell lines induces apoptosis by increasing sensitivity to TNF-related apoptosis-inducing ligand (TRAIL) through its multi-pronged action on STAT3, DRP1 and DR5. The Boolean framework is used to detect additional drug intervention points in the pathway that could amplify the action of Cryptotanshinone.

Integrative Network Modeling Highlights the Crucial Roles of Rho-GDI Signaling Pathway in the Progression of non-Small Cell Lung Cancer.

Non-small cell lung cancer (NSCLC) is the most prevalent form of lung cancer and a leading cause of cancer-related deaths worldwide. Using an integrative approach, we analyzed a publicly available merged NSCLC transcriptome dataset using machine learning, protein-protein interaction (PPI) networks and bayesian modeling to pinpoint key cellular factors and pathways likely to be involved with the onset and progression of NSCLC. First, we generated multiple prediction models using various machine learning classifiers to classify NSCLC and healthy cohorts. Our models achieved prediction accuracies ranging from 0.83 to 1.0, with XGBoost emerging as the best performer. Next, using functional enrichment analysis (and gene co-expression network analysis with WGCNA) of the machine learning feature-selected genes, we determined that genes involved in Rho GTPase signaling that modulate actin stability and cytoskeleton were likely to be crucial in NSCLC. We further assembled a PPI network for the feature-selected genes that was partitioned using Markov clustering to detect protein complexes functionally relevant to NSCLC. Finally, we modeled the perturbations in RhoGDI signaling using a bayesian network; our simulations suggest that aberrations in ARHGEF19 and/or RAC2 gene activities contributed to impaired MAPK signaling and disrupted actin and cytoskeleton organization and were arguably key contributors to the onset of tumorigenesis in NSCLC. We hypothesize that targeted measures to restore aberrant ARHGEF19 and/or RAC2 functions could conceivably rescue the cancerous phenotype in NSCLC. Our findings offer promising avenues for early predictive biomarker discovery, targeted therapeutic intervention and improved clinical outcomes in NSCLC.

Anti-tumor effects of cryptotanshinone (C19H20O3) in human osteosarcoma cell lines.

Osteosarcoma is the most prevalent malignant bone tumor and occurs most commonly in the adolescent and young adult population. Despite the recent advances in surgeries and chemotherapy, the overall survival in patients with resectable metastases is around 20%. This challenge in osteosarcoma is often attributed to the drastic differences in the tumorigenic profiles and mutations among patients. With diverse mutations and multiple oncogenes, it is necessary to identify the therapies that can attack various mutations and simultaneously have minor side-effects. In this paper, we constructed the osteosarcoma pathway from literature and modeled it using ordinary differential equations. We then simulated this network for every possible gene mutation and their combinations and ranked different drug combinations based on their efficacy to drive a mutated osteosarcoma network towards cell death. Our theoretical results predict that drug combinations with Cryptotanshinone (C19H20O3), a traditional Chinese herb derivative, have the best overall performance. Specifically, Cryptotanshinone in combination with Temsirolimus inhibit the JAK/STAT, MAPK/ERK, and PI3K/Akt/mTOR pathways and induce cell death in tumor cells. We corroborated our theoretical predictions using wet-lab experiments on SaOS2, 143B, G292, and HU03N1 human osteosarcoma cell lines, thereby demonstrating the potency of Cryptotanshinone in fighting osteosarcoma.

Drug target ranking for glioblastoma multiforme.

BACKGROUND: Glioblastoma Multiforme, an aggressive primary brain tumor, has a poor prognosis and no effective standard of care treatments. Most patients undergoing radiotherapy, along with Temozolomide chemotherapy, develop resistance to the drug, and recurrence of the tumor is a common issue after the treatment. We propose to model the pathways active in Glioblastoma using Boolean network techniques. The network captures the genetic interactions and possible mutations that are involved in the development of the brain tumor. The model is used to predict the theoretical efficacies of drugs for the treatment of cancer. RESULTS: We use the Boolean network to rank the critical intervention points in the pathway to predict an effective therapeutic strategy for Glioblastoma. Drug repurposing helps to identify non-cancer drugs that could be effective in cancer treatment. We predict the effectiveness of drug combinations of anti-cancer and non-cancer drugs for Glioblastoma. CONCLUSIONS: Given the genetic profile of a GBM tumor, the Boolean model can predict the most effective targets for treatment. We also identified two-drug combinations that could be more effective in killing GBM cells than conventional chemotherapeutic agents. The non-cancer drug Aspirin could potentially increase the cytotoxicity of TMZ in GBM patients.

Detecting drought regulators using stochastic inference in Bayesian networks.

Drought is a natural hazard that affects crops by inducing water stress. Water stress, induced by drought accounts for more loss in crop yield than all the other causes combined. With the increasing frequency and intensity of droughts worldwide, it is essential to develop drought-resistant crops to ensure food security. In this paper, we model multiple drought signaling pathways in Arabidopsis using Bayesian networks to identify potential regulators of drought-responsive reporter genes. Genetically intervening at these regulators can help develop drought-resistant crops. We create the Bayesian network model from the biological literature and determine its parameters from publicly available data. We conduct inference on this model using a stochastic simulation technique known as likelihood weighting to determine the best regulators of drought-responsive reporter genes. Our analysis reveals that activating MYC2 or inhibiting ATAF1 are the best single node intervention strategies to regulate the drought-responsive reporter genes. Additionally, we observe simultaneously activating MYC2 and inhibiting ATAF1 is a better strategy. The Bayesian network model indicated that MYC2 and ATAF1 are possible regulators of the drought response. Validation experiments showed that ATAF1 negatively regulated the drought response. Thus intervening at ATAF1 has the potential to create drought-resistant crops.

Comprehensive live-cell imaging analysis of cryptotanshinone and synergistic drug-screening effects in various human and canine cancer cell lines.

BACKGROUND: Several studies have highlighted both the extreme anticancer effects of Cryptotanshinone (CT), a Stat3 crippling component from Salvia miltiorrhiza, as well as other STAT3 inhibitors to fight cancer. METHODS: Data presented in this experiment incorporates 2 years of in vitro studies applying a comprehensive live-cell drug-screening analysis of human and canine cancer cells exposed to CT at 20 µM concentration, as well as to other drug combinations. As previously observed in other studies, dogs are natural cancer models, given to their similarity in cancer genetics, epidemiology and disease progression compared to humans. RESULTS: Results obtained from several types of human and canine cancer cells exposed to CT and varied drug combinations, verified CT efficacy at combating cancer by achieving an extremely high percentage of apoptosis within 24 hours of drug exposure. CONCLUSIONS: CT anticancer efficacy in various human and canine cancer cell lines denotes its ability to interact across different biological processes and cancer regulatory cell networks, driving inhibition of cancer cell survival.

Targeting oncogenic mutations in colorectal cancer using cryptotanshinone.

Colorectal cancer (CRC) is one of the most prevalent types of cancer in the world and ranks second in cancer deaths in the US. Despite the recent improvements in screening and treatment, the number of deaths associated with CRC is still very significant. The complexities involved in CRC therapy stem from multiple oncogenic mutations and crosstalk between abnormal pathways. This calls for using advanced molecular genetics to understand the underlying pathway interactions responsible for this cancer. In this paper, we construct the CRC pathway from the literature and using an existing public dataset on healthy vs tumor colon cells, we identify the genes and pathways that are mutated and are possibly responsible for the disease progression. We then introduce drugs in the CRC pathway, and using a boolean modeling technique, we deduce the drug combinations that produce maximum cell death. Our theoretical simulations demonstrate the effectiveness of Cryptotanshinone, a traditional Chinese herb derivative, achieved by targeting critical oncogenic mutations and enhancing cell death. Finally, we validate our theoretical results using wet lab experiments on HT29 and HCT116 human colorectal carcinoma cell lines.

Adaptive intervention in probabilistic boolean networks.

MOTIVATION: A basic problem of translational systems biology is to utilize gene regulatory networks as a vehicle to design therapeutic intervention strategies to beneficially alter network and, therefore, cellular dynamics. One strain of research has this problem from the perspective of control theory via the design of optimal Markov chain decision processes, mainly in the framework of probabilistic Boolean networks (PBNs). Full optimization assumes that the network is accurately modeled and, to the extent that model inference is inaccurate, which can be expected for gene regulatory networks owing to the combination of model complexity and a paucity of time-course data, the designed intervention strategy may perform poorly. We desire intervention strategies that do not assume accurate full-model inference. RESULTS: This article demonstrates the feasibility of applying on-line adaptive control to improve intervention performance in genetic regulatory networks modeled by PBNs. It shows via simulations that when the network is modeled by a member of a known family of PBNs, an adaptive design can yield improved performance in terms of the average cost. Two algorithms are presented, one better suited for instantaneously random PBNs and the other better suited for context-sensitive PBNs with low switching probability between the constituent BNs.