Multi-tissue profiling of oxylipins reveal a conserved up-regulation of epoxide:diol ratio that associates with white adipose tissue inflammation and liver steatosis in obesity

BACKGROUND

Obesity drives maladaptive changes in the white adipose tissue (WAT) which can progressively cause insulin resistance, type 2 diabetes mellitus (T2DM) and metabolic dysfunction-associated liver disease (MASLD). Obesity-mediated loss of WAT homeostasis can trigger liver steatosis through dysregulated lipid pathways such as those related to polyunsaturated fatty acid (PUFA)-derived oxylipins. However, the exact relationship between oxylipins and metabolic syndrome remains elusive and cross-tissue dynamics of oxylipins are ill-defined.

METHODS

We quantified PUFA-related oxylipin species in the omental WAT, liver biopsies and plasma of 88 patients undergoing bariatric surgery (female N = 79) and 9 patients (female N = 4) undergoing upper gastrointestinal surgery, using UPLC-MS/MS. We integrated oxylipin abundance with WAT phenotypes (adipogenesis, adipocyte hypertrophy, macrophage infiltration, type I and VI collagen remodelling) and the severity of MASLD (steatosis, inflammation, fibrosis) quantified in each biopsy. The integrative analysis was subjected to (i) adjustment for known risk factors and, (ii) control for potential drug-effects through UPLC-MS/MS analysis of metformin-treated fat explants ex vivo.

FINDINGS

We reveal a generalized down-regulation of cytochrome P450 (CYP)-derived diols during obesity conserved between the WAT and plasma. Notably, epoxide:diol ratio, indicative of soluble epoxide hydrolyse (sEH) activity, increases with WAT inflammation/fibrosis, hepatic steatosis and T2DM. Increased 12,13-EpOME:DiHOME in WAT and liver is a marker of worsening metabolic syndrome in patients with obesity.

INTERPRETATION

These findings suggest a dampened sEH activity and a possible role of fatty acid diols during metabolic syndrome in major metabolic organs such as WAT and liver. They also have implications in view of the clinical trials based on sEH inhibition for metabolic syndrome.

FUNDING

