Joseph Cursons

Predicting expression: the complementary power of histone modification and transcription factor binding data

BackgroundTranscription factors (TFs) and histone modifications (HMs) play critical roles in gene expression by regulating mRNA transcription. Modelling frameworks have been developed to integrate high-throughput omics data, with the aim of elucidating the regulatory logic that results from the interactions of DNA, TFs and HMs. These models have yielded an unexpected and poorly understood result: that TFs and HMs are statistically redundant in explaining mRNA transcript abundance at a genome-wide level.ResultsWe constructed predictive models of gene expression by integrating RNA-sequencing, TF and HM chromatin immunoprecipitation sequencing and DNase I hypersensitivity data for two mammalian cell types. All models identified genome-wide statistical redundancy both within and between TFs and HMs, as previously reported. To investigate potential explanations, groups of genes were constructed for ontology-classified biological processes. Predictive models were constructed for each process to explore the distribution of statistical redundancy. We found significant variation in the predictive capacity of TFs and HMs across these processes and demonstrated the predictive power of HMs to be inversely proportional to process enrichment for housekeeping genes.ConclusionsIt is well established that the roles played by TFs and HMs are not functionally redundant. Instead, we attribute the statistical redundancy reported in this and previous genome-wide modelling studies to the heterogeneous distribution of HMs across chromatin domains. Furthermore, we conclude that statistical redundancy between individual TFs can be readily explained by nucleosome-mediated cooperative binding. This could possibly help the cell confer regulatory robustness by rejecting signalling noise and allowing control via multiple pathways.

Toward Community Standards and Software for Whole-Cell Modeling

Objective: Whole-cell (WC) modeling is a promising tool for biological research, bioengineering, and medicine. However, substantial work remains to create accurate comprehensive models of complex cells. Methods: We organized the 2015 Whole-Cell Modeling Summer School to teach WC modeling and evaluate the need for new WC modeling standards and software by recoding a recently published WC model in the Systems Biology Markup Language. Results: Our analysis revealed several challenges to representing WC models using the current standards. Conclusion: We, therefore, propose several new WC modeling standards, software, and databases. Significance: We anticipate that these new standards and software will enable more comprehensive models.

Hierarchical bond graph modelling of biochemical networks

The bond graph approach to modelling biochemical networks is extended to allow hierarchical construction of complex models from simpler components. This is made possible by representing the simpler components as thermodynamically open systems exchanging mass and energy via ports. A key feature of this approach is that the resultant models are robustly thermodynamically compliant: the thermodynamic compliance is not dependent on precise numerical values of parameters. Moreover, the models are reusable owing to the well-defined interface provided by the energy ports. To extract bond graph model parameters from parameters found in the literature, general and compact formulae are developed to relate free-energy constants and equilibrium constants. The existence and uniqueness of solutions is considered in terms of fundamental properties of stoichiometric matrices. The approach is illustrated by building a hierarchical bond graph model of glycogenolysis in skeletal muscle.

NAIL, a software toolset for inferring, analyzing and visualizing regulatory networks

Summary: The wide variety of published approaches for the problem of regulatory network inference makes using multiple inference algorithms complex and time-consuming. Network Analysis and Inference Library (NAIL) is a set of software tools to simplify the range of computational activities involved in regulatory network inference. It uses a modular approach to connect different network inference algorithms to the same visualization and network-based analyses. NAIL is technologyindependent and includes an interface layer to allow easy integration of components into other applications. Availability and implementation: NAIL is implemented in MATLAB, runs on Windows, Linux and OSX, and is available from SourceForge at https://sourceforge.net/projects/nailsystemsbiology/ for all researchers

The complexities and caveats of lineage tracing in the mammary gland

