Improved risk stratification and prognosis prediction in sepsis is a critical unmet need. Clinical severity scores and available assays such as blood lactate reflect global illness severity with ...suboptimal performance, and do not specifically reveal the underlying dysregulation of sepsis. Here, we present prognostic models for 30-day mortality generated independently by three scientific groups by using 12 discovery cohorts containing transcriptomic data collected from primarily community-onset sepsis patients. Predictive performance is validated in five cohorts of community-onset sepsis patients in which the models show summary AUROCs ranging from 0.765-0.89. Similar performance is observed in four cohorts of hospital-acquired sepsis. Combining the new gene-expression-based prognostic models with prior clinical severity scores leads to significant improvement in prediction of 30-day mortality as measured via AUROC and net reclassification improvement index These models provide an opportunity to develop molecular bedside tests that may improve risk stratification and mortality prediction in patients with sepsis.
Effective targeted therapy for sepsis requires an understanding of the heterogeneity in the individual host response to infection. We investigated this heterogeneity by defining interindividual ...variation in the transcriptome of patients with sepsis and related this to outcome and genetic diversity.
We assayed peripheral blood leucocyte global gene expression for a prospective discovery cohort of 265 adult patients admitted to UK intensive care units with sepsis due to community-acquired pneumonia and evidence of organ dysfunction. We then validated our findings in a replication cohort consisting of a further 106 patients. We mapped genomic determinants of variation in gene transcription between patients as expression quantitative trait loci (eQTL).
We discovered that following admission to intensive care, transcriptomic analysis of peripheral blood leucocytes defines two distinct sepsis response signatures (SRS1 and SRS2). The presence of SRS1 (detected in 108 41% patients in discovery cohort) identifies individuals with an immunosuppressed phenotype that included features of endotoxin tolerance, T-cell exhaustion, and downregulation of human leucocyte antigen (HLA) class II. SRS1 was associated with higher 14 day mortality than was SRS2 (discovery cohort hazard ratio (HR) 2·4, 95% CI 1·3–4·5, p=0·005; validation cohort HR 2·8, 95% CI 1·5–5·1, p=0·0007). We found that a predictive set of seven genes enabled the classification of patients as SRS1 or SRS2. We identified cis-acting and trans-acting eQTL for key immune and metabolic response genes and sepsis response networks. Sepsis eQTL were enriched in endotoxin-induced epigenetic marks and modulated the individual host response to sepsis, including effects specific to SRS group. We identified regulatory genetic variants involving key mediators of gene networks implicated in the hypoxic response and the switch to glycolysis that occurs in sepsis, including HIF1α and mTOR, and mediators of endotoxin tolerance, T-cell activation, and viral defence.
Our integrated genomics approach advances understanding of heterogeneity in sepsis by defining subgroups of patients with different immune response states and prognoses, as well as revealing the role of underlying genetic variation. Our findings provide new insights into the pathogenesis of sepsis and create opportunities for a precision medicine approach to enable targeted therapeutic intervention to improve sepsis outcomes.
European Commission, Medical Research Council (UK), and the Wellcome Trust.
Although alterations in myeloid cells have been observed in COVID-19, the specific underlying mechanisms are not completely understood. Here, we examine the function of classical CD14
monocytes in ...patients with mild and moderate COVID-19 during the acute phase of infection and in healthy individuals. Monocytes from COVID-19 patients display altered expression of cell surface receptors and a dysfunctional metabolic profile that distinguish them from healthy monocytes. Secondary pathogen sensing ex vivo leads to defects in pro-inflammatory cytokine and type-I IFN production in moderate COVID-19 cases, together with defects in glycolysis. COVID-19 monocytes switch their gene expression profile from canonical innate immune to pro-thrombotic signatures and are functionally pro-thrombotic, both at baseline and following ex vivo stimulation with SARS-CoV-2. Transcriptionally, COVID-19 monocytes are characterized by enrichment of pathways involved in hemostasis, immunothrombosis, platelet aggregation and other accessory pathways to platelet activation and clot formation. These results identify a potential mechanism by which monocyte dysfunction may contribute to COVID-19 pathology.
Despite significant progress in annotating the genome with experimental methods, much of the regulatory noncoding genome remains poorly defined. Here we assert that regulatory elements may be ...characterized by leveraging local epigenomic signatures where specific transcription factors (TFs) are bound. To link these two features, we introduce IMPACT, a genome annotation strategy that identifies regulatory elements defined by cell-state-specific TF binding profiles, learned from 515 chromatin and sequence annotations. We validate IMPACT using multiple compelling applications. First, IMPACT distinguishes between bound and unbound TF motif sites with high accuracy (average AUPRC 0.81, SE 0.07; across 8 tested TFs) and outperforms state-of-the-art TF binding prediction methods, MocapG, MocapS, and Virtual ChIP-seq. Second, in eight tested cell types, RNA polymerase II IMPACT annotations capture more cis-eQTL variation than sequence-based annotations, such as promoters and TSS windows (25% average increase in enrichment). Third, integration with rheumatoid arthritis (RA) summary statistics from European (N = 38,242) and East Asian (N = 22,515) populations revealed that the top 5% of CD4+ Treg IMPACT regulatory elements capture 85.7% of RA h2, the most comprehensive explanation for RA h2 to date. In comparison, the average RA h2 captured by compared CD4+ T histone marks is 42.3% and by CD4+ T specifically expressed gene sets is 36.4%. Lastly, we find that IMPACT may be used in many different cell types to identify complex trait associated regulatory elements.
Sepsis arises from diverse and incompletely understood dysregulated host response processes following infection that leads to life-threatening organ dysfunction. Here we showed that neutrophils and ...emergency granulopoiesis drove a maladaptive response during sepsis. We generated a whole-blood single-cell multiomic atlas (272,993 cells, n = 39 individuals) of the sepsis immune response that identified populations of immunosuppressive mature and immature neutrophils. In co-culture, CD66b
sepsis neutrophils inhibited proliferation and activation of CD4
T cells. Single-cell multiomic mapping of circulating hematopoietic stem and progenitor cells (HSPCs) (29,366 cells, n = 27) indicated altered granulopoiesis in patients with sepsis. These features were enriched in a patient subset with poor outcome and a specific sepsis response signature that displayed higher frequencies of IL1R2
immature neutrophils, epigenetic and transcriptomic signatures of emergency granulopoiesis in HSPCs and STAT3-mediated gene regulation across different infectious etiologies and syndromes. Our findings offer potential therapeutic targets and opportunities for stratified medicine in severe infection.
Heterogeneity in the septic response has hindered efforts to understand pathophysiology and develop targeted therapies. Source of infection, with different causative organisms and temporal changes, ...might influence this heterogeneity.
To investigate individual and temporal variations in the transcriptomic response to sepsis due to fecal peritonitis, and to compare these with the same parameters in community-acquired pneumonia.
We performed genome-wide gene expression profiling in peripheral blood leukocytes of adult patients admitted to intensive care with sepsis due to fecal peritonitis (n = 117) or community-acquired pneumonia (n = 126), and of control subjects without sepsis (n = 10).
A substantial portion of the transcribed genome (18%) was differentially expressed compared with that of control subjects, independent of source of infection, with eukaryotic initiation factor 2 signaling being the most enriched canonical pathway. We identified two sepsis response signature (SRS) subgroups in fecal peritonitis associated with early mortality (P = 0.01; hazard ratio, 4.78). We defined gene sets predictive of SRS group, and serial sampling demonstrated that subgroup membership is dynamic during intensive care unit admission. We found that SRS is the major predictor of transcriptomic variation; a small number of genes (n = 263) were differentially regulated according to the source of infection, enriched for IFN signaling and antigen presentation. We define temporal changes in gene expression from disease onset involving phagosome formation as well as natural killer cell and IL-3 signaling.
The majority of the sepsis transcriptomic response is independent of the source of infection and includes signatures reflecting immune response state and prognosis. A modest number of genes show evidence of specificity. Our findings highlight opportunities for patient stratification and precision medicine in sepsis.
The apelin receptor, a G protein-coupled receptor, has emerged as a key regulator of cardiovascular development, physiology, and disease. However, there is a lack of suitable human in vitro models to ...investigate the apelinergic system in cardiovascular cell types. For the first time we have used human embryonic stem cell-derived cardiomyocytes (hESC-CMs) and a novel inducible knockdown system to examine the role of the apelin receptor in both cardiomyocyte development and to determine the consequences of loss of apelin receptor function as a model of disease.
Expression of the apelin receptor and its ligands in hESCs and hESC-CMs was determined. hESCs carrying a tetracycline-inducible short hairpin RNA targeting the apelin receptor were generated using the sOPTiKD system. Phenotypic assays characterized the consequences of either apelin receptor knockdown before hESC-CM differentiation (early knockdown) or in 3D engineered heart tissues as a disease model (late knockdown). hESC-CMs expressed the apelin signalling system at a similar level to the adult heart. Early apelin receptor knockdown decreased cardiomyocyte differentiation efficiency and prolonged voltage sensing, associated with asynchronous contraction. Late apelin receptor knockdown had detrimental consequences on 3D engineered heart tissue contractile properties, decreasing contractility and increasing stiffness.
We have successfully knocked down the apelin receptor, using an inducible system, to demonstrate a key role in hESC-CM differentiation. Knockdown in 3D engineered heart tissues recapitulated the phenotype of apelin receptor down-regulation in a failing heart, providing a potential platform for modelling heart failure and testing novel therapeutic strategies.
Abstract High blood pressure in the portal vein, portal hypertension (PH), is the final common pathway in liver cirrhosis regardless of aetiology. Complications from PH are the major cause of ...morbidity and mortality in these patients. Current drug therapy to reduce portal pressure is mainly limited to β-adrenergic receptor blockade but approximately 40% of patients do not respond. Our aim was to use microarray to measure the expression of ∼20,800 genes in portal vein from patients with PH undergoing transplantation for liver cirrhosis (PH, n=12) versus healthy vessels (control, n=9) to identify potential drug targets to improve therapy. Expression of 9,964 genes above background was detected in portal vein samples. Comparing PH veins versus control (adjusted P-value < 0.05, fold change > 1.5) identified 548 up-regulated genes and 1,996 down-regulated genes. The 2,544 differentially expressed genes were subjected to pathway analysis. We identified 49 significantly enriched pathways. The endothelin pathway was ranked the tenth most significant, the only vasoconstrictive pathway to be identified. ET-1 gene (EDN1) was significantly up-regulated, consistent with elevated levels of ET-1 peptide previously measured in PH and cirrhosis. ETA receptor gene (EDNRA) was significantly down-regulated, consistent with an adaptive response to increased peptide levels in the portal vein but there was no change in the ETB gene (EDNRB). The results provide further support for evaluating the efficacy of ETA receptor antagonists as a potential therapy in addition to β-blockers in patients with PH and cirrhosis.
Understanding how genetic variants influence disease risk and complex traits (variant-to-function) is one of the major challenges in human genetics. Here we present a model-driven framework to ...leverage human genome-scale metabolic networks to define how genetic variants affect biochemical reaction fluxes across major human tissues, including skeletal muscle, adipose, liver, brain and heart. As proof of concept, we build personalised organ-specific metabolic flux models for 524,615 individuals of the INTERVAL and UK Biobank cohorts and perform a fluxome-wide association study (FWAS) to identify 4312 associations between personalised flux values and the concentration of metabolites in blood. Furthermore, we apply FWAS to identify 92 metabolic fluxes associated with the risk of developing coronary artery disease, many of which are linked to processes previously described to play in role in the disease. Our work demonstrates that genetically personalised metabolic models can elucidate the downstream effects of genetic variants on biochemical reactions involved in common human diseases.
The use of omic modalities to dissect the molecular underpinnings of common diseases and traits is becoming increasingly common. But multi-omic traits can be genetically predicted, which enables ...highly cost-effective and powerful analyses for studies that do not have multi-omics
. Here we examine a large cohort (the INTERVAL study
; n = 50,000 participants) with extensive multi-omic data for plasma proteomics (SomaScan, n = 3,175; Olink, n = 4,822), plasma metabolomics (Metabolon HD4, n = 8,153), serum metabolomics (Nightingale, n = 37,359) and whole-blood Illumina RNA sequencing (n = 4,136), and use machine learning to train genetic scores for 17,227 molecular traits, including 10,521 that reach Bonferroni-adjusted significance. We evaluate the performance of genetic scores through external validation across cohorts of individuals of European, Asian and African American ancestries. In addition, we show the utility of these multi-omic genetic scores by quantifying the genetic control of biological pathways and by generating a synthetic multi-omic dataset of the UK Biobank
to identify disease associations using a phenome-wide scan. We highlight a series of biological insights with regard to genetic mechanisms in metabolism and canonical pathway associations with disease; for example, JAK-STAT signalling and coronary atherosclerosis. Finally, we develop a portal ( https://www.omicspred.org/ ) to facilitate public access to all genetic scores and validation results, as well as to serve as a platform for future extensions and enhancements of multi-omic genetic scores.