Multiomic analyses uncover immunological signatures in acute and chronic coronary syndromes.

Acute and chronic coronary syndromes (ACS and CCS) are leading causes of mortality. Inflammation is considered a key pathogenic driver of these diseases, but the underlying immune states and their clinical implications remain poorly understood. Multiomic factor analysis (MOFA) allows unsupervised data exploration across multiple data types, identifying major axes of variation and associating these with underlying molecular processes. We hypothesized that applying MOFA to multiomic data obtained from blood might uncover hidden sources of variance and provide pathophysiological insights linked to clinical needs. Here we compile a longitudinal multiomic dataset of the systemic immune landscape in both ACS and CCS (n = 62 patients in total, n = 15 women and n = 47 men) and validate this in an external cohort (n = 55 patients in total, n = 11 women and n = 44 men). MOFA reveals multicellular immune signatures characterized by distinct monocyte, natural killer and T cell substates and immune-communication pathways that explain a large proportion of inter-patient variance. We also identify specific factors that reflect disease state or associate with treatment outcome in ACS as measured using left ventricular ejection fraction. Hence, this study provides proof-of-concept evidence for the ability of MOFA to uncover multicellular immune programs in cardiovascular disease, opening new directions for mechanistic, biomarker and therapeutic studies.

Resident and recruited macrophages differentially contribute to cardiac healing after myocardial ischemia.

Cardiac macrophages are heterogenous in phenotype and functions, which has been associated with differences in their ontogeny. Despite extensive research, our understanding of the precise role of different subsets of macrophages in ischemia/reperfusion (I/R) injury remains incomplete. We here investigated macrophage lineages and ablated tissue macrophages in homeostasis and after I/R injury in a CSF1R-dependent manner. Genomic deletion of a fms-intronic regulatory element (FIRE) in the Csf1r locus resulted in specific absence of resident homeostatic and antigen-presenting macrophages, without affecting the recruitment of monocyte-derived macrophages to the infarcted heart. Specific absence of homeostatic, monocyte-independent macrophages altered the immune cell crosstalk in response to injury and induced proinflammatory neutrophil polarization, resulting in impaired cardiac remodeling without influencing infarct size. In contrast, continuous CSF1R inhibition led to depletion of both resident and recruited macrophage populations. This augmented adverse remodeling after I/R and led to an increased infarct size and deterioration of cardiac function. In summary, resident macrophages orchestrate inflammatory responses improving cardiac remodeling, while recruited macrophages determine infarct size after I/R injury. These findings attribute distinct beneficial effects to different macrophage populations in the context of myocardial infarction.

Peripheral priming induces plastic transcriptomic and proteomic responses in circulating neutrophils required for pathogen containment.

Neutrophils rapidly respond to inflammation and infection, but to which degree their functional trajectories after mobilization from the bone marrow are shaped within the circulation remains vague. Experimental limitations have so far hampered neutrophil research in human disease. Here, using innovative fixation and single-cell-based toolsets, we profile human and murine neutrophil transcriptomes and proteomes during steady state and bacterial infection. We find that peripheral priming of circulating neutrophils leads to dynamic shifts dominated by conserved up-regulation of antimicrobial genes across neutrophil substates, facilitating pathogen containment. We show the TLR4/NF-κB signaling-dependent up-regulation of canonical neutrophil activation markers like CD177/NB-1 during acute inflammation, resulting in functional shifts in vivo. Blocking de novo RNA synthesis in circulating neutrophils abrogates these plastic shifts and prevents the adaptation of antibacterial neutrophil programs by up-regulation of distinct effector molecules upon infection. These data underline transcriptional plasticity as a relevant mechanism of functional neutrophil reprogramming during acute infection to foster bacterial containment within the circulation.

Augusta: From RNA-Seq to gene regulatory networks and Boolean models.

