hCoCena: A toolbox for network-based co-expression analysis and horizontal integration of transcriptomic datasets

As the number and complexity of transcriptomic datasets increase, there is a rising demand for accessible and user-friendly analysis tools. Here, we present hCoCena (horizontal construction of co-expression networks and analysis), a toolbox facilitating the analysis of a single dataset, as well as the joint analysis of multiple datasets. We describe steps for workspace setup, formatting tables, data processing, and network integration. We then detail procedures for gene clustering, gene set enrichment analysis, and transcription factor enrichment analysis. For complete details on the use and execution of this protocol, please refer to Oestreich et al.1.

Molecular mechanisms and treatment responses of pulmonary fibrosis in severe COVID-19

BACKGROUND

Coronavirus disease 2019 (COVID-19) patients can develop pulmonary fibrosis (PF), which is associated with impaired outcome. We assessed specific leukocytic transcriptome profiles associated with PF and the influence of early dexamethasone (DEXA) treatment on the clinical course of PF in critically ill COVID-19 patients.

METHODS

We performed a pre-post design study in 191 COVID-19 patients admitted to the Intensive Care Unit (ICU) spanning two treatment cohorts: the pre-DEXA- (n = 67) and the DEXA-cohort (n = 124). PF was identified based on radiological findings, worsening of ventilatory parameters and elevated circulating PIIINP levels. Longitudinal transcriptome profiles of 52 pre-DEXA patients were determined using RNA sequencing. Effects of prednisone treatment on clinical fibrosis parameters and outcomes were analyzed between PF- and no-PF-patients within both cohorts.

RESULTS

Transcriptome analyses revealed upregulation of inflammatory, coagulation and neutrophil extracellular trap-related pathways in PF-patients compared to no-PF patients. Key genes involved included PADI4, PDE4D, MMP8, CRISP3, and BCL2L15. Enrichment of several identified pathways was associated with impaired survival in a external cohort of patients with idiopathic pulmonary fibrosis. Following prednisone treatment, PF-related profiles reverted towards those observed in the no-PF-group. Likewise, PIIINP levels decreased significantly following prednisone treatment. PF incidence was 28% and 25% in the pre-DEXA- and DEXA-cohort, respectively (p = 0.61). ICU length-of-stay (pre-DEXA: 42 [29-49] vs. 18 [13-27] days, p < 0.001; DEXA: 42 [28-57] vs. 13 [7-24] days, p < 0.001) and mortality (pre-DEXA: 47% vs. 15%, p = 0.009; DEXA: 61% vs. 19%, p < 0.001) were higher in the PF-groups compared to the no-PF-groups within both cohorts. Early dexamethasone therapy did not influence these outcomes.

CONCLUSIONS

ICU patients with COVID-19 who develop PF exhibit upregulated coagulation, inflammation, and neutrophil extracellular trap-related pathways as well as prolonged ICU length-of-stay and mortality. This study indicates that early dexamethasone treatment neither influences the incidence or clinical course of PF, nor clinical outcomes.

Identification of the novel FOXP3-dependent Treg cell transcription factor MEOX1 by high-dimensional analysis of human CD4+ T cells

CD4+ T cells play a central role in the adaptive immune response through their capacity to activate, support and control other immune cells. Although these cells have become the focus of intense research, a comprehensive understanding of the underlying regulatory networks that orchestrate CD4+ T cell function and activation is still incomplete. Here, we analyzed a large transcriptomic dataset consisting of 48 different human CD4+ T cell conditions. By performing reverse network engineering, we identified six common denominators of CD4+ T cell functionality (CREB1, E2F3, AHR, STAT1, NFAT5 and NFATC3). Moreover, we also analyzed condition-specific genes which led us to the identification of the transcription factor MEOX1 in Treg cells. Expression of MEOX1 was comparable to FOXP3 in Treg cells and can be upregulated by IL-2. Epigenetic analyses revealed a permissive epigenetic landscape for MEOX1 solely in Treg cells. Knockdown of MEOX1 in Treg cells revealed a profound impact on downstream gene expression programs and Treg cell suppressive capacity. These findings in the context of CD4+ T cells contribute to a better understanding of the transcriptional networks and biological mechanisms controlling CD4+ T cell functionality, which opens new avenues for future therapeutic strategies.

Systemic alterations in neutrophils and their precursors in early-stage chronic obstructive pulmonary disease

Systemic inflammation is established as part of late-stage severe lung disease, but molecular, functional, and phenotypic changes in peripheral immune cells in early disease stages remain ill defined. Chronic obstructive pulmonary disease (COPD) is a major respiratory disease characterized by small-airway inflammation, emphysema, and severe breathing difficulties. Using single-cell analyses we demonstrate that blood neutrophils are already increased in early-stage COPD, and changes in molecular and functional neutrophil states correlate with lung function decline. Assessing neutrophils and their bone marrow precursors in a murine cigarette smoke exposure model identified similar molecular changes in blood neutrophils and precursor populations that also occur in the blood and lung. Our study shows that systemic molecular alterations in neutrophils and their precursors are part of early-stage COPD, a finding to be further explored for potential therapeutic targets and biomarkers for early diagnosis and patient stratification.

Circuit-selective cell-autonomous regulation of inhibition in pyramidal neurons by Ste20-like kinase

Maintaining an appropriate balance between excitation and inhibition is critical for neuronal information processing. Cortical neurons can cell-autonomously adjust the inhibition they receive to individual levels of excitatory input, but the underlying mechanisms are unclear. We describe that Ste20-like kinase (SLK) mediates cell-autonomous regulation of excitation-inhibition balance in the thalamocortical feedforward circuit, but not in the feedback circuit. This effect is due to regulation of inhibition originating from parvalbumin-expressing interneurons, while inhibition via somatostatin-expressing interneurons is unaffected. Computational modeling shows that this mechanism promotes stable excitatory-inhibitory ratios across pyramidal cells and ensures robust and sparse coding. Patch-clamp RNA sequencing yields genes differentially regulated by SLK knockdown, as well as genes associated with excitation-inhibition balance participating in transsynaptic communication and cytoskeletal dynamics. These data identify a mechanism for cell-autonomous regulation of a specific inhibitory circuit that is critical to ensure that a majority of cortical pyramidal cells participate in information coding.

Human variation in population-wide gene expression data predicts gene perturbation phenotype

Population-scale datasets of healthy individuals capture genetic and environmental factors influencing gene expression. The expression variance of a gene of interest (GOI) can be exploited to set up a quasi loss- or gain-of-function "in population" experiment. We describe here an approach, huva (human variation), taking advantage of population-scale multi-layered data to infer gene function and relationships between phenotypes and expression. Within a reference dataset, huva derives two experimental groups with LOW or HIGH expression of the GOI, enabling the subsequent comparison of their transcriptional profile and functional parameters. We demonstrate that this approach robustly identifies the phenotypic relevance of a GOI allowing the stratification of genes according to biological functions, and we generalize this concept to almost 16,000 genes in the human transcriptome. Additionally, we describe how huva predicts monocytes to be the major cell type in the pathophysiology of STAT1 mutations, evidence validated in a clinical cohort.

BCL7A-containing SWI/SNF/BAF complexes modulate mitochondrial bioenergetics during neural progenitor differentiation

