Expression Atlas update: gene and protein expression in multiple species.

The EMBL-EBI Expression Atlas is an added value knowledge base that enables researchers to answer the question of where (tissue, organism part, developmental stage, cell type) and under which conditions (disease, treatment, gender, etc) a gene or protein of interest is expressed. Expression Atlas brings together data from >4500 expression studies from >65 different species, across different conditions and tissues. It makes these data freely available in an easy to visualise form, after expert curation to accurately represent the intended experimental design, re-analysed via standardised pipelines that rely on open-source community developed tools. Each study's metadata are annotated using ontologies. The data are re-analyzed with the aim of reproducing the original conclusions of the underlying experiments. Expression Atlas is currently divided into Bulk Expression Atlas and Single Cell Expression Atlas. Expression Atlas contains data from differential studies (microarray and bulk RNA-Seq) and baseline studies (bulk RNA-Seq and proteomics), whereas Single Cell Expression Atlas is currently dedicated to Single Cell RNA-Sequencing (scRNA-Seq) studies. The resource has been in continuous development since 2009 and it is available at https://www.ebi.ac.uk/gxa.

A compendium of uniformly processed human gene expression and splicing quantitative trait loci.

Many gene expression quantitative trait locus (eQTL) studies have published their summary statistics, which can be used to gain insight into complex human traits by downstream analyses, such as fine mapping and co-localization. However, technical differences between these datasets are a barrier to their widespread use. Consequently, target genes for most genome-wide association study (GWAS) signals have still not been identified. In the present study, we present the eQTL Catalogue ( https://www.ebi.ac.uk/eqtl ), a resource of quality-controlled, uniformly re-computed gene expression and splicing QTLs from 21 studies. We find that, for matching cell types and tissues, the eQTL effect sizes are highly reproducible between studies. Although most QTLs were shared between most bulk tissues, we identified a greater diversity of cell-type-specific QTLs from purified cell types, a subset of which also manifested as new disease co-localizations. Our summary statistics are freely available to enable the systematic interpretation of human GWAS associations across many cell types and tissues.

A single-cell RNA-sequencing training and analysis suite using the Galaxy framework.

BACKGROUND: The vast ecosystem of single-cell RNA-sequencing tools has until recently been plagued by an excess of diverging analysis strategies, inconsistent file formats, and compatibility issues between different software suites. The uptake of 10x Genomics datasets has begun to calm this diversity, and the bioinformatics community leans once more towards the large computing requirements and the statistically driven methods needed to process and understand these ever-growing datasets. RESULTS: Here we outline several Galaxy workflows and learning resources for single-cell RNA-sequencing, with the aim of providing a comprehensive analysis environment paired with a thorough user learning experience that bridges the knowledge gap between the computational methods and the underlying cell biology. The Galaxy reproducible bioinformatics framework provides tools, workflows, and trainings that not only enable users to perform 1-click 10x preprocessing but also empower them to demultiplex raw sequencing from custom tagged and full-length sequencing protocols. The downstream analysis supports a range of high-quality interoperable suites separated into common stages of analysis: inspection, filtering, normalization, confounder removal, and clustering. The teaching resources cover concepts from computer science to cell biology. Access to all resources is provided at the singlecell.usegalaxy.eu portal. CONCLUSIONS: The reproducible and training-oriented Galaxy framework provides a sustainable high-performance computing environment for users to run flexible analyses on both 10x and alternative platforms. The tutorials from the Galaxy Training Network along with the frequent training workshops hosted by the Galaxy community provide a means for users to learn, publish, and teach single-cell RNA-sequencing analysis.

Constitutively Active SMAD2/3 Are Broad-Scope Potentiators of Transcription-Factor-Mediated Cellular Reprogramming.

Reprogramming of cellular identity using exogenous expression of transcription factors (TFs) is a powerful and exciting tool for tissue engineering, disease modeling, and regenerative medicine. However, generation of desired cell types using this approach is often plagued by inefficiency, slow conversion, and an inability to produce mature functional cells. Here, we show that expression of constitutively active SMAD2/3 significantly improves the efficiency of induced pluripotent stem cell (iPSC) generation by the Yamanaka factors. Mechanistically, SMAD3 interacts with reprogramming factors and co-activators and co-occupies OCT4 target loci during reprogramming. Unexpectedly, active SMAD2/3 also markedly enhances three other TF-mediated direct reprogramming conversions, from B cells to macrophages, myoblasts to adipocytes, and human fibroblasts to neurons, highlighting broad and general roles for SMAD2/3 as cell-reprogramming potentiators. Our results suggest that co-expression of active SMAD2/3 could enhance multiple types of TF-based cell identity conversion and therefore be a powerful tool for cellular engineering.

Urinary peptidomics in a rodent model of diabetic nephropathy highlights epidermal growth factor as a biomarker for renal deterioration in patients with type 2 diabetes.

Many diabetic patients suffer from declining renal function without developing albuminuria. To identify alternative biomarkers for diabetic nephropathy (DN) we performed urinary peptidomic analysis in a rodent model in which hyperglycemia and hypertension synergize to promote renal pathologic changes consistent with human DN. We identified 297 increased and 15 decreased peptides in the urine of rats with DN compared with controls, including peptides derived from proteins associated with DN and novel candidate biomarkers. We confirmed by ELISA that one of the parent proteins, urinary epidermal growth factor (uEGF), was more than 2-fold reduced in rats with DN in comparison with controls. To assess the clinical utility of uEGF we examined renal outcomes in 642 participants from the Edinburgh Type 2 Diabetes Study who were normoalbuminuric and had preserved renal function at baseline. After adjustment for established renal risk factors, a lower uEGF to creatinine ratio was associated with new-onset estimated glomerular filtration rate less than 60 ml/min per 1.73m(2) (odds ratio 0.48; 95% confidence interval, 0.26-0.90), rapid (over 5% per annum) decline in renal function (odds ratio 0.44; 95% confidence interval, 0.27-0.72) or the composite of both outcomes (odds ratio 0.38; 95% confidence interval, 0.24-0.62). Thus, the utility of a low uEGF to creatinine ratio as a biomarker of progressive decline in renal function in normoalbuminuric patients should be assessed in additional populations.

Glucocorticoids accelerate maturation of the heme pathway in fetal liver through effects on transcription and DNA methylation.

Glucocorticoids are widely used in threatened preterm labor to promote maturation in many organ systems in preterm babies and have significant beneficial effects on morbidity and mortality. We performed transcriptional profiling in fetal liver in a rat model of prenatal glucocorticoid exposure and identified marked gene expression changes in heme biosynthesis, utilization, and degradation pathways in late gestation. These changes in gene expression associated with alterations in DNA methylation and with a reduction in hepatic heme concentration. There were no persistent differences in gene expression, DNA methylation, or heme concentrations at 4 weeks of age, suggesting that these are transient effects. Our findings are consistent with glucocorticoid-induced accelerated maturation of the haematopoietic system and support the hypothesis that glucocorticoids can drive changes in gene expression in association with alterations in DNA methylation.

Fibrinogen production is enhanced in an in-vitro model of non-alcoholic fatty liver disease: an isolated risk factor for cardiovascular events?

BACKGROUND: Cardiovascular disease (CVD) remains the major cause of excess mortality in patients with non-alcoholic fatty liver disease (NAFLD). The aim of this study was to investigate the individual contribution of NAFLD to CVD risk factors in the absence of pathogenic influences from other comorbidities often found in NAFLD patients, by using an established in-vitro model of hepatic steatosis. METHODS: Histopathological events in non-alcoholic fatty liver disease were recapitulated by focused metabolic nutrient overload of hepatoblastoma C3A cells, using oleate-treated-cells and untreated controls for comparison. Microarray and proteomic data from cell culture experiments were integrated into a custom-built systems biology database and proteogenomics analysis performed. Candidate genes with significant dysregulation and concomitant changes in protein abundance were identified and STRING association and enrichment analysis performed to identify putative pathogenic pathways. RESULTS: The search strategy yielded 3 candidate genes that were specifically and significantly up-regulated in nutrient-overloaded cells compared to untreated controls: fibrinogen alpha chain (2.2 fold), fibrinogen beta chain (2.3 fold) and fibrinogen gamma chain (2.1 fold) (all rank products pfp <0.05). Fibrinogen alpha and gamma chain also demonstrated significant concomitant increases in protein abundance (3.8-fold and 2.0-fold, respectively, p <0.05). CONCLUSIONS: In-vitro modelling of NAFLD and reactive oxygen species formation in nutrient overloaded C3A cells, in the absence of pathogenic influences from other comorbidities, suggests that NAFLD is an isolated determinant of CVD. Nutrient overload-induced up-regulation of all three fibrinogen component subunits of the coagulation cascade provides a possible mechanism to explain the excess CVD mortality observed in NAFLD patients.

Genome-wide mRNA expression profiling in vastus lateralis of COPD patients with low and normal fat free mass index and healthy controls.

BACKGROUND: Chronic Obstructive Pulmonary Disease (COPD) has significant systemic effects beyond the lungs amongst which muscle wasting is a prominent contributor to exercise limitation and an independent predictor of morbidity and mortality. The molecular mechanisms leading to skeletal muscle dysfunction/wasting are not fully understood and are likely to be multi-factorial. The need to develop therapeutic strategies aimed at improving skeletal muscle dysfunction/wasting requires a better understanding of the molecular mechanisms responsible for these abnormalities. Microarrays are powerful tools that allow the investigation of the expression of thousands of genes, virtually the whole genome, simultaneously. We aim at identifying genes and molecular pathways involved in skeletal muscle wasting in COPD. METHODS: We assessed and compared the vastus lateralis transcriptome of COPD patients with low fat free mass index (FFMI) as a surrogate of muscle mass (COPDL) (FEV1 30 ± 3.6%pred, FFMI 15 ± 0.2 Kg.m(-2)) with patients with COPD and normal FFMI (COPDN) (FEV1 44 ± 5.8%pred, FFMI 19 ± 0.5 Kg.m(-2)) and a group of age and sex matched healthy controls (C) (FEV1 95 ± 3.9%pred, FFMI 20 ± 0.8 Kg.m(-2)) using Agilent Human Whole Genome 4x44K microarrays. The altered expression of several of these genes was confirmed by real time TaqMan PCR. Protein levels of P21 were assessed by immunoblotting. RESULTS: A subset of 42 genes was differentially expressed in COPDL in comparison to both COPDN and C (PFP < 0.05; -1.5 ≥ FC ≥ 1.5). The altered expression of several of these genes was confirmed by real time TaqMan PCR and correlated with different functional and structural muscle parameters. Five of these genes (CDKN1A, GADD45A, PMP22, BEX2, CGREF1, CYR61), were associated with cell cycle arrest and growth regulation and had been previously identified in studies relating muscle wasting and ageing. Protein levels of CDKN1A, a recognized marker of premature ageing/cell cycle arrest, were also found to be increased in COPDL. CONCLUSIONS: This study provides evidence of differentially expressed genes in peripheral muscle in COPD patients corresponding to relevant biological processes associated with skeletal muscle wasting and provides potential targets for future therapeutic interventions to prevent loss of muscle function and mass in COPD.

Tight blood glycaemic and blood pressure control in experimental diabetic nephropathy reduces extracellular matrix production without regression of fibrosis.

AIMS: Regression of albuminuria and renal fibrosis occurs in patients with diabetic nephropathy (DN) following tight control of blood glucose and blood pressure, however the pathways that promote regression remain poorly understood and we wished to characterize these using a rodent model. METHODS: Diabetes was induced with streptozotocin in Cyp1a1mRen2 rats and hypertension was generated by inducing renin transgene expression with dietary indole-3-carbinol (I-3-C) for 28 weeks. At this point an 'injury cohort' was culled, while in a 'reversal cohort' glycaemia was tightly controlled using insulin implants and blood pressure normalized by withdrawing dietary I-3-C for a further 8 weeks. Pathways activated during and following reversal of diabetes and hypertension were assessed by microarray profiling. RESULTS: Tight control of blood glucose and blood pressure reduced albuminuria and renal hypertrophy, but had no impact on renal fibrosis. 85 genes were up-regulated specifically during the injury phase, including genes encoding multiple myofibroblast and extracellular matrix (ECM) proteins. Conversely, 314 genes remained persistently elevated during reversal including genes linked to innate/adaptive immunity, phagocytosis, lysosomal processing and degradative metalloproteinases (MMPs). Despite increased MMP gene expression, MMP activity was suppressed during both injury and reversal, in association with up-regulation of tissue inhibitor of metalloproteinase-1 (TIMP-1) protein. Physical separation of the TIMP-1/MMP complexes during zymography of tissue homogenate restored MMP activity. CONCLUSION: Normalization of blood glucose and pressure ameliorates albuminuria and inhibits excess ECM production, however persistent TIMP-1 expression hinders attempts at ECM remodelling. Therapies which counteract the action of TIMPs may accelerate scar resolution.

The value of immunophenotyping hepatocellular adenomas: consecutive resections at one UK centre.

AIMS: To determine the utility of immunophenotyping for classification of hepatocellular adenomas resected at one Scottish centre. METHODS AND RESULTS: This study comprised a retrospective review and immunophenotyping of consecutive resected benign hepatocellular tumours. Fifty-five patients (seven men) had 64 adenomas and 26 focal nodular hyperplasias (FNHs) resected. Map-like glutamine synthetase (GS) staining was specific for FNH. Immunophenotyping changed the morphological typing for three adenomas and resolved 16 of 18 unclassified or equivocal cases, revealing GS positivity in these (seven) and four others. Steatotic/liver fatty acid binding protein-deficient adenomas were the commonest type in women (12/29 women, 41%) but were absent from men. Where one of multiple adenomas was morphologically unclassified, there was still a shared immunophenotype. Diffuse CD34 positivity correlated with GS positivity or unclassified status (P < 0.0001). Supervised cluster analysis identified morphological discriminants for FNH and predictors of adenoma type and their insensitivity in predicting GS status. Forty per cent of men and 7% of women with adenomas had a specific adenoma risk, including danazol and portal venopathies. Inflammatory adenomas were associated with metabolic syndrome, steatosis, or alcohol (P = 0.053). Four patients showed carcinoma ex-adenoma. CONCLUSIONS: The distribution of adenoma types in this population matches that in others, and immunoprofiling is required for accurate typing. Carcinoma ex-adenoma is uncommon and fits the published risk profile (large size and GS-positive).

Differential Ly-6C expression identifies the recruited macrophage phenotype, which orchestrates the regression of murine liver fibrosis.

Although macrophages are widely recognized to have a profibrotic role in inflammation, we have used a highly tractable CCl(4)-induced model of reversible hepatic fibrosis to identify and characterize the macrophage phenotype responsible for tissue remodeling: the hitherto elusive restorative macrophage. This CD11B(hi) F4/80(int) Ly-6C(lo) macrophage subset was most abundant in livers during maximal fibrosis resolution and represented the principle matrix metalloproteinase (MMP) -expressing subset. Depletion of this population in CD11B promoter-diphtheria toxin receptor (CD11B-DTR) transgenic mice caused a failure of scar remodeling. Adoptive transfer and in situ labeling experiments showed that these restorative macrophages derive from recruited Ly-6C(hi) monocytes, a common origin with profibrotic Ly-6C(hi) macrophages, indicative of a phenotypic switch in vivo conferring proresolution properties. Microarray profiling of the Ly-6C(lo) subset, compared with Ly-6C(hi) macrophages, showed a phenotype outside the M1/M2 classification, with increased expression of MMPs, growth factors, and phagocytosis-related genes, including Mmp9, Mmp12, insulin-like growth factor 1 (Igf1), and Glycoprotein (transmembrane) nmb (Gpnmb). Confocal microscopy confirmed the postphagocytic nature of restorative macrophages. Furthermore, the restorative macrophage phenotype was recapitulated in vitro by the phagocytosis of cellular debris with associated activation of the ERK signaling cascade. Critically, induced phagocytic behavior in vivo, through administration of liposomes, increased restorative macrophage number and accelerated fibrosis resolution, offering a therapeutic strategy to this orphan pathological process.

Hyperglycemia and renin-dependent hypertension synergize to model diabetic nephropathy.

Rodent models exhibit only the earliest features of human diabetic nephropathy, which limits our ability to investigate new therapies. Hypertension is a prerequisite for advanced diabetic nephropathy in humans, so its rarity in typical rodent models may partly explain their resistance to nephropathy. Here, we used the Cyp1a1mRen2 rat, in which the murine renin-2 gene is incorporated under the Cytochrome P4501a1 promoter. In this transgenic strain, administration of low-dose dietary indole-3-carbinol induces moderate hypertension. In the absence of hypertension, streptozotocin-induced diabetes resulted in a 14-fold increase in albuminuria but only mild changes in histology and gene expression despite 28 weeks of marked hyperglycemia. In the presence of induced hypertension, hyperglycemia resulted in a 500-fold increase in albuminuria, marked glomerulosclerosis and tubulointerstitial fibrosis, and induction of many of the same pathways that are upregulated in the tubulointerstitium in human diabetic nephropathy. In conclusion, although induction of diabetes alone in rodents has limited utility to model human diabetic nephropathy, renin-dependent hypertension and hyperglycemia synergize to recapitulate many of the clinical, histological, and gene expression changes observed in humans.

Underestimated and under-recognized: the late consequences of acute coronary syndrome (GRACE UK-Belgian Study).

AIM: To define the long-term outcome of patients presenting with acute coronary syndrome [ST-segment elevation myocardial infarction (STEMI), and non-STEMI and unstable angina acute coronary syndrome (ACS) without biomarker elevation] and to test the hypothesis that the GRACE (Global Registry of Acute Coronary Events) risk score predicts mortality and death/MI at 5 years. METHODS AND RESULTS: In the GRACE long-term study, UK and Belgian centres prospectively recruited and followed ACS patients for a median of 5 years (1797 days). Primary outcome events: deaths, cardiovascular deaths (CVDs) and MIs. Secondary events: stroke and re-hospitalization for ACS. There were 736 deaths, 19.8% (482 CVDs, 13%) and 347 (9.3%) MIs (>24 h), 261 strokes (7.7%), and 452 (17%) subsequent revascularizations. Rehospitalization was common: average 1.6 per patient; 31.2% had >1 admission, 9.2% had 5+ admissions. These events were despite high rates of guideline indicated therapies. The GRACE score was highly predictive of all-cause death, CVD, and CVD/MI at 5 years (death: χ(2) likelihood ratio 632; Wald 709.9, P< 0.0001, C-statistic 0.77; for CVD C-statistic 0.75, P < 0.0001; CVD/MI C-statistic 0.70, P < 0.0001). Compared with the low-risk GRACE stratum (ESC Guideline criteria), those with intermediate [hazard ratio (HR) 2.14, 95% CI 1.63, 2.81] and those with high-risk (HR 6.36, 95% CI 4.95, 8.16) had two- and six-fold higher risk of later death (Cox proportional hazard). A landmark analysis after 6 months confirmed that the GRACE score predicted long-term death (χ(2) likelihood ratio 265.4; Wald 289.5, P < 0.0001). Although in-hospital rates of death and MI are higher following STEMI, the cumulative rates of death (and CVD) were not different, by class of ACS, over the duration of follow-up (Wilcoxon = 1.5597, df = 1, P = 0.21). At 5 years after STEMI 269/1403 (19%) died; after non-STEMI 262/1170 (22%) after unstable angina (UA) 149/850 (17%). Two-thirds (68%) of STEMI deaths occurred after initial hospital discharge, but this was 86% for non-STEMI and 97% for UA. CONCLUSION: The GRACE risk score predicts early and 5 year death and CVD/MI. Five year morbidity and mortality are as high in patients following non-ST MI and UA as seen following STEMI. Their morbidity burden is high (MI, stroke, readmissions) and the substantial late mortality in non-STE ACS is under-recognized. The findings highlight the importance of pursuing novel approaches to diminish long-term risk.

In silico structure-function analysis of pathological variation in the HSD11B2 gene sequence.

11beta-Hydroxysteroid dehydrogenase type 2 (11betaHSD2) is a short-chain dehydrogenase/reductase (SDR) responsible for inactivating cortisol and preventing its binding to the mineralocorticoid receptor (MR). Nonfunctional mutations in HSD11B2, the gene encoding 11betaHSD2, cause the hypertensive syndrome of apparent mineralocorticoid excess (AME). Like other such Mendelian disorders, AME is rare but has nevertheless helped to illuminate principles fundamental to the regulation of blood pressure. Furthermore, polymorphisms in HSD11B2 have been associated with salt sensitivity, a major risk factor for cardiovascular mortality. It is therefore highly likely that sequence variation in HSD11B2, having subtle functional ramifications, will affect blood pressure in the wider population. In this study, a three-dimensional homology model of 11betaHSD2 was created and used to hypothesize the functional consequences in terms of protein structure of published mutations in HSD11B2. This approach underscored the strong genotype-phenotype correlation of AME: severe forms of the disease, associated with little in vivo enzyme activity, arise from mutations occurring in invariant alignment positions. These were predicted to exert gross structural changes in the protein. In contrast, those mutations causing a mild clinical phenotype were in less conserved regions of the protein that were predicted to be relatively more tolerant to substitution. Finally, a number of pathogenic mutations are shown to be associated with regions predicted to participate in dimer formation, and in protein stabilization, which may therefore suggest molecular mechanisms of disease.

Automated seeding of specialised wiki knowledgebases with BioKb.

BACKGROUND: Wiki technology has become a ubiquitous mechanism for dissemination of information, and places strong emphasis on collaboration. We aimed to leverage wiki technology to allow small groups of researchers to collaborate around a specific domain, for example a biological pathway. Automatically gathered seed data could be modified by the group and enriched with domain specific information. RESULTS: We describe a software system, BioKb, implemented as a plugin for the TWiki engine, and designed to facilitate construction of a field-specific wiki containing collaborative and automatically generated content. Features of this system include: query of publicly available resources such as KEGG, iHOP and MeSH, to generate 'seed' content for topics; simple definition of structure for topics of different types via an administration page; and interactive incorporation of relevant PubMed references. An exemplar is shown for the use of this system, in the creation of the RAASWiki knowledgebase on the renin-angiotensin-aldosterone system (RAAS). RAASWiki has been seeded with data by use of BioKb, and will be the subject of ongoing development into an extensive knowledgebase on the RAAS. CONCLUSION: The BioKb system is available from http://www.bioinf.mvm.ed.ac.uk/twiki/bin/view/TWiki/BioKbPlugin as a plugin for the TWiki engine.

The contrasting properties of conservation and correlated phylogeny in protein functional residue prediction.

BACKGROUND: Amino acids responsible for structure, core function or specificity may be inferred from multiple protein sequence alignments where a limited set of residue types are tolerated. The rise in available protein sequences continues to increase the power of techniques based on this principle. RESULTS: A new algorithm, SMERFS, for predicting protein functional sites from multiple sequences alignments was compared to 14 conservation measures and to the MINER algorithm. Validation was performed on an automatically generated dataset of 1457 families derived from the protein interactions database SNAPPI-DB, and a smaller manually curated set of 148 families. The best performing measure overall was Williamson property entropy, with ROC0.1 scores of 0.0087 and 0.0114 for domain and small molecule contact prediction, respectively. The Lancet method performed worse than random on protein-protein interaction site prediction (ROC0.1 score of 0.0008). The SMERFS algorithm gave similar accuracy to the phylogenetic tree-based MINER algorithm but was superior to Williamson in prediction of non-catalytic transient complex interfaces. SMERFS predicts sites that are significantly more solvent accessible compared to Williamson. CONCLUSION: Williamson property entropy is the the best performing of 14 conservation measures examined. The difference in performance of SMERFS relative to Williamson in manually defined complexes was dependent on complex type. The best choice of analysis method is therefore dependent on the system of interest. Additional computation employed by Miner in calculation of phylogenetic trees did not produce improved results over SMERFS. SMERFS performance was improved by use of windows over alignment columns, illustrating the necessity of considering the local environment of positions when assessing their functional significance.