Latent inter-organ mechanism of idiopathic pulmonary fibrosis unveiled by a generative computational approach.

Idiopathic pulmonary fibrosis (IPF) is a chronic and progressive disease characterized by complex lung pathogenesis affecting approximately three million people worldwide. While the molecular and cellular details of the IPF mechanism is emerging, our current understanding is centered around the lung itself. On the other hand, many human diseases are the products of complex multi-organ interactions. Hence, we postulate that a dysfunctional crosstalk of the lung with other organs plays a causative role in the onset, progression and/or complications of IPF. In this study, we employed a generative computational approach to identify such inter-organ mechanism of IPF. This approach found unexpected molecular relatedness of IPF to neoplasm, diabetes, Alzheimer's disease, obesity, atherosclerosis, and arteriosclerosis. Furthermore, as a potential mechanism underlying this relatedness, we uncovered a putative molecular crosstalk system across the lung and the liver. In this inter-organ system, a secreted protein, kininogen 1, from hepatocytes in the liver interacts with its receptor, bradykinin receptor B1 in the lung. This ligand-receptor interaction across the liver and the lung leads to the activation of calmodulin pathways in the lung, leading to the activation of interleukin 6 and phosphoenolpyruvate carboxykinase 1 pathway across these organs. Importantly, we retrospectively identified several pre-clinical and clinical evidence supporting this inter-organ mechanism of IPF. In conclusion, such feedforward and feedback loop system across the lung and the liver provides a unique opportunity for the development of the treatment and/or diagnosis of IPF. Furthermore, the result illustrates a generative computational framework for machine-mediated synthesis of mechanisms that facilitates and complements the traditional experimental approaches in biomedical sciences.

Latent disease similarities and therapeutic repurposing possibilities uncovered by multi-modal generative topic modeling of human diseases.

Motivation: Human diseases are characterized by multiple features such as their pathophysiological, molecular and genetic changes. The rapid expansion of such multi-modal disease-omics space provides an opportunity to re-classify diverse human diseases and to uncover their latent molecular similarities, which could be exploited to repurpose a therapeutic-target for one disease to another. Results: Herein, we probe this underexplored space by soft-clustering 6955 human diseases by multi-modal generative topic modeling. Focusing on chronic kidney disease and myocardial infarction, two most life-threatening diseases, unveiled are their previously underrecognized molecular similarities to neoplasia and mental/neurological-disorders, and 69 repurposable therapeutic-targets for these diseases. Using an edit-distance-based pathway-classifier, we also find molecular pathways by which these targets could elicit their clinical effects. Importantly, for the 17 targets, the evidence for their therapeutic usefulness is retrospectively found in the pre-clinical and clinical space, illustrating the effectiveness of the method, and suggesting its broader applications across diverse human diseases. Availability and implementation: The code reported in this article is available at: https://github.com/skozawa170301ktx/MultiModalDiseaseModeling. Supplementary information: Supplementary data are available at Bioinformatics Advances online.

Predicting Human Clinical Outcomes Using Mouse Multi-Organ Transcriptome.

Approximately 90% of pre-clinically validated drugs fail in clinical trials owing to unanticipated clinical outcomes, costing over several hundred million US dollars per drug. Despite such critical importance, translating pre-clinical data to clinical outcomes remain a major challenge. Herein, we designed a modality-independent and unbiased approach to predict clinical outcomes of drugs. The approach exploits their multi-organ transcriptome patterns induced in mice and a unique mouse-transcriptome database "humanized" by machine learning algorithms and human clinical outcome datasets. The cross-validation with small-molecule, antibody, and peptide drugs shows effective and efficient identification of the previously known outcomes of 5,519 adverse events and 11,312 therapeutic indications. In addition, the approach is adaptable to deducing potential molecular mechanisms underlying these outcomes. Furthermore, the approach identifies previously unsuspected repositioning targets. These results, together with the fact that it requires no prior structural or mechanistic information of drugs, illustrate its versatile applications to drug development process.