Computational models of gene regulations help to understand regulatory mechanisms and are extensively used in a wide range of areas, e.g., biotechnology or medicine, with significant benefits. Unfortunately, there are only a few computational gene regulatory models of whole genomes allowing static and dynamic analysis due to the lack of sophisticated tools for their reconstruction. Here, we describe Augusta, an open-source Python package for Gene Regulatory Network (GRN) and Boolean Network (BN) inference from the high-throughput gene expression data. Augusta can reconstruct genome-wide models suitable for static and dynamic analyses. Augusta uses a unique approach where the first estimation of a GRN inferred from expression data is further refined by predicting transcription factor binding motifs in promoters of regulated genes and by incorporating verified interactions obtained from databases. Moreover, a refined GRN is transformed into a draft BN by searching in the curated model database and setting logical rules to incoming edges of target genes, which can be further manually edited as the model is provided in the SBML file format. The approach is applicable even if information about the organism under study is not available in the databases, which is typically the case for non-model organisms including most microbes. Augusta can be operated from the command line and, thus, is easy to use for automated prediction of models for various genomes. The Augusta package is freely available at github.com/JanaMus/Augusta. Documentation and tutorials are available at augusta.readthedocs.io.

Deletion of the transcription factors Hsf1, Msn2 and Msn4 in yeast uncovers transcriptional reprogramming in response to proteotoxic stress.

The response to proteotoxic stresses such as heat shock allows organisms to maintain protein homeostasis under changing environmental conditions. We asked what happens if an organism can no longer react to cytosolic proteotoxic stress. To test this, we deleted or depleted, either individually or in combination, the stress-responsive transcription factors Msn2, Msn4, and Hsf1 in Saccharomyces cerevisiae. Our study reveals a combination of survival strategies, which together protect essential proteins. Msn2 and 4 broadly reprogram transcription, triggering the response to oxidative stress, as well as biosynthesis of the protective sugar trehalose and glycolytic enzymes, while Hsf1 mainly induces the synthesis of molecular chaperones and reverses the transcriptional response upon prolonged mild heat stress (adaptation).

Mural cell-derived chemokines provide a protective niche to safeguard vascular macrophages and limit chronic inflammation.

Maladaptive, non-resolving inflammation contributes to chronic inflammatory diseases such as atherosclerosis. Because macrophages remove necrotic cells, defective macrophage programs can promote chronic inflammation with persistent tissue injury. Here, we investigated the mechanisms sustaining vascular macrophages. Intravital imaging revealed a spatiotemporal macrophage niche across vascular beds alongside mural cells (MCs)-pericytes and smooth muscle cells. Single-cell transcriptomics, co-culture, and genetic deletion experiments revealed MC-derived expression of the chemokines CCL2 and MIF, which actively preserved macrophage survival and their homeostatic functions. In atherosclerosis, this positioned macrophages in viable plaque areas, away from the necrotic core, and maintained a homeostatic macrophage phenotype. Disruption of this MC-macrophage unit via MC-specific deletion of these chemokines triggered detrimental macrophage relocalizing, exacerbated plaque necrosis, inflammation, and atheroprogression. In line, CCL2 inhibition at advanced stages of atherosclerosis showed detrimental effects. This work presents a MC-driven safeguard toward maintaining the homeostatic vascular macrophage niche.

Human genetic defects in SRP19 and SRPRA cause severe congenital neutropenia with distinctive proteome changes.

The mechanisms of coordinated changes in proteome composition and their relevance for the differentiation of neutrophil granulocytes are not well studied. Here, we discover 2 novel human genetic defects in signal recognition particle receptor alpha (SRPRA) and SRP19, constituents of the mammalian cotranslational targeting machinery, and characterize their roles in neutrophil granulocyte differentiation. We systematically study the proteome of neutrophil granulocytes from patients with variants in the SRP genes, HAX1, and ELANE, and identify global as well as specific proteome aberrations. Using in vitro differentiation of human induced pluripotent stem cells and in vivo zebrafish models, we study the effects of SRP deficiency on neutrophil granulocyte development. In a heterologous cell-based inducible protein expression system, we validate the effects conferred by SRP dysfunction for selected proteins that we identified in our proteome screen. Thus, SRP-dependent protein processing, intracellular trafficking, and homeostasis are critically important for the differentiation of neutrophil granulocytes.