Mammalian SWI/SNF/BAF chromatin remodeling complexes influence cell lineage determination. While the contribution of these complexes to neural progenitor cell (NPC) proliferation and differentiation has been reported, little is known about the transcriptional profiles that determine neurogenesis or gliogenesis. Here, we report that BCL7A is a modulator of the SWI/SNF/BAF complex that stimulates the genome-wide occupancy of the ATPase subunit BRG1. We demonstrate that BCL7A is dispensable for SWI/SNF/BAF complex integrity, whereas it is essential to regulate Notch/Wnt pathway signaling and mitochondrial bioenergetics in differentiating NPCs. Pharmacological stimulation of Wnt signaling restores mitochondrial respiration and attenuates the defective neurogenic patterns observed in NPCs lacking BCL7A. Consistently, treatment with an enhancer of mitochondrial biogenesis, pioglitazone, partially restores mitochondrial respiration and stimulates neuronal differentiation of BCL7A-deficient NPCs. Using conditional BCL7A knockout mice, we reveal that BCL7A expression in NPCs and postmitotic neurons is required for neuronal plasticity and supports behavioral and cognitive performance. Together, our findings define the specific contribution of BCL7A-containing SWI/SNF/BAF complexes to mitochondria-driven NPC commitment, thereby providing a better understanding of the cell-intrinsic transcriptional processes that connect metabolism, neuronal morphogenesis, and cognitive flexibility.

Decoding mechanism of action and sensitivity to drug candidates from integrated transcriptome and chromatin state

Omics-based technologies are driving major advances in precision medicine, but efforts are still required to consolidate their use in drug discovery. In this work, we exemplify the use of multi-omics to support the development of 3-chloropiperidines, a new class of candidate anticancer agents. Combined analyses of transcriptome and chromatin accessibility elucidated the mechanisms underlying sensitivity to test agents. Furthermore, we implemented a new versatile strategy for the integration of RNA- and ATAC-seq (Assay for Transposase-Accessible Chromatin) data, able to accelerate and extend the standalone analyses of distinct omic layers. This platform guided the construction of a perturbation-informed basal signature predicting cancer cell lines' sensitivity and to further direct compound development against specific tumor types. Overall, this approach offers a scalable pipeline to support the early phases of drug discovery, understanding of mechanisms, and potentially inform the positioning of therapeutics in the clinic.

hCoCena: Horizontal integration and analysis of transcriptomics datasets

Motivation

Transcriptome-based gene co-expression analysis has become a standard procedure for structured and contextualized understanding and comparison of different conditions and phenotypes. Since large study designs with a broad variety of conditions are costly and laborious, extensive comparisons are hindered when utilizing only a single dataset. Thus, there is an increased need for tools that allow the integration of multiple transcriptomic datasets with subsequent joint analysis, which can provide a more systematic understanding of gene co-expression and co-functionality within and across conditions. To make such an integrative analysis accessible to a wide spectrum of users with differing levels of programming expertise it is essential to provide user-friendliness and customizability as well as thorough documentation.

Results

This article introduces horizontal CoCena (hCoCena: horizontal construction of co-expression networks and analysis), an R-package for network-based co-expression analysis that allows the analysis of a single transcriptomic dataset as well as the joint analysis of multiple datasets. With hCoCena, we provide a freely available, user-friendly and adaptable tool for integrative multi-study or single-study transcriptomics analyses alongside extensive comparisons to other existing tools.

Availability and implementation

The hCoCena R-package is provided together with R Markdowns that implement an exemplary analysis workflow including extensive documentation and detailed descriptions of data structures and objects. Such efforts not only make the tool easy to use but also enable the seamless integration of user-written scripts and functions into the workflow, creating a tool that provides a clear design while remaining flexible and highly customizable. The package and additional information including an extensive Wiki are freely available on GitHub: https://github.com/MarieOestreich/hCoCena. The version at the time of writing has been added to Zenodo under the following link: https://doi.org/10.5281/zenodo.6911782.

Supplementary information

Supplementary data are available at Bioinformatics online.

Alveolar macrophages in early stage COPD show functional deviations with properties of impaired immune activation

Despite its high prevalence, the cellular and molecular mechanisms of chronic obstructive pulmonary disease (COPD) are far from being understood. Here, we determine disease-related changes in cellular and molecular compositions within the alveolar space and peripheral blood of a cohort of COPD patients and controls. Myeloid cells were the largest cellular compartment in the alveolar space with invading monocytes and proliferating macrophages elevated in COPD. Modeling cell-to-cell communication, signaling pathway usage, and transcription factor binding predicts TGF-β1 to be a major upstream regulator of transcriptional changes in alveolar macrophages of COPD patients. Functionally, macrophages in COPD showed reduced antigen presentation capacity, accumulation of cholesteryl ester, reduced cellular chemotaxis, and mitochondrial dysfunction, reminiscent of impaired immune activation.

Does hand mirror morphology a useful diagnostic tool for the treatment decision of acute undifferentiated leukemia at admission?

It is difficult how to manage acute undifferentiated leukemia in daily practice, but generally, hand mirror morphology provides ease to treat these patients. Thirty-nine years old male patient was admitted to with the complaints of echymosis, and pain at his left buttock due to an intramuscular injection for the treatment of previously diagnosed of the lower respiratory infection. Peripheral blood smear revealed >%50 blasts cells with a moderate nuclear: cytoplasmic ratio and one or more nucleoli. The blast cells showed a hand-mirror morphology and not harboring auer rods. According to the flow cytometric analysis the blastic cells do not represent to be originated from myeloid or lymphoid origin, because the cells harboring both of two cell lineages. AML-like therapy was commenced based on the positive myeloid markers including CD117 and CD135. Even though hand mirror morphology of the blasts usually demonstrates the lymphoid origin and the patients are treated as ALL like therapy, myeloid blasts rarely represents the same morphology, as was in our patient.

Mature neutrophils and a NF-κB-to-IFN transition determine the unifying disease recovery dynamics in COVID-19

Disease recovery dynamics are often difficult to assess, as patients display heterogeneous recovery courses. To model recovery dynamics, exemplified by severe COVID-19, we apply a computational scheme on longitudinally sampled blood transcriptomes, generating recovery states, which we then link to cellular and molecular mechanisms, presenting a framework for studying the kinetics of recovery compared with non-recovery over time and long-term effects of the disease. Specifically, a decrease in mature neutrophils is the strongest cellular effect during recovery, with direct implications on disease outcome. Furthermore, we present strong indications for global regulatory changes in gene programs, decoupled from cell compositional changes, including an early rise in T cell activation and differentiation, resulting in immune rebalancing between interferon and NF-κB activity and restoration of cell homeostasis. Overall, we present a clinically relevant computational framework for modeling disease recovery, paving the way for future studies of the recovery dynamics in other diseases and tissues.

Impaired neurogenesis alters brain biomechanics in a neuroprogenitor-based genetic subtype of congenital hydrocephalus

Hydrocephalus, characterized by cerebral ventricular dilatation, is routinely attributed to primary defects in cerebrospinal fluid (CSF) homeostasis. This fosters CSF shunting as the leading reason for brain surgery in children despite considerable disease heterogeneity. In this study, by integrating human brain transcriptomics with whole-exome sequencing of 483 patients with congenital hydrocephalus (CH), we found convergence of CH risk genes in embryonic neuroepithelial stem cells. Of all CH risk genes, TRIM71/lin-41 harbors the most de novo mutations and is most specifically expressed in neuroepithelial cells. Mice harboring neuroepithelial cell-specific Trim71 deletion or CH-specific Trim71 mutation exhibit prenatal hydrocephalus. CH mutations disrupt TRIM71 binding to its RNA targets, causing premature neuroepithelial cell differentiation and reduced neurogenesis. Cortical hypoplasia leads to a hypercompliant cortex and secondary ventricular enlargement without primary defects in CSF circulation. These data highlight the importance of precisely regulated neuroepithelial cell fate for normal brain-CSF biomechanics and support a clinically relevant neuroprogenitor-based paradigm of CH.