The Body-wide Transcriptome Landscape of Disease Models.

Virtually all diseases affect multiple organs. However, our knowledge of the body-wide effects remains limited. Here, we report the body-wide transcriptome landscape across 13-23 organs of mouse models of myocardial infarction, diabetes, kidney diseases, cancer, and pre-mature aging. Using such datasets, we find (1) differential gene expression in diverse organs across all models; (2) skin as a disease-sensor organ represented by disease-specific activities of putative gene-expression network; (3) a bone-skin cross talk mediated by a bone-derived hormone, FGF23, in response to dysregulated phosphate homeostasis, a known risk-factor for kidney diseases; (4) candidates for the signature activities of many more putative inter-organ cross talk for diseases; and (5) a cross-species map illustrating organ-to-organ and model-to-disease relationships between human and mouse. These findings demonstrate the usefulness and the potential of such body-wide datasets encompassing mouse models of diverse disease types as a resource in biological and medical sciences. Furthermore, the findings described herein could be exploited for designing disease diagnosis and treatment.

Retinal Angiogenesis Regulates Astrocytic Differentiation in Neonatal Mouse Retinas by Oxygen Dependent Mechanisms.

In mice, retinal vascular and astrocyte networks begin to develop at birth, expanding radially from the optic nerve head (ONH) towards the retinal periphery. The retinal vasculature grows towards the periphery ahead of differentiated astrocytes, but behind astrocytic progenitor cells (APCs) and immature astrocytes. Endothelial cell specific Vegfr-2 disruption in newborn mice not only blocked retinal vascular development but also suppressed astrocytic differentiation, reducing the abundance of differentiated astrocytes while causing the accumulation of precursors. By contrast, retinal astrocytic differentiation was accelerated by the exposure of wild-type newborn mice to hyperoxia for 24 hours, or by APC specific deficiency in hypoxia inducible factor (HIF)-2α, an oxygen labile transcription factor. These findings reveal a novel function of the retinal vasculature, and imply that in normal neonatal mice, oxygen from the retinal circulation may promote astrocytic differentiation, in part by triggering oxygen dependent HIF-2α degradation in astrocytic precursors.

Re-evaluating the functional landscape of the cardiovascular system during development.

The cardiovascular system facilitates body-wide distribution of oxygen, a vital process for the development and survival of virtually all vertebrates. However, the zebrafish, a vertebrate model organism, appears to form organs and survive mid-larval periods without a functional cardiovascular system. Despite such dispensability, it is the first organ to develop. Such enigma prompted us to hypothesize other cardiovascular functions that are important for developmental and/or physiological processes. Hence, systematic cellular ablations and functional perturbations were performed on the zebrafish cardiovascular system to gain comprehensive and body-wide understanding of such functions and to elucidate the underlying mechanisms. This approach identifies a set of organ-specific genes, each implicated for important functions. The study also unveils distinct cardiovascular mechanisms, each differentially regulating their expressions in organ-specific and oxygen-independent manners. Such mechanisms are mediated by organ-vessel interactions, circulation-dependent signals, and circulation-independent beating-heart-derived signals. A comprehensive and body-wide functional landscape of the cardiovascular system reported herein may provide clues as to why it is the first organ to develop. Furthermore, these data could serve as a resource for the study of organ development and function.

Real-time prediction of cell division timing in developing zebrafish embryo.

Combination of live-imaging and live-manipulation of developing embryos in vivo provides a useful tool to study developmental processes. Identification and selection of target cells for an in vivo live-manipulation are generally performed by experience- and knowledge-based decision-making of the observer. Computer-assisted live-prediction method would be an additional approach to facilitate the identification and selection of the appropriate target cells. Herein we report such a method using developing zebrafish embryos. We choose V2 neural progenitor cells in developing zebrafish embryo as their successive shape changes can be visualized in real-time in vivo. We developed a relatively simple mathematical method of describing cellular geometry of V2 cells to predict cell division-timing based on their successively changing shapes in vivo. Using quantitatively measured 4D live-imaging data, features of V2 cell-shape at each time point prior to division were extracted and a statistical model capturing the successive changes of the V2 cell-shape was developed. By applying sequential Bayesian inference method to the model, we successfully predicted division-timing of randomly selected individual V2 cells while the cell behavior was being live-imaged. This system could assist pre-selecting target cells desirable for real-time manipulation-thus, presenting a new opportunity for in vivo experimental systems.