Wellcome Trust (PS3431_WMIH); Duke-NUS (Intramural Goh Cardiovascular Research Award (Duke-NUS-GCR/2022/0020); National Medical Research Council (OFLCG22may-0011); National Institute of Environmental Health Sciences (Z01 ES025034); NIHR Imperial Biomedical Research Centre.

AI is a viable alternative to high throughput screening: a 318-target study

High throughput screening (HTS) is routinely used to identify bioactive small molecules. This requires physical compounds, which limits coverage of accessible chemical space. Computational approaches combined with vast on-demand chemical libraries can access far greater chemical space, provided that the predictive accuracy is sufficient to identify useful molecules. Through the largest and most diverse virtual HTS campaign reported to date, comprising 318 individual projects, we demonstrate that our AtomNet® convolutional neural network successfully finds novel hits across every major therapeutic area and protein class. We address historical limitations of computational screening by demonstrating success for target proteins without known binders, high-quality X-ray crystal structures, or manual cherry-picking of compounds. We show that the molecules selected by the AtomNet® model are novel drug-like scaffolds rather than minor modifications to known bioactive compounds. Our empirical results suggest that computational methods can substantially replace HTS as the first step of small-molecule drug discovery.

Comparative Analysis of Single-Cell RNA Sequencing Methods with and without Sample Multiplexing

Single-cell RNA sequencing (scRNA-seq) has emerged as a powerful technique for investigating biological heterogeneity at the single-cell level in human systems and model organisms. Recent advances in scRNA-seq have enabled the pooling of cells from multiple samples into single libraries, thereby increasing sample throughput while reducing technical batch effects, library preparation time, and the overall cost. However, a comparative analysis of scRNA-seq methods with and without sample multiplexing is lacking. In this study, we benchmarked methods from two representative platforms: Parse Biosciences (Parse; with sample multiplexing) and 10x Genomics (10x; without sample multiplexing). By using peripheral blood mononuclear cells (PBMCs) obtained from two healthy individuals, we demonstrate that demultiplexed scRNA-seq data obtained from Parse showed similar cell type frequencies compared to 10x data where samples were not multiplexed. Despite relatively lower cell capture affecting library preparation, Parse can detect rare cell types (e.g., plasmablasts and dendritic cells) which is likely due to its relatively higher sensitivity in gene detection. Moreover, a comparative analysis of transcript quantification between the two platforms revealed platform-specific distributions of gene length and GC content. These results offer guidance for researchers in designing high-throughput scRNA-seq studies.

WWP2 Regulates Renal Fibrosis and the Metabolic Reprogramming of Profibrotic Myofibroblasts

Key Points

Background

Renal fibrosis is a common pathologic end point in CKD that is challenging to reverse, and myofibroblasts are responsible for the accumulation of a fibrillar collagen–rich extracellular matrix. Recent studies have unveiled myofibroblasts' diversity in proliferative and fibrotic characteristics, which are linked to different metabolic states. We previously demonstrated the regulation of extracellular matrix genes and tissue fibrosis by WWP2, a multifunctional E3 ubiquitin–protein ligase. Here, we investigate WWP2 in renal fibrosis and in the metabolic reprograming of myofibroblasts in CKD.

Methods

We used kidney samples from patients with CKD and WWP2-null kidney disease mice models and leveraged single-cell RNA sequencing analysis to detail the cell-specific regulation of WWP2 in fibrotic kidneys. Experiments in primary cultured myofibroblasts by bulk-RNA sequencing, chromatin immunoprecipitation sequencing, metabolomics, and cellular metabolism assays were used to study the metabolic regulation of WWP2 and its downstream signaling.

Results

The tubulointerstitial expression of WWP2 was associated with fibrotic progression in patients with CKD and in murine kidney disease models. WWP2 deficiency promoted myofibroblast proliferation and halted profibrotic activation, reducing the severity of renal fibrosis in vivo. In renal myofibroblasts, WWP2 deficiency increased fatty acid oxidation and activated the pentose phosphate pathway, boosting mitochondrial respiration at the expense of glycolysis. WWP2 suppressed the transcription of peroxisome proliferator-activated receptor gamma coactivator 1-alpha (PGC-1α), a metabolic mediator of fibrotic response, and pharmacologic inhibition of PGC-1α partially abrogated the protective effects of WWP2 deficiency on myofibroblasts.

Conclusions

WWP2 regulates the metabolic reprogramming of profibrotic myofibroblasts by a WWP2-PGC-1α axis, and WWP2 deficiency protects against renal fibrosis in CKD.

The cholesterol biosynthesis enzyme FAXDC2 couples Wnt/β-catenin to RTK/MAPK signaling

Wnts, cholesterol, and MAPK signaling are essential for development and adult homeostasis. Here, we report that fatty acid hydroxylase domain containing 2 (FAXDC2), a previously uncharacterized enzyme, functions as a methyl sterol oxidase catalyzing C4 demethylation in the Kandutsch-Russell branch of the cholesterol biosynthesis pathway. FAXDC2, a paralog of MSMO1, regulated the abundance of the specific C4-methyl sterols lophenol and dihydro-T-MAS. Highlighting its clinical relevance, FAXDC2 was repressed in Wnt/β-catenin-high cancer xenografts, in a mouse genetic model of Wnt activation, and in human colorectal cancers. Moreover, in primary human colorectal cancers, the sterol lophenol, regulated by FAXDC2, accumulated in the cancerous tissues and not in adjacent normal tissues. FAXDC2 linked Wnts to RTK/MAPK signaling. Wnt inhibition drove increased recycling of RTKs and activation of the MAPK pathway, and this required FAXDC2. Blocking Wnt signaling in Wnt-high cancers caused both differentiation and senescence; and this was prevented by knockout of FAXDC2. Our data show the integration of 3 ancient pathways, Wnts, cholesterol synthesis, and RTK/MAPK signaling, in cellular proliferation and differentiation.

Systems level identification of a matrisome-associated macrophage polarisation state in multi-organ fibrosis

Tissue fibrosis affects multiple organs and involves a master-regulatory role of macrophages which respond to an initial inflammatory insult common in all forms of fibrosis. The recently unravelled multi-organ heterogeneity of macrophages in healthy and fibrotic human disease suggests that macrophages expressing osteopontin (SPP1) associate with lung and liver fibrosis. However, the conservation of this SPP1+ macrophage population across different tissues and its specificity to fibrotic diseases with different etiologies remain unclear. Integrating 15 single-cell RNA-sequencing datasets to profile 235,930 tissue macrophages from healthy and fibrotic heart, lung, liver, kidney, skin, and endometrium, we extended the association of SPP1+ macrophages with fibrosis to all these tissues. We also identified a subpopulation expressing matrisome-associated genes (e.g., matrix metalloproteinases and their tissue inhibitors), functionally enriched for ECM remodelling and cell metabolism, representative of a matrisome-associated macrophage (MAM) polarisation state within SPP1+ macrophages. Importantly, the MAM polarisation state follows a differentiation trajectory from SPP1+ macrophages and is associated with a core set of regulon activity. SPP1+ macrophages without the MAM polarisation state (SPP1+MAM-) show a positive association with ageing lung in mice and humans. These results suggest an advanced and conserved polarisation state of SPP1+ macrophages in fibrotic tissues resulting from prolonged inflammatory cues within each tissue microenvironment.

Pluripotent stem cell-derived committed cardiac progenitors remuscularize damaged ischemic hearts and improve their function in pigs

Ischemic heart disease, which is often associated with irreversibly damaged heart muscle, is a major global health burden. Here, we report the potential of stem cell-derived committed cardiac progenitors (CCPs) have in regenerative cardiology. Human pluripotent embryonic stem cells were differentiated to CCPs on a laminin 521 + 221 matrix, characterized with bulk and single-cell RNA sequencing, and transplanted into infarcted pig hearts. CCPs differentiated for eleven days expressed a set of genes showing higher expression than cells differentiated for seven days. Functional heart studies revealed significant improvement in left ventricular ejection fraction at four and twelve weeks following transplantation. We also observed significant improvements in ventricular wall thickness and a reduction in infarction size after CCP transplantation (p-value < 0.05). Immunohistology analyses revealed in vivo maturation of the CCPs into cardiomyocytes (CM). We observed temporary episodes of ventricular tachyarrhythmia (VT) in four pigs and persistent VT in one pig, but the remaining five pigs exhibited normal sinus rhythm. Importantly, all pigs survived without the formation of any tumors or VT-related abnormalities. We conclude that pluripotent stem cell-derived CCPs constitute a promising possibility for myocardial infarction treatment and that they may positively impact regenerative cardiology.

Photoreceptor laminin drives differentiation of human pluripotent stem cells to photoreceptor progenitors that partially restore retina function

Blindness caused by advanced stages of inherited retinal diseases and age-related macular degeneration are characterized by photoreceptor loss. Cell therapy involving replacement with functional photoreceptor-like cells generated from human pluripotent stem cells holds great promise. Here, we generated a human recombinant retina-specific laminin isoform, LN523, and demonstrated the role in promoting the differentiation of human embryonic stem cells into photoreceptor progenitors. This chemically defined and xenogen-free method enables reproducible production of photoreceptor progenitors within 32 days. We observed that the transplantation into rd10 mice were able to protect the host photoreceptor outer nuclear layer (ONL) up to 2 weeks post transplantation as measured by full-field electroretinogram. At 4 weeks post transplantation, the engrafted cells were found to survive, mature, and associate with the host's rod bipolar cells. Visual behavioral assessment using the water maze swimming test demonstrated visual improvement in the cell-transplanted rodents. At 20 weeks post transplantation, the maturing engrafted cells were able to replace the loss of host ONL by extensive association with host bipolar cells and synapses. Post-transplanted rabbit model also provided congruent evidence for synaptic connectivity with the degenerated host retina. The results may pave the way for the development of stem cell-based therapeutics for retina degeneration.

Targeting endogenous kidney regeneration using anti-IL11 therapy in acute and chronic models of kidney disease

The kidney has large regenerative capacity, but this is compromised when kidney damage is excessive and renal tubular epithelial cells (TECs) undergo SNAI1-driven growth arrest. Here we investigate the role of IL11 in TECs, kidney injury and renal repair. IL11 stimulation of TECs induces ERK- and p90RSK-mediated GSK3β inactivation, SNAI1 upregulation and pro-inflammatory gene expression. Mice with acute kidney injury upregulate IL11 in TECs leading to SNAI1 expression and kidney dysfunction, which is not seen in Il11 deleted mice or in mice administered a neutralizing IL11 antibody in either preemptive or treatment modes. In acute kidney injury, anti-TGFβ reduces renal fibrosis but exacerbates inflammation and tubule damage whereas anti-IL11 reduces all pathologies. Mice with TEC-specific deletion of Il11ra1 have reduced pathogenic signaling and are protected from renal injury-induced inflammation, fibrosis, and failure. In a model of chronic kidney disease, anti-IL11 therapy promotes TEC proliferation and parenchymal regeneration, reverses fibroinflammation and restores renal mass and function. These data highlight IL11-induced mesenchymal transition of injured TECs as an important renal pathology and suggest IL11 as a therapeutic target for restoring stalled endogenous regeneration in the diseased kidney.

The E3 ubiquitin ligase WWP2 regulates pro-fibrogenic monocyte infiltration and activity in heart fibrosis

Non-ischemic cardiomyopathy (NICM) can cause left ventricular dysfunction through interstitial fibrosis, which corresponds to the failure of cardiac tissue remodeling. Recent evidence implicates monocytes/macrophages in the etiopathology of cardiac fibrosis, but giving their heterogeneity and the antagonizing roles of macrophage subtypes in fibrosis, targeting these cells has been challenging. Here we focus on WWP2, an E3 ubiquitin ligase that acts as a positive genetic regulator of human and murine cardiac fibrosis, and show that myeloid specific deletion of WWP2 reduces cardiac fibrosis in hypertension-induced NICM. By using single cell RNA sequencing analysis of immune cells in the same model, we establish the functional heterogeneity of macrophages and define an early pro-fibrogenic phase of NICM that is driven by Ccl5-expressing Ly6c^high monocytes. Among cardiac macrophage subtypes, WWP2 dysfunction primarily affects Ly6c^high monocytes via modulating Ccl5, and consequentially macrophage infiltration and activation, which contributes to reduced myofibroblast trans-differentiation. WWP2 interacts with transcription factor IRF7, promoting its non-degradative mono-ubiquitination, nuclear translocation and transcriptional activity, leading to upregulation of Ccl5 at transcriptional level. We identify a pro-fibrogenic macrophage subtype in non-ischemic cardiomyopathy, and demonstrate that WWP2 is a key regulator of IRF7-mediated Ccl5/Ly6c^high monocyte axis in heart fibrosis.

PPMS: A framework to Profile Primary MicroRNAs from Single-cell RNA-sequencing datasets

MOTIVATION

Single-cell/nuclei RNA-sequencing (scRNA-seq) technologies can simultaneously quantify gene expression in thousands of cells across the genome. However, the majority of the noncoding RNAs, such as microRNAs (miRNAs), cannot currently be profiled at the same scale. MiRNAs are a class of small noncoding RNAs and play an important role in gene regulation. MiRNAs originate from the processing of primary transcripts, known as primary-microRNAs (pri-miRNAs). The pri-miRNA transcripts, independent of their cognate miRNAs, can also function as long noncoding RNAs, code for micropeptides or even interact with DNA, acting like enhancers. Therefore, it is apparent that the significance of scRNA-seq pri-miRNA profiling expands beyond using pri-miRNA as proxies of mature miRNAs. However, there are no computational methods that allow profiling and quantification of pri-miRNAs at the single-cell-type resolution.

RESULTS

We have developed a simple yet effective computational framework to profile pri-MiRNAs from single-cell RNA-sequencing datasets (PPMS). Based on user input, PPMS can profile pri-miRNAs at cell-type resolution. PPMS can be applied to both newly produced and publicly available datasets obtained via single cell or single-nuclei RNA-seq. It allows users to (i) investigate the distribution of pri-miRNAs across cell types and cell states and (ii) establish a relationship between the number of cells/reads sequenced and the detection of pri-miRNAs. Here, to demonstrate its efficacy, we have applied PPMS to publicly available scRNA-seq data generated from (i) individual chambers (ventricles and atria) of the human heart, (ii) human pluripotent stem cells during their differentiation into cardiomyocytes (the heart beating cells) and (iii) hiPSCs-derived cardiomyocytes infected with severe acute respiratory syndrome coronavirus 2.

A novel network pharmacology approach for leukaemia differentiation therapy using Mogrify®

Acute myeloid leukaemia (AML) is a rapidly fatal blood cancer that is characterised by the accumulation of immature myeloid cells in the blood and bone marrow as a result of blocked differentiation. Methods which identify master transcriptional regulators of AML subtype-specific leukaemia cell states and their combinations could be critical for discovering novel differentiation-inducing therapies. In this proof-of-concept study, we demonstrate a novel utility of the Mogrify® algorithm in identifying combinations of transcription factors (TFs) and drugs, which recapitulate granulocytic differentiation of the NB4 acute promyelocytic leukaemia (APL) cell line, using two different approaches. In the first approach, Connectivity Map (CMAP) analysis of these TFs and their target networks outperformed standard approaches, retrieving ATRA as the top hit. We identify dimaprit and mebendazole as a drug combination which induces myeloid differentiation. In the second approach, we show that genetic manipulation of specific Mogrify®-identified TFs (MYC and IRF1) leads to co-operative induction of APL differentiation, as does pharmacological targeting of these TFs using currently available compounds. We also show that loss of IRF1 blunts ATRA-mediated differentiation, and that MYC represses IRF1 expression through recruitment of PML-RARα, the driver fusion oncoprotein in APL, to the IRF1 promoter. Finally, we demonstrate that these drug combinations can also induce differentiation of primary patient-derived APL cells, and highlight the potential of targeting MYC and IRF1 in high-risk APL. Thus, these results suggest that Mogrify® could be used for drug discovery or repositioning in leukaemia differentiation therapy for other subtypes of leukaemia or cancers.

Immunolipidomics Reveals a Globoside Network During the Resolution of Pro-Inflammatory Response in Human Macrophages

Toll-like receptor 4 (TLR4)-mediated changes in macrophages reshape intracellular lipid pools to coordinate an effective innate immune response. Although this has been previously well-studied in different model systems, it remains incompletely understood in primary human macrophages. Here we report time-dependent lipidomic and transcriptomic responses to lipopolysaccharide (LPS) in primary human macrophages from healthy donors. We grouped the variation of ~200 individual lipid species measured by LC-MS/MS into eight temporal clusters. Among all other lipids, glycosphingolipids (glycoSP) and cholesteryl esters (CE) showed a sharp increase during the resolution phase (between 8h or 16h post LPS). GlycoSP, belonging to the globoside family (Gb3 and Gb4), showed the greatest inter-individual variability among all lipids quantified. Integrative network analysis between GlycoSP/CE levels and genome-wide transcripts, identified Gb4 d18:1/16:0 and CE 20:4 association with subnetworks enriched for T cell receptor signaling (PDCD1, CD86, PTPRC, CD247, IFNG) and DC-SIGN signaling (RAF1, CD209), respectively. Our findings reveal Gb3 and Gb4 globosides as sphingolipids associated with the resolution phase of inflammatory response in human macrophages.

Single-cell transcriptomics and surface epitope detection in human brain epileptic lesions identifies pro-inflammatory signaling

Epileptogenic triggers are multifactorial and not well understood. Here we aimed to address the hypothesis that inappropriate pro-inflammatory mechanisms contribute to the pathogenesis of refractory epilepsy (non-responsiveness to antiepileptic drugs) in human patients. We used single-cell cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) to reveal the immunotranscriptome of surgically resected epileptic lesion tissues. Our approach uncovered a pro-inflammatory microenvironment, including extensive activation of microglia and infiltration of other pro-inflammatory immune cells. These findings were supported by ligand–receptor (LR) interactome analysis, which demonstrated potential mechanisms of infiltration and evidence of direct physical interactions between microglia and T cells. Together, these data provide insight into the immune microenvironment in epileptic tissue, which may aid the development of new therapeutics.