Severe COVID-19 Shares a Common Neutrophil Activation Signature with Other Acute Inflammatory States

Severe COVID-19 patients present a clinical and laboratory overlap with other hyperinflammatory conditions such as hemophagocytic lymphohistiocytosis (HLH). However, the underlying mechanisms of these conditions remain to be explored. Here, we investigated the transcriptome of 1596 individuals, including patients with COVID-19 in comparison to healthy controls, other acute inflammatory states (HLH, multisystem inflammatory syndrome in children [MIS-C], Kawasaki disease [KD]), and different respiratory infections (seasonal coronavirus, influenza, bacterial pneumonia). We observed that COVID-19 and HLH share immunological pathways (cytokine/chemokine signaling and neutrophil-mediated immune responses), including gene signatures that stratify COVID-19 patients admitted to the intensive care unit (ICU) and COVID-19_nonICU patients. Of note, among the common differentially expressed genes (DEG), there is a cluster of neutrophil-associated genes that reflects a generalized hyperinflammatory state since it is also dysregulated in patients with KD and bacterial pneumonia. These genes are dysregulated at the protein level across several COVID-19 studies and form an interconnected network with differentially expressed plasma proteins that point to neutrophil hyperactivation in COVID-19 patients admitted to the intensive care unit. scRNAseq analysis indicated that these genes are specifically upregulated across different leukocyte populations, including lymphocyte subsets and immature neutrophils. Artificial intelligence modeling confirmed the strong association of these genes with COVID-19 severity. Thus, our work indicates putative therapeutic pathways for intervention.

The clinicopathological features and survival of Castleman disease: a multicenter Turkish study

OBJECTIVE

In this study, we aimed to investigate the clinicopathological features and survival of CD, which is quite rare and has many unknowns.

PATIENTS AND METHODS

This study was conducted by retrospectively evaluating patients diagnosed with CD in six different centers in Turkey.

RESULTS

The median age of 33 patients included in the study was 49 and 51.5% (n = 17) of these patients were women. 18 (54.5%) patients were in the hyaline vascular subtype and most of the patients were UCD (n = 20, 60.6%). The most common involvement region was head and neck (n = 19, 57.5%). The UCD group was younger than the MCD group (p=0.027). Visceral lymph node involvement was higher in MCD than in UCD (p=0.001). Similarly, it was observed that there was more hepatomegaly (p=0.035) and splenomegaly (p=0.013) in the MCD group. During the median 19.5 months follow-up period, there were no patients who died.

CONCLUSIONS

It was observed that UCD and MCD are different clinical entities. Promising survival times can be achieved with surgical and systemic treatments in both subtypes of this extremely rare disease. However, this result should be supported by well-designed prospective comprehensive studies.

Creld2 function during unfolded protein response is essential for liver metabolism homeostasis

The unfolded protein response (UPR) is associated with hepatic metabolic function, yet it is not well understood how endoplasmic reticulum (ER) disturbance might influence metabolic homeostasis. Here, we describe the physiological function of Cysteine-rich with EGF-like domains 2 (Creld2), previously characterized as a downstream target of the ER-stress signal transducer Atf6. To this end, we generated Creld2-deficient mice and induced UPR by injection of tunicamycin. Creld2 augments protein folding and creates an interlink between the UPR axes through its interaction with proteins involved in the cellular stress response. Thereby, Creld2 promotes tolerance to ER stress and recovery from acute stress. Creld2-deficiency leads to a dysregulated UPR and causes the development of hepatic steatosis during ER stress conditions. Moreover, Creld2-dependent enhancement of the UPR assists in the regulation of energy expenditure. Furthermore, we observed a sex dimorphism in human and mouse livers with only male patients showing an accumulation of CRELD2 protein during the progression from non-alcoholic fatty liver disease to non-alcoholic steatohepatitis and only male Creld2-deficient mice developing hepatic steatosis upon aging. These results reveal a Creld2 function at the intersection between UPR and metabolic homeostasis and suggest a mechanism in which chronic ER stress underlies fatty liver disease in males.

Alveolar macrophage transcriptomic profiling in COPD shows major lipid metabolism changes

Background

Immune cells play a major role in the pathogenesis of COPD. Changes in the distribution and cellular functions of major immune cells, such as alveolar macrophages (AMs) and neutrophils are well known; however, their transcriptional reprogramming and contribution to the pathophysiology of COPD are still not fully understood.

Method

To determine changes in transcriptional reprogramming and lipid metabolism in the major immune cell type within bronchoalveolar lavage fluid, we analysed whole transcriptomes and lipidomes of sorted CD45⁺Lin⁻HLA-DR⁺CD66b⁻Autofluorescence^hi AMs from controls and COPD patients.

Results

We observed global transcriptional reprogramming featuring a spectrum of activation states, including pro- and anti-inflammatory signatures. We further detected significant changes between COPD patients and controls in genes involved in lipid metabolism, such as fatty acid biosynthesis in GOLD2 patients. Based on these findings, assessment of a total of 202 lipid species in sorted AMs revealed changes of cholesteryl esters, monoacylglycerols and phospholipids in a disease grade-dependent manner.

Conclusions

Transcriptome and lipidome profiling of COPD AMs revealed GOLD grade-dependent changes, such as in cholesterol metabolism and interferon-α and γ responses.

AMs from COPD patients undergo GOLD grade-specific transcriptional reprogramming and acquire a complex activation profile. Among the observed changes are gene programmes involved in lipid metabolism that translate into alterations in the AM lipidome.https://bit.ly/3sYAqgd

Cannabinoid receptor 2 is necessary to induce toll-like receptor-mediated microglial activation

The tight regulation of microglia activity is key for precise responses to potential threats, while uncontrolled and exacerbated microglial activity is neurotoxic. Microglial toll-like receptors (TLRs) are indispensable for sensing different types of assaults and triggering an innate immune response. Cannabinoid receptor 2 (CB2) signaling is a key pathway to control microglial homeostasis and activation, and its activation is connected to changes in microglial activity. We aimed to investigate how CB2 signaling impacts TLR-mediated microglial activation. Here, we demonstrate that deletion of CB2 causes a dampened transcriptional response to prototypic TLR ligands in microglia. Loss of CB2 results in distinct microglial gene expression profiles, morphology, and activation. We show that the CB2-mediated attenuation of TLR-induced microglial activation is mainly p38 MAPK-dependent. Taken together, we demonstrate that CB2 expression and signaling are necessary to fine-tune TLR-induced activation programs in microglia.

Creld1 regulates myocardial development and function