Memory of cell shape biases stochastic fate decision-making despite mitotic rounding.

Cell shape influences function, and the current model suggests that such shape effect is transient. However, cells dynamically change their shapes, thus, the critical question is whether shape information remains influential on future cell function even after the original shape is lost. We address this question by integrating experimental and computational approaches. Quantitative live imaging of asymmetric cell-fate decision-making and their live shape manipulation demonstrates that cellular eccentricity of progenitor cell indeed biases stochastic fate decisions of daughter cells despite mitotic rounding. Modelling and simulation indicates that polarized localization of Delta protein instructs by the progenitor eccentricity is an origin of the bias. Simulation with varying parameters predicts that diffusion rate and abundance of Delta molecules quantitatively influence the bias. These predictions are experimentally validated by physical and genetic methods, showing that cells exploit a mechanism reported herein to influence their future fates based on their past shape despite dynamic shape changes.

Whole-mount single molecule FISH method for zebrafish embryo.

Noise in gene expression renders cells more adaptable to changing environment by imposing phenotypic and functional heterogeneity on genetically identical individual cells. Hence, quantitative measurement of noise in gene expression is essential for the study of biological processes in cells. Currently, there are two complementary methods for quantitatively measuring noise in gene expression at the single cell level: single molecule FISH (smFISH) and single cell qRT-PCR (or single cell RNA-seq). While smFISH has been developed for culture cells, tissue sections and whole-mount invertebrate organisms, the method has not been reported for whole-mount vertebrate organisms. Here, we report an smFISH method that is suitable for whole-mount zebrafish embryo, a popular vertebrate model organism for the studies of development, physiology and disease. We show the detection of individual transcripts for several cell-type specific and ubiquitously expressed genes at the single cell level in whole-mount zebrafish embryo. We also demonstrate that the method can be adapted to detect two different genes in individual cells simultaneously. The whole-mount smFISH method described in this report is expected to facilitate the study of noise in gene expression and its role in zebrafish, a vertebrate animal model relevant to human biology.

Neurons limit angiogenesis by titrating VEGF in retina.

Vascular and nervous systems, two major networks in mammalian bodies, show a high degree of anatomical parallelism and functional crosstalk. During development, neurons guide and attract blood vessels, and consequently this parallelism is established. Here, we identified a noncanonical neurovascular interaction in eye development and disease. VEGFR2, a critical endothelial receptor for VEGF, was more abundantly expressed in retinal neurons than in endothelial cells, including endothelial tip cells. Genetic deletion of VEGFR2 in neurons caused misdirected angiogenesis toward neurons, resulting in abnormally increased vascular density around neurons. Further genetic experiments revealed that this misdirected angiogenesis was attributable to an excessive amount of VEGF protein around neurons caused by insufficient engulfment of VEGF by VEGFR2-deficient neurons. Moreover, absence of neuronal VEGFR2 caused misdirected regenerative angiogenesis in ischemic retinopathy. Thus, this study revealed neurovascular crosstalk and unprecedented cellular regulation of VEGF: retinal neurons titrate VEGF to limit neuronal vascularization. PAPERFLICK:

Synthetic cell division system: controlling equal vs. unequal divisions by design.

Cell division is one of the most fundamental and evolutionarily conserved biological processes. Here, we report a synthetic system where we can control by design equal vs. unequal divisions. We synthesized a micro-scale inverse amphipathic droplet of which division is triggered by the increase of surface to volume ratio. Using this system, we succeeded in selectively inducing equal vs. unequal divisions of the droplet cells by adjusting the temperature or the viscosity of the solvent outside the droplet cell accordingly. Our synthetic division system may provide a platform for further development to a system where intracellular contents of the parent droplet cell could be divided into various ratios between the two daughter droplet cells to control their functions and fates.