Lineage tracing is increasingly being utilised to probe different cell types that exist within the mammary gland. Whilst this technique is powerful for tracking cells in vivo and dissecting the roles of different cellular subsets in development, homeostasis and oncogenesis, there are important caveats associated with lineage tracing strategies. Here we highlight key parameters of particular relevance for the mammary gland. These include tissue preparation for whole-mount imaging, whereby the inclusion of enzymatic digestion can drastically alter tissue architecture and cell morphology, and therefore should be avoided. Other factors include the scoring of clones in three dimensions versus two dimensions, the timing of induction, and the marked variability in labelling efficiency that is evident not only between different mouse models harbouring a similar gene promoter but also within a given strain and even within a single mammary gland. Thus, it becomes crucial to visualise extensive areas of ductal tissue and to consider the intricacies of the methodology for lineage tracing studies on normal mammary development and on potential ‘cells of origin’ of cancer.

Regulation of ERK-MAPK signaling in human epidermis

BackgroundThe skin is largely comprised of keratinocytes within the interfollicular epidermis. Over approximately two weeks these cells differentiate and traverse the thickness of the skin. The stage of differentiation is therefore reflected in the positions of cells within the tissue, providing a convenient axis along which to study the signaling events that occur in situ during keratinocyte terminal differentiation, over this extended two-week timescale. The canonical ERK-MAPK signaling cascade (Raf-1, MEK-1/2 and ERK-1/2) has been implicated in controlling diverse cellular behaviors, including proliferation and differentiation. While the molecular interactions involved in signal transduction through this cascade have been well characterized in cell culture experiments, our understanding of how this sequence of events unfolds to determine cell fate within a homeostatic tissue environment has not been fully characterized.MethodsWe measured the abundance of total and phosphorylated ERK-MAPK signaling proteins within interfollicular keratinocytes in transverse cross-sections of human epidermis using immunofluorescence microscopy. To investigate these data we developed a mathematical model of the signaling cascade using a normalized-Hill differential equation formalism.ResultsThese data show coordinated variation in the abundance of phosphorylated ERK-MAPK components across the epidermis. Statistical analysis of these data shows that associations between phosphorylated ERK-MAPK components which correspond to canonical molecular interactions are dependent upon spatial position within the epidermis. The model demonstrates that the spatial profile of activation for ERK-MAPK signaling components across the epidermis may be maintained in a cell-autonomous fashion by an underlying spatial gradient in calcium signaling.ConclusionsOur data demonstrate an extended phospho-protein profile of ERK-MAPK signaling cascade components across the epidermis in situ, and statistical associations in these data indicate canonical ERK-MAPK interactions underlie this spatial profile of ERK-MAPK activation. Using mathematical modelling we have demonstrated that spatially varying calcium signaling components across the epidermis may be sufficient to maintain the spatial profile of ERK-MAPK signaling cascade components in a cell-autonomous manner. These findings may have significant implications for the wide range of cancer drugs which therapeutically target ERK-MAPK signaling components.

A Transcriptional Program for Detecting TGFbeta-induced EMT in Cancer

Most cancer deaths are due to metastasis, and epithelial-to-mesenchymal transition (EMT) plays a central role in driving cancer cell metastasis. EMT is induced by different stimuli, leading to different signaling patterns and therapeutic responses. TGFβ is one of the best-studied drivers of EMT, and many drugs are available to target this signaling pathway. A comprehensive bioinformatics approach was employed to derive a signature for TGFβ-induced EMT which can be used to score TGFβ-driven EMT in cells and clinical specimens. Considering this signature in pan-cancer cell and tumor datasets, a number of cell lines (including basal B breast cancer and cancers of the central nervous system) show evidence for TGFβ-driven EMT and carry a low mutational burden across the TGFβ signaling pathway. Furthermore, significant variation is observed in the response of high scoring cell lines to some common cancer drugs. Finally, this signature was applied to pan-cancer data from The Cancer Genome Atlas to identify tumor types with evidence of TGFβ-induced EMT. Tumor types with high scores showed significantly lower survival rates than those with low scores and also carry a lower mutational burden in the TGFβ pathway. The current transcriptomic signature demonstrates reproducible results across independent cell line and cancer datasets and identifies samples with strong mesenchymal phenotypes likely to be driven by TGFβ. Implications: The TGFβ-induced EMT signature may be useful to identify patients with mesenchymal-like tumors who could benefit from targeted therapeutics to inhibit promesenchymal TGFβ signaling and disrupt the metastatic cascade. Mol Cancer Res; 15(5); 619–31. ©2017 AACR.