CRELD1 (Cysteine-Rich with EGF-Like Domains 1) is a risk gene for non-syndromic atrioventricular septal defects in human patients. In a mouse model, Creld1 has been shown to be essential for heart development, particularly in septum and valve formation. However, due to the embryonic lethality of global Creld1 knockout (KO) mice, its cell type-specific function during peri- and postnatal stages remains unknown. Here, we generated conditional Creld1 KO mice lacking Creld1 either in the endocardium (KOTie2) or the myocardium (KOMyHC). Using a combination of cardiac phenotyping, histology, immunohistochemistry, RNA-sequencing, and flow cytometry, we demonstrate that Creld1 function in the endocardium is dispensable for heart development. Lack of myocardial Creld1 causes extracellular matrix remodeling and trabeculation defects by modulation of the Notch1 signaling pathway. Hence, KOMyHC mice die early postnatally due to myocardial hypoplasia. Our results reveal that Creld1 not only controls the formation of septa and valves at an early stage during heart development, but also cardiac maturation and function at a later stage. These findings underline the central role of Creld1 in mammalian heart development and function.

Neue Krankheiten mit Bluttranskriptomik entschlüsseln

The COVID-19 pandemic is leading to increasing numbers of patients all over the world. Reports on a dysregulated immune system in the severe cases calls for a better characterization of the ongoing changes. To dissect COVID-19-driven immune host responses, we profiled whole blood transcriptomes enabling a data-driven stratification based on molecular phenotype. This analysis allowed prediction of patient subgroup-specific drug candidates targeting the dysregulated systemic immune response of the host.

Swarm Learning for decentralized and confidential clinical machine learning

Fast and reliable detection of patients with severe and heterogeneous illnesses is a major goal of precision medicine^1,2. Patients with leukaemia can be identified using machine learning on the basis of their blood transcriptomes³. However, there is an increasing divide between what is technically possible and what is allowed, because of privacy legislation^4,5. Here, to facilitate the integration of any medical data from any data owner worldwide without violating privacy laws, we introduce Swarm Learning—a decentralized machine-learning approach that unites edge computing, blockchain-based peer-to-peer networking and coordination while maintaining confidentiality without the need for a central coordinator, thereby going beyond federated learning. To illustrate the feasibility of using Swarm Learning to develop disease classifiers using distributed data, we chose four use cases of heterogeneous diseases (COVID-19, tuberculosis, leukaemia and lung pathologies). With more than 16,400 blood transcriptomes derived from 127 clinical studies with non-uniform distributions of cases and controls and substantial study biases, as well as more than 95,000 chest X-ray images, we show that Swarm Learning classifiers outperform those developed at individual sites. In addition, Swarm Learning completely fulfils local confidentiality regulations by design. We believe that this approach will notably accelerate the introduction of precision medicine.

Swarm Learning is a decentralized machine learning approach that outperforms classifiers developed at individual sites for COVID-19 and other diseases while preserving confidentiality and privacy.

The stem cell-specific protein TRIM71 inhibits maturation and activity of the pro-differentiation miRNA let-7 via two independent molecular mechanisms

The stem cell-specific RNA-binding protein TRIM71/LIN-41 was the first identified target of the pro-differentiation and tumor suppressor miRNA let-7. TRIM71 has essential functions in embryonic development and a proposed oncogenic role in several cancer types, such as hepatocellular carcinoma. Here, we show that TRIM71 regulates let-7 expression and activity via two independent mechanisms. On the one hand, TRIM71 enhances pre-let-7 degradation through its direct interaction with LIN28 and TUT4, thereby inhibiting let-7 maturation and indirectly promoting the stabilization of let-7 targets. On the other hand, TRIM71 represses the activity of mature let-7 via its RNA-dependent interaction with the RNA-Induced Silencing Complex (RISC) effector protein AGO2. We found that TRIM71 directly binds and stabilizes let-7 targets, suggesting that let-7 activity inhibition occurs on active RISCs. MiRNA enrichment analysis of several transcriptomic datasets from mouse embryonic stem cells and human hepatocellular carcinoma cells suggests that these let-7 regulatory mechanisms shape transcriptomic changes during developmental and oncogenic processes. Altogether, our work reveals a novel role for TRIM71 as a miRNA repressor and sheds light on a dual mechanism of let-7 regulation.

The gut microbiota is associated with the small intestinal paracellular permeability and the development of the immune system in healthy children during the first two years of life

BACKGROUND

The intestinal barrier plays an important role in the defense against infections, and nutritional, endocrine, and immune functions. The gut microbiota playing an important role in development of the gastrointestinal tract can impact intestinal permeability and immunity during early life, but data concerning this problem are scarce.

METHODS

We analyzed the microbiota in fecal samples (101 samples in total) collected longitudinally over 24 months from 21 newborns to investigate whether the markers of small intestinal paracellular permeability (zonulin) and immune system development (calprotectin) are linked to the gut microbiota. The results were validated using data from an independent cohort that included the calprotectin and gut microbiota in children during the first year of life.

RESULTS

Zonulin levels tended to increase for up to 6 months after childbirth and stabilize thereafter remaining at a high level while calprotectin concentration was high after childbirth and began to decline from 6 months of life. The gut microbiota composition and the related metabolic potentials changed during the first 2 years of life and were correlated with zonulin and calprotectin levels. Faecal calprotectin correlated inversely with alpha diversity (Shannon index, r = - 0.30, FDR P (Q) = 0.039). It also correlated with seven taxa; i.a. negatively with Ruminococcaceae (r = - 0.34, Q = 0.046), and Clostridiales (r = - 0.34, Q = 0.048) and positively with Staphylococcus (r = 0.38, Q = 0.023) and Staphylococcaceae (r = 0.35, Q = 0.04), whereas zonulin correlated with 19 taxa; i.a. with Bacillales (r = - 0.52, Q = 0.0004), Clostridiales (r = 0.48, Q = 0.001) and the Ruminococcus (torques group) (r = 0.40, Q = 0.026). When time intervals were considered only changes in abundance of the Ruminococcus (torques group) were associated with changes in calprotectin (β = 2.94, SE = 0.8, Q = 0.015). The dynamics of stool calprotectin was negatively associated with changes in two MetaCyc pathways: pyruvate fermentation to butanoate (β = - 4.54, SE = 1.08, Q = 0.028) and Clostridium acetobutylicum fermentation (β = - 4.48, SE = 1.16, Q = 0.026).

CONCLUSIONS

The small intestinal paracellular permeability, immune system-related markers and gut microbiota change dynamically during the first 2 years of life. The Ruminococcus (torques group) seems to be especially involved in controlling paracellular permeability. Staphylococcus, Staphylococcaceae, Ruminococcaceae, and Clostridiales, may be potential biomarkers of the immune system. Despite observed correlations their clear causation and health consequences were not proven. Mechanistic studies are required.

Disease severity-specific neutrophil signatures in blood transcriptomes stratify COVID-19 patients

BACKGROUND

The SARS-CoV-2 pandemic is currently leading to increasing numbers of COVID-19 patients all over the world. Clinical presentations range from asymptomatic, mild respiratory tract infection, to severe cases with acute respiratory distress syndrome, respiratory failure, and death. Reports on a dysregulated immune system in the severe cases call for a better characterization and understanding of the changes in the immune system.

METHODS

In order to dissect COVID-19-driven immune host responses, we performed RNA-seq of whole blood cell transcriptomes and granulocyte preparations from mild and severe COVID-19 patients and analyzed the data using a combination of conventional and data-driven co-expression analysis. Additionally, publicly available data was used to show the distinction from COVID-19 to other diseases. Reverse drug target prediction was used to identify known or novel drug candidates based on finding from data-driven findings.

RESULTS