An anthelmintic drug, pyrvinium pamoate, thwarts fibrosis and ameliorates myocardial contractile dysfunction in a mouse model of myocardial infarction.

Metabolic adaptation to limited supplies of oxygen and nutrients plays a pivotal role in health and disease. Heart attack results from insufficient delivery of oxygen and nutrients to the heart, where cardiomyocytes die and cardiac fibroblasts proliferate--the latter causing scar formation, which impedes regeneration and impairs contractility of the heart. We postulated that cardiac fibroblasts survive metabolic stress by adapting their intracellular metabolism to low oxygen and nutrients, and impeding this metabolic adaptation would thwart their survival and facilitate the repair of scarred heart. Herein, we show that an anthelmintic drug, Pyrvinium pamoate, which has been previously shown to compromise cancer cell survival under glucose starvation condition, also disables cardiac fibroblast survival specifically under glucose deficient condition. Furthermore, Pyrvinium pamoate reduces scar formation and improves cardiac contractility in a mouse model of myocardial infarction. As Pyrvinium pamoate is an FDA-approved drug, our results suggest a therapeutic use of this or other related drugs to repair scarred heart and possibly other organs.

Identification of inter-organ vascular network: vessels bridging between organs.

Development and homeostasis of organs and whole body is critically dependent on the circulatory system. In particular, the circulatory system, the railways shuttling oxygen and nutrients among various organs, is indispensible for inter-organ humoral communication. Since the modern view of the anatomy and mechanics of the circulatory system was established in 17(th) century, it has been assumed that humoral factors are carried to and from organs via vascular branches of the central arteries and veins running along the body axis. Over the past few decades, major advances have been made in understanding molecular and cellular mechanisms underlying the vascularization of organs. However, very little is known about how each organ is linked by vasculature (i.e., inter-organ vascular networks). In fact, the exact anatomy of inter-organ vascular networks has remained obscure. Herein, we report the identification of four distinct vessels, V1(LP), V2(LP), V3(LP) and V4(LP), that bridge between two organs, liver and pancreas in developing zebrafish. We found that these inter-organ vessels can be classified into two types: direct and indirect types. The direct type vessels are those that bridge between two organs via single distinct vessel, to which V1(LP) and V2(LP) vessels belong. The indirect type bridges between two organs via separate branches that emanate from a stem vessel, and V3(LP) and V4(LP) vessels belong to this type. Our finding of V1(LP), V2(LP), V3(LP) and V4(LP) vessels provides the proof of the existence of inter-organ vascular networks. These and other yet-to-be-discovered inter-organ vascular networks may facilitate the direct exchange of humoral factors that are necessary for the coordinated growth, differentiation and homeostasis of the connected organs. It is also possible that the inter-organ vessels serve as tracks for their connected organs to follow during their growth to establish their relative positions and size differences.

Endothelial-derived angiocrine signals induce and sustain regenerative lung alveolarization.

To identify pathways involved in adult lung regeneration, we employ a unilateral pneumonectomy (PNX) model that promotes regenerative alveolarization in the remaining intact lung. We show that PNX stimulates pulmonary capillary endothelial cells (PCECs) to produce angiocrine growth factors that induce proliferation of epithelial progenitor cells supporting alveologenesis. Endothelial cells trigger expansion of cocultured epithelial cells, forming three-dimensional angiospheres reminiscent of alveolar-capillary sacs. After PNX, endothelial-specific inducible genetic ablation of Vegfr2 and Fgfr1 in mice inhibits production of MMP14, impairing alveolarization. MMP14 promotes expansion of epithelial progenitor cells by unmasking cryptic EGF-like ectodomains that activate the EGF receptor (EGFR). Consistent with this, neutralization of MMP14 impairs EGFR-mediated alveolar regeneration, whereas administration of EGF or intravascular transplantation of MMP14(+) PCECs into pneumonectomized Vegfr2/Fgfr1-deficient mice restores alveologenesis and lung inspiratory volume and compliance function. VEGFR2 and FGFR1 activation in PCECs therefore increases MMP14-dependent bioavailability of EGFR ligands to initiate and sustain alveologenesis.