Single-cell analysis of immune cells from surgically resected human epileptic brain tissues showed heterogeneity and pro-inflammatory signaling in microglia and evidence for direct interaction of microglia with T cells.

Cross-species integration of single-cell RNA-seq resolved alveolar-epithelial transitional states in idiopathic pulmonary fibrosis

Single-cell transcriptomics analyses of the fibrotic lung uncovered two cell states critical to lung injury recovery in the alveolar epithelium-a reparative transitional cell state in the mouse and a disease-specific cell state (KRT5^-/KRT17⁺) in human idiopathic pulmonary fibrosis (IPF). The murine transitional cell state lies between the differentiation from type 2 (AT2) to type 1 pneumocyte (AT1), and the human KRT5^-/KRT17⁺ cell state may arise from the dysregulation of this differentiation process. We review major findings of single-cell transcriptomics analyses of the fibrotic lung and reanalyzed data from seven single-cell RNA sequencing studies of human and murine models of IPF, focusing on the alveolar epithelium. Our comparative and cross-species single-cell transcriptomics analyses allowed us to further delineate the differentiation trajectories from AT2 to AT1 and AT2 to the KRT5^-/KRT17⁺ cell state. We observed AT1 cells in human IPF retain the transcriptional signature of the murine transitional cell state. Using pseudotime analysis, we recapitulated the differentiation trajectories from AT2 to AT1 and from AT2 to KRT5^-/KRT17⁺ cell state in multiple human IPF studies. We further delineated transcriptional programs underlying cell-state transitions and determined the molecular phenotypes at terminal differentiation. We hypothesize that in addition to the reactivation of developmental programs (SOX4, SOX9), senescence (TP63, SOX4) and the Notch pathway (HES1) are predicted to steer intermediate progenitors to the KRT5^-/KRT17⁺ cell state. Our analyses suggest that activation of SMAD3 later in the differentiation process may explain the fibrotic molecular phenotype typical of KRT5^-/KRT17⁺ cells.

Brugada syndrome genetics is associated with phenotype severity

Aims 

Brugada syndrome (BrS) is associated with an increased risk of sudden cardiac death due to ventricular tachycardia/fibrillation (VT/VF) in young, otherwise healthy individuals. Despite SCN5A being the most commonly known mutated gene to date, the genotype–phenotype relationship is poorly understood and remains uncertain. This study aimed to elucidate the genotype–phenotype correlation in BrS.

Methods and results