Here, we profiled whole blood transcriptomes of 39 COVID-19 patients and 10 control donors enabling a data-driven stratification based on molecular phenotype. Neutrophil activation-associated signatures were prominently enriched in severe patient groups, which was corroborated in whole blood transcriptomes from an independent second cohort of 30 as well as in granulocyte samples from a third cohort of 16 COVID-19 patients (44 samples). Comparison of COVID-19 blood transcriptomes with those of a collection of over 3100 samples derived from 12 different viral infections, inflammatory diseases, and independent control samples revealed highly specific transcriptome signatures for COVID-19. Further, stratified transcriptomes predicted patient subgroup-specific drug candidates targeting the dysregulated systemic immune response of the host.

CONCLUSIONS

Our study provides novel insights in the distinct molecular subgroups or phenotypes that are not simply explained by clinical parameters. We show that whole blood transcriptomes are extremely informative for COVID-19 since they capture granulocytes which are major drivers of disease severity.

Disease severity-specific neutrophil signatures in blood transcriptomes stratify COVID-19 patients

BACKGROUND

The SARS-CoV-2 pandemic is currently leading to increasing numbers of COVID-19 patients all over the world. Clinical presentations range from asymptomatic, mild respiratory tract infection, to severe cases with acute respiratory distress syndrome, respiratory failure, and death. Reports on a dysregulated immune system in the severe cases call for a better characterization and understanding of the changes in the immune system.

METHODS

In order to dissect COVID-19-driven immune host responses, we performed RNA-seq of whole blood cell transcriptomes and granulocyte preparations from mild and severe COVID-19 patients and analyzed the data using a combination of conventional and data-driven co-expression analysis. Additionally, publicly available data was used to show the distinction from COVID-19 to other diseases. Reverse drug target prediction was used to identify known or novel drug candidates based on finding from data-driven findings.

RESULTS

Here, we profiled whole blood transcriptomes of 39 COVID-19 patients and 10 control donors enabling a data-driven stratification based on molecular phenotype. Neutrophil activation-associated signatures were prominently enriched in severe patient groups, which was corroborated in whole blood transcriptomes from an independent second cohort of 30 as well as in granulocyte samples from a third cohort of 16 COVID-19 patients (44 samples). Comparison of COVID-19 blood transcriptomes with those of a collection of over 3100 samples derived from 12 different viral infections, inflammatory diseases, and independent control samples revealed highly specific transcriptome signatures for COVID-19. Further, stratified transcriptomes predicted patient subgroup-specific drug candidates targeting the dysregulated systemic immune response of the host.

CONCLUSIONS

Our study provides novel insights in the distinct molecular subgroups or phenotypes that are not simply explained by clinical parameters. We show that whole blood transcriptomes are extremely informative for COVID-19 since they capture granulocytes which are major drivers of disease severity.

S100A8 and S100A9 Are Important for Postnatal Development of Gut Microbiota and Immune System in Mice and Infants

BACKGROUND & AIMS

After birth, the immune system matures via interactions with microbes in the gut. The S100 calcium binding proteins S100A8 and S100A9, and their extracellular complex form, S100A8-A9, are found in high amounts in human breast milk. We studied levels of S100A8-A9 in fecal samples (also called fecal calprotectin) from newborns and during infancy, and their effects on development of the intestinal microbiota and mucosal immune system.

METHODS

We collected stool samples (n = 517) from full-term (n = 72) and preterm infants (n = 49) at different timepoints over the first year of life (days 1, 3, 10, 30, 90, 180, and 360). We measured levels of S100A8-A9 by enzyme-linked immunosorbent assay and analyzed fecal microbiomes by 16S sRNA gene sequencing. We also obtained small and large intestine biopsies from 8 adults and 10 newborn infants without inflammatory bowel diseases (controls) and 8 infants with necrotizing enterocolitis and measured levels of S100A8 by immunofluorescence microscopy. Children were followed for 2.5 years and anthropometric data and medical information on infections were collected. We performed studies with newborn C57BL/6J wild-type and S100a9^-/- mice (which also lack S100A8). Some mice were fed or given intraperitoneal injections of S100A8 or subcutaneous injections of Staphylococcus aureus. Blood and intestine, mesenterial and celiac lymph nodes were collected; cells and cytokines were measured by flow cytometry and studied in cell culture assays. Colon contents from mice were analyzed by culture-based microbiology assays.

RESULTS

Loss of S100A8 and S100A9 in mice altered the phenotypes of colonic lamina propria macrophages, compared with wild-type mice. Intestinal tissues from neonatal S100-knockout mice had reduced levels of CX3CR1 protein, and Il10 and Tgfb1 mRNAs, compared with wild-type mice, and fewer T-regulatory cells. S100-knockout mice weighed 21% more than wild-type mice at age 8 weeks and a higher proportion developed fatal sepsis during the neonatal period. S100-knockout mice had alterations in their fecal microbiomes, with higher abundance of Enterobacteriaceae. Feeding mice S100 at birth prevented the expansion of Enterobacteriaceae, increased numbers of T-regulatory cells and levels of CX3CR1 protein and Il10 mRNA in intestine tissues, and reduced body weight and death from neonatal sepsis. Fecal samples from term infants, but not preterm infants, had significantly higher levels of S100A8-A9 during the first 3 months of life than fecal samples from adults; levels decreased to adult levels after weaning. Fecal samples from infants born by cesarean delivery had lower levels of S100A8-A9 than from infants born by vaginal delivery. S100 proteins were expressed by lamina propria macrophages in intestinal tissues from infants, at higher levels than in intestinal tissues from adults. High fecal levels of S100 proteins, from 30 days to 1 year of age, were associated with higher abundance of Actinobacteria and Bifidobacteriaceae, and lower abundance of Gammaproteobacteria-particularly opportunistic Enterobacteriaceae. A low level of S100 proteins in infants' fecal samples associated with development of sepsis and obesity by age 2 years.

CONCLUSION

S100A8 and S100A9 regulate development of the intestinal microbiota and immune system in neonates. Nutritional supplementation with these proteins might aide in development of preterm infants and prevent microbiota-associated disorders in later years.

Longitudinal Multi-omics Analyses Identify Responses of Megakaryocytes, Erythroid Cells, and Plasmablasts as Hallmarks of Severe COVID-19

Temporal resolution of cellular features associated with a severe COVID-19 disease trajectory is needed for understanding skewed immune responses and defining predictors of outcome. Here, we performed a longitudinal multi-omics study using a two-center cohort of 14 patients. We analyzed the bulk transcriptome, bulk DNA methylome, and single-cell transcriptome (>358,000 cells, including BCR profiles) of peripheral blood samples harvested from up to 5 time points. Validation was performed in two independent cohorts of COVID-19 patients. Severe COVID-19 was characterized by an increase of proliferating, metabolically hyperactive plasmablasts. Coinciding with critical illness, we also identified an expansion of interferon-activated circulating megakaryocytes and increased erythropoiesis with features of hypoxic signaling. Megakaryocyte- and erythroid-cell-derived co-expression modules were predictive of fatal disease outcome. The study demonstrates broad cellular effects of SARS-CoV-2 infection beyond adaptive immune cells and provides an entry point toward developing biomarkers and targeted treatments of patients with COVID-19.

•

SARS-CoV2 infection elicits dynamic changes of circulating cells in the blood

•

Severe COVID-19 is characterized by increased metabolically active plasmablasts

•

Elevation of IFN-activated megakaryocytes and erythroid cells in severe COVID-19

•

Cell-type-specific expression signatures are associated with a fatal COVID-19 outcome