Analysis of proteome adaptation reveals a key role of the bacterial envelope in starvation survival.

Bacteria reorganize their physiology upon entry to stationary phase. What part of this reorganization improves starvation survival is a difficult question because the change in physiology includes a global reorganization of the proteome, envelope, and metabolism of the cell. In this work, we used several trade-offs between fast growth and long survival to statistically score over 2,000 Escherichia coli proteins for their global correlation with death rate. The combined ranking allowed us to narrow down the set of proteins that positively correlate with survival and validate the causal role of a subset of proteins. Remarkably, we found that important survival genes are related to the cell envelope, i.e., periplasm and outer membrane, because the maintenance of envelope integrity of E. coli plays a crucial role during starvation. Our results uncover a new protective feature of the outer membrane that adds to the growing evidence that the outer membrane is not only a barrier that prevents abiotic substances from reaching the cytoplasm but also essential for bacterial proliferation and survival.

The Amino Acid Homoarginine Inhibits Atherogenesis by Modulating T-Cell Function.

BACKGROUND: Amino acid metabolism is crucial for inflammatory processes during atherogenesis. The endogenous amino acid homoarginine is a robust biomarker for cardiovascular outcome and mortality with high levels being protective. However, the underlying mechanisms remain elusive. We investigated the effect of homoarginine supplementation on atherosclerotic plaque development with a particular focus on inflammation. METHODS: Female ApoE-deficient mice were supplemented with homoarginine (14 mg/L) in drinking water starting 2 weeks before and continuing throughout a 6-week period of Western-type diet feeding. Control mice received normal drinking water. Immunohistochemistry and flow cytometry were used for plaque- and immunological phenotyping. T cells were characterized using mass spectrometry-based proteomics, by functional in vitro approaches, for example, proliferation and migration/chemotaxis assays as well as by super-resolution microscopy. RESULTS: Homoarginine supplementation led to a 2-fold increase in circulating homoarginine concentrations. Homoarginine-treated mice exhibited reduced atherosclerosis in the aortic root and brachiocephalic trunk. A substantial decrease in CD3+ T cells in the atherosclerotic lesions suggested a T-cell-related effect of homoarginine supplementation, which was mainly attributed to CD4+ T cells. Macrophages, dendritic cells, and B cells were not affected. CD4+ T-cell proteomics and subsequent pathway analysis together with in vitro studies demonstrated that homoarginine profoundly modulated the spatial organization of the T-cell actin cytoskeleton and increased filopodia formation via inhibition of Myh9 (myosin heavy chain 9). Further mechanistic studies revealed an inhibition of T-cell proliferation as well as a striking impairment of the migratory capacities of T cells in response to relevant chemokines by homoarginine, all of which likely contribute to its atheroprotective effects. CONCLUSIONS: Our study unravels a novel mechanism by which the amino acid homoarginine reduces atherosclerosis, establishing that homoarginine modulates the T-cell cytoskeleton and thereby mitigates T-cell functions important during atherogenesis. These findings provide a molecular explanation for the beneficial effects of homoarginine in atherosclerotic cardiovascular disease.

Protective immune trajectories in early viral containment of non-pneumonic SARS-CoV-2 infection.

The antiviral immune response to SARS-CoV-2 infection can limit viral spread and prevent development of pneumonic COVID-19. However, the protective immunological response associated with successful viral containment in the upper airways remains unclear. Here, we combine a multi-omics approach with longitudinal sampling to reveal temporally resolved protective immune signatures in non-pneumonic and ambulatory SARS-CoV-2 infected patients and associate specific immune trajectories with upper airway viral containment. We see a distinct systemic rather than local immune state associated with viral containment, characterized by interferon stimulated gene (ISG) upregulation across circulating immune cell subsets in non-pneumonic SARS-CoV2 infection. We report reduced cytotoxic potential of Natural Killer (NK) and T cells, and an immune-modulatory monocyte phenotype associated with protective immunity in COVID-19. Together, we show protective immune trajectories in SARS-CoV2 infection, which have important implications for patient prognosis and the development of immunomodulatory therapies.