Brugada syndrome probands deemed at high risk of future arrhythmic events underwent genetic testing and phenotype characterization by the means of epicardial arrhythmogenic substrate (AS) mapping, and were divided into two groups according to the presence or absence of SCN5A mutation. Two-hundred probands (160 males, 80%; mean age 42.6 ± 12.2 years) were included in this study. Patients harbouring SCN5A mutations exhibited a spontaneous type 1 pattern and experienced aborted cardiac arrest or spontaneous VT/VF more frequently than the other subjects. SCN5A-positive patients exhibited a larger epicardial AS area, more prolonged electrograms and more frequently observed non-invasive late potentials. The presence of an SCN5A mutation explained >26% of the variation in the epicardial AS area and was the strongest predictor of a large epicardial area.

Conclusion 

In BrS, the genetic background is the main determinant for the extent of the electrophysiological abnormalities. SCN5A mutation carriers exhibit more pronounced epicardial electrical abnormalities and a more aggressive clinical presentation. These results contribute to the understanding of the genetic determinants of the BrS phenotypic expression and provide possible explanations for the varying degrees of disease expression.

WNT inhibition creates a BRCA-like state in Wnt-addicted cancer

Wnt signaling maintains diverse adult stem cell compartments and is implicated in chemotherapy resistance in cancer. PORCN inhibitors that block Wnt secretion have proven effective in Wnt-addicted preclinical cancer models and are in clinical trials. In a survey for potential combination therapies, we found that Wnt inhibition synergizes with the PARP inhibitor olaparib in Wnt-addicted cancers. Mechanistically, we find that multiple genes in the homologous recombination and Fanconi anemia repair pathways, including BRCA1, FANCD2, and RAD51, are dependent on Wnt/β-catenin signaling in Wnt-high cancers, and treatment with a PORCN inhibitor creates a BRCA-like state. This coherent regulation of DNA repair genes occurs in part via a Wnt/β-catenin/MYBL2 axis. Importantly, this pathway also functions in intestinal crypts, where high expression of BRCA and Fanconi anemia genes is seen in intestinal stem cells, with further upregulation in Wnt-high APCmin mutant polyps. Our findings suggest a general paradigm that Wnt/β-catenin signaling enhances DNA repair in stem cells and cancers to maintain genomic integrity. Conversely, interventions that block Wnt signaling may sensitize cancers to radiation and other DNA damaging agents.

Implication of sestrin3 in epilepsy and its comorbidities

Epilepsy is a serious neurological disorder affecting about 1% of the population worldwide. Epilepsy may arise as a result of acquired brain injury, or as a consequence of genetic predisposition. To date, genome-wide association studies and exome sequencing approaches have provided limited insights into the mechanisms of acquired brain injury. We have previously reported a pro-epileptic gene network, which is conserved across species, encoding inflammatory processes and positively regulated by sestrin3 (SESN3). In this study, we investigated the phenotype of SESN3 knock-out rats in terms of susceptibility to seizures and observed a significant delay in status epilepticus onset in SESN3 knock-out compared to control rats. This finding confirms previous in vitro and in vivo evidence indicating that SESN3 may favour occurrence and/or severity of seizures. We also analysed the phenotype of SESN3 knock-out rats for common comorbidities of epilepsy, i.e., anxiety, depression and cognitive impairment. SESN3 knock-out rats proved less anxious compared to control rats in a selection of behavioural tests. Taken together, the present results suggest that SESN3 may regulate mechanisms involved in the pathogenesis of epilepsy and its comorbidities.

Widespread Repression of Gene Expression in Cancer by a Wnt/β-Catenin/MAPK Pathway

Aberrant Wnt signaling drives a number of cancers through regulation of diverse downstream pathways. Wnt/β-catenin signaling achieves this in part by increasing the expression of proto-oncogenes such as MYC and cyclins. However, global assessment of the Wnt-regulated transcriptome in vivo in genetically distinct cancers demonstrates that Wnt signaling suppresses the expression of as many genes as it activates. In this study, we examined the set of genes that are upregulated upon inhibition of Wnt signaling in Wnt-addicted pancreatic and colorectal cancer models. Decreasing Wnt signaling led to a marked increase in gene expression by activating ERK and JNK; these changes in gene expression could be mitigated in part by concurrent inhibition of MEK. These findings demonstrate that increased Wnt signaling in cancer represses MAPK activity, preventing RAS-mediated senescence while allowing cancer cells to proliferate. These results shift the paradigm from Wnt/β-catenin primarily as an activator of transcription to a more nuanced view where Wnt/β-catenin signaling drives both widespread gene repression and activation. SIGNIFICANCE: These findings show that Wnt/β-catenin signaling causes widespread gene repression via inhibition of MAPK signaling, thus fine tuning the RAS-MAPK pathway to optimize proliferation in cancer.

Wnt-regulated lncRNA discovery enhanced by in vivo identification and CRISPRi functional validation

BACKGROUND

Wnt signaling is an evolutionarily conserved developmental pathway that is frequently hyperactivated in cancer. While multiple protein-coding genes regulated by Wnt signaling are known, the functional lncRNAs regulated by Wnt signaling have not been systematically characterized.

METHODS

We comprehensively mapped Wnt-regulated lncRNAs from an orthotopic Wnt-addicted pancreatic cancer model and examined the response of lncRNAs to Wnt inhibition between in vivo and in vitro cancer models. We further annotated and characterized these Wnt-regulated lncRNAs using existing genomic classifications (using data from FANTOM5) in the context of Wnt signaling and inferred their role in cancer pathogenesis (using GWAS and expression data from the TCGA). To functionally validate Wnt-regulated lncRNAs, we performed CRISPRi screens to assess their role in cancer cell proliferation both in vivo and in vitro.

RESULTS

We identified 3633 lncRNAs, of which 1503 were regulated by Wnt signaling in an orthotopic Wnt-addicted pancreatic cancer model. These lncRNAs were much more sensitive to changes in Wnt signaling in xenografts than in cultured cells. Our analysis suggested that Wnt signaling inhibition could influence the co-expression relationship of Wnt-regulated lncRNAs and their eQTL-linked protein-coding genes. Wnt-regulated lncRNAs were also implicated in specific gene networks involved in distinct biological processes that contribute to the pathogenesis of cancers. Consistent with previous genome-wide lncRNA CRISPRi screens, around 1% (13/1503) of the Wnt-regulated lncRNAs were found to modify cancer cell growth in vitro. This included CCAT1 and LINC00263, previously reported to regulate cancer growth. Using an in vivo CRISPRi screen, we doubled the discovery rate, identifying twice as many Wnt-regulated lncRNAs (25/1503) that had a functional effect on cancer cell growth.

CONCLUSIONS

Our study demonstrates the value of studying lncRNA functions in vivo, provides a valuable resource of lncRNAs regulated by Wnt signaling, and establishes a framework for systematic discovery of functional lncRNAs.

A trans-eQTL network regulates osteoclast multinucleation and bone mass

Functional characterisation of cell-type-specific regulatory networks is key to establish a causal link between genetic variation and phenotype. The osteoclast offers a unique model for interrogating the contribution of co-regulated genes to in vivo phenotype as its multinucleation and resorption activities determine quantifiable skeletal traits. Here we took advantage of a trans-regulated gene network (MMnet, macrophage multinucleation network) which we found to be significantly enriched for GWAS variants associated with bone-related phenotypes. We found that the network hub gene Bcat1 and seven other co-regulated MMnet genes out of 13, regulate bone function. Specifically, global (Pik3cb^-/-, Atp8b2^+/-, Igsf8^-/-, Eml1^-/-, Appl2^-/-, Deptor^-/-) and myeloid-specific Slc40a1 knockout mice displayed abnormal bone phenotypes. We report opposing effects of MMnet genes on bone mass in mice and osteoclast multinucleation/resorption in humans with strong correlation between the two. These results identify MMnet as a functionally conserved network that regulates osteoclast multinucleation and bone mass.