Transcriptional Signature Derived from Murine Tumor-Associated Macrophages Correlates with Poor Outcome in Breast Cancer Patients

Tumor-associated macrophages (TAMs) are frequently the most abundant immune cells in cancers and are associated with poor survival. Here, we generated TAM molecular signatures from K14cre;Cdh1^flox/flox;Trp53^flox/flox (KEP) and MMTV-NeuT (NeuT) transgenic mice that resemble human invasive lobular carcinoma (ILC) and HER2⁺ tumors, respectively. Determination of TAM-specific signatures requires comparison with healthy mammary tissue macrophages to avoid overestimation of gene expression differences. TAMs from the two models feature a distinct transcriptomic profile, suggesting that the cancer subtype dictates their phenotype. The KEP-derived signature reliably correlates with poor overall survival in ILC but not in triple-negative breast cancer patients, indicating that translation of murine TAM signatures to patients is cancer subtype dependent. Collectively, we show that a transgenic mouse tumor model can yield a TAM signature relevant for human breast cancer outcome prognosis and provide a generalizable strategy for determining and applying immune cell signatures provided the murine model reflects the human disease.

•

Murine TAM signatures prognosticate outcomes in corresponding cancer patients

•

TAM signatures are robust when they are compared with healthy tissue macrophages

•

TAM transcriptome is dictated by tissue and tumor subtype-related signals

•

Murine TAM signatures can be translated only when a suitable model is chosen

Severe COVID-19 Is Marked by a Dysregulated Myeloid Cell Compartment

Coronavirus disease 2019 (COVID-19) is a mild to moderate respiratory tract infection, however, a subset of patients progress to severe disease and respiratory failure. The mechanism of protective immunity in mild forms and the pathogenesis of severe COVID-19 associated with increased neutrophil counts and dysregulated immune responses remain unclear. In a dual-center, two-cohort study, we combined single-cell RNA-sequencing and single-cell proteomics of whole-blood and peripheral-blood mononuclear cells to determine changes in immune cell composition and activation in mild versus severe COVID-19 (242 samples from 109 individuals) over time. HLA-DRhiCD11chi inflammatory monocytes with an interferon-stimulated gene signature were elevated in mild COVID-19. Severe COVID-19 was marked by occurrence of neutrophil precursors, as evidence of emergency myelopoiesis, dysfunctional mature neutrophils, and HLA-DRlo monocytes. Our study provides detailed insights into the systemic immune response to SARS-CoV-2 infection and reveals profound alterations in the myeloid cell compartment associated with severe COVID-19.

CD163 expression defines specific, IRF8-dependent, immune-modulatory macrophages in the bone marrow

BACKGROUND

Scavenger receptor CD163 is exclusively expressed on monocytes/macrophages and is widely used as a marker for alternatively activated macrophages. However, the role of CD163 is not yet clear.

OBJECTIVES

We sought to examine the function of CD163 in steady-state as well as in sterile and infectious inflammation.

METHODS

Expression of CD163 was analyzed under normal and inflammatory conditions in mice. Functional relevance of CD163 was investigated in models of inflammation in wild-type and CD163^-/- mice.

RESULTS

We describe a subpopulation of bone marrow-resident macrophages (BMRMs) characterized by a high expression of CD163 and functionally distinct from classical bone marrow-derived macrophages. Development of CD163⁺ BMRMs is strictly dependent on IFN regulatory factor-8. CD163⁺ BMRMs show a specific transcriptome and cytokine secretion pattern demonstrating a specific immunomodulatory profile of these cells. Accordingly, CD163^-/- mice show a stronger inflammation in allergic contact dermatitis, indicating a regulatory role of CD163. However, CD163^-/- mice are highly susceptible to S aureus infections, demonstrating the relevance of CD163 for antimicrobial defense as well.

CONCLUSIONS

Our data indicate that anti-inflammatory and immunosuppressive mechanisms are not necessarily associated with a decreased antimicrobial activity. In contrast, our data define a novel macrophage population that controls overwhelming inflammation on one hand but is also necessary for an effective control of infections on the other hand.

Tumor endothelial cell up-regulation of IDO1 is an immunosuppressive feed-back mechanism that reduces the response to CD40-stimulating immunotherapy

CD40-stimulating immunotherapy can elicit potent anti-tumor responses by activating dendritic cells and enhancing T-cell priming. Tumor vessels orchestrate T-cell recruitment during immune response, but the effect of CD40-stimulating immunotherapy on tumor endothelial cells has not been evaluated. Here, we have investigated how tumor endothelial cells transcriptionally respond to CD40-stimulating immunotherapy by isolating tumor endothelial cells from agonistic CD40 mAb- or isotype-treated mice bearing B16-F10 melanoma, and performing RNA-sequencing. Gene set enrichment analysis revealed that agonistic CD40 mAb therapy increased interferon (IFN)-related responses in tumor endothelial cells, including up-regulation of the immunosuppressive enzyme Indoleamine 2, 3-Dioxygenase 1 (IDO1). IDO1 was predominantly expressed in endothelial cells within the tumor microenvironment, and its expression in tumor endothelium was positively correlated to T-cell infiltration and to increased intratumoral expression of IFNγ. In vitro, endothelial cells up-regulated IDO1 in response to T-cell-derived IFNγ, but not in response to CD40-stimulation. Combining agonistic CD40 mAb therapy with the IDO1 inhibitor epacadostat delayed tumor growth in B16-F10 melanoma, associated with increased activation of tumor-infiltrating T-cells. Hereby, we show that the tumor endothelial cells up-regulate IDO1 upon CD40-stimulating immunotherapy in response to increased IFNγ-secretion by T-cells, revealing a novel immunosuppressive feedback mechanism whereby tumor vessels limit T-cell activation.

Cxcr4 distinguishes HSC-derived monocytes from microglia and reveals monocyte immune responses to experimental stroke

Monocyte-derived and tissue-resident macrophages are ontogenetically distinct components of the innate immune system. Assessment of their respective functions in pathology is complicated by changes to the macrophage phenotype during inflammation. Here we find that Cxcr4-CreER enables permanent genetic labeling of hematopoietic stem cells (HSCs) and distinguishes HSC-derived monocytes from microglia and other tissue-resident macrophages. By combining Cxcr4-CreER-mediated lineage tracing with Cxcr4 inhibition or conditional Cxcr4 ablation in photothrombotic stroke, we find that Cxcr4 promotes initial monocyte infiltration and subsequent territorial restriction of monocyte-derived macrophages to infarct tissue. After transient focal ischemia, Cxcr4 deficiency reduces monocyte infiltration and blunts the expression of pattern recognition and defense response genes in monocyte-derived macrophages. This is associated with an altered microglial response and deteriorated outcomes. Thus, Cxcr4 is essential for an innate-immune-system-mediated defense response after cerebral ischemia. We further propose Cxcr4-CreER as a universal tool to study functions of HSC-derived cells.

Modeling population heterogeneity from microbial communities to immune response in cells

Heterogeneity is universally observed in all natural systems and across multiple scales. Understanding population heterogeneity is an intriguing and attractive topic of research in different disciplines, including microbiology and immunology. Microbes and mammalian immune cells present obviously rather different system-specific biological features. Nevertheless, as typically occurs in science, similar methods can be used to study both types of cells. This is particularly true for mathematical modeling, in which key features of a system are translated into algorithms to challenge our mechanistic understanding of the underlying biology. In this review, we first present a broad overview of the experimental developments that allowed observing heterogeneity at the single cell level. We then highlight how this "data revolution" requires the parallel advancement of algorithms and computing infrastructure for data processing and analysis, and finally present representative examples of computational models of population heterogeneity, from microbial communities to immune response in cells.

