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.

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.

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.

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.

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.

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.

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.

Robust Recurrent CNV Detection in the Presence of Inter-Subject Variability.

The study of recurrent copy number variations (CNVs) plays an important role in understanding the onset and evolution of complex diseases such as cancer. Array-based comparative genomic hybridization (aCGH) is a widely used microarray based technology for identifying CNVs. However, due to high noise levels and inter-sample variability, detecting recurrent CNVs from aCGH data remains a challenging topic. This paper proposes a novel method for identification of the recurrent CNVs. In the proposed method, the noisy aCGH data is modeled as the superposition of three matrices: a full-rank matrix of weighted piece-wise generating signals accounting for the clean aCGH data, a Gaussian noise matrix to model the inherent experimentation errors and other sources of error, and a sparse matrix to capture the sparse inter-sample (sample-specific) variations. We demonstrated the ability of our method to separate accurately recurrent CNVs from sample-specific variations and noise in both simulated (artificial) data and real data. The proposed method produced more accurate results than current state-of-the-art methods used in recurrent CNV detection and exhibited robustness to noise and sample-specific variations.

In Silico Design and Experimental Validation of Combination Therapy for Pancreatic Cancer.

The number of deaths associated with Pancreatic Cancer has been on the rise in the United States making it an especially dreaded disease. The overall prognosis for pancreatic cancer patients continues to be grim because of the complexity of the disease at the molecular level involving the potential activation/inactivation of several diverse signaling pathways. In this paper, we first model the aberrant signaling in pancreatic cancer using a multi-fault Boolean Network. Thereafter, we theoretically evaluate the efficacy of different drug combinations by simulating this boolean network with drugs at the relevant intervention points and arrive at the most effective drug(s) to achieve cell death. The simulation results indicate that drug combinations containing Cryptotanshinone, a traditional Chinese herb derivative, result in considerably enhanced cell death. These in silico results are validated using wet lab experiments we carried out on Human Pancreatic Cancer (HPAC) cell lines.

A Gaussian Mixture-Model Exploiting Pathway Knowledge for Dissecting Cancer Heterogeneity.

In this work, we develop a systematic approach for applying pathway knowledge to a multivariate Gaussian mixture model for dissecting a heterogeneous cancer tissue. The downstream transcription factors are selected as observables from available partial pathway knowledge in such a way that the subpopulations produce some differential behavior in response to the drugs selected in the upstream. For each subpopulation, each unique (drug, observable) pair is considered as a unique dimension of a multivariate Gaussian distribution. Expectation-maximization (EM) algorithm with hill-climbing is then used to rank the most probable estimates of the mixture composition based on the log-likelihood value. A major contribution of this work is to examine the efficacy of the EM based approach in estimating the composition of experimental mixture sets from cell-by-cell measurements collected on a dynamic cell imaging platform. Towards this end, we apply the algorithm on hourly data collected for two different mixture compositions of A2058, HCT116, and SW480 cell lines for three scenarios: untreated, Lapatinib-treated, and Temsirolimus-treated. Additionally, we show how this methodology can provide a basis for comparing the killing rate of different drugs for a heterogeneous cancer tissue. This obviously has important implications for designing efficient drugs for treating heterogeneous malignant tumors.

Cryptotanshinone Induces Cell Death in Lung Cancer by Targeting Aberrant Feedback Loops.

Signaling pathways oversee highly efficient cellular mechanisms such as growth, division, and death. These processes are controlled by robust negative feedback loops that inhibit receptor-mediated growth factor pathways. Specifically, the ERK, the AKT, and the S6K feedback loops attenuate signaling via growth factor receptors and other kinase receptors to regulate cell growth. Irregularity in any of these supervised processes can lead to uncontrolled cell proliferation and possibly Cancer. These irregularities primarily occur as mutated genes, and an exhaustive search of the perfect drug combination by performing experiments can be both costly and complex. Hence, in this paper, we model the Lung Cancer pathway as a Modified Boolean Network that incorporates feedback. By simulating this network, we theoretically predict the drug combinations that achieve the desired goal for the majority of mutations. Our theoretical analysis identifies Cryptotanshinone, a traditional Chinese herb derivative, as a potent drug component in the fight against cancer. We validated these theoretical results using multiple wet lab experiments carried out on H2073 and SW900 lung cancer cell lines.

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.

Bayesian Inference Identifies Combination Therapeutic Targets in Breast Cancer.

OBJECTIVE: Breast cancer is the second leading cause of cancer death among US women; hence, identifying potential drug targets is an ever increasing need. In this paper, we integrate existing biological information with graphical models to deduce the significant nodes in the breast cancer signaling pathway. METHODS: We make use of biological information from the literature to develop a Bayesian network. Using the relevant gene expression data we estimate the parameters of this network. Then, using a message passing algorithm, we infer the network. The inferred network is used to quantitatively rank different interventions for achieving a desired phenotypic outcome. The particular phenotype considered here is the induction of apoptosis. RESULTS: Theoretical analysis pinpoints to the role of Cryptotanshinone, a compound found in traditional Chinese herbs, as a potent modulator for bringing about cell death in the treatment of cancer. CONCLUSION: Using a mathematical framework, we showed that the combination therapy of mTOR and STAT3 genes yields the best apoptosis in breast cancer. SIGNIFICANCE: The computational results we arrived at are consistent with the experimental results that we obtained using Cryptotanshinone on MCF-7 breast cancer cell lines and also by the past results of others from the literature, thereby demonstrating the effectiveness of our model.

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.

Understanding the Bioinformatics Challenges of Integrating Genomics into Healthcare.

Genomic data is paving the way towards personalized healthcare. By unveiling genetic disease-contributing factors, genomic data can aid in the detection, diagnosis, and treatment of a wide range of complex diseases. Integrating genomic data into healthcare is riddled with a wide range of challenges spanning social, ethical, legal, educational, economic, and technical aspects. Bioinformatics is a core integration aspect presenting an overwhelming number of unaddressed challenges. In this paper we tackle the fundamental bioinformatics integration concerns including: genomic data generation, storage, representation, and utilization in conjunction with clinical data. We divide the bioinformatics challenges into a series of seven intertwined integration aspects spanning the areas of informatics, knowledge management, and communication. For each aspect, we provide a detailed discussion of the current research directions, outstanding challenges, and possible resolutions. This paper seeks to help narrow the gap between the genomic applications, which are being predominantly utilized in research settings, and the clinical adoption of these applications.

Examining De Novo Transcriptome Assemblies via a Quality Assessment Pipeline.

New de novo transcriptome assembly and annotation methods provide an incredible opportunity to study the transcriptome of organisms that lack an assembled and annotated genome. There are currently a number of de novo transcriptome assembly methods, but it has been difficult to evaluate the quality of these assemblies. In order to assess the quality of the transcriptome assemblies, we composed a workflow of multiple quality check measurements that in combination provide a clear evaluation of the assembly performance. We presented novel transcriptome assemblies and functional annotations for Pacific Whiteleg Shrimp (Litopenaeus vannamei ), a mariculture species with great national and international interest, and no solid transcriptome/genome reference. We examined Pacific Whiteleg transcriptome assemblies via multiple metrics, and provide an improved gene annotation. Our investigations show that assessing the quality of an assembly purely based on the assembler's statistical measurements can be misleading; we propose a hybrid approach that consists of statistical quality checks and further biological-based evaluations.

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.