In Vivo Generation of Post-infarct Human Cardiac Muscle by Laminin-Promoted Cardiovascular Progenitors

Regeneration of injured human heart muscle is limited and an unmet clinical need. There are no methods for the reproducible generation of clinical-quality stem cell-derived cardiovascular progenitors (CVPs). We identified laminin-221 (LN-221) as the most likely expressed cardiac laminin. We produced it as human recombinant protein and showed that LN-221 promotes differentiation of pluripotent human embryonic stem cells (hESCs) toward cardiomyocyte lineage and downregulates pluripotency and teratoma-associated genes. We developed a chemically defined, xeno-free laminin-based differentiation protocol to generate CVPs. We show high reproducibility of the differentiation protocol using time-course bulk RNA sequencing developed from different hESC lines. Single-cell RNA sequencing of CVPs derived from hESC lines supported reproducibility and identified three main progenitor subpopulations. These CVPs were transplanted into myocardial infarction mice, where heart function was measured by echocardiogram and human heart muscle bundle formation was identified histologically. This method may provide clinical-quality cells for use in regenerative cardiology.

Impulsivity is a heritable trait in rodents and associated with a novel quantitative trait locus on chromosome 1

Impulsivity describes the tendency to act prematurely without appropriate foresight and is symptomatic of a number of neuropsychiatric disorders. Although a number of genes for impulsivity have been identified, no study to date has carried out an unbiased, genome-wide approach to identify genetic markers associated with impulsivity in experimental animals. Herein we report a linkage study of a six-generational pedigree of adult rats phenotyped for one dimension of impulsivity, namely premature responding on the five-choice serial reaction time task, combined with genome wide sequencing and transcriptome analysis to identify candidate genes associated with the expression of the impulsivity trait. Premature responding was found to be heritable (h² = 13-16%), with significant linkage (LOD 5.2) identified on chromosome 1. Fine mapping of this locus identified a number of polymorphic candidate genes, however only one, beta haemoglobin, was differentially expressed in both the founder strain and F6 generation. These findings provide novel insights into the genetic substrates and putative neurobiological mechanisms of impulsivity with broader translational relevance for impulsivity-related disorders in humans.

Whole-Genome Sequencing of Finnish Type 1 Diabetic Siblings Discordant for Kidney Disease Reveals DNA Variants associated with Diabetic Nephropathy

BACKGROUND

Several genetic susceptibility loci associated with diabetic nephropathy have been documented, but no causative variants implying novel pathogenetic mechanisms have been elucidated.

METHODS

We carried out whole-genome sequencing of a discovery cohort of Finnish siblings with type 1 diabetes who were discordant for the presence (case) or absence (control) of diabetic nephropathy. Controls had diabetes without complications for 15-37 years. We analyzed and annotated variants at genome, gene, and single-nucleotide variant levels. We then replicated the associated variants, genes, and regions in a replication cohort from the Finnish Diabetic Nephropathy study that included 3531 unrelated Finns with type 1 diabetes.

RESULTS

We observed protein-altering variants and an enrichment of variants in regions associated with the presence or absence of diabetic nephropathy. The replication cohort confirmed variants in both regulatory and protein-coding regions. We also observed that diabetic nephropathy-associated variants, when clustered at the gene level, are enriched in a core protein-interaction network representing proteins essential for podocyte function. These genes include protein kinases (protein kinase C isoforms ε and ι) and protein tyrosine kinase 2.

CONCLUSIONS

Our comprehensive analysis of a diabetic nephropathy cohort of siblings with type 1 diabetes who were discordant for kidney disease points to variants and genes that are potentially causative or protective for diabetic nephropathy. This includes variants in two isoforms of the protein kinase C family not previously linked to diabetic nephropathy, adding support to previous hypotheses that the protein kinase C family members play a role in diabetic nephropathy and might be attractive therapeutic targets.

Broad regulation of gene isoform expression by Wnt signaling in cancer

Differential gene isoform expression is a ubiquitous mechanism to enhance proteome diversity and maintain cell homeostasis. Mechanisms such as splicing that drive gene isoform variability are highly dynamic and responsive to changes in cell signaling pathways. Wnt/β-catenin signaling has profound effects on cell activity and cell fate and is known to modify several splicing events by altering the expression of individual splicing factors. However, a global assessment of how extensively Wnt signaling regulates splicing and other mechanisms that determine mRNA isoform composition in cancer is lacking. We used deep time-resolved RNA-seq in two independent in vivo Wnt-addicted tumor models during treatment with the potent Wnt inhibitor ETC-159 and examined Wnt regulated splicing events and splicing regulators. We found 1025 genes that underwent Wnt regulated variable exon usage leading to isoform expression changes. This was accompanied by extensive Wnt regulated changes in the expression of splicing regulators. Many of these Wnt regulated events were conserved in multiple human cancers, and many were linked to previously defined cancer-associated splicing quantitative trait loci. This suggests that the Wnt regulated splicing events are components of fundamental oncogenic processes. These findings demonstrate the wide-ranging effects of Wnt signaling on the isoform composition of the cell and provides an extensive resource of expression changes of splicing regulators and gene isoforms regulated by Wnt signaling.

A single-cell atlas of entorhinal cortex from individuals with Alzheimer's disease reveals cell-type-specific gene expression regulation

There is currently little information available about how individual cell types contribute to Alzheimer's disease. Here we applied single-nucleus RNA sequencing to entorhinal cortex samples from control and Alzheimer's disease brains (n = 6 per group), yielding a total of 13,214 high-quality nuclei. We detail cell-type-specific gene expression patterns, unveiling how transcriptional changes in specific cell subpopulations are associated with Alzheimer's disease. We report that the Alzheimer's disease risk gene APOE is specifically repressed in Alzheimer's disease oligodendrocyte progenitor cells and astrocyte subpopulations and upregulated in an Alzheimer's disease-specific microglial subopulation. Integrating transcription factor regulatory modules with Alzheimer's disease risk loci revealed drivers of cell-type-specific state transitions towards Alzheimer's disease. For example, transcription factor EB, a master regulator of lysosomal function, regulates multiple disease genes in a specific Alzheimer's disease astrocyte subpopulation. These results provide insights into the coordinated control of Alzheimer's disease risk genes and their cell-type-specific contribution to disease susceptibility. These results are available at http://adsn.ddnetbio.com.

Systems Genetics Approaches in Rat Identify Novel Genes and Gene Networks Associated With Cardiac Conduction

Background

Electrocardiographic (ECG) parameters are regarded as intermediate phenotypes of cardiac arrhythmias. Insight into the genetic underpinnings of these parameters is expected to contribute to the understanding of cardiac arrhythmia mechanisms. Here we used HXB/BXH recombinant inbred rat strains to uncover genetic loci and candidate genes modulating ECG parameters.

Methods and Results

