Single nuclei transcriptomics of the in situ human limbal stem cell niche

The corneal epithelium acts as a barrier to pathogens entering the eye; corneal epithelial cells are continuously renewed by uni-potent, quiescent limbal stem cells (LSCs) located at the limbus, where the cornea transitions to conjunctiva. There has yet to be a consensus on LSC markers and their transcriptome profile is not fully understood, which may be due to using cadaveric tissue without an intact stem cell niche for transcriptomics. In this study, we addressed this problem by using single nuclei RNA sequencing (snRNAseq) on healthy human limbal tissue that was immediately snap-frozen after excision from patients undergoing cataract surgery. We identified the quiescent LSCs as a sub-population of corneal epithelial cells with a low level of total transcript counts. Moreover, TP63, KRT15, CXCL14, and ITGβ4 were found to be highly expressed in LSCs and transiently amplifying cells (TACs), which constitute the corneal epithelial progenitor populations at the limbus. The surface markers SLC6A6 and ITGβ4 could be used to enrich human corneal epithelial cell progenitors, which were also found to specifically express the putative limbal progenitor cell markers MMP10 and AC093496.1.

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.

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.

High-dose corticosteroid pulse therapy increases the survival rate in COVID-19 patients at risk of hyper-inflammatory response

OBJECTIVE

Test whether high dose corticosteroid pulse therapy (HDCPT) with either methylprednisolone or dexamethasone is associated with increased survival in COVID-19 patients at risk of hyper-inflammatory response. Provide some initial diagnostic criteria using laboratory markers to stratify these patients.

METHODS

This is a prospective observational study, 318 met the inclusion criteria. 64 patients (20.1%) were treated with HDCPT by using at least 1.5mg/kg/24h of methylprednisolone or dexamethasone equivalent. A multivariate Cox regression (controlling for co-morbidities and other therapies) was carried out to determine whether HDCPT (among other interventions) was associated with decreased mortality. We also carried out a 30-day time course analysis of laboratory markers between survivors and non-survivors, to identify potential markers for patient stratification.

RESULTS

HDCPT showed a statistically significant decrease in mortality (HR = 0.087 [95% CI 0.021-0.36]; P < 0.001). 30-day time course analysis of laboratory marker tests showed marked differences in pro-inflammatory markers between survivors and non-survivors. As diagnostic criteria to define the patients at risk of developing a COVID-19 hyper-inflammatory response, we propose the following parameters (IL-6 > = 40 pg/ml, and/or two of the following: C-reactive protein > = 100 mg/L, D-dimer > = 1000 ng/ml, ferritin > = 500 ng/ml and lactate dehydrogenase > = 300 U/L).

CONCLUSIONS

HDCPT can be an effective intervention to increase COVID-19 survival rates in patients at risk of developing a COVID-19 hyper-inflammatory response, laboratory marker tests can be used to stratify these patients who should be given HDCPT. This study is not a randomized clinical trial (RCT). Future RCTs should be carried out to confirm the efficacy of HDCPT to increase the survival rates of COVID-19.

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.

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.

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.

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.

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.

Epoxygenase inactivation exacerbates diet and aging-associated metabolic dysfunction resulting from impaired adipogenesis

OBJECTIVE

When molecular drivers of healthy adipogenesis are perturbed, this can cause hepatic steatosis. The role of arachidonic acid (AA) and its downstream enzymatic cascades, such as cyclooxygenase, in adipogenesis is well established. The exact contribution of the P450 epoxygenase pathway, however, remains to be established. Enzymes belonging to this pathway are mainly encoded by the CYP2J locus which shows extensive allelic expansion in mice. Here we aimed to establish the role of endogenous epoxygenase during adipogenesis under homeostatic and metabolic stress conditions.

METHODS

We took advantage of the simpler genetic architecture of the Cyp2j locus in the rat and used a Cyp2j4 (orthologue of human CYP2J2) knockout rat in two models of metabolic dysfunction: physiological aging and cafeteria diet (CAF). The phenotyping of Cyp2j4-/- rats under CAF was integrated with proteomics (LC-MS/MS) and lipidomics (LC-MS) analyses in the liver and the adipose tissue.