Inductive angiocrine signals from sinusoidal endothelium are required for liver regeneration.

During embryogenesis, endothelial cells induce organogenesis before the development of circulation. These findings suggest that endothelial cells not only form passive conduits to deliver nutrients and oxygen, but also establish an instructive vascular niche, which through elaboration of paracrine trophogens stimulates organ regeneration, in a manner similar to endothelial-cell-derived angiocrine factors that support haematopoiesis. However, the precise mechanism by which tissue-specific subsets of endothelial cells promote organogenesis in adults is unknown. Here we demonstrate that liver sinusoidal endothelial cells (LSECs) constitute a unique population of phenotypically and functionally defined VEGFR3(+)CD34(-)VEGFR2(+)VE-cadherin(+)FactorVIII(+)CD45(-) endothelial cells, which through the release of angiocrine trophogens initiate and sustain liver regeneration induced by 70% partial hepatectomy. After partial hepatectomy, residual liver vasculature remains intact without experiencing hypoxia or structural damage, which allows study of physiological liver regeneration. Using this model, we show that inducible genetic ablation of vascular endothelial growth factor (VEGF)-A receptor-2 (VEGFR2) in the LSECs impairs the initial burst of hepatocyte proliferation (days 1-3 after partial hepatectomy) and subsequent reconstitution of the hepatovascular mass (days 4-8 after partial hepatectomy) by inhibiting upregulation of the endothelial-cell-specific transcription factor Id1. Accordingly, Id1-deficient mice also manifest defects throughout liver regeneration, owing to diminished expression of LSEC-derived angiocrine factors, including hepatocyte growth factor (HGF) and Wnt2. Notably, in in vitro co-cultures, VEGFR2-Id1 activation in LSECs stimulates hepatocyte proliferation. Indeed, intrasplenic transplantation of Id1(+/+) or Id1(-/-) LSECs transduced with Wnt2 and HGF (Id1(-/-)Wnt2(+)HGF(+) LSECs) re-establishes an inductive vascular niche in the liver sinusoids of the Id1(-/-) mice, initiating and restoring hepatovascular regeneration. Therefore, in the early phases of physiological liver regeneration, VEGFR2-Id1-mediated inductive angiogenesis in LSECs through release of angiocrine factors Wnt2 and HGF provokes hepatic proliferation. Subsequently, VEGFR2-Id1-dependent proliferative angiogenesis reconstitutes liver mass. Therapeutic co-transplantation of inductive VEGFR2(+)Id1(+)Wnt2(+)HGF(+) LSECs with hepatocytes provides an effective strategy to achieve durable liver regeneration.

Engraftment and reconstitution of hematopoiesis is dependent on VEGFR2-mediated regeneration of sinusoidal endothelial cells.

Myelosuppression damages the bone marrow (BM) vascular niche, but it is unclear how regeneration of bone marrow vessels contributes to engraftment of transplanted hematopoietic stem and progenitor cells (HSPCs) and restoration of hematopoiesis. We found that chemotherapy and sublethal irradiation induced minor regression of BM sinusoidal endothelial cells (SECs), while lethal irradiation induced severe regression of SECs and required BM transplantation (BMT) for regeneration. Within the BM, VEGFR2 expression specifically demarcated a continuous network of arterioles and SECs, with arterioles uniquely expressing Sca1 and SECs uniquely expressing VEGFR3. Conditional deletion of VEGFR2 in adult mice blocked regeneration of SECs in sublethally irradiated animals and prevented hematopoietic reconstitution. Similarly, inhibition of VEGFR2 signaling in lethally irradiated wild-type mice rescued with BMT severely impaired SEC reconstruction and prevented engraftment and reconstitution of HSPCs. Therefore, regeneration of SECs via VEGFR2 signaling is essential for engraftment of HSPCs and restoration of hematopoiesis.