Differential Gene Expression in Circulating CD14+ Monocytes Indicates the Prognosis of Critically Ill Patients with Sepsis

Critical illness and sepsis are characterized by drastic changes in the systemic innate immune response, particularly involving monocytes. The exact monocyte activation profile during sepsis, however, has remained obscure. Therefore, we prospectively analyzed the gene expression profile of circulating CD14⁺ monocytes from healthy volunteers (n = 54) and intensive care unit (ICU) patients (n = 76), of which n = 36 had sepsis. RNA sequencing of selected samples revealed that monocytes from septic ICU patients display a peculiar activation pattern, which resembles characteristic functional stages of monocyte-derived macrophages and is distinct from controls or non-sepsis ICU patients. Focusing on 55 highly variable genes selected for further investigation, arachidonate 5-lipoxygenase-activating protein (ALOX5AP) was highly upregulated in monocytes of ICU patients and only normalized during 7 days in the ICU in non-sepsis patients. Strikingly, low monocytic guanine nucleotide exchange factor 10-like protein (ARHGEF10L) mRNA expression was associated with the disease severity and mortality of ICU patients. Collectively, our comprehensive analysis of circulating monocytes in critically ill patients revealed a distinct activation pattern, particularly in ICU patients with sepsis. The association with disease severity, the longitudinal recovery or lack thereof during the ICU stay, and the association with prognosis indicate the clinical relevance of monocytic gene expression profiles during sepsis.

Scaling Genomics Data Processing with Memory-Driven Computing to Accelerate Computational Biology

Research is increasingly becoming data-driven, and natural sciences are not an exception. In both biology and medicine, we are observing an exponential growth of structured data collections from experiments and population studies, enabling us to gain novel insights that would otherwise not be possible. However, these growing data sets pose a challenge for existing compute infrastructures since data is outgrowing limits within compute. In this work, we present the application of a novel approach, Memory-Driven Computing (MDC), in the life sciences. MDC proposes a data-centric approach that has been designed for growing data sizes and provides a composable infrastructure for changing workloads. In particular, we show how a typical pipeline for genomics data processing can be accelerated, and application modifications required to exploit this novel architecture. Furthermore, we demonstrate how the isolated evaluation of individual tasks misses significant overheads of typical pipelines in genomics data processing.

Shiny-Seq: advanced guided transcriptome analysis

Objective

A comprehensive analysis of RNA-Seq data uses a wide range of different tools and algorithms, which are normally limited to R users only. While several tools and advanced analysis pipelines are available, some require programming skills and others lack the support for many important features that enable a more comprehensive data analysis. There is thus, a need for a guided and easy to use comprehensive RNA-Seq data platform, which integrates the state of the art analysis workflow.

Results