RESULTS

We report that Cyp2j4 deletion causes adipocyte dysfunction under metabolic challenges. This is characterized by (i) down-regulation of white adipose tissue (WAT) PPARγ and C/EBPα, (ii) adipocyte hypertrophy, (iii) extracellular matrix remodeling, and (iv) alternative usage of AA pathway. Specifically, in Cyp2j4-/- rats treated with a cafeteria diet, the dysfunctional adipogenesis is accompanied by exacerbated weight gain, hepatic lipid accumulation, and dysregulated gluconeogenesis.

CONCLUSION

These results suggest that AA epoxygenases are essential regulators of healthy adipogenesis. Our results uncover their synergistic role in fine-tuning AA pathway in obesity-mediated hepatic steatosis.

Systems biology of the IMIDIA biobank from organ donors and pancreatectomised patients defines a novel transcriptomic signature of islets from individuals with type 2 diabetes

AIMS/HYPOTHESIS

Pancreatic islet beta cell failure causes type 2 diabetes in humans. To identify transcriptomic changes in type 2 diabetic islets, the Innovative Medicines Initiative for Diabetes: Improving beta-cell function and identification of diagnostic biomarkers for treatment monitoring in Diabetes (IMIDIA) consortium ( www.imidia.org ) established a comprehensive, unique multicentre biobank of human islets and pancreas tissues from organ donors and metabolically phenotyped pancreatectomised patients (PPP).

METHODS

Affymetrix microarrays were used to assess the islet transcriptome of islets isolated either by enzymatic digestion from 103 organ donors (OD), including 84 non-diabetic and 19 type 2 diabetic individuals, or by laser capture microdissection (LCM) from surgical specimens of 103 PPP, including 32 non-diabetic, 36 with type 2 diabetes, 15 with impaired glucose tolerance (IGT) and 20 with recent-onset diabetes (<1 year), conceivably secondary to the pancreatic disorder leading to surgery (type 3c diabetes). Bioinformatics tools were used to (1) compare the islet transcriptome of type 2 diabetic vs non-diabetic OD and PPP as well as vs IGT and type 3c diabetes within the PPP group; and (2) identify transcription factors driving gene co-expression modules correlated with insulin secretion ex vivo and glucose tolerance in vivo. Selected genes of interest were validated for their expression and function in beta cells.

RESULTS

Comparative transcriptomic analysis identified 19 genes differentially expressed (false discovery rate ≤0.05, fold change ≥1.5) in type 2 diabetic vs non-diabetic islets from OD and PPP. Nine out of these 19 dysregulated genes were not previously reported to be dysregulated in type 2 diabetic islets. Signature genes included TMEM37, which inhibited Ca²⁺-influx and insulin secretion in beta cells, and ARG2 and PPP1R1A, which promoted insulin secretion. Systems biology approaches identified HNF1A, PDX1 and REST as drivers of gene co-expression modules correlated with impaired insulin secretion or glucose tolerance, and 14 out of 19 differentially expressed type 2 diabetic islet signature genes were enriched in these modules. None of these signature genes was significantly dysregulated in islets of PPP with impaired glucose tolerance or type 3c diabetes.

CONCLUSIONS/INTERPRETATION

These studies enabled the stringent definition of a novel transcriptomic signature of type 2 diabetic islets, regardless of islet source and isolation procedure. Lack of this signature in islets from PPP with IGT or type 3c diabetes indicates differences possibly due to peculiarities of these hyperglycaemic conditions and/or a role for duration and severity of hyperglycaemia. Alternatively, these transcriptomic changes capture, but may not precede, beta cell failure.

Changes in macrophage transcriptome associate with systemic sclerosis and mediate GSDMA contribution to disease risk

OBJECTIVES

Several common and rare risk variants have been reported for systemic sclerosis (SSc), but the effector cell(s) mediating the function of these genetic variants remains to be elucidated. While innate immune cells have been proposed as the critical targets to interfere with the disease process underlying SSc, no studies have comprehensively established their effector role. Here we investigated the contribution of monocyte-derived macrophages (MDMs) in mediating genetic susceptibility to SSc.

METHODS