iBRET Screen of the ABCD1 Peroxisomal Network and Mutation-Induced Network Perturbations.

Mapping the network of proteins provides a powerful means to investigate the function of disease genes and to unravel the molecular basis of phenotypes. We present an automated informatics-aided and bioluminescence resonance energy transfer-based approach (iBRET) enabling high-confidence detection of protein-protein interactions in living mammalian cells. A screen of the ABCD1 protein, which is affected in X-linked adrenoleukodystrophy (X-ALD), against an organelle library of peroxisomal proteins demonstrated applicability of iBRET for large-scale experiments. We identified novel protein-protein interactions for ABCD1 (with ALDH3A2, DAO, ECI2, FAR1, PEX10, PEX13, PEX5, PXMP2, and PIPOX), mapped its position within the peroxisomal protein-protein interaction network, and determined that pathogenic missense variants in ABCD1 alter the interaction with selected binding partners. These findings provide mechanistic insights into pathophysiology of X-ALD and may foster the identification of new disease modifiers.

Self-sustaining IL-8 loops drive a prothrombotic neutrophil phenotype in severe COVID-19.

Neutrophils provide a critical line of defense in immune responses to various pathogens, inflicting self-damage upon transition to a hyperactivated, procoagulant state. Recent work has highlighted proinflammatory neutrophil phenotypes contributing to lung injury and acute respiratory distress syndrome (ARDS) in patients with coronavirus disease 2019 (COVID-19). Here, we use state-of-the art mass spectrometry-based proteomics and transcriptomic and correlative analyses as well as functional in vitro and in vivo studies to dissect how neutrophils contribute to the progression to severe COVID-19. We identify a reinforcing loop of both systemic and neutrophil intrinsic IL-8 (CXCL8/IL-8) dysregulation, which initiates and perpetuates neutrophil-driven immunopathology. This positive feedback loop of systemic and neutrophil autocrine IL-8 production leads to an activated, prothrombotic neutrophil phenotype characterized by degranulation and neutrophil extracellular trap (NET) formation. In severe COVID-19, neutrophils directly initiate the coagulation and complement cascade, highlighting a link to the immunothrombotic state observed in these patients. Targeting the IL-8-CXCR-1/-2 axis interferes with this vicious cycle and attenuates neutrophil activation, degranulation, NETosis, and IL-8 release. Finally, we show that blocking IL-8-like signaling reduces severe acute respiratory distress syndrome of coronavirus 2 (SARS-CoV-2) spike protein-induced, human ACE2-dependent pulmonary microthrombosis in mice. In summary, our data provide comprehensive insights into the activation mechanisms of neutrophils in COVID-19 and uncover a self-sustaining neutrophil-IL-8 axis as a promising therapeutic target in severe SARS-CoV-2 infection.

Toward a gold standard for benchmarking gene set enrichment analysis.

MOTIVATION: Although gene set enrichment analysis has become an integral part of high-throughput gene expression data analysis, the assessment of enrichment methods remains rudimentary and ad hoc. In the absence of suitable gold standards, evaluations are commonly restricted to selected datasets and biological reasoning on the relevance of resulting enriched gene sets. RESULTS: We develop an extensible framework for reproducible benchmarking of enrichment methods based on defined criteria for applicability, gene set prioritization and detection of relevant processes. This framework incorporates a curated compendium of 75 expression datasets investigating 42 human diseases. The compendium features microarray and RNA-seq measurements, and each dataset is associated with a precompiled GO/KEGG relevance ranking for the corresponding disease under investigation. We perform a comprehensive assessment of 10 major enrichment methods, identifying significant differences in runtime and applicability to RNA-seq data, fraction of enriched gene sets depending on the null hypothesis tested and recovery of the predefined relevance rankings. We make practical recommendations on how methods originally developed for microarray data can efficiently be applied to RNA-seq data, how to interpret results depending on the type of gene set test conducted and which methods are best suited to effectively prioritize gene sets with high phenotype relevance. AVAILABILITY: http://bioconductor.org/packages/GSEABenchmarkeR. CONTACT: ludwig.geistlinger@sph.cuny.edu.