We present the tool Shiny-Seq, which provides a guided and easy to use comprehensive RNA-Seq data analysis pipeline. It has many features such as batch effect estimation and removal, quality check with several visualization options, enrichment analysis with multiple biological databases, identification of patterns using advanced methods such as weighted gene co-expression network analysis, summarizing analysis as power point presentation and all results as tables via a one-click feature. The source code is published on GitHub (https://github.com/schultzelab/Shiny-Seq) and licensed under GPLv3. Shiny-Seq is written in R using the Shiny framework. In addition, the application is hosted on a public website hosted by the shinyapps.io server (https://schultzelab.shinyapps.io/Shiny-Seq/) and as a Docker image https://hub.docker.com/r/makaho/shiny-seq.

Membrane Cholesterol Efflux Drives Tumor-Associated Macrophage Reprogramming and Tumor Progression

Macrophages possess intrinsic tumoricidal activity, yet tumor-associated macrophages (TAMs) rapidly adopt an alternative phenotype within the tumor microenvironment that is marked by tumor-promoting immunosuppressive and trophic functions. The mechanisms that promote such TAM polarization remain poorly understood, but once identified, they may represent important therapeutic targets to block the tumor-promoting functions of TAMs and restore their anti-tumor potential. Here, we have characterized TAMs in a mouse model of metastatic ovarian cancer. We show that ovarian cancer cells promote membrane-cholesterol efflux and depletion of lipid rafts from macrophages. Increased cholesterol efflux promoted IL-4-mediated reprogramming, including inhibition of IFNγ-induced gene expression. Genetic deletion of ABC transporters, which mediate cholesterol efflux, reverts the tumor-promoting functions of TAMs and reduces tumor progression. These studies reveal an unexpected role for membrane-cholesterol efflux in driving TAM-mediated tumor progression while pointing to a potentially novel anti-tumor therapeutic strategy.

Enzymatic Activity of HPGD in Treg Cells Suppresses Tconv Cells to Maintain Adipose Tissue Homeostasis and Prevent Metabolic Dysfunction

Regulatory T cells (Treg cells) are important for preventing autoimmunity and maintaining tissue homeostasis, but whether Treg cells can adopt tissue- or immune-context-specific suppressive mechanisms is unclear. Here, we found that the enzyme hydroxyprostaglandin dehydrogenase (HPGD), which catabolizes prostaglandin E₂ (PGE₂) into the metabolite 15-keto PGE₂, was highly expressed in Treg cells, particularly those in visceral adipose tissue (VAT). Nuclear receptor peroxisome proliferator-activated receptor-γ (PPARγ)-induced HPGD expression in VAT Treg cells, and consequential Treg-cell-mediated generation of 15-keto PGE₂ suppressed conventional T cell activation and proliferation. Conditional deletion of Hpgd in mouse Treg cells resulted in the accumulation of functionally impaired Treg cells specifically in VAT, causing local inflammation and systemic insulin resistance. Consistent with this mechanism, humans with type 2 diabetes showed decreased HPGD expression in Treg cells. These data indicate that HPGD-mediated suppression is a tissue- and context-dependent suppressive mechanism used by Treg cells to maintain adipose tissue homeostasis.

Therapeutic targeting of macrophages enhances chemotherapy efficacy by unleashing type I interferon response

Recent studies have revealed a role for macrophages and neutrophils in limiting chemotherapy efficacy; however, the mechanisms underlying the therapeutic benefit of myeloid-targeting agents in combination with chemotherapy are incompletely understood. Here, we show that targeting tumour-associated macrophages by colony-stimulating factor-1 receptor (CSF-1R) blockade in the K14cre;Cdh1^F/F;Trp53^F/F transgenic mouse model for breast cancer stimulates intratumoural type I interferon (IFN) signalling, which enhances the anticancer efficacy of platinum-based chemotherapeutics. Notably, anti-CSF-1R treatment also increased intratumoural expression of type I IFN-stimulated genes in patients with cancer, confirming that CSF-1R blockade is a powerful strategy to trigger an intratumoural type I IFN response. By inducing an inflamed, type I IFN-enriched tumour microenvironment and by further targeting immunosuppressive neutrophils during cisplatin therapy, antitumour immunity was activated in this poorly immunogenic breast cancer mouse model. These data illustrate the importance of breaching multiple layers of immunosuppression during cytotoxic therapy to successfully engage antitumour immunity in breast cancer.

Transcriptional profiling of human microglia reveals grey-white matter heterogeneity and multiple sclerosis-associated changes

Here we report the transcriptional profile of human microglia, isolated from normal-appearing grey matter (GM) and white matter (WM) of multiple sclerosis (MS) and non-neurological control donors, to find possible early changes related to MS pathology. Microglia show a clear region-specific profile, indicated by higher expression of type-I interferon genes in GM and higher expression of NF-κB pathway genes in WM. Transcriptional changes in MS microglia also differ between GM and WM. MS WM microglia show increased lipid metabolism gene expression, which relates to MS pathology since active MS lesion-derived microglial nuclei show similar altered gene expression. Microglia from MS GM show increased expression of genes associated with glycolysis and iron homeostasis, possibly reflecting microglia reacting to iron depositions. Except for ADGRG1/GPR56, expression of homeostatic genes, such as P2RY12 and TMEM119, is unaltered in normal-appearing MS tissue, demonstrating overall preservation of microglia homeostatic functions in the initiation phase of MS.

Expression of the Phosphatase Ppef2 Controls Survival and Function of CD8+ Dendritic Cells

Apoptotic cell death of Dendritic cells (DCs) is critical for immune homeostasis. Although intrinsic mechanisms controlling DC death have not been fully characterized up to now, experimentally enforced inhibition of DC-death causes various autoimmune diseases in model systems. We have generated mice deficient for Protein Phosphatase with EF-Hands 2 (Ppef2), which is selectively expressed in CD8⁺ DCs, but not in other related DC subtypes such as tissue CD103⁺ DCs. Ppef2 is down-regulated rapidly upon maturation of DCs by toll-like receptor stimuli, but not upon triggering of CD40. Ppef2-deficient CD8⁺ DCs accumulate the pro-apoptotic Bcl-2-like protein 11 (Bim) and show increased apoptosis and reduced competitve repopulation capacities. Furthermore, Ppef2^−/− CD8⁺ DCs have strongly diminished antigen presentation capacities in vivo, as CD8⁺ T cells primed by Ppef2^−/− CD8⁺ DCs undergo reduced expansion. In conclusion, our data suggests that Ppef2 is crucial to support survival of immature CD8⁺ DCs, while Ppef2 down-regulation during DC-maturation limits T cell responses.

CD83 expression is essential for Treg cell differentiation and stability

Foxp3-positive regulatory T cells (Tregs) are crucial for the maintenance of immune homeostasis and keep immune responses in check. Upon activation, Tregs are transferred into an effector state expressing transcripts essential for their suppressive activity, migration, and survival. However, it is not completely understood how different intrinsic and environmental factors control differentiation. Here, we present for the first time to our knowledge data suggesting that Treg-intrinsic expression of CD83 is essential for Treg differentiation upon activation. Interestingly, mice with Treg-intrinsic CD83 deficiency are characterized by a proinflammatory phenotype. Furthermore, the loss of CD83 expression by Tregs leads to the downregulation of Treg-specific differentiation markers and the induction of an inflammatory profile. In addition, Treg-specific conditional knockout mice showed aggravated autoimmunity and an impaired resolution of inflammation. Altogether, our results show that CD83 expression in Tregs is an essential factor for the development and function of effector Tregs upon activation. Since Tregs play a crucial role in the maintenance of immune tolerance and thus prevention of autoimmune disorders, our findings are also clinically relevant.

Exposure to the gut microbiota drives distinct methylome and transcriptome changes in intestinal epithelial cells during postnatal development

BACKGROUND

The interplay of epigenetic processes and the intestinal microbiota may play an important role in intestinal development and homeostasis. Previous studies have established that the microbiota regulates a large proportion of the intestinal epithelial transcriptome in the adult host, but microbial effects on DNA methylation and gene expression during early postnatal development are still poorly understood. Here, we sought to investigate the microbial effects on DNA methylation and the transcriptome of intestinal epithelial cells (IECs) during postnatal development.

METHODS

We collected IECs from the small intestine of each of five 1-, 4- and 12 to 16-week-old mice representing the infant, juvenile, and adult states, raised either in the presence or absence of a microbiota. The DNA methylation profile was determined using reduced representation bisulfite sequencing (RRBS) and the epithelial transcriptome by RNA sequencing using paired samples from each individual mouse to analyze the link between microbiota, gene expression, and DNA methylation.

RESULTS

We found that microbiota-dependent and -independent processes act together to shape the postnatal development of the transcriptome and DNA methylation signatures of IECs. The bacterial effect on the transcriptome increased over time, whereas most microbiota-dependent DNA methylation differences were detected already early after birth. Microbiota-responsive transcripts could be attributed to stage-specific cellular programs during postnatal development and regulated gene sets involved primarily immune pathways and metabolic processes. Integrated analysis of the methylome and transcriptome data identified 126 genomic loci at which coupled differential DNA methylation and RNA transcription were associated with the presence of intestinal microbiota. We validated a subset of differentially expressed and methylated genes in an independent mouse cohort, indicating the existence of microbiota-dependent "functional" methylation sites which may impact on long-term gene expression signatures in IECs.

CONCLUSIONS

Our study represents the first genome-wide analysis of microbiota-mediated effects on maturation of DNA methylation signatures and the transcriptional program of IECs after birth. It indicates that the gut microbiota dynamically modulates large portions of the epithelial transcriptome during postnatal development, but targets only a subset of microbially responsive genes through their DNA methylation status.

Human glioblastoma-associated microglia/monocytes express a distinct RNA profile compared to human control and murine samples

Glioblastoma (GBM) is the most aggressive brain tumor in adults. It is strongly infiltrated by microglia and peripheral monocytes that support tumor growth. In the present study we used RNA sequencing to compare the expression profile of CD11b(+) human glioblastoma-associated microglia/monocytes (hGAMs) to CD11b(+) microglia isolated from non-tumor samples. Hierarchical clustering and principal component analysis showed a clear separation of the two sample groups and we identified 334 significantly regulated genes in hGAMs. In comparison to human control microglia hGAMs upregulated genes associated with mitotic cell cycle, cell migration, cell adhesion, and extracellular matrix organization. We validated the expression of several genes associated with extracellular matrix organization in samples of human control microglia, hGAMs, and the hGAMs-depleted fraction via qPCR. The comparison to murine GAMs (mGAMs) showed that both cell populations share a significant fraction of upregulated transcripts compared with their respective controls. These genes were mostly related to mitotic cell cycle. However, in contrast to murine cells, human GAMs did not upregulate genes associated to immune activation. Comparison of human and murine GAMs expression data to several data sets of in vitro-activated human macrophages and murine microglia showed that, in contrast to mGAMs, hGAMs share a smaller overlap to these data sets in general and in particular to cells activated by proinflammatory stimulation with LPS + INFγ or TNFα. Our findings provide new insights into the biology of human glioblastoma-associated microglia/monocytes and give detailed information about the validity of murine experimental models. GLIA 2016 GLIA 2016;64:1416-1436.

Bioinformatic Assessment of Macrophage Activation by the Innate Immune System

Weighted gene co-expression network analysis (WGCNA) allows for the identification and characterization of cell type-specific gene modules in complex transcriptome datasets. Here, we use a microarray dataset of human macrophages comprising 29 conditions and 299 samples generated by differentiation of CD14⁺ monocytes into macrophages followed by in vitro stimulations to identify stimulation-specific gene modules. These gene modules can be used for experimental validation, as well as further bioinformatic analysis to determine key pathways or upstream transcription factors.