We carried out RNA sequencing and genome-wide genotyping in MDMs from 57 patients with SSc and 15 controls. Our differential expression and expression quantitative trait locus (eQTL) analysis in SSc was further integrated with epigenetic, expression and eQTL data from skin, monocytes, neutrophils and lymphocytes.

RESULTS

We identified 602 genes upregulated and downregulated in SSc macrophages that were significantly enriched for genes previously implicated in SSc susceptibility (P=5×10^-4), and 270 cis-regulated genes in MDMs. Among these, GSDMA was reported to carry an SSc risk variant (rs3894194) regulating expression of neighbouring genes in blood. We show that GSDMA is upregulated in SSc MDMs (P=8.4×10^-4) but not in the skin, and is a significant eQTL in SSc macrophages and lipopolysaccharide/interferon gamma (IFNγ)-stimulated monocytes. Furthermore, we identify an SSc macrophage transcriptome signature characterised by upregulation of glycolysis, hypoxia and mTOR signalling and a downregulation of IFNγ response pathways.

CONCLUSIONS

Our data further establish the link between macrophages and SSc, and suggest that the contribution of the rs3894194 risk variant to SSc susceptibility can be mediated by GSDMA expression in macrophages.

IL-11 is a crucial determinant of cardiovascular fibrosis

Fibrosis is a common pathology in cardiovascular disease. In the heart, fibrosis causes mechanical and electrical dysfunction and in the kidney, it predicts the onset of renal failure. Transforming growth factor β1 (TGFβ1) is the principal pro-fibrotic factor, but its inhibition is associated with side effects due to its pleiotropic roles. We hypothesized that downstream effectors of TGFβ1 in fibroblasts could be attractive therapeutic targets and lack upstream toxicity. Here we show, using integrated imaging-genomics analyses of primary human fibroblasts, that upregulation of interleukin-11 (IL-11) is the dominant transcriptional response to TGFβ1 exposure and required for its pro-fibrotic effect. IL-11 and its receptor (IL11RA) are expressed specifically in fibroblasts, in which they drive non-canonical, ERK-dependent autocrine signalling that is required for fibrogenic protein synthesis. In mice, fibroblast-specific Il11 transgene expression or Il-11 injection causes heart and kidney fibrosis and organ failure, whereas genetic deletion of Il11ra1 protects against disease. Therefore, inhibition of IL-11 prevents fibroblast activation across organs and species in response to a range of important pro-fibrotic stimuli. These results reveal a central role of IL-11 in fibrosis and we propose that inhibition of IL-11 is a potential therapeutic strategy to treat fibrotic diseases.

Natural genetic variation of the cardiac transcriptome in non-diseased donors and patients with dilated cardiomyopathy

BACKGROUND

Genetic variation is an important determinant of RNA transcription and splicing, which in turn contributes to variation in human traits, including cardiovascular diseases.

RESULTS

Here we report the first in-depth survey of heart transcriptome variation using RNA-sequencing in 97 patients with dilated cardiomyopathy and 108 non-diseased controls. We reveal extensive differences of gene expression and splicing between dilated cardiomyopathy patients and controls, affecting known as well as novel dilated cardiomyopathy genes. Moreover, we show a widespread effect of genetic variation on the regulation of transcription, isoform usage, and allele-specific expression. Systematic annotation of genome-wide association SNPs identifies 60 functional candidate genes for heart phenotypes, representing 20% of all published heart genome-wide association loci. Focusing on the dilated cardiomyopathy phenotype we found that eQTL variants are also enriched for dilated cardiomyopathy genome-wide association signals in two independent cohorts.

CONCLUSIONS

RNA transcription, splicing, and allele-specific expression are each important determinants of the dilated cardiomyopathy phenotype and are controlled by genetic factors. Our results represent a powerful resource for the field of cardiovascular genetics.

Systems Genetics as a Tool to Identify Master Genetic Regulators in Complex Disease