RR interval, PR interval, QRS duration, and QTc interval were measured from ECGs obtained in 6 male rats from each of the 29 available HXB/BXH recombinant inbred strains. Genes at loci displaying significant quantitative trait loci (QTL) effects were prioritized by assessing the presence of protein‐altering variants, and by assessment of cis expression QTL (eQTL) effects and correlation of transcript abundance to the respective trait in the heart. Cardiac RNA‐seq data were additionally used to generate gene co‐expression networks. QTL analysis of ECG parameters identified 2 QTL for PR interval, respectively, on chromosomes 10 and 17. At the chromosome 10 QTL, cis‐eQTL effects were identified for Acbd4, Cd300lg, Fam171a2, and Arhgap27; the transcript abundance in the heart of these 4 genes was correlated with PR interval. At the chromosome 17 QTL, a cis‐eQTL was uncovered for Nhlrc1 candidate gene; the transcript abundance of this gene was also correlated with PR interval. Co‐expression analysis furthermore identified 50 gene networks, 6 of which were correlated with PR interval or QRS duration, both parameters of cardiac conduction.

Conclusions

These newly identified genetic loci and gene networks associated with the ECG parameters of cardiac conduction provide a starting point for future studies with the potential of identifying novel mechanisms underlying cardiac electrical function.

Transcriptional Networks of Microglia in Alzheimer's Disease and Insights into Pathogenesis

Microglia, the main immune cells of the central nervous system, are increasingly implicated in Alzheimer's disease (AD). Manifold transcriptomic studies in the brain have not only highlighted microglia's role in AD pathogenesis, but also mapped crucial pathological processes and identified new therapeutic targets. An important component of many of these transcriptomic studies is the investigation of gene expression networks in AD brain, which has provided important new insights into how coordinated gene regulatory programs in microglia (and other cell types) underlie AD pathogenesis. Given the rapid technological advancements in transcriptional profiling, spanning from microarrays to single-cell RNA sequencing (scRNA-seq), tools used for mapping gene expression networks have evolved to keep pace with the unique features of each transcriptomic platform. In this article, we review the trajectory of transcriptomic network analyses in AD from brain to microglia, highlighting the corresponding methodological developments. Lastly, we discuss examples of how transcriptional network analysis provides new insights into AD mechanisms and pathogenesis.

P6309Impaired secretion of clusterin in pericardial fluid of diabetics, a deleterious outcome for the cardiac micro-vasculature

Background

Extracellular vesicles are vital mediators of cell-to-cell communications. We previously demonstrated that in the adult heart, exosomes released in the pericardial fluid (PF) have the functional capability to induce cell survival and angiogenesis in recipient endothelial cells. It has previously been shown that PF exosomes from patients with an existing myocardial infarction (MI) contain clusterin (CLU). Remarkably, secreted CLU is cardio-protective as evidenced by enhanced arteriogenesis and reduced apoptosis in animal models of MI. We hypothesize that type 2 Diabetes Mellitus (T2DM) along with ischemic heart disease (IHD) alters PF EV content impairing its reparative angiogenic potential.

Aim

To investigate the role of T2DM in modulating human PF exosomal content and its vascular action.

Methods and results

PF samples were collected from 3 different groups of patients (N=3): 1) Ischemic heart disease (IHD) with or 2) without T2DM, and 3) non-ischaemic, non-diabetic controls operated for mitral valve repair.

The impact of PF exosomes on human coronary microvascular ECs (HCMECs) was evaluated by cell based functional assays for measuring apoptosis, proliferation and angiogenesis. Interestingly, unlike PF exosomes from control patients, PF exosomes from DM patients with IHD promoted EC apoptosis and impaired angiogenesis.

To investigate the effect of T2DM on PF content, we performed high throughput proteomic and metabolomic analysis of whole PF and PF EVs. We identified proteins and metabolites that were differentially expressed under T2DM condition by using the R package Limma. Employing a network approach, protein and metabolic data were integrated by using our newly developed inhouse R package Metabosignal. This networks approach revealed an interesting interaction circuit involving protein CLU. The T2DM-associated reduction of CLU level in the PF, identified from the “omics” analyses was confirmed by ELISA (n=8 patients per group).

To examine the effects of DM and ischemia on CLU intracellular expression and secretion, HCMECs were exposed to high D-glucose (HG, 25 mM) and hypoxia (1%O2). In line with our PF data, HCMECs exposed to HG and hypoxia showed lower levels of both CLU mRNA and secreted CLU (ELISA). Importantly, recombinant CLU (rCLU) treatment on HCMECs rescued their angiogenic potential, induced cell proliferation and protected them from apoptosis under both HG and hypoxia conditions. In addition, silencing of endogenous levels of CLU in HCMECs impaired angiogenesis, suggesting that CLU might be an inherent component of the angiogenesis machinery in ECs.

Conclusions

This study suggests that DM reduces the level of secreted CLU in the PF, depriving myocardial micro-vessels of this protective and regenerative mediator. Restoring diminished CLU levels in the diabetic heart could be a possible therapeutic approach for contrasting diabetic micro-angiopathy.

Acknowledgement/Funding

British Heart Foundation Program Grant

Author Correction: WWP2 regulates pathological cardiac fibrosis by modulating SMAD2 signaling

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

miR-15a/-16 Inhibit Angiogenesis by Targeting the Tie2 Coding Sequence: Therapeutic Potential of a miR-15a/16 Decoy System in Limb Ischemia

MicroRNA-15a (miR-15a) and miR-16, which are transcribed from the miR-15a/miR-16-1 cluster, inhibit post-ischemic angiogenesis. MicroRNA (miRNA) binding to mRNA coding sequences (CDSs) is a newly emerging mechanism of gene expression regulation. We aimed to (1) identify new mediators of the anti-angiogenic action of miR-15a and -16, (2) develop an adenovirus (Ad)-based miR-15a/16 decoy system carrying a luciferase reporter (Luc) to both sense and inhibit miR-15a/16 activity, and (3) investigate Ad.Luc-Decoy-15a/16 therapeutic potential in a mouse limb ischemia (LI) model. LI increased miR-15a and -16 expression in mouse muscular endothelial cells (ECs). The miRNAs also increased in cultured human umbilical vein ECs (HUVECs) exposed to serum starvation, but not hypoxia. Using bioinformatic tools and luciferase activity assays, we characterized miR-15a and -16 binding to Tie2 CDS. In HUVECs, miR-15a or -16 overexpression reduced Tie2 at the protein, but not the mRNA, level. Conversely, miR-15a or -16 inhibition improved angiogenesis in a Tie2-dependent manner. Local Ad.Luc-Decoy-15a/16 delivery increased Tie2 levels in ischemic skeletal muscle and improved post-LI angiogenesis and perfusion recovery, with reduced toe necrosis. Bioluminescent imaging (in vivo imaging system [IVIS]) provided evidence that the Ad.Luc-Decoy-15a/16 system responds to miR-15a/16 increases. In conclusion, we have provided novel mechanistic evidence of the therapeutic potential of local miR-15a/16 inhibition in LI.

Integrated systems-genetic analyses reveal a network target for delaying glioma progression

Objective

To identify a convergent, multitarget proliferation characteristic for astrocytoma transformation that could be targeted for therapy discovery.

Methods

Using an integrated functional genomics approach, we prioritized networks associated with astrocytoma progression using the following criteria: differential co‐expression between grade II and grade III IDH1‐mutated and 1p/19q euploid astrocytomas, preferential enrichment for genetic risk to cancer, association with patient survival and sample‐level genomic features. Drugs targeting the identified multitarget network characteristic for astrocytoma transformation were computationally predicted using drug transcriptional perturbation data and validated using primary human astrocytoma cells.

Results

A single network, M2, consisting of 177 genes, was associated with glioma progression on the basis of the above criteria. Functionally, M2 encoded physically interacting proteins regulating cell cycle processes and analysis of genome‐wide gene‐regulatory interactions using mutual information and DNA–protein interactions revealed the known regulators of cell cycle processes FoxM1, B‐Myb, and E2F2 as key regulators of M2. These results suggest functional disruption of M2 via gene mutation or altered expression as a convergent pathway regulating astrocytoma transformation. By considering M2 as a multitarget drug target regulating astrocytoma transformation, we identified several drugs that are predicted to restore M2 expression in anaplastic astrocytoma toward its low‐grade profile and of these, we validated the known antiproliferative drug resveratrol as down‐regulating multiple nodes of M2 including at nanomolar concentrations achievable in human cerebrospinal fluid by oral dosing.

