tRigon: an R package and Shiny App for integrative (path-)omics data analysis.

BACKGROUND: Pathomics facilitates automated, reproducible and precise histopathology analysis and morphological phenotyping. Similar to molecular omics, pathomics datasets are high-dimensional, but also face large outlier variability and inherent data missingness, making quick and comprehensible data analysis challenging. To facilitate pathomics data analysis and interpretation as well as support a broad implementation we developed tRigon (Toolbox foR InteGrative (path-)Omics data aNalysis), a Shiny application for fast, comprehensive and reproducible pathomics analysis. RESULTS: tRigon is available via the CRAN repository ( https://cran.r-project.org/web/packages/tRigon ) with its source code available on GitLab ( https://git-ce.rwth-aachen.de/labooratory-ai/trigon ). The tRigon package can be installed locally and its application can be executed from the R console via the command 'tRigon::run_tRigon()'. Alternatively, the application is hosted online and can be accessed at https://labooratory.shinyapps.io/tRigon . We show fast computation of small, medium and large datasets in a low- and high-performance hardware setting, indicating broad applicability of tRigon. CONCLUSIONS: tRigon allows researchers without coding abilities to perform exploratory feature analyses of pathomics and non-pathomics datasets on their own using a variety of hardware.

Lineage tracing reveals transient phenotypic adaptation of tubular cells during acute kidney injury.

Tubular injury is the hallmark of acute kidney injury (AKI) with a tremendous impact on patients and health-care systems. During injury, any differentiated proximal tubular cell (PT) may transition into a specific injured phenotype, so-called "scattered tubular cell" (STC)-phenotype. To understand the fate of this specific phenotype, we generated transgenic mice allowing inducible, reversible, and irreversible tagging of these cells in a murine AKI model, the unilateral ischemia-reperfusion injury (IRI). For lineage tracing, we analyzed the kidneys using single-cell profiling during disease development at various time points. Labeled cells, which we defined by established endogenous markers, already appeared 8 h after injury and showed a distinct expression set of genes. We show that STCs re-differentiate back into fully differentiated PTs upon the resolution of the injury. In summary, we show the dynamics of the phenotypic transition of PTs during injury, revealing a reversible transcriptional program as an adaptive response during disease.

Proinflammatory phenotype of iPS cell-derived JAK2 V617F megakaryocytes induces fibrosis in 3D in vitro bone marrow niche.

The myeloproliferative disease polycythemia vera (PV) driven by the JAK2 V617F mutation can transform into myelofibrosis (post-PV-MF). It remains an open question how JAK2 V617F in hematopoietic stem cells induces MF. Megakaryocytes are major players in murine PV models but are difficult to study in the human setting. We generated induced pluripotent stem cells (iPSCs) from JAK2 V617F PV patients and differentiated them into megakaryocytes. In differentiation assays, JAK2 V617F iPSCs recapitulated the pathognomonic skewed megakaryocytic and erythroid differentiation. JAK2 V617F iPSCs had a TPO-independent and increased propensity to differentiate into megakaryocytes. RNA sequencing of JAK2 V617F iPSC-derived megakaryocytes reflected a proinflammatory, profibrotic phenotype and decreased ribosome biogenesis. In three-dimensional (3D) coculture, JAK2 V617F megakaryocytes induced a profibrotic phenotype through direct cell contact, which was reversed by the JAK2 inhibitor ruxolitinib. The 3D coculture system opens the perspective for further disease modeling and drug discovery.

Detection of PatIent-Level distances from single cell genomics and pathomics data with Optimal Transport (PILOT).

Although clinical applications represent the next challenge in single-cell genomics and digital pathology, we still lack computational methods to analyze single-cell or pathomics data to find sample-level trajectories or clusters associated with diseases. This remains challenging as single-cell/pathomics data are multi-scale, i.e., a sample is represented by clusters of cells/structures, and samples cannot be easily compared with each other. Here we propose PatIent Level analysis with Optimal Transport (PILOT). PILOT uses optimal transport to compute the Wasserstein distance between two individual single-cell samples. This allows us to perform unsupervised analysis at the sample level and uncover trajectories or cellular clusters associated with disease progression. We evaluate PILOT and competing approaches in single-cell genomics or pathomics studies involving various human diseases with up to 600 samples/patients and millions of cells or tissue structures. Our results demonstrate that PILOT detects disease-associated samples from large and complex single-cell or pathomics data. Moreover, PILOT provides a statistical approach to find changes in cell populations, gene expression, and tissue structures related to the trajectories or clusters supporting interpretation of predictions.

Non-canonical Hedgehog signaling mediates profibrotic hematopoiesis-stroma crosstalk in myeloproliferative neoplasms.

The role of hematopoietic Hedgehog signaling in myeloproliferative neoplasms (MPNs) remains incompletely understood despite data suggesting that Hedgehog (Hh) pathway inhibitors have therapeutic activity in patients. We aim to systematically interrogate the role of canonical vs. non-canonical Hh signaling in MPNs. We show that Gli1 protein levels in patient peripheral blood mononuclear cells (PBMCs) mark fibrotic progression and that, in murine MPN models, absence of hematopoietic Gli1, but not Gli2 or Smo, significantly reduces MPN phenotype and fibrosis, indicating that GLI1 in the MPN clone can be activated in a non-canonical fashion. Additionally, we establish that hematopoietic Gli1 has a significant effect on stromal cells, mediated through a druggable MIF-CD74 axis. These data highlight the complex interplay between alterations in the MPN clone and activation of stromal cells and indicate that Gli1 represents a promising therapeutic target in MPNs, particularly that Hh signaling is dispensable for normal hematopoiesis.

Enhanced fatty acid oxidation through metformin and baicalin as therapy for COVID-19 and associated inflammatory states in lung and kidney.

Progressive respiratory failure is the primary cause of death in the coronavirus disease 2019 (COVID-19) pandemic. It is the final outcome of the acute respiratory distress syndrome (ARDS), characterized by an initial exacerbated inflammatory response, metabolic derangement and ultimate tissue scarring. A positive balance of cellular energy may result crucial for the recovery of clinical COVID-19. Hence, we asked if two key pathways involved in cellular energy generation, AMP-activated protein kinase (AMPK)/acetyl-CoA carboxylase (ACC) signaling and fatty acid oxidation (FAO) could be beneficial. We tested the drugs metformin (AMPK activator) and baicalin (CPT1A activator) in different experimental models mimicking COVID-19 associated inflammation in lung and kidney. We also studied two different cohorts of COVID-19 patients that had been previously treated with metformin. These drugs ameliorated lung damage in an ARDS animal model, while activation of AMPK/ACC signaling increased mitochondrial function and decreased TGF-ß-induced fibrosis, apoptosis and inflammation markers in lung epithelial cells. Similar results were observed with two indole derivatives, IND6 and IND8 with AMPK activating capacity. Consistently, a reduced time of hospitalization and need of intensive care was observed in COVID-19 patients previously exposed to metformin. Baicalin also mitigated the activation of pro-inflammatory bone marrow-derived macrophages (BMDMs) and reduced kidney fibrosis in two animal models of kidney injury, another key target of COVID-19. In human epithelial lung and kidney cells, both drugs improved mitochondrial function and prevented TGF-ß-induced renal epithelial cell dedifferentiation. Our results support that favoring cellular energy production through enhanced FAO may prove useful in the prevention of COVID-19-induced lung and renal damage.

Mapping cardiac remodeling in chronic kidney disease.

Patients with advanced chronic kidney disease (CKD) mostly die from sudden cardiac death and recurrent heart failure. The mechanisms of cardiac remodeling are largely unclear. To dissect molecular and cellular mechanisms of cardiac remodeling in CKD in an unbiased fashion, we performed left ventricular single-nuclear RNA sequencing in two mouse models of CKD. Our data showed a hypertrophic response trajectory of cardiomyocytes with stress signaling and metabolic changes driven by soluble uremia-related factors. We mapped fibroblast to myofibroblast differentiation in this process and identified notable changes in the cardiac vasculature, suggesting inflammation and dysfunction. An integrated analysis of cardiac cellular responses to uremic toxins pointed toward endothelin-1 and methylglyoxal being involved in capillary dysfunction and TNFα driving cardiomyocyte hypertrophy in CKD, which was validated in vitro and in vivo. TNFα inhibition in vivo ameliorated the cardiac phenotype in CKD. Thus, interventional approaches directed against uremic toxins, such as TNFα, hold promise to ameliorate cardiac remodeling in CKD.

EphrinA5 regulates cell motility by modulating Snhg15/DNA triplex-dependent targeting of DNMT1 to the Ncam1 promoter.

Cell-cell communication is mediated by membrane receptors and their ligands, such as the Eph/ephrin system, orchestrating cell migration during development and in diverse cancer types. Epigenetic mechanisms are key for integrating external "signals", e.g., from neighboring cells, into the transcriptome in health and disease. Previously, we reported ephrinA5 to trigger transcriptional changes of lncRNAs and protein-coding genes in cerebellar granule cells, a cell model for medulloblastoma. LncRNAs represent important adaptors for epigenetic writers through which they regulate gene expression. Here, we investigate a lncRNA-mediated targeting of DNMT1 to specific gene loci by the combined power of in silico modeling of RNA/DNA interactions and wet lab approaches, in the context of the clinically relevant use case of ephrinA5-dependent regulation of cellular motility of cerebellar granule cells. We provide evidence that Snhg15, a cancer-related lncRNA, recruits DNMT1 to the Ncam1 promoter through RNA/DNA triplex structure formation and the interaction with DNMT1. This mediates DNA methylation-dependent silencing of Ncam1, being abolished by ephrinA5 stimulation-triggered reduction of Snhg15 expression. Hence, we here propose a triple helix recognition mechanism, underlying cell motility regulation via lncRNA-targeted DNA methylation in a clinically relevant context.

Targeted DNA Methylation Analysis Facilitates Leukocyte Counts in Dried Blood Samples.

BACKGROUND: Cell-type specific DNA methylation (DNAm) can be employed to determine the numbers of leukocyte subsets in blood. In contrast to conventional methods for leukocyte counts, which are based on cellular morphology or surface marker protein expression, the cellular deconvolution based on DNAm levels is applicable for frozen or dried blood. Here, we further enhanced targeted DNAm assays for leukocyte counts in clinical application. METHODS: DNAm profiles of 40 different studies were compiled to identify CG dinucleotides (CpGs) with cell-type specific DNAm using a computational framework, CimpleG. DNAm levels at these CpGs were then measured with digital droplet PCR in venous blood from 160 healthy donors and 150 patients with various hematological disorders. Deconvolution was further validated with venous blood (n = 75) and capillary blood (n = 31) that was dried on Whatman paper or on Mitra microsampling devices. RESULTS: In venous blood, automated cell counting or flow cytometry correlated well with epigenetic estimates of relative leukocyte counts for granulocytes (r = 0.95), lymphocytes (r = 0.97), monocytes (r = 0.82), CD4 T cells (r = 0.84), CD8 T cells (r = 0.94), B cells (r = 0.96), and NK cells (r = 0.72). Similar correlations and precisions were achieved for dried blood samples. Spike-in with a reference plasmid enabled accurate epigenetic estimation of absolute leukocyte counts from dried blood samples, correlating with conventional venous (r = 0.86) and capillary (r = 0.80) blood measurements. CONCLUSIONS: The advanced selection of cell-type specific CpGs and utilization of digital droplet PCR analysis provided accurate epigenetic blood counts. Analysis of dried blood facilitates self-sampling with a finger prick, thereby enabling easier accessibility to testing.

Myelofibrosis at diagnosis is associated with the failure of treatment-free remission in CML patients.

The management of patients with chronic myeloid leukemia (CML) has been revolutionized by the introduction of tyrosine kinase inhibitors (TKIs), which induce deep molecular responses so that treatment can eventually be discontinued, leading to treatment-free remission (TFR) in a subset of patients. Unfortunately, leukemic stem cells (LSCs) often persist and a fraction of these can again expand in about half of patients that attempt TKI discontinuation. In this study, we show that presence of myelofibrosis (MF) at the time of diagnosis is a factor associating with TFR failure. Fibrotic transformation is governed by the action of several cytokines, and interestingly, some of them have also been described to support LSC persistence. At the cellular level, these could be produced by both malignant cells and by components of the bone marrow (BM) niche, including megakaryocytes (MKs) and mesenchymal stromal cells (MSCs). In our cohort of 57 patients, around 40% presented with MF at diagnosis and the number of blasts in the peripheral blood and BM was significantly elevated in patients with higher grade of MF. Employing a CML transgenic mouse model, we could observe higher levels of alpha-smooth muscle actin (α-SMA) in the BM when compared to control mice. Short-term treatment with the TKI nilotinib, efficiently reduced spleen weight and BCR::ABL1 mRNA levels, while α-SMA expression was only partially reduced. Interestingly, the number of MKs was increased in the spleen of CML mice and elevated in both BM and spleen upon nilotinib treatment. Analysis of human CML-vs healthy donor (HD)-derived MSCs showed an altered expression of gene signatures reflecting fibrosis as well as hematopoietic support, thus suggesting MSCs as a potential player in these two processes. Finally, in our cohort, 12 patients qualified for TKI discontinuation, and here we observed that all patients who failed TFR had BM fibrosis at diagnosis, whereas this was only the case in 25% of patients with achieved TFR, further supporting the link between fibrosis and LSC persistence.

Response of Soybean Cultivars to Oxidative Stress caused by Meloidogyne javanica.

The present study aimed to investigate the response of soybean cultivars with different susceptibility levels to the root-knot nematode Meloidogyne javanica at varied time intervals by analyzing the initial plant-nematode interaction using antioxidant enzymes as oxidative stress markers. A 4 × 4 × 2 factorial method with 5 repetitions was used to analyze 4 soybean cultivars at 4 different collection times-6, 12, 24, and 48 h-with and without M. javanica inoculation. The parameters evaluated were the activities of antioxidant enzymes phenol peroxidase (POX) and ascorbate peroxidase (APX); the concentrations of hydrogen peroxide (H2O2) and malondialdehyde (MDA); and the number of M. javanica juveniles penetrated into each plant. H2O2 concentration varied among the cultivars with and without inoculation and at different collection times as indicated by MDA concentration and POX and APX activities, demonstrating a rapid response of the host to an infection by M. javanica. Oxidative stress caused by M. javanica did not vary among the soybean cultivars regardless of their susceptibility level; however, the antioxidant enzymes POX and APX responded according to the susceptibility level of the cultivars.

CimpleG: finding simple CpG methylation signatures.

DNA methylation signatures are usually based on multivariate approaches that require hundreds of sites for predictions. Here, we propose a computational framework named CimpleG for the detection of small CpG methylation signatures used for cell-type classification and deconvolution. We show that CimpleG is both time efficient and performs as well as top performing methods for cell-type classification of blood cells and other somatic cells, while basing its prediction on a single DNA methylation site per cell type. Altogether, CimpleG provides a complete computational framework for the delineation of DNAm signatures and cellular deconvolution.

Avaliação da cinética on da frequência cardíaca em pacientes asmáticos e em indivíduos saudáveis durante o Endurance Shuttle Walk Test / Evaluation of heart rate on kinetics in asthmatic patients and in healthy individuals during the endurance shuttle walk test

Introdução: O comportamento cinético da frequência cardíaca (FC) na transição do repouso para o exercício nos informa sobre a integridade do sistema nervoso autônomo. Recuperações mais lentas associam-se ao risco de mortalidade por eventos cardiovasculares, tornando-se imprescindível sua avaliação. Objetivo: Avaliar e comparar a resposta da cinética on da FC em pacientes asmáticos e indivíduos saudáveis durante o Endurance Shuttle Walk Test (ESWT). Métodos: Trata-se de um estudo prospectivo, transversal e controlado, com 14 adultos asmáticos e 8 controles saudáveis. Os indivíduos realizaram as seguintes avaliações: Teste de função pulmonar, Variabilidade da Frequência Cardíaca (VFC) e Incremental Shuttle Walk Test e ESWT. Resultados: O grupo asmático apresentou um atraso da cinética on da FC na transição do repouso para o teste, e uma correlação negativa moderada (r=-0,60; p<0,05) entre a distância percorrida (m) e o tempo de resposta (TRM) cinética on da FC. Conclusão: Os pacientes asmáticos apresentaram um atraso da cinética "on", quando comparados ao grupo de indivíduos saudáveis, demonstrando ser um importante marcador na avaliação da performance física.

Introduction: The kinetic behavior of heart rate (HR) in the transition from rest to exercise, as this assessment informs us about the integrity of the autonomic nervous system. Slower recoveries are associated with the risk of mortality from cardiovascular events, making their evaluation, essential. Objective: To evaluate and compare the HR on kinetics response in asthmatic patients and healthy individuals during the Endurance Shuttle Walk Test (ESWT). Methods: This is a prospective, cross-sectional, controlled study with 14 asthmatic adults and 8 healthy controls. Subjects performed the following assessments: Pulmonary Function Test, Heart Rate Variability (HRV) and Incremental Shuttle Walk Test and ESWT. Results: The asthmatic group showed a delay in the HR on kinetics in the transition from rest to the test, and a moderate negative correlation (r=-0.60; p<0.05) between the distance covered (m) and the response time (TRM) kinetics on from FC. Conclusion: Asthmatic patients showed a delay in "on" kinetics, in comparison to the group of healthy individuals, proving to be an important marker in physical performance assessments.

M cell maturation and cDC activation determine the onset of adaptive immune priming in the neonatal Peyer's patch.

Early-life immune development is critical to long-term host health. However, the mechanisms that determine the pace of postnatal immune maturation are not fully resolved. Here, we analyzed mononuclear phagocytes (MNPs) in small intestinal Peyer's patches (PPs), the primary inductive site of intestinal immunity. Conventional type 1 and 2 dendritic cells (cDC1 and cDC2) and RORgt+ antigen-presenting cells (RORgt+ APC) exhibited significant age-dependent changes in subset composition, tissue distribution, and reduced cell maturation, subsequently resulting in a lack in CD4+ T cell priming during the postnatal period. Microbial cues contributed but could not fully explain the discrepancies in MNP maturation. Type I interferon (IFN) accelerated MNP maturation but IFN signaling did not represent the physiological stimulus. Instead, follicle-associated epithelium (FAE) M cell differentiation was required and sufficient to drive postweaning PP MNP maturation. Together, our results highlight the role of FAE M cell differentiation and MNP maturation in postnatal immune development.

Mast cell deficiency prevents BCR::ABL1 induced splenomegaly and cytokine elevation in a CML mouse model.

The persistence of leukemic stem cells (LSCs) represents a problem in the therapy of chronic myeloid leukemia (CML). Hence, it is of utmost importance to explore the underlying mechanisms to develop new therapeutic approaches to cure CML. Using the genetically engineered ScltTA/TRE-BCR::ABL1 mouse model for chronic phase CML, we previously demonstrated that the loss of the docking protein GAB2 counteracts the infiltration of mast cells (MCs) in the bone marrow (BM) of BCR::ABL1 positive mice. Here, we show for the first time that BCR::ABL1 drives the cytokine independent expansion of BM derived MCs and sensitizes them for FcεRI triggered degranulation. Importantly, we demonstrate that genetic mast cell deficiency conferred by the Cpa3Cre allele prevents BCR::ABL1 induced splenomegaly and impairs the production of pro-inflammatory cytokines. Furthermore, we show in CML patients that splenomegaly is associated with high BM MC counts and that upregulation of pro-inflammatory cytokines in patient serum samples correlates with tryptase levels. Finally, MC-associated transcripts were elevated in human CML BM samples. Thus, our study identifies MCs as essential contributors to disease progression and suggests considering them as an additional target in CML therapy. Mast cells play a key role in the pro-inflammatory tumor microenvironment of the bone marrow. Shown is a cartoon summarizing our results from the mouse model. BCR::ABL1 transformed MCs, as part of the malignant clone, are essential for the elevation of pro-inflammatory cytokines, known to be important in disease initiation and progression.

RGT: a toolbox for the integrative analysis of high throughput regulatory genomics data.

BACKGROUND: Massive amounts of data are produced by combining next-generation sequencing with complex biochemistry techniques to characterize regulatory genomics profiles, such as protein-DNA interaction and chromatin accessibility. Interpretation of such high-throughput data typically requires different computation methods. However, existing tools are usually developed for a specific task, which makes it challenging to analyze the data in an integrative manner. RESULTS: We here describe the Regulatory Genomics Toolbox (RGT), a computational library for the integrative analysis of regulatory genomics data. RGT provides different functionalities to handle genomic signals and regions. Based on that, we developed several tools to perform distinct downstream analyses, including the prediction of transcription factor binding sites using ATAC-seq data, identification of differential peaks from ChIP-seq data, and detection of triple helix mediated RNA and DNA interactions, visualization, and finding an association between distinct regulatory factors. CONCLUSION: We present here RGT; a framework to facilitate the customization of computational methods to analyze genomic data for specific regulatory genomics problems. RGT is a comprehensive and flexible Python package for analyzing high throughput regulatory genomics data and is available at: https://github.com/CostaLab/reg-gen . The documentation is available at: https://reg-gen.readthedocs.io.

KIT D816V Mast Cells Derived from Induced Pluripotent Stem Cells Recapitulate Systemic Mastocytosis Transcriptional Profile.

Mast cells (MCs) represent a population of hematopoietic cells with a key role in innate and adaptive immunity and are well known for their detrimental role in allergic responses. Yet, MCs occur in low abundance, which hampers their detailed molecular analysis. Here, we capitalized on the potential of induced pluripotent stem (iPS) cells to give rise to all cells in the body and established a novel and robust protocol for human iPS cell differentiation toward MCs. Relying on a panel of systemic mastocytosis (SM) patient-specific iPS cell lines carrying the KIT D816V mutation, we generated functional MCs that recapitulate SM disease features: increased number of MCs, abnormal maturation kinetics and activated phenotype, CD25 and CD30 surface expression and a transcriptional signature characterized by upregulated expression of innate and inflammatory response genes. Therefore, human iPS cell-derived MCs are a reliable, inexhaustible, and close-to-human tool for disease modeling and pharmacological screening to explore novel MC therapeutics.

A lncRNA identifies Irf8 enhancer element in negative feedback control of dendritic cell differentiation.

Transcription factors play a determining role in lineage commitment and cell differentiation. Interferon regulatory factor 8 (IRF8) is a lineage determining transcription factor in hematopoiesis and master regulator of dendritic cells (DC), an important immune cell for immunity and tolerance. IRF8 is prominently upregulated in DC development by autoactivation and controls both DC differentiation and function. However, it is unclear how Irf8 autoactivation is controlled and eventually limited. Here, we identified a novel long non-coding RNA transcribed from the +32 kb enhancer downstream of Irf8 transcription start site and expressed specifically in mouse plasmacytoid DC (pDC), referred to as lncIrf8. The lncIrf8 locus interacts with the lrf8 promoter and shows differential epigenetic signatures in pDC versus classical DC type 1 (cDC1). Interestingly, a sequence element of the lncIrf8 promoter, but not lncIrf8 itself, is crucial for mouse pDC and cDC1 differentiation, and this sequence element confers feedback inhibition of Irf8 expression. Taken together, in DC development Irf8 autoactivation is first initiated by flanking enhancers and then second controlled by feedback inhibition through the lncIrf8 promoter element in the +32 kb enhancer. Our work reveals a previously unrecognized negative feedback loop of Irf8 that orchestrates its own expression and thereby controls DC differentiation.