Systems genetics stems from systems biology and similarly employs integrative modeling approaches to describe the perturbations and phenotypic effects observed in a complex system. However, in the case of systems genetics the main source of perturbation is naturally occurring genetic variation, which can be analyzed at the systems-level to explain the observed variation in phenotypic traits. In contrast with conventional single-variant association approaches, the success of systems genetics has been in the identification of gene networks and molecular pathways that underlie complex disease. In addition, systems genetics has proven useful in the discovery of master trans-acting genetic regulators of functional networks and pathways, which in many cases revealed unexpected gene targets for disease. Here we detail the central components of a fully integrated systems genetics approach to complex disease, starting from assessment of genetic and gene expression variation, linking DNA sequence variation to mRNA (expression QTL mapping), gene regulatory network analysis and mapping the genetic control of regulatory networks. By summarizing a few illustrative (and successful) examples, we highlight how different data-modeling strategies can be effectively integrated in a systems genetics study.

Rare and common epilepsies converge on a shared gene regulatory network providing opportunities for novel antiepileptic drug discovery

BACKGROUND

The relationship between monogenic and polygenic forms of epilepsy is poorly understood and the extent to which the genetic and acquired epilepsies share common pathways is unclear. Here, we use an integrated systems-level analysis of brain gene expression data to identify molecular networks disrupted in epilepsy.

RESULTS

We identified a co-expression network of 320 genes (M30), which is significantly enriched for non-synonymous de novo mutations ascertained from patients with monogenic epilepsy and for common variants associated with polygenic epilepsy. The genes in the M30 network are expressed widely in the human brain under tight developmental control and encode physically interacting proteins involved in synaptic processes. The most highly connected proteins within the M30 network were preferentially disrupted by deleterious de novo mutations for monogenic epilepsy, in line with the centrality-lethality hypothesis. Analysis of M30 expression revealed consistent downregulation in the epileptic brain in heterogeneous forms of epilepsy including human temporal lobe epilepsy, a mouse model of acquired temporal lobe epilepsy, and a mouse model of monogenic Dravet (SCN1A) disease. These results suggest functional disruption of M30 via gene mutation or altered expression as a convergent mechanism regulating susceptibility to epilepsy broadly. Using the large collection of drug-induced gene expression data from Connectivity Map, several drugs were predicted to preferentially restore the downregulation of M30 in epilepsy toward health, most notably valproic acid, whose effect on M30 expression was replicated in neurons.

CONCLUSIONS

Taken together, our results suggest targeting the expression of M30 as a potential new therapeutic strategy in epilepsy.

From integrative genomics to systems genetics in the rat to link genotypes to phenotypes

Complementary to traditional gene mapping approaches used to identify the hereditary components of complex diseases, integrative genomics and systems genetics have emerged as powerful strategies to decipher the key genetic drivers of molecular pathways that underlie disease. Broadly speaking, integrative genomics aims to link cellular-level traits (such as mRNA expression) to the genome to identify their genetic determinants. With the characterization of several cellular-level traits within the same system, the integrative genomics approach evolved into a more comprehensive study design, called systems genetics, which aims to unravel the complex biological networks and pathways involved in disease, and in turn map their genetic control points. The first fully integrated systems genetics study was carried out in rats, and the results, which revealed conserved trans-acting genetic regulation of a pro-inflammatory network relevant to type 1 diabetes, were translated to humans. Many studies using different organisms subsequently stemmed from this example. The aim of this Review is to describe the most recent advances in the fields of integrative genomics and systems genetics applied in the rat, with a focus on studies of complex diseases ranging from inflammatory to cardiometabolic disorders. We aim to provide the genetics community with a comprehensive insight into how the systems genetics approach came to life, starting from the first integrative genomics strategies [such as expression quantitative trait loci (eQTLs) mapping] and concluding with the most sophisticated gene network-based analyses in multiple systems and disease states. Although not limited to studies that have been directly translated to humans, we will focus particularly on the successful investigations in the rat that have led to primary discoveries of genes and pathways relevant to human disease.

MT-HESS: an efficient Bayesian approach for simultaneous association detection in OMICS datasets, with application to eQTL mapping in multiple tissues

MOTIVATION