Interpretation

Our results identify M2 as a multitarget network characteristic for astrocytoma progression and encourage M2‐based drug screening to identify new compounds for preventing glioma transformation.

WWP2 regulates pathological cardiac fibrosis by modulating SMAD2 signaling

Cardiac fibrosis is a final common pathology in inherited and acquired heart diseases that causes cardiac electrical and pump failure. Here, we use systems genetics to identify a pro-fibrotic gene network in the diseased heart and show that this network is regulated by the E3 ubiquitin ligase WWP2, specifically by the WWP2-N terminal isoform. Importantly, the WWP2-regulated pro-fibrotic gene network is conserved across different cardiac diseases characterized by fibrosis: human and murine dilated cardiomyopathy and repaired tetralogy of Fallot. Transgenic mice lacking the N-terminal region of the WWP2 protein show improved cardiac function and reduced myocardial fibrosis in response to pressure overload or myocardial infarction. In primary cardiac fibroblasts, WWP2 positively regulates the expression of pro-fibrotic markers and extracellular matrix genes. TGFβ1 stimulation promotes nuclear translocation of the WWP2 isoforms containing the N-terminal region and their interaction with SMAD2. WWP2 mediates the TGFβ1-induced nucleocytoplasmic shuttling and transcriptional activity of SMAD2.

Abstract 130: Development of Bioinspired Synthetic Exosomes With Proangiogenic Potential

Exosomes with variable microRNA cargos are released from different cell types and stimulate angiogenesis in animal models. Therefore, exosomes are currently considered for their potential to represent a safer and easier biological-based alternative to stem cell therapies. However, the limited amount of exosomes that can be reproducibly prepared from cultured stem cells and the complexity of the cargo of endogenous exosomes make them difficult to develop as off-the-shelf “pharmaceutical products” and limit their clinical potential. These issues could be surpassed by the use of artificial (i.e. nanoparticles mimicking the natural product) exosomes carrying a therapeutic cargo. In consideration of the possibility that not all cargo components of exosomes are required for their therapeutic functions, we reasoned that bioinspired artificial exosomes (AEs) incorporating key therapeutically relevant molecules and devoid of potentially negative or indifferent factors could further improve their therapeutic potential. Therefore, we developed and tested AEs containing exosomal proangiogenic miRNAs as a “proof-of-concept” for therapeutic angiogenesis. In order to achieve this, we: (1) characterized the common microRNA cargo of endogenous proangiogenic exosomes (from stem cells, pericytes and pericardial fluid) using bioinformatics, (2) exploited this knowledge to develop off-the-shelf artificial exosomes (AEs) with proangiogenic capacities, (3) validated the angiogenic potential of the bioinspired AEs. Bioinformatics analyses integrating data of miRNA arrays and panels of proangiogenic exosomes confirmed the enrichment of let-7 family in these exosomes. After testing the angiogenic potential of different members of let-7 family, we produced AEs containing let-7b (the most potent in the family) by microfluidic micromixing. The AEs were uptaken by human ECs cultured under hypoxic conditions, without causing toxicity. let-7b-AEs transferred functional let-7b, thus decreased the expression of validated targets of let-7b in recipient cells. let-7b-AEs improved EC survival, proliferation and angiogenesis in vitro and in vivo . These data suggest the therapeutic potential of bioinspired AEs containing let-7b.

The State of Systems Genetics in 2017

Cell Systems invited 16 experts to share their views on the field of systems genetics. In questions repeated in the headings, we asked them to define systems genetics, highlight its relevance to researchers outside the field, discuss what makes a strong systems genetics paper, and paint a picture of where the field is heading in the coming years. Their responses, ordered by the journal but otherwise unedited, make it clear that deciphering genotype to phenotype relationships is a central challenge of systems genetics and will require understanding how networks and higher-order properties of biological systems underlie complex traits. In addition, our experts illuminate the applications and relevance of systems genetics to human disease, the gut microbiome, development of tools that connect the global research community, sustainability, drug discovery, patient-specific disease and network models, and personalized treatments. Finally, a table of suggested reading provides a sample of influential work in the field.

Sleeve gastrectomy rapidly enhances islet function independently of body weight

Bariatric surgeries including vertical sleeve gastrectomy (VSG) ameliorate obesity and diabetes. Weight loss and accompanying increases to insulin sensitivity contribute to improved glycemia after surgery; however, studies in humans also suggest weight-independent actions of bariatric procedures to lower blood glucose, possibly by improving insulin secretion. To evaluate this hypothesis, we compared VSG-operated mice with pair-fed, sham-surgical controls (PF-Sham) 2 weeks after surgery. This paradigm yielded similar postoperative body weight and insulin sensitivity between VSG and calorically restricted PF-Sham animals. However, VSG improved glucose tolerance and markedly enhanced insulin secretion during oral nutrient and i.p. glucose challenges compared with controls. Islets from VSG mice displayed a unique transcriptional signature enriched for genes involved in Ca2+ signaling and insulin secretion pathways. This finding suggests that bariatric surgery leads to intrinsic changes within the islet that alter function. Indeed, islets isolated from VSG mice had increased glucose-stimulated insulin secretion and a left-shifted glucose sensitivity curve compared with islets from PF-Sham mice. Isolated islets from VSG animals showed corresponding increases in the pulse duration of glucose-stimulated Ca2+ oscillations. Together, these findings demonstrate a weight-independent improvement in glycemic control following VSG, which is, in part, driven by improved insulin secretion and associated with substantial changes in islet gene expression. These results support a model in which β cells play a key role in the adaptation to bariatric surgery and the improved glucose tolerance that is typical of these procedures.

Optimisation of laboratory methods for whole transcriptomic RNA analyses in human left ventricular biopsies and blood samples of clinical relevance

This study aimed to optimise techniques for whole transcriptome and small RNA analyses on clinical tissue samples from patients with cardiovascular disease. Clinical samples often represent a particular challenge to extracting RNA of sufficient quality for robust RNA sequencing analysis, and due to availability, it is rarely possible to optimise techniques on the samples themselves. Therefore, we have used equivalent samples from pigs undergoing cardiopulmonary bypass surgery to test different protocols for optimal RNA extraction, and then validated the protocols in human samples. Here we present an assessment of the quality and quantity of RNA obtained using a variety of commercially-available RNA extraction kits on both left ventricular biopsies and blood plasma. RNA extraction from these samples presents different difficulties; left ventricular biopsies are small and fibrous, while blood plasma has a low RNA content. We have validated our optimised extraction techniques on human clinical samples collected as part of the ARCADIA (Association of non-coding RNAs with Coronary Artery Disease and type 2 Diabetes) cohort study, resulting in successful whole transcriptome and small RNA sequencing of human left ventricular tissue.

Ethnicity-Specific Skeletal Muscle Transcriptional Signatures and Their Relevance to Insulin Resistance in Singapore

CONTEXT

Insulin resistance (IR) and obesity differ among ethnic groups in Singapore, with the Malays more obese yet less IR than Asian-Indians. However, the molecular basis underlying these differences is not clear.

OBJECTIVE

As the skeletal muscle (SM) is metabolically relevant to IR, we investigated molecular pathways in SM that are associated with ethnic differences in IR, obesity, and related traits.

DESIGN, SETTING, AND MAIN OUTCOME MEASURES