Vascular neutrophilic inflammation and immunothrombosis distinguish severe COVID-19 from influenza pneumonia.

OBJECTIVE: Infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) can lead to severe pneumonia, but also thrombotic complications and non-pulmonary organ failure. Recent studies suggest intravascular neutrophil activation and subsequent immune cell-triggered immunothrombosis as a central pathomechanism linking the heterogenous clinical picture of coronavirus disease 2019 (COVID-19). We sought to study whether immunothrombosis is a pathognomonic factor in COVID-19 or a general feature of (viral) pneumonia, as well as to better understand its upstream regulation. APPROACH AND RESULTS: By comparing histopathological specimens of SARS-CoV-2 with influenza-affected lungs, we show that vascular neutrophil recruitment, NETosis, and subsequent immunothrombosis are typical features of severe COVID-19, but less prominent in influenza pneumonia. Activated neutrophils were typically found in physical association with monocytes. To explore this further, we combined clinical data of COVID-19 cases with comprehensive immune cell phenotyping and bronchoalveolar lavage fluid scRNA-seq data. We show that a HLADRlow CD9low monocyte population expands in severe COVID-19, which releases neutrophil chemokines in the lungs, and might in turn explain neutrophil expansion and pulmonary recruitment in the late stages of severe COVID-19. CONCLUSIONS: Our data underline an innate immune cell axis causing vascular inflammation and immunothrombosis in severe SARS-CoV-2 infection.

SEQU-INTO: Early detection of impurities, contamination and off-targets (ICOs) in long read/MinION sequencing.

The MinION sequencer by Oxford Nanopore Technologies turns DNA and RNA sequencing into a routine task in biology laboratories or in field research. For downstream analysis it is required to have a sufficient amount of target reads. Especially prokaryotic or bacteriophagic sequencing samples can contain a significant amount of off-target sequences in the processed sample, stemming from human DNA/RNA contamination, insufficient rRNA depletion, or remaining DNA/RNA from other organisms (e.g. host organism from bacteriophage cultivation). Such impurity, contamination and off-targets (ICOs) block read capacity, requiring to sequence deeper. In comparison to second-generation sequencing, MinION sequencing allows to reuse its chip after a (partial) run. This allows further usage of the same chip with more sample, even after adjusting the library preparation to reduce ICOs. The earlier a sample's ICOs are detected, the better the sequencing chip can be conserved for future use. Here we present sequ-into, a low-resource and user-friendly cross-platform tool to detect ICO sequences from a predefined ICO database in samples early during a MinION sequencing run. The data provided by sequ-into empowers the user to quickly take action to preserve sample material and chip capacity. sequ-into is available from https://github.com/mjoppich/sequ-into.

Integrative functional genomics decodes herpes simplex virus 1.

The predicted 80 open reading frames (ORFs) of herpes simplex virus 1 (HSV-1) have been intensively studied for decades. Here, we unravel the complete viral transcriptome and translatome during lytic infection with base-pair resolution by computational integration of multi-omics data. We identify a total of 201 transcripts and 284 ORFs including all known and 46 novel large ORFs. This includes a so far unknown ORF in the locus deleted in the FDA-approved oncolytic virus Imlygic. Multiple transcript isoforms expressed from individual gene loci explain translation of the vast majority of ORFs as well as N-terminal extensions (NTEs) and truncations. We show that NTEs with non-canonical start codons govern the subcellular protein localization and packaging of key viral regulators and structural proteins. We extend the current nomenclature to include all viral gene products and provide a genome browser that visualizes all the obtained data from whole genome to single-nucleotide resolution.

CNVRanger: association analysis of CNVs with gene expression and quantitative phenotypes.

SUMMARY: Copy number variation (CNV) is a major type of structural genomic variation that is increasingly studied across different species for association with diseases and production traits. Established protocols for experimental detection and computational inference of CNVs from SNP array and next-generation sequencing data are available. We present the CNVRanger R/Bioconductor package which implements a comprehensive toolbox for structured downstream analysis of CNVs. This includes functionality for summarizing individual CNV calls across a population, assessing overlap with functional genomic regions, and genome-wide association analysis with gene expression and quantitative phenotypes. AVAILABILITY AND IMPLEMENTATION: http://bioconductor.org/packages/CNVRanger.