Inference of an in situ epidermal intracellular signaling cascade

The stratified architecture of the epidermis makes it an ideal system in which to investigate intracellular signaling pathways within the context of a native tissue. We have applied quantitative imaging protocols to investigate the expression of 13 total-proteins and 4 phosphorylated-proteins across human epidermis. Plasma membrane, nuclear and/or cytoplasmic protein expression levels were measured along the gradient of keratinocyte differentiation. Dynamic Bayesian network techniques were used to infer conditional dependencies between the expression levels of target molecules and construct an associated cascade topology. The resulting networks were compared against a canonical network to investigate the extent to which known biochemical interactions could be recapitulated in situ. Biochemical evidence from the literature supported the majority (71–86%) of inferred network edges, however, overall coverage of the canonical network was relatively low (12–31%). Identified edges may represent key signaling pathway interactions which occur during keratinocyte differentiation. Inferred networks were ranked by model likelihood given the data and the top five were used to construct a consensus network. Several edges were present within this consensus network yet absent from the canonical network. These edges may represent putative interactions which occur in human epidermis.

Systems analysis identifies miR-29b regulation of invasiveness in melanoma.

BackgroundIn many cancers, microRNAs (miRs) contribute to metastatic progression by modulating phenotypic reprogramming processes such as epithelial-mesenchymal plasticity. This can be driven by miRs targeting multiple mRNA transcripts, inducing regulated changes across large sets of genes. The miR-target databases TargetScan and DIANA-microT predict putative relationships by examining sequence complementarity between miRs and mRNAs. However, it remains a challenge to identify which miR-mRNA interactions are active at endogenous expression levels, and of biological consequence.MethodsWe developed a workflow to integrate TargetScan and DIANA-microT predictions into the analysis of data-driven associations calculated from transcript abundance (RNASeq) data, specifically the mutual information and Pearson’s correlation metrics. We use this workflow to identify putative relationships of miR-mediated mRNA repression with strong support from both lines of evidence. Applying this approach systematically to a large, published collection of unique melanoma cell lines – the Ludwig Melbourne melanoma (LM-MEL) cell line panel – we identified putative miR-mRNA interactions that may contribute to invasiveness. This guided the selection of interactions of interest for further in vitro validation studies.ResultsSeveral miR-mRNA regulatory relationships supported by TargetScan and DIANA-microT demonstrated differential activity across cell lines of varying matrigel invasiveness. Strong negative statistical associations for these putative regulatory relationships were consistent with target mRNA inhibition by the miR, and suggest that differential activity of such miR-mRNA relationships contribute to differences in melanoma invasiveness. Many of these relationships were reflected across the skin cutaneous melanoma TCGA dataset, indicating that these observations also show graded activity across clinical samples. Several of these miRs are implicated in cancer progression (miR-211, -340, -125b, −221, and -29b). The specific role for miR-29b-3p in melanoma has not been well studied. We experimentally validated the predicted miR-29b-3p regulation of LAMC1 and PPIC and LASP1, and show that dysregulation of miR-29b-3p or these mRNA targets can influence cellular invasiveness in vitro.ConclusionsThis analytic strategy provides a comprehensive, systems-level approach to identify miR-mRNA regulation in high-throughput cancer data, identifies novel putative interactions with functional phenotypic relevance, and can be used to direct experimental resources for subsequent experimental validation.Computational scripts are available: http://github.com/uomsystemsbiology/LMMEL-miR-miner

Determining the Significance of Protein Network Features and Attributes Using Permutation Testing

Network analysis methods are increasing in popularity. An approach commonly applied to analyze proteomics data involves the use of protein-protein interaction (PPI) networks to explore the systems-level cooperation between proteins identified in a study. In this context, protein interaction networks can be used alongside the statistical analysis of proteomics data and traditional functional enrichment or pathway enrichment analyses. In network analysis it is possible to adjust for some of the complexities that arise due to the known, explicit interdependence between the measured quantities, in particular, differences in the number of interactions between proteins. Here we describe a method for calculating robust empirical p-values for protein interaction networks. We also provide a worked example with python code demonstrating the implementation of this methodology.