We integrated transcriptomic, genomic, and phenotypic analyses in 156 healthy subjects representing three major ethnicities in the Singapore Adult Metabolism Study.

PATIENTS

This study contains Chinese (n = 63), Malay (n = 51), and Asian-Indian (n = 42) men, aged 21 to 40 years, without systemic diseases.

RESULTS

We found remarkable diversity in the SM transcriptome among the three ethnicities, with >8000 differentially expressed genes (40% of all genes expressed in SM). Comparison with blood transcriptome from a separate Singaporean cohort showed that >95% of SM expression differences among ethnicities were unique to SM. We identified a network of 46 genes that were specifically downregulated in Malays, suggesting dysregulation of components of cellular respiration in SM of Malay individuals. We also report 28 differentially expressed gene clusters, four of which were also enriched for genes that were found in genome-wide association studies of metabolic traits and disease and correlated with variation in IR, obesity, and related traits.

CONCLUSION

We identified extensive gene-expression changes in SM among the three Singaporean ethnicities and report specific genes and molecular pathways that might underpin and explain the differences in IR among these ethnic groups.

Systems genetics identifies a macrophage cholesterol network associated with physiological wound healing

Among other cells, macrophages regulate the inflammatory and reparative phases during wound healing but genetic determinants and detailed molecular pathways that modulate these processes are not fully elucidated. Here, we took advantage of normal variation in wound healing in 1,378 genetically outbred mice, and carried out macrophage RNA-sequencing profiling of mice with extreme wound healing phenotypes (i.e., slow and fast healers, n = 146 in total). The resulting macrophage coexpression networks were genetically mapped and led to the identification of a unique module under strong trans-acting genetic control by the Runx2 locus. This macrophage-mediated healing network was specifically enriched for cholesterol and fatty acid biosynthetic processes. Pharmacological blockage of fatty acid synthesis with cerulenin resulted in delayed wound healing in vivo, and increased macrophage infiltration in the wounded skin, suggesting the persistence of an unresolved inflammation. We show how naturally occurring sequence variation controls transcriptional networks in macrophages, which in turn regulate specific metabolic pathways that could be targeted in wound healing.

MWASTools: an R/bioconductor package for metabolome-wide association studies

Summary

MWASTools is an R package designed to provide an integrated pipeline to analyse metabonomic data in large-scale epidemiological studies. Key functionalities of our package include: quality control analysis; metabolome-wide association analysis using various models (partial correlations, generalized linear models); visualization of statistical outcomes; metabolite assignment using statistical total correlation spectroscopy (STOCSY); and biological interpretation of metabolome-wide association studies results.

Availability and implementation

The MWASTools R package is implemented in R (version  > =3.4) and is available from Bioconductor: https://bioconductor.org/packages/MWASTools/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Temporal dynamics of Wnt-dependent transcriptome reveal an oncogenic Wnt/MYC/ribosome axis

Activating mutations in the Wnt pathway drive a variety of cancers, but the specific targets and pathways activated by Wnt ligands are not fully understood. To bridge this knowledge gap, we performed a comprehensive time-course analysis of Wnt-dependent signaling pathways in an orthotopic model of Wnt-addicted pancreatic cancer, using a porcupine (PORCN) inhibitor currently in clinical trials, and validated key results in additional Wnt-addicted models. The temporal analysis of the drug-perturbed transcriptome demonstrated direct and indirect regulation of more than 3,500 Wnt-activated genes (23% of the transcriptome). Regulation was both via Wnt/β-catenin and through the modulation of protein abundance of important transcription factors, including MYC, via Wnt-dependent stabilization of proteins (Wnt/STOP). Our study identifies a central role of Wnt/β-catenin and Wnt/STOP signaling in controlling ribosome biogenesis, a key driver of cancer proliferation.

Biologically relevant laminin as chemically defined and fully human platform for human epidermal keratinocyte culture

The current expansion of autologous human keratinocytes to resurface severe wound defects still relies on murine feeder layer and calf serum in the cell culture system. Through our characterization efforts of the human skin basement membrane and murine feeder layer 3T3-J2, we identified two biologically relevant recombinant laminins-LN-511 and LN-421- as potential candidates to replace the murine feeder. Herein, we report a completely xeno-free and defined culture system utilizing these laminins which enables robust expansion of adult human skin keratinocytes. We demonstrate that our laminin system is comparable to the 3T3-J2 co-culture system in terms of basal markers' profile, colony-forming efficiency and the ability to form normal stratified epidermal structure in both in vitro and in vivo models. These results show that the proposed system may not only provide safer keratinocyte use in the clinics, but also facilitate the broader use of other cultured human epithelial cells in regenerative medicine.

A systems-level framework for drug discovery identifies Csf1R as an anti-epileptic drug target

The identification of drug targets is highly challenging, particularly for diseases of the brain. To address this problem, we developed and experimentally validated a general computational framework for drug target discovery that combines gene regulatory information with causal reasoning (“Causal Reasoning Analytical Framework for Target discovery”—CRAFT). Using a systems genetics approach and starting from gene expression data from the target tissue, CRAFT provides a predictive framework for identifying cell membrane receptors with a direction-specified influence over disease-related gene expression profiles. As proof of concept, we applied CRAFT to epilepsy and predicted the tyrosine kinase receptor Csf1R as a potential therapeutic target. The predicted effect of Csf1R blockade in attenuating epilepsy seizures was validated in three pre-clinical models of epilepsy. These results highlight CRAFT as a systems-level framework for target discovery and suggest Csf1R blockade as a novel therapeutic strategy in epilepsy. CRAFT is applicable to disease settings other than epilepsy.

The identification of new drug targets is highly challenging, particularly for diseases of the brain. This study describes a general computational gene regulatory framework called CRAFT for drug target discovery, and the authors use CRAFT to identify the microglial membrane receptor Csf1R as a potential therapeutic target for epilepsy.

Control of inducible gene expression links cohesin to hematopoietic progenitor self-renewal and differentiation

Cohesin is important for 3D genome organization. Nevertheless, even the complete removal of cohesin has surprisingly little impact on steady-state gene transcription and enhancer activity. Here we show that cohesin is required for the core transcriptional response of primary macrophages to microbial signals, and for inducible enhancer activity that underpins inflammatory gene expression. Consistent with a role for inflammatory signals in promoting myeloid differentiation of hematopoietic stem and progenitor cells (HPSCs), cohesin mutations in HSPCs led to reduced inflammatory gene expression and increased resistance to differentiation-inducing inflammatory stimuli. These findings uncover an unexpected dependence of inducible gene expression on cohesin, link cohesin with myeloid differentiation, and may help explain the prevalence of cohesin mutations in human acute myeloid leukemia.

Common signalling pathways in macrophage and osteoclast multinucleation

Macrophage cell fusion and multinucleation are fundamental processes in the formation of multinucleated giant cells (MGCs) in chronic inflammatory disease and osteoclasts in the regulation of bone mass. However, this basic cell phenomenon is poorly understood despite its pathophysiological relevance. Granulomas containing multinucleated giant cells are seen in a wide variety of complex inflammatory disorders, as well as in infectious diseases. Dysregulation of osteoclastic bone resorption underlies the pathogenesis of osteoporosis and malignant osteolytic bone disease. Recent reports have shown that the formation of multinucleated giant cells and osteoclast fusion display a common molecular signature, suggesting shared genetic determinants. In this Review, we describe the background of cell-cell fusion and the similar origin of macrophages and osteoclasts. We specifically focus on the common pathways involved in osteoclast and MGC fusion. We also highlight potential approaches that could help to unravel the core mechanisms underlying bone and granulomatous disorders in humans.