Analysing the joint association between a large set of responses and predictors is a fundamental statistical task in integrative genomics, exemplified by numerous expression Quantitative Trait Loci (eQTL) studies. Of particular interest are the so-called ': hotspots ': , important genetic variants that regulate the expression of many genes. Recently, attention has focussed on whether eQTLs are common to several tissues, cell-types or, more generally, conditions or whether they are specific to a particular condition.

RESULTS

We have implemented MT-HESS, a Bayesian hierarchical model that analyses the association between a large set of predictors, e.g. SNPs, and many responses, e.g. gene expression, in multiple tissues, cells or conditions. Our Bayesian sparse regression algorithm goes beyond ': one-at-a-time ': association tests between SNPs and responses and uses a fully multivariate model search across all linear combinations of SNPs, coupled with a model of the correlation between condition/tissue-specific responses. In addition, we use a hierarchical structure to leverage shared information across different genes, thus improving the detection of hotspots. We show the increase of power resulting from our new approach in an extensive simulation study. Our analysis of two case studies highlights new hotspots that would remain undetected by standard approaches and shows how greater prediction power can be achieved when several tissues are jointly considered.

AVAILABILITY AND IMPLEMENTATION

C[Formula: see text] source code and documentation including compilation instructions are available under GNU licence at http://www.mrc-bsu.cam.ac.uk/software/.

Multi-tissue analysis of co-expression networks by higher-order generalized singular value decomposition identifies functionally coherent transcriptional modules

Recent high-throughput efforts such as ENCODE have generated a large body of genome-scale transcriptional data in multiple conditions (e.g., cell-types and disease states). Leveraging these data is especially important for network-based approaches to human disease, for instance to identify coherent transcriptional modules (subnetworks) that can inform functional disease mechanisms and pathological pathways. Yet, genome-scale network analysis across conditions is significantly hampered by the paucity of robust and computationally-efficient methods. Building on the Higher-Order Generalized Singular Value Decomposition, we introduce a new algorithmic approach for efficient, parameter-free and reproducible identification of network-modules simultaneously across multiple conditions. Our method can accommodate weighted (and unweighted) networks of any size and can similarly use co-expression or raw gene expression input data, without hinging upon the definition and stability of the correlation used to assess gene co-expression. In simulation studies, we demonstrated distinctive advantages of our method over existing methods, which was able to recover accurately both common and condition-specific network-modules without entailing ad-hoc input parameters as required by other approaches. We applied our method to genome-scale and multi-tissue transcriptomic datasets from rats (microarray-based) and humans (mRNA-sequencing-based) and identified several common and tissue-specific subnetworks with functional significance, which were not detected by other methods. In humans we recapitulated the crosstalk between cell-cycle progression and cell-extracellular matrix interactions processes in ventricular zones during neocortex expansion and further, we uncovered pathways related to development of later cognitive functions in the cortical plate of the developing brain which were previously unappreciated. Analyses of seven rat tissues identified a multi-tissue subnetwork of co-expressed heat shock protein (Hsp) and cardiomyopathy genes (Bag3, Cryab, Kras, Emd, Plec), which was significantly replicated using separate failing heart and liver gene expression datasets in humans, thus revealing a conserved functional role for Hsp genes in cardiovascular disease.

Transcriptional network analysis for the regulation of left ventricular hypertrophy and microvascular remodeling

Hypertension and cardiomyopathies share maladaptive changes of cardiac morphology, eventually leading to heart failure. These include left ventricular hypertrophy (LVH), myocardial fibrosis, and structural remodeling of coronary microcirculation, which is the morphologic hallmark of coronary microvascular dysfunction. To pinpoint the complex molecular mechanisms and pathways underlying LVH-associated cardiac remodeling independent of blood pressure effects, we employed gene network approaches to the rat heart. We used the Spontaneously Hypertensive Rat model showing many features of human hypertensive cardiomyopathy, for which we collected histological and histomorphometric data of the heart and coronary vasculature, and genome-wide cardiac gene expression. Here, we provide a large catalogue of gene co-expression networks in the heart that are significantly associated with quantitative variation in LVH, microvascular remodeling, and fibrosis-related traits. Many of these networks were significantly conserved to human idiopathic and/or ischemic cardiomyopathy patients, suggesting a potential role for these co-expressed genes in human heart disease.