The Heat Shock Response in Yeast Maintains Protein Homeostasis by Chaperoning and Replenishing Proteins.

Life is resilient because living systems are able to respond to elevated temperatures with an ancient gene expression program called the heat shock response (HSR). In yeast, the transcription of hundreds of genes is upregulated at stress temperatures. Besides stress protection conferred by chaperones, the function of the majority of the upregulated genes under stress has remained enigmatic. We show that those genes are required to directly counterbalance increased protein turnover at stress temperatures and to maintain the metabolism. This anaplerotic reaction together with molecular chaperones allows yeast to efficiently buffer proteotoxic stress. When the capacity of this system is exhausted at extreme temperatures, aggregation processes stop translation and growth pauses. The emerging concept is that the HSR is modular with distinct programs dependent on the severity of the stress.

From command-line bioinformatics to bioGUI.

Bioinformatics is a highly interdisciplinary field providing (bioinformatics) applications for scientists from many disciplines. Installing and starting applications on the command-line (CL) is inconvenient and/or inefficient for many scientists. Nonetheless, most methods are implemented with a command-line interface only. Providing a graphical user interface (GUI) for bioinformatics applications is one step toward routinely making CL-only applications available to more scientists and, thus, toward a more effective interdisciplinary work. With our bioGUI framework we address two main problems of using CL bioinformatics applications: First, many tools work on UNIX-systems only, while many scientists use Microsoft Windows. Second, scientists refrain from using CL tools which, however, could well support them in their research. With bioGUI install modules and templates, installing and using CL tools is made possible for most scientists-even on Windows, due to bioGUI's support for Windows Subsystem for Linux. In addition, bioGUI templates can easily be created, making the bioGUI framework highly rewarding for developers. From the bioGUI repository it is possible to download, install and use bioinformatics tools with just a few clicks.

A transcriptional response of Clostridium beijerinckii NRRL B-598 to a butanol shock.

BACKGROUND: One of the main obstacles preventing solventogenic clostridia from achieving higher yields in biofuel production is the toxicity of produced solvents. Unfortunately, regulatory mechanisms responsible for the shock response are poorly described on the transcriptomic level. Although the strain Clostridium beijerinckii NRRL B-598, a promising butanol producer, has been studied under different conditions in the past, its transcriptional response to a shock caused by butanol in the cultivation medium remains unknown. RESULTS: In this paper, we present a transcriptional response of the strain during a butanol challenge, caused by the addition of butanol to the cultivation medium at the very end of the acidogenic phase, using RNA-Seq. We resequenced and reassembled the genome sequence of the strain and prepared novel genome and gene ontology annotation to provide the most accurate results. When compared to samples under standard cultivation conditions, samples gathered during butanol shock represented a well-distinguished group. Using reference samples gathered directly before the addition of butanol, we identified genes that were differentially expressed in butanol challenge samples. We determined clusters of 293 down-regulated and 301 up-regulated genes whose expression was affected by the cultivation conditions. Enriched term "RNA binding" among down-regulated genes corresponded to the downturn of translation and the cluster contained a group of small acid-soluble spore proteins. This explained phenotype of the culture that had not sporulated. On the other hand, up-regulated genes were characterized by the term "protein binding" which corresponded to activation of heat-shock proteins that were identified within this cluster. CONCLUSIONS: We provided an overall transcriptional response of the strain C. beijerinckii NRRL B-598 to butanol shock, supplemented by auxiliary technologies, including high-pressure liquid chromatography and flow cytometry, to capture the corresponding phenotypic response. We identified genes whose regulation was affected by the addition of butanol to the cultivation medium and inferred related molecular functions that were significantly influenced. Additionally, using high-quality genome assembly and custom-made gene ontology annotation, we demonstrated that this settled terminology, widely used for the analysis of model organisms, could also be applied to non-model organisms and for research in the field of biofuels.