Motivation: DNA enrichment followed by sequencing is a versatile tool in molecular biology, with a wide variety of applications including genome-wide analysis of epigenetic marks and mechanisms. A common requirement of these diverse applications is a comparison of read coverage between experimental conditions. The amount of samples generated for such comparisons ranges from few replicates to hundreds of samples per condition for epigenome-wide association studies. Consequently, there is an urgent need for software that allows for fast and simple processing and comparison of sequencing data derived from enriched DNA. Results: Here, we present a major update of the R/Bioconductor package MEDIPS, which allows for an arbitrary number of replicates per group and integrates sophisticated statistical methods for the detection of differential coverage between experimental conditions. Our approach can be applied to a diversity of quantitative sequencing data. In addition, our update adds novel functionality to MEDIPS, including correlation analysis between samples, and takes advantage of Bioconductor’s annotation databases to facilitate annotation of specific genomic regions. Availability and implementation: The latest version of MEDIPS is available as version 1.12.0 and part of Bioconductor 2.13. The package comes with a manual containing detailed description of its functionality and is available at http://www.bioconductor.org. Contact:lienhard@molgen.mpg.de Supplementary information: Supplementary data are available at Bioinformatics online.

Recent developments in immuno-oncology demonstrate that not only cancer cells, but also the tumor microenvironment can guide precision medicine. A comprehensive and in-depth characterization of the tumor microenvironment is challenging since its cell populations are diverse and can be important even if scarce. To identify clinically relevant microenvironmental and cancer features, we applied single-cell RNA sequencing to ten human lung adenocarcinomas and ten normal control tissues. Our analyses revealed heterogeneous carcinoma cell transcriptomes reflecting histological grade and oncogenic pathway activities, and two distinct microenvironmental patterns. The immune-activated CP²E microenvironment was composed of cancer-associated myofibroblasts, proinflammatory monocyte-derived macrophages, plasmacytoid dendritic cells and exhausted CD8+ T cells, and was prognostically unfavorable. In contrast, the inert N³MC microenvironment was characterized by normal-like myofibroblasts, non-inflammatory monocyte-derived macrophages, NK cells, myeloid dendritic cells and conventional T cells, and was associated with a favorable prognosis. Microenvironmental marker genes and signatures identified in single-cell profiles had progonostic value in bulk tumor profiles. In summary, single-cell RNA profiling of lung adenocarcinoma provides additional prognostic information based on the microenvironment, and may help to predict therapy response and to reveal possible target cell populations for future therapeutic approaches.

Viruses manipulate cellular metabolism and macromolecule recycling processes like autophagy. Dysregulated metabolism might lead to excessive inflammatory and autoimmune responses as observed in severe and long COVID-19 patients. Here we show that SARS-CoV-2 modulates cellular metabolism and reduces autophagy. Accordingly, compound-driven induction of autophagy limits SARS-CoV-2 propagation. In detail, SARS-CoV-2-infected cells show accumulation of key metabolites, activation of autophagy inhibitors (AKT1, SKP2) and reduction of proteins responsible for autophagy initiation (AMPK, TSC2, ULK1), membrane nucleation, and phagophore formation (BECN1, VPS34, ATG14), as well as autophagosome-lysosome fusion (BECN1, ATG14 oligomers). Consequently, phagophore-incorporated autophagy markers LC3B-II and P62 accumulate, which we confirm in a hamster model and lung samples of COVID-19 patients. Single-nucleus and single-cell sequencing of patient-derived lung and mucosal samples show differential transcriptional regulation of autophagy and immune genes depending on cell type, disease duration, and SARS-CoV-2 replication levels. Targeting of autophagic pathways by exogenous administration of the polyamines spermidine and spermine, the selective AKT1 inhibitor MK-2206, and the BECN1-stabilizing anthelmintic drug niclosamide inhibit SARS-CoV-2 propagation in vitro with IC₅₀ values of 136.7, 7.67, 0.11, and 0.13 μM, respectively. Autophagy-inducing compounds reduce SARS-CoV-2 propagation in primary human lung cells and intestinal organoids emphasizing their potential as treatment options against COVID-19.

Wound healing of the skin is a crucial regenerative process in adult mammals. We examined wound healing in conditional mutant mice, in which the c-Met gene that encodes the receptor of hepatocyte growth factor/scatter factor was mutated in the epidermis by cre recombinase. c-Met–deficient keratinocytes were unable to contribute to the reepithelialization of skin wounds. In conditional c-Met mutant mice, wound closure was slightly attenuated, but occurred exclusively by a few (5%) keratinocytes that had escaped recombination. This demonstrates that the wound process selected and amplified residual cells that express a functional c-Met receptor. We also cultured primary keratinocytes from the skin of conditional c-Met mutant mice and examined them in scratch wound assays. Again, closure of scratch wounds occurred by the few remaining c-Met–positive cells. Our data show that c-Met signaling not only controls cell growth and migration during embryogenesis but is also essential for the generation of the hyperproliferative epithelium in skin wounds, and thus for a fundamental regenerative process in the adult.

In colorectal cancer, oncogenic mutations transform a hierarchically organized and homeostatic epithelium into invasive cancer tissue lacking visible organization. We sought to define transcriptional states of colorectal cancer cells and signals controlling their development by performing single‐cell transcriptome analysis of tumors and matched non‐cancerous tissues of twelve colorectal cancer patients. We defined patient‐overarching colorectal cancer cell clusters characterized by differential activities of oncogenic signaling pathways such as mitogen‐activated protein kinase and oncogenic traits such as replication stress. RNA metabolic labeling and assessment of RNA velocity in patient‐derived organoids revealed developmental trajectories of colorectal cancer cells organized along a mitogen‐activated protein kinase activity gradient. This was in contrast to normal colon organoid cells developing along graded Wnt activity. Experimental targeting of EGFR‐BRAF‐MEK in cancer organoids affected signaling and gene expression contingent on predictive KRAS/BRAF mutations and induced cell plasticity overriding default developmental trajectories. Our results highlight directional cancer cell development as a driver of non‐genetic cancer cell heterogeneity and re‐routing of trajectories as a response to targeted therapy.

Gene expression profiling of β-catenin, Cripto and Wnt3 mutant mouse embryos has been used to characterise the genetic networks that regulate early embryonic development. We have defined genes whose expression is regulated by β-catenin during formation of the anteroposterior axis and the mesoderm, and have identified Cripto, which encodes a Nodal co-receptor, as a primary target of β-catenin signals both in embryogenesis as well as in colon carcinoma cell lines and tissues. We have also defined groups of genes regulated by Wnt3/β-catenin signalling during primitive streak and mesoderm formation. Our data assign a key role to β-catenin upstream of two distinct gene expression programs during anteroposterior axis and mesoderm formation.