Secreted Frizzled-related protein 2 is a procollagen C proteinase enhancer with a role in fibrosis associated with myocardial infarction.

Secreted Frizzled-related proteins (sFRPs) have emerged as key regulators of a wide range of developmental and disease processes. Most of the known functions of mammalian sFRPs have been attributed to their ability to antagonize Wnt signalling. Recently however, Xenopus laevis and zebrafish sFRP, Sizzled, was shown to function as an antagonist of Chordin processing by Tolloid-like metalloproteinases. This has led to the proposal that sFRPs may function as evolutionarily conserved antagonists of chordinase activities of this class of proteinases. In contrast to this proposal, we show here that the mammalian sFRP, sFRP2, does not affect Chordin processing, but instead, can serve as a direct enhancer of procollagen C proteinase activity of Tolloid-like metalloproteinases. We also show that the level of fibrosis, in which procollagen processing by Tolloid-like proteinases has a rate-limiting role, is markedly reduced in Sfrp2-null mice subjected to myocardial infarction. Importantly, this reduced level of fibrosis is accompanied by significantly improved cardiac function. This study thus uncovers a function for sFRP2 and a potential therapeutic application for sFRP2 antagonism in controlling fibrosis in the infarcted heart.

On the mechanics of cardiac function of Drosophila embryo.

The heart is a vital organ that provides essential circulation throughout the body. Malfunction of cardiac pumping, thus, leads to serious and most of the times, to fatal diseases. Mechanics of cardiac pumping is a complex process, and many experimental and theoretical approaches have been undertaken to understand this process. We have taken advantage of the simplicity of the embryonic heart of an invertebrate, Drosophila melanogaster, to understand the fundamental mechanics of the beating heart. We applied a live imaging technique to the beating embryonic heart combined with analytical imaging tools to study the dynamic mechanics of the pumping. Furthermore, we have identified one mutant line that exhibits aberrant pumping mechanics. The Drosophila embryonic heart consists of only 104 cardiac cells forming a simple straight tube that can be easily accessed for real-time imaging. Therefore, combined with the wealth of available genetic tools, the embryonic Drosophila heart may serve as a powerful model system for studies of human heart diseases, such as arrhythmia and congenital heart diseases. We, furthermore, believe our mechanistic data provides important information that is useful for our further understanding of the design of biological structure and function and for engineering the pumps for medical uses.

The BTB-kelch protein KLHL6 is involved in B-lymphocyte antigen receptor signaling and germinal center formation.

BTB-kelch proteins can elicit diverse biological functions but very little is known about the physiological role of these proteins in vivo. Kelch-like protein 6 (KLHL6) is a BTB-kelch protein with a lymphoid tissue-restricted expression pattern. In the B-lymphocyte lineage, KLHL6 is expressed throughout ontogeny, and KLHL6 expression is strongly upregulated in germinal center (GC) B cells. To analyze the role of KLHL6 in vivo, we have generated mouse mutants of KLHL6. Development of pro- and pre-B cells was normal but numbers of subsequent stages, transitional 1 and 2, and mature B cells were reduced in KLHL6-deficient mice. The antigen-dependent GC reaction was blunted (smaller GCs, reduced B-cell expansion, and reduced memory antibody response) in the absence of KLHL6. Comparison of mutants with global loss of KLHL6 to mutants lacking KLHL6 specifically in B cells demonstrated a B-cell-intrinsic requirement for KLHL6 expression. Finally, B-cell antigen receptor (BCR) cross-linking was less sensitive in KLHL6-deficient B cells compared to wild-type B cells as measured by proliferation, Ca2+ response, and activation of phospholipase Cgamma2. Our results strongly point to a role for KLHL6 in BCR signal transduction and formation of the full germinal center response.