Novel stochastic and entropy-based Expectation-Maximisation algorithm for transcription factor binding site motif discovery

The discovery of transcription factor binding site (TFBS) motifs remains an important and challenging problem in computational biology. This thesis presents MITSU, a novel algorithm for TFBS motif discovery which exploits stochastic methods as a means of both overcoming optimality limitations in current algorithms and as a framework for incorporating relevant prior knowledge in order to improve results. The current state of the TFBS motif discovery field is surveyed, with a focus on probabilistic algorithms that typically take the promoter regions of coregulated genes as input. A case is made for an approach based on the stochastic Expectation- Maximisation (sEM) algorithm; its position amongst existing probabilistic algorithms for motif discovery is shown. The algorithm developed in this thesis is unique amongst existing motif discovery algorithms in that it combines the sEM algorithm with a derived data set which leads to an improved approximation to the likelihood function. This likelihood function is unconstrained with regard to the distribution of motif occurrences within the input dataset. MITSU also incorporates a novel heuristic to automatically determine TFBS motif width. This heuristic, known as MCOIN, is shown to outperform current methods for determining motif width. MITSU is implemented in Java and an executable is available for download. MITSU is evaluated quantitatively using realistic synthetic data and several collections of previously characterised prokaryotic TFBS motifs. The evaluation demonstrates that MITSU improves on a deterministic EM-based motif discovery algorithm and an alternative sEM-based algorithm, in terms of previously established metrics. The ability of the sEM algorithm to escape stable fixed points of the EM algorithm, which trap deterministic motif discovery algorithms and the ability of MITSU to discover multiple motif occurrences within a single input sequence are also demonstrated. MITSU is validated using previously characterised Alphaproteobacterial motifs, before being applied to motif discovery in uncharacterised Alphaproteobacterial data. A number of novel results from this analysis are presented and motivate two extensions of MITSU: a strategy for the discovery of multiple different motifs within a single dataset and a higher order Markov background model. The effects of incorporating these extensions within MITSU are evaluated quantitatively using previously characterised prokaryotic TFBS motifs and demonstrated using Alphaproteobacterial motifs. Finally, an information-theoretic measure of motif palindromicity is presented and its advantages over existing approaches for discovering palindromic motifs discussed

Protocol to analyze and validate transcriptomic changes in PDGFRβ-KO mesenchymal stem cell osteogenic potential in the mouse embryo

Mesenchymal stem/stromal cells (MSCs) can differentiate into osteoblasts under appropriate conditions. PDGFRβ signaling controls MSC osteogenic potential both transcriptomically and in culture. Here, we present a “computer to the bench” protocol to analyze changes in MSC osteogenic potential at transcriptomic and cellular level in the absence of PDGFRβ. We detail the preparation of cells from mouse embryos, the analysis of transcriptomic changes from single-cell RNA-sequencing data, the procedure for MSC derivation and culture, and an osteogenic assay for functional validation. For complete details on the use and execution of this protocol, please refer to Sá da Bandeira et al. (2022).(1

Failures of nerve regeneration caused by aging or chronic denervation are rescued by restoring Schwann cell c-Jun.

After nerve injury, myelin and Remak Schwann cells reprogram to repair cells specialized for regeneration. Normally providing strong regenerative support, these cells fail in aging animals, and during chronic denervation that results from slow axon growth. This impairs axonal regeneration and causes significant clinical problems. In mice, we find that repair cells express reduced c-Jun protein as regenerative support provided by these cells declines during aging and chronic denervation. In both cases, genetically restoring Schwann cell c-Jun levels restores regeneration to control levels. We identify potential gene candidates mediating this effect and implicate Shh in the control of Schwann cell c-Jun levels. This establishes that a common mechanism, reduced c-Jun in Schwann cells, regulates success and failure of nerve repair both during aging and chronic denervation. This provides a molecular framework for addressing important clinical problems, suggesting molecular pathways that can be targeted to promote repair in the PNS

Runx1+ vascular smooth muscle cells are essential for hematopoietic stem and progenitor cell development in vivo

Hematopoietic stem cells (HSCs) produce all essential cellular components of the blood. Stromal cell lines supporting HSCs follow a vascular smooth muscle cell (vSMC) differentiation pathway, suggesting that some hematopoiesis-supporting cells originate from vSMC precursors. These pericyte-like precursors were recently identified in the aorta-gonad-mesonephros (AGM) region; however, their role in the hematopoietic development in vivo remains unknown. Here, we identify a subpopulation of NG2 +Runx1 + perivascular cells that display a sclerotome-derived vSMC transcriptomic profile. We show that deleting Runx1 in NG2 + cells impairs the hematopoietic development in vivo and causes transcriptional changes in pericytes/vSMCs, endothelial cells and hematopoietic cells in the murine AGM. Importantly, this deletion leads also to a significant reduction of HSC reconstitution potential in the bone marrow in vivo. This defect is developmental, as NG2 +Runx1 + cells were not detected in the adult bone marrow, demonstrating the existence of a specialised pericyte population in the HSC-generating niche, unique to the embryo. </p

Brain matters: unveiling the distinct contributions of region, age, and sex to glia diversity and CNS function

The myelinated white matter tracts of the central nervous system (CNS) are essential for fast transmission of electrical impulses and are often differentially affected in human neurodegenerative diseases across CNS region, age and sex. We hypothesize that this selective vulnerability is underpinned by physiological variation in white matter glia. Using single nucleus RNA sequencing of human post-mortem white matter samples from the brain, cerebellum and spinal cord and subsequent tissue-based validation we found substantial glial heterogeneity with tissue region: we identified region-specific oligodendrocyte precursor cells (OPCs) that retain developmental origin markers into adulthood, distinguishing them from mouse OPCs. Region-specific OPCs give rise to similar oligodendrocyte populations, however spinal cord oligodendrocytes exhibit markers such as SKAP2 which are associated with increased myelin production and we found a spinal cord selective population particularly equipped for producing long and thick myelin sheaths based on the expression of genes/proteins such as HCN2. Spinal cord microglia exhibit a more activated phenotype compared to brain microglia, suggesting that the spinal cord is a more pro-inflammatory environment, a difference that intensifies with age. Astrocyte gene expression correlates strongly with CNS region, however, astrocytes do not show a more activated state with region or age. Across all glia, sex differences are subtle but the consistent increased expression of protein-folding genes in male donors hints at pathways that may contribute to sex differences in disease susceptibility. These findings are essential to consider for understanding selective CNS pathologies and developing tailored therapeutic strategies

TWEAK/Fn14 signalling promotes cholangiocarcinoma niche formation and progression.

Background & Aims: Cholangiocarcinoma (CCA) is a cancer of the hepatic bile ducts that is rarely resectable and is associated with poor prognosis. Tumour necrosis factor-like weak inducer of apoptosis (TWEAK) is known to signal via its receptor fibroblast growth factor-inducible 14 (Fn14) and induce cholangiocyte and myofibroblast proliferation in liver injury. We aimed to characterise its role in CCA. Methods: The expression of the TWEAK ligand and Fn14 receptor was assessed immunohistochemically and by bulk RNA and single cell transcriptomics of human liver tissue. Spatiotemporal dynamics of pathway regulation were comprehensively analysed in rat and mouse models of thioacetamide (TAA)-mediated CCA. Flow cytometry, qPCR and proteomic analyses of CCA cell lines and conditioned medium experiments with primary macrophages were performed to evaluate the downstream functions of TWEAK/Fn14. In vivo pathway manipulation was assessed via TWEAK overexpression in NICD/AKT-induced CCA or genetic Fn14 knockout during TAA-mediated carcinogenesis. Results: Our data reveal TWEAK and Fn14 overexpression in multiple human CCA cohorts, and Fn14 upregulation in early TAA-induced carcinogenesis. TWEAK regulated the secretion of factors from CC-SW-1 and SNU-1079 CCA cells, inducing polarisation of proinflammatory CD206+ macrophages. Pharmacological blocking of the TWEAK downstream target chemokine monocyte chemoattractant protein 1 (MCP-1 or CCL2) significantly reduced CCA xenograft growth, while TWEAK overexpression drove cancer-associated fibroblast proliferation and collagen deposition in the tumour niche. Genetic Fn14 ablation significantly reduced inflammatory, fibrogenic and ductular responses during carcinogenic TAA-mediated injury. Conclusion: These novel data provide evidence for the action of TWEAK/Fn14 on macrophage recruitment and phenotype, and cancer-associated fibroblast proliferation in CCA. Targeting TWEAK/Fn14 and its downstream signals may provide a means to inhibit CCA niche development and tumour growth. Lay summary: Cholangiocarcinoma is an aggressive, chemotherapy-resistant liver cancer. Interactions between tumour cells and cells that form a supportive environment for the tumour to grow are a source of this aggressiveness and resistance to chemotherapy. Herein, we describe interactions between tumour cells and their supportive environment via a chemical messenger, TWEAK and its receptor Fn14. TWEAK/Fn14 alters the recruitment and type of immune cells in tumours, increases the growth of cancer-associated fibroblasts in the tumour environment, and is a potential target to reduce tumour formation