Cell-free DNA methylation-based preeclampsia prediction: A journey to improve maternal health

Presymptomatic prediction of preeclampsia (PE) is crucial to enable early prophylactic treatment. Current screening tools are either complex or lack predictive value. We recently demonstrated that cell-free DNA methylation can be leveraged to predict early-onset PE in 57% at a 10% false positive rate. Importantly, this minimally invasive screening test can be implemented as an add-on to current widespread noninvasive prenatal aneuploidy screening. Here, we highlight the pitfalls and promising prospects of this research.

Evaluating offspring Genomic and Epigenomic alterations after prenatal exposure to Cancer treatment In Pregnancy (GE-CIP): a multicentric observational study

INTRODUCTION

Around 1 in 1000-2000 pregnancies are affected by a cancer diagnosis. Previous studies have shown that chemotherapy during pregnancy has reassuring cognitive and cardiac neonatal outcomes, and hence has been proposed as standard of care. However, although these children perform within normal ranges for their age, subtle differences have been identified. Given that chemotherapeutic compounds can cross the placenta, the possibility that prenatal chemotherapy exposure mutates the offspring's genome and/or epigenome, with potential deleterious effects later in life, urges to be investigated.

METHODS AND ANALYSES

This multicentric observational study aims to collect cord blood, meconium and neonatal buccal cells at birth, as well as peripheral blood, buccal cells and urine from infants when 6, 18 and/or 36 months of age. Using bulk and single-cell approaches, we will compare samples from chemotherapy-treated pregnant patients with cancer, pregnant patients with cancer not treated with chemotherapy and healthy pregnant women. Potential chemotherapy-related newborn genomic and/or epigenomic alterations, such as single nucleotide variants, copy number variants and DNA-methylation alterations, will be identified in mononuclear and epithelial cells, isolated from blood, buccal swabs and urine. DNA from maternal peripheral blood and paternal buccal cells will be used to determine de novo somatic mutations in the neonatal blood and epithelial cells. Additionally, the accumulated exposure of the fetus, and biological effective dose of alkylating agents, will be assessed in meconium and cord blood via mass spectrometry approaches.

ETHICS AND DISSEMINATION

The Ethics Committee Research of UZ/KU Leuven (EC Research) and the Medical Ethical Review Committee of University Medical Center Amsterdam have approved the study. Results of this study will be disseminated via presentations at (inter)national conferences, through peer-reviewed, open-access publications, via social media platforms aimed to inform patients and healthcare workers, and through the website of the International Network on Cancer, Infertility and Pregnancy (www.cancerinpregnancy.org).

Cell type signatures in cell-free DNA fragmentation profiles reveal disease biology

Circulating cell-free DNA (cfDNA) fragments have characteristics that are specific to the cell types that release them. Current methods for cfDNA deconvolution typically use disease tailored marker selection in a limited number of bulk tissues or cell lines. Here, we utilize single cell transcriptome data as a comprehensive cellular reference set for disease-agnostic cfDNA cell-of-origin analysis. We correlate cfDNA-inferred nucleosome spacing with gene expression to rank the relative contribution of over 490 cell types to plasma cfDNA. In 744 healthy individuals and patients, we uncover cell type signatures in support of emerging disease paradigms in oncology and prenatal care. We train predictive models that can differentiate patients with colorectal cancer (84.7%), early-stage breast cancer (90.1%), multiple myeloma (AUC 95.0%), and preeclampsia (88.3%) from matched controls. Importantly, our approach performs well in ultra-low coverage cfDNA datasets and can be readily transferred to diverse clinical settings for the expansion of liquid biopsy.

Epigenetic remodelling under hypoxia

Hypoxia is intrinsic to tumours and contributes to malignancy and metastasis while hindering the efficiency of existing treatments. Epigenetic mechanisms play a crucial role in the regulation of hypoxic cancer cell programs, both in the initial phases of sensing the decrease in oxygen levels and during adaptation to chronic lack of oxygen. During the latter, the epigenetic regulation of tumour biology intersects with hypoxia-sensitive transcription factors in a complex network of gene regulation that also involves metabolic reprogramming. Here, we review the current literature on the epigenetic control of gene programs in hypoxic cancer cells. We highlight common themes and features of such epigenetic remodelling and discuss their relevance for the development of therapeutic strategies.

The serine synthesis pathway drives osteoclast differentiation through epigenetic regulation of NFATc1 expression

Bone-resorbing osteoclasts are vital for postnatal bone health, as increased differentiation or activity results in skeletal pathologies such as osteoporosis. The metabolism of mature osteoclasts differs from their progenitor cells, but whether the observed metabolic changes are secondary to the altered cell state or actively drive the process of cell differentiation is unknown. Here, we show that transient activation of the serine synthesis pathway (SSP) is essential for osteoclastogenesis, as deletion of the rate-limiting enzyme phosphoglycerate dehydrogenase in osteoclast progenitors impairs their differentiation and results in increased bone mass. In addition, pharmacological phosphoglycerate dehydrogenase inhibition abrogated bone loss in a mouse model of postmenopausal osteoporosis by blocking bone resorption. Mechanistically, SSP-derived α-ketoglutarate is necessary for histone demethylases that remove repressive histone methylation marks at the nuclear factor of activated T cells, cytoplasmic 1 (Nfatc1) gene locus, thereby inducing NFATc1 expression and consequent osteoclast maturation. Taken together, this study reveals a metabolic-epigenetic coupling mechanism that directs osteoclast differentiation and suggests that the SSP can be therapeutically targeted to prevent osteoporotic bone loss.

A Metabolic Gene Survey Pinpoints Fucosylation as a Key Pathway Underlying the Suppressive Function of Regulatory T Cells in Cancer

Forkhead box P3 (Foxp3)-expressing regulatory T cells (Treg) are the guardians of controlled immune reactions and prevent the development of autoimmune diseases. However, in the tumor context, their increased number suppresses antitumor immune responses, indicating the importance of understanding the mechanisms behind their function and stability. Metabolic reprogramming can affect Foxp3 regulation and, therefore, Treg suppressive function and fitness. Here, we performed a metabolic CRISPR/Cas9 screen and pinpointed novel candidate positive and negative metabolic regulators of Foxp3. Among the positive regulators, we revealed that targeting the GDP-fucose transporter Slc35c1, and more broadly fucosylation (Fuco), in Tregs compromises their proliferation and suppressive function both in vitro and in vivo, leading to alteration of the tumor microenvironment and impaired tumor progression and protumoral immune responses. Pharmacologic inhibition of Fuco dampened tumor immunosuppression mostly by targeting Tregs, thus resulting in reduced tumor growth. In order to substantiate these findings in humans, tumoral Tregs from patients with colorectal cancer were clustered on the basis of the expression of Fuco-related genes. FucoLOW Tregs were found to exhibit a more immunogenic profile compared with FucoHIGH Tregs. Furthermore, an enrichment of a FucoLOW signature, mainly derived from Tregs, correlated with better prognosis and response to immune checkpoint blockade in melanoma patients. In conclusion, Slc35c1-dependent Fuco is able to regulate the suppressive function of Tregs, and measuring its expression in Tregs might pave the way towards a useful biomarker model for patients with cancer. See related Spotlight by Silveria and DuPage, p. 1570.

GiRAFR improves gRNA detection and annotation in single-cell CRISPR screens

Novel methods that combine single cell RNA-seq with CRISPR screens enable high-throughput characterization of transcriptional changes caused by genetic perturbations. Dedicated software is however lacking to annotate CRISPR guide RNA (gRNA) libraries and associate them with single cell transcriptomes. Here, we describe a CRISPR droplet sequencing (CROP-seq) dataset. During analysis, we observed that the most commonly used method fails to detect mutant gRNAs. We therefore developed a python tool to identify and characterize intact and mutant gRNAs, called GiRAFR. We show that mutant gRNAs are dysfunctional, and failure to detect and annotate them leads to an inflated estimate of the number of untransformed cells, attenuated downregulation of target genes, as well as an underestimated multiplet frequency. These findings are mirrored in publicly available datasets, where we find that up to 35% of cells are transduced with a mutant gRNA. Applying GiRAFR hence stands to improve the annotation and quality of single cell CRISPR screens.

Cell-free DNA methylome analysis for early preeclampsia prediction

Preeclampsia (PE) is a leading cause for peripartal morbidity, especially if developing early in gestation. To enable prophylaxis in the prevention of PE, pregnancies at risk of PE must be identified early-in the first trimester. To identify at-risk pregnancies we profiled methylomes of plasma-derived, cell-free DNA from 498 pregnant women, of whom about one-third developed early-onset PE. We detected DNA methylation differences between control and PE pregnancies that enabled risk stratification at PE diagnosis but also presymptomatically, at around 12 weeks of gestation (range 9-14 weeks). The first-trimester risk prediction model was validated in an external cohort collected from two centers (area under the curve (AUC) = 0.75) and integrated with routinely available maternal risk factors (AUC = 0.85). The combined risk score correctly predicted 72% of patients with early-onset PE at 80% specificity. These preliminary results suggest that cell-free DNA methylation profiling is a promising tool for presymptomatic PE risk assessment, and has the potential to improve treatment and follow-up in the obstetric clinic.

Methods and applications for single-cell and spatial multi-omics

The joint analysis of the genome, epigenome, transcriptome, proteome and/or metabolome from single cells is transforming our understanding of cell biology in health and disease. In less than a decade, the field has seen tremendous technological revolutions that enable crucial new insights into the interplay between intracellular and intercellular molecular mechanisms that govern development, physiology and pathogenesis. In this Review, we highlight advances in the fast-developing field of single-cell and spatial multi-omics technologies (also known as multimodal omics approaches), and the computational strategies needed to integrate information across these molecular layers. We demonstrate their impact on fundamental cell biology and translational research, discuss current challenges and provide an outlook to the future.

Metabolic profiling stratifies colorectal cancer and reveals adenosylhomocysteinase as a therapeutic target

The genomic landscape of colorectal cancer (CRC) is shaped by inactivating mutations in tumour suppressors such as APC, and oncogenic mutations such as mutant KRAS. Here we used genetically engineered mouse models, and multimodal mass spectrometry-based metabolomics to study the impact of common genetic drivers of CRC on the metabolic landscape of the intestine. We show that untargeted metabolic profiling can be applied to stratify intestinal tissues according to underlying genetic alterations, and use mass spectrometry imaging to identify tumour, stromal and normal adjacent tissues. By identifying ions that drive variation between normal and transformed tissues, we found dysregulation of the methionine cycle to be a hallmark of APC-deficient CRC. Loss of Apc in the mouse intestine was found to be sufficient to drive expression of one of its enzymes, adenosylhomocysteinase (AHCY), which was also found to be transcriptionally upregulated in human CRC. Targeting of AHCY function impaired growth of APC-deficient organoids in vitro, and prevented the characteristic hyperproliferative/crypt progenitor phenotype driven by acute deletion of Apc in vivo, even in the context of mutant Kras. Finally, pharmacological inhibition of AHCY reduced intestinal tumour burden in ApcMin/+ mice indicating its potential as a metabolic drug target in CRC.

Heterogeneity of Repolarization and Cell-Cell Variability of Cardiomyocyte Remodeling Within the Myocardial Infarction Border Zone Contribute to Arrhythmia Susceptibility

BACKGROUND

After myocardial infarction, the infarct border zone (BZ) is the dominant source of life-threatening arrhythmias, where fibrosis and abnormal repolarization create a substrate for reentry. We examined whether repolarization abnormalities are heterogeneous within the BZ in vivo and could be related to heterogeneous cardiomyocyte remodeling.

METHODS

Myocardial infarction was induced in domestic pigs by 120-minute ischemia-reperfusion injury. After 1 month, remodeling was assessed by magnetic resonance imaging, and electroanatomical mapping was performed to determine the spatial distribution of activation-recovery intervals. Cardiomyocytes were isolated and tissue samples collected from the BZ and remote regions. Optical recording allowed assessment of action potential duration (di-8-Anepps, stimulation at 1 Hz, 37 °C) of large cardiomyocyte populations while gene expression in cardiomyocytes was determined by single nuclear RNA sequencing.

RESULTS

In vivo, activation-recovery intervals in the BZ tended to be longer than in remote with increased spatial heterogeneity evidenced by a greater local SD (3.5±1.3 ms versus remote: 2.0±0.5 ms, P=0.036, n_pigs=5). Increased activation-recovery interval heterogeneity correlated with enhanced arrhythmia susceptibility. Cellular population studies (n_cells=635-862 cells per region) demonstrated greater heterogeneity of action potential duration in the BZ (SD, 105.9±17.0 ms versus remote: 73.9±8.6 ms; P=0.001; n_pigs=6), which correlated with heterogeneity of activation-recovery interval in vivo. Cell-cell gene expression heterogeneity in the BZ was evidenced by increased Euclidean distances between nuclei of the BZ (12.1 [9.2-15.0] versus 10.6 [7.5-11.6] in remote; P<0.0001). Differentially expressed genes characterizing BZ cardiomyocyte remodeling included hypertrophy-related and ion channel-related genes with high cell-cell variability of expression. These gene expression changes were driven by stress-responsive TFs (transcription factors). In addition, heterogeneity of left ventricular wall thickness was greater in the BZ than in remote.

CONCLUSIONS

Heterogeneous cardiomyocyte remodeling in the BZ is driven by uniquely altered gene expression, related to heterogeneity in the local microenvironment, and translates to heterogeneous repolarization and arrhythmia vulnerability in vivo.

Lupus, DNA Methylation, and Air Pollution: A Malicious Triad

The pathogenesis of systemic lupus erythematosus (SLE) remains elusive to this day; however, genetic, epigenetic, and environmental factors have been implicated to be involved in disease pathogenesis. Recently, it was demonstrated that in systemic lupus erythematosus (SLE) patients, interferon-regulated genes are hypomethylated in naïve CD4⁺ T cells, CD19⁺ B lymphocytes, and CD14⁺ monocytes. This suggests that interferon-regulated genes may have been epigenetically poised in SLE patients for rapid expression upon stimulation by different environmental factors. Additionally, environmental studies have identified DNA (hypo)methylation changes as a potential mechanism of environmentally induced health effects in utero, during childhood and in adults. Finally, epidemiologic studies have firmly established air pollution as a crucial SLE risk factor, as studies showed an association between fine particulate matter (PM_2.5) and traditional SLE biomarkers related to disease flare, hospital admissions, and an increased SLEDAI score. In this review, the relationship between aberrant epigenetic regulation, the environment, and the development of SLE will be discussed.

The pulmonary vasculature in lethal COVID-19 and idiopathic pulmonary fibrosis at single cell resolution

Aims

SARS-CoV-2 infection causes COVID-19, which in severe cases evokes life-threatening acute respiratory distress syndrome (ARDS). Transcriptome signatures and the functional relevance of non-vascular cell types (e.g. immune and epithelial cells) in COVID-19 are becoming increasingly evident. However, despite its known contribution to vascular inflammation, recruitment/invasion of immune cells, vascular leakage and perturbed hemostasis in the lungs of severe COVID-19 patients, an in-depth interrogation of the endothelial cell (EC) compartment in lethal COVID-19 is lacking. Moreover, progressive fibrotic lung disease represents one of the complications of COVID-19 pneumonia and ARDS. Analogous features between idiopathic pulmonary fibrosis (IPF) and COVID-19 suggest partial similarities in their pathophysiology, yet, a head-to-head comparison of pulmonary cell transcriptomes between both conditions has not been implemented to date.

Methods and Results

We performed single nucleus RNA-seq (snRNA-seq) on frozen lungs from 7 deceased COVID-19 patients, 6 IPF explant lungs and 12 controls. The vascular fraction, comprising 38,794 nuclei, could be subclustered into 14 distinct EC subtypes. Non-vascular cell types, comprising 137,746 nuclei, were subclustered and used for EC-interactome analyses. Pulmonary ECs of deceased COVID-19 patients showed an enrichment of genes involved in cellular stress, as well as signatures suggestive of dampened immunomodulation and impaired vessel wall integrity. In addition, increased abundance of a population of systemic capillary and venous ECs was identified in COVID-19 and IPF. COVID-19 systemic ECs closely resembled their IPF counterparts, and a set of 30 genes was found congruently enriched in systemic ECs across studies. Receptor-ligand interaction analysis of ECs with non-vascular cell types in the pulmonary micro-environment revealed numerous previously unknown interactions specifically enriched/depleted in COVID-19 and/or IPF.

Conclusions

This study uncovered novel insights into the abundance, expression patterns and interactomes of EC subtypes in COVID-19 and IPF, relevant for future investigations into the progression and treatment of both lethal conditions.

Translational perspective

While assessing clinical and molecular characteristics of severe and lethal COVID-19 cases, the vasculature’s undeniable role in disease progression has been widely acknowledged. COVID-19 lung pathology moreover shares certain clinical features with late-stage IPF – yet an in-depth interrogation and direct comparison of the endothelium at single-cell level in both conditions is still lacking. By comparing the transcriptomes of ECs from lungs of deceased COVID-19 patients to those from IPF explant and control lungs, we gathered key insights the heterogeneous composition and potential roles of ECs in both lethal diseases, which may serve as a foundation for development of novel therapeutics.

Skeletal progenitors preserve proliferation and self-renewal upon inhibition of mitochondrial respiration by rerouting the TCA cycle

A functional electron transport chain (ETC) is crucial for supporting bioenergetics and biosynthesis. Accordingly, ETC inhibition decreases proliferation in cancer cells but does not seem to impair stem cell proliferation. However, it remains unclear how stem cells metabolically adapt. In this study, we show that pharmacological inhibition of complex III of the ETC in skeletal stem and progenitor cells induces glycolysis side pathways and reroutes the tricarboxylic acid (TCA) cycle to regenerate NAD⁺ and preserve cell proliferation. These metabolic changes also culminate in increased succinate and 2-hydroxyglutarate levels that inhibit Ten-eleven translocation (TET) DNA demethylase activity, thereby preserving self-renewal and multilineage potential. Mechanistically, mitochondrial malate dehydrogenase and reverse succinate dehydrogenase activity proved to be essential for the metabolic rewiring in response to ETC inhibition. Together, these data show that the metabolic plasticity of skeletal stem and progenitor cells allows them to bypass ETC blockade and preserve their self-renewal.

•

Skeletal stem/progenitor cells can proliferate upon electron transport chain blockade

•

Succinate dehydrogenase is reversed with fumarate functioning as electron acceptor

•

Pyruvate and aspartate are critical for NAD⁺ regeneration and proliferation

•

Metabolic changes prevent DNA demethylation and preserve self-renewal

More of less: Novel multi-ome profiling of single human neurons

Epigenetic modifications to DNA and chromatin interact to influence gene expression and cellular phenotypes, but defining these omics layers in complex tissues is a daunting task. In this issue of Cell Genomics, Luo et al. describe a novel single-cell multi-omic method, simultaneously profiling transcriptome, DNA methylome, and chromatin accessibility, to shed light on human neurons.

Novel next-generation sequencing-based methodologies to characterize the mutational consequences of (prenatal) chemotherapy exposure in noncancerous tissue

PURPOSE OF REVIEW

Although chemotherapeutics are considered as genotoxins for decades, their exact mutagenic impact on the genome of cancerous and normal cells of cancer patients was unknown for a long time. However, this knowledge is necessary to understand the long-term side effects of chemotherapy. A particular condition represents pregnant cancer patients being treated with chemotherapy. Since certain chemotherapeutics can cross the placenta, concerns exist about possible mutational effects on the fetus' genome with potential long-term health consequences.

RECENT FINDINGS

Recent advances of next-generation sequencing (NGS) techniques have opened possibilities to explore the exact mutational footprint of chemotherapies in healthy tissue from treated cancer patients. However, the ultra-low frequency of chemotherapy-induced mutations, introduction of technical artefacts, and inaccessibility of normal tissue has posed important limitations. This review discusses five state-of-the-art approaches that were recently designed to overcome these drawbacks.

SUMMARY

Results of the latest investigations give valuable insights into the genome-wide genotoxicity profile of frequently applied chemotherapies, with most of these drugs being associated with a signature of random base substitutions and small indels. Though these findings still might be limited to extrapolate to healthy tissue, they pave the way for studies on the origin of long-term chemotherapy-related adverse health effects.

Correction: Integrative analysis of genomic variants reveals new associations of candidate haploinsufficient genes with congenital heart disease

[This corrects the article DOI: 10.1371/journal.pgen.1009679.].

Integrative analysis of genomic variants reveals new associations of candidate haploinsufficient genes with congenital heart disease

Numerous genetic studies have established a role for rare genomic variants in Congenital Heart Disease (CHD) at the copy number variation (CNV) and de novo variant (DNV) level. To identify novel haploinsufficient CHD disease genes, we performed an integrative analysis of CNVs and DNVs identified in probands with CHD including cases with sporadic thoracic aortic aneurysm. We assembled CNV data from 7,958 cases and 14,082 controls and performed a gene-wise analysis of the burden of rare genomic deletions in cases versus controls. In addition, we performed variation rate testing for DNVs identified in 2,489 parent-offspring trios. Our analysis revealed 21 genes which were significantly affected by rare CNVs and/or DNVs in probands. Fourteen of these genes have previously been associated with CHD while the remaining genes (FEZ1, MYO16, ARID1B, NALCN, WAC, KDM5B and WHSC1) have only been associated in small cases series or show new associations with CHD. In addition, a systems level analysis revealed affected protein-protein interaction networks involved in Notch signaling pathway, heart morphogenesis, DNA repair and cilia/centrosome function. Taken together, this approach highlights the importance of re-analyzing existing datasets to strengthen disease association and identify novel disease genes and pathways.

Regulatory Dynamics of Tet1 and Oct4 Resolve Stages of Global DNA Demethylation and Transcriptomic Changes in Reprogramming

Reprogramming somatic cells into induced pluripotent stem cells (iPSCs) involves the reactivation of endogenous pluripotency genes and global DNA demethylation, but temporal resolution of these events using existing markers is limited. Here, we generate murine transgenic lines harboring reporters for the 5-methylcytosine dioxygenase Tet1 and for Oct4. By monitoring dual reporter fluorescence during pluripotency entry, we identify a sequential order of Tet1 and Oct4 activation by proximal and distal regulatory elements. Full Tet1 activation marks an intermediate stage that accompanies predominantly repression of somatic genes, preceding full Oct4 activation, and distinguishes two waves of global DNA demethylation that target distinct genomic features but are uncoupled from transcriptional changes. Tet1 knockout shows that TET1 contributes to both waves of demethylation and activates germline regulatory genes in reprogramming intermediates but is dispensable for Oct4 reactivation. Our dual reporter system for time-resolving pluripotency entry thus refines the molecular roadmap of iPSC maturation.

Overcoming immunotherapy resistance in non-small cell lung cancer (NSCLC) - novel approaches and future outlook

Immunotherapy (IO) has revolutionized the therapy landscape of non-small cell lung cancer (NSCLC), significantly prolonging the overall survival (OS) of advanced stage patients. Over the recent years IO therapy has been broadly integrated into the first-line setting of non-oncogene driven NSCLC, either in combination with chemotherapy, or in selected patients with PD-L1high expression as monotherapy. Still, a significant proportion of patients suffer from disease progression. A better understanding of resistance mechanisms depicts a central goal to avoid or overcome IO resistance and to improve patient outcome.We here review major cellular and molecular pathways within the tumor microenvironment (TME) that may impact the evolution of IO resistance. We summarize upcoming treatment options after IO resistance including novel IO targets (e.g. RIG-I, STING) as well as interesting combinational approaches such as IO combined with anti-angiogenic agents or metabolic targets (e.g. IDO-1, adenosine signaling, arginase). By discussing the fundamental mode of action of IO within the TME, we aim to understand and manage IO resistance and to seed new ideas for effective therapeutic IO concepts.

DNA methylation repels binding of hypoxia-inducible transcription factors to maintain tumor immunotolerance

BACKGROUND

Hypoxia is pervasive in cancer and other diseases. Cells sense and adapt to hypoxia by activating hypoxia-inducible transcription factors (HIFs), but it is still an outstanding question why cell types differ in their transcriptional response to hypoxia.

RESULTS

We report that HIFs fail to bind CpG dinucleotides that are methylated in their consensus binding sequence, both in in vitro biochemical binding assays and in vivo studies of differentially methylated isogenic cell lines. Based on in silico structural modeling, we show that 5-methylcytosine indeed causes steric hindrance in the HIF binding pocket. A model wherein cell-type-specific methylation landscapes, as laid down by the differential expression and binding of other transcription factors under normoxia, control cell-type-specific hypoxia responses is observed. We also discover ectopic HIF binding sites in repeat regions which are normally methylated. Genetic and pharmacological DNA demethylation, but also cancer-associated DNA hypomethylation, expose these binding sites, inducing HIF-dependent expression of cryptic transcripts. In line with such cryptic transcripts being more prone to cause double-stranded RNA and viral mimicry, we observe low DNA methylation and high cryptic transcript expression in tumors with high immune checkpoint expression, but not in tumors with low immune checkpoint expression, where they would compromise tumor immunotolerance. In a low-immunogenic tumor model, DNA demethylation upregulates cryptic transcript expression in a HIF-dependent manner, causing immune activation and reducing tumor growth.

CONCLUSIONS

Our data elucidate the mechanism underlying cell-type-specific responses to hypoxia and suggest DNA methylation and hypoxia to underlie tumor immunotolerance.

IGF1R is an entry receptor for respiratory syncytial virus

Pneumonia resulting from infection is one of the leading causes of death worldwide. Pulmonary infection by the respiratory syncytial virus (RSV) is a large burden on human health, for which there are few therapeutic options1. RSV targets ciliated epithelial cells in the airways, but how viruses such as RSV interact with receptors on these cells is not understood. Nucleolin is an entry coreceptor for RSV2 and also mediates the cellular entry of influenza, the parainfluenza virus, some enteroviruses and the bacterium that causes tularaemia3,4. Here we show a mechanism of RSV entry into cells in which outside-in signalling, involving binding of the prefusion RSV-F glycoprotein with the insulin-like growth factor-1 receptor, triggers the activation of protein kinase C zeta (PKCζ). This cellular signalling cascade recruits nucleolin from the nuclei of cells to the plasma membrane, where it also binds to RSV-F on virions. We find that inhibiting PKCζ activation prevents the trafficking of nucleolin to RSV particles on airway organoid cultures, and reduces viral replication and pathology in RSV-infected mice. These findings reveal a mechanism of virus entry in which receptor engagement and signal transduction bring the coreceptor to viral particles at the cell surface, and could form the basis of new therapeutics to treat RSV infection.

A pan-cancer blueprint of the heterogeneous tumor microenvironment revealed by single-cell profiling

The stromal compartment of the tumor microenvironment consists of a heterogeneous set of tissue-resident and tumor-infiltrating cells, which are profoundly moulded by cancer cells. An outstanding question is to what extent this heterogeneity is similar between cancers affecting different organs. Here, we profile 233,591 single cells from patients with lung, colorectal, ovary and breast cancer (n = 36) and construct a pan-cancer blueprint of stromal cell heterogeneity using different single-cell RNA and protein-based technologies. We identify 68 stromal cell populations, of which 46 are shared between cancer types and 22 are unique. We also characterise each population phenotypically by highlighting its marker genes, transcription factors, metabolic activities and tissue-specific expression differences. Resident cell types are characterised by substantial tissue specificity, while tumor-infiltrating cell types are largely shared across cancer types. Finally, by applying the blueprint to melanoma tumors treated with checkpoint immunotherapy and identifying a naïve CD4⁺ T-cell phenotype predictive of response to checkpoint immunotherapy, we illustrate how it can serve as a guide to interpret scRNA-seq data. In conclusion, by providing a comprehensive blueprint through an interactive web server, we generate the first panoramic view on the shared complexity of stromal cells in different cancers.

A new protocol for single-cell RNA-seq reveals stochastic gene expression during lag phase in budding yeast

Current methods for single-cell RNA sequencing (scRNA-seq) of yeast cells do not match the throughput and relative simplicity of the state-of-the-art techniques that are available for mammalian cells. In this study, we report how 10x Genomics' droplet-based single-cell RNA sequencing technology can be modified to allow analysis of yeast cells. The protocol, which is based on in-droplet spheroplasting of the cells, yields an order-of-magnitude higher throughput in comparison to existing methods. After extensive validation of the method, we demonstrate its use by studying the dynamics of the response of isogenic yeast populations to a shift in carbon source, revealing the heterogeneity and underlying molecular processes during this shift. The method we describe opens new avenues for studies focusing on yeast cells, as well as other cells with a degradable cell wall.

Single-Cell RNA Sequencing Maps Endothelial Metabolic Plasticity in Pathological Angiogenesis

Endothelial cell (EC) metabolism is an emerging target for anti-angiogenic therapy in tumor angiogenesis and choroidal neovascularization (CNV), but little is known about individual EC metabolic transcriptomes. By single-cell RNA sequencing 28,337 murine choroidal ECs (CECs) and sprouting CNV-ECs, we constructed a taxonomy to characterize their heterogeneity. Comparison with murine lung tumor ECs (TECs) revealed congruent marker gene expression by distinct EC phenotypes across tissues and diseases, suggesting similar angiogenic mechanisms. Trajectory inference predicted that differentiation of venous to angiogenic ECs was accompanied by metabolic transcriptome plasticity. ECs displayed metabolic transcriptome heterogeneity during cell-cycle progression and in quiescence. Hypothesizing that conserved genes are important, we used an integrated analysis, based on congruent transcriptome analysis, CEC-tailored genome-scale metabolic modeling, and gene expression meta-analysis in cross-species datasets, followed by in vitro and in vivo validation, to identify SQLE and ALDH18A1 as previously unknown metabolic angiogenic targets.

DNA methylation-driven EMT is a common mechanism of resistance to various therapeutic agents in cancer

Background

Overcoming therapeutic resistance is one of the major hurdles in cancer care. One mechanism contributing to therapeutic resistance is a process in which epithelial cells switch to a mesenchymal state (epithelial-to-mesenchymal transition or EMT). The precise mechanisms driving EMT-mediated therapeutic resistance have, however, not been elucidated.

Results

Here, we study ten cell line pairs, for which parental cell lines were made resistant to either a targeted or chemotherapy-based treatment. First, we show by miRNA-200 overexpression that treatment resistance is driven by EMT. Next, we demonstrate that DNA methylation changes occur within each cell line pair and show that exposure to 5-azacytidine or knock down of DNA methyltransferases (DNMTs), both of which globally demethylate cells, result in EMT reversal and increased therapeutic sensitivity. This suggests DNA methylation to causally underlie EMT and treatment resistance. We also observe significant overlap in methylation profiles between resistant lines, suggesting a common epigenetic mechanism to cause resistance to therapy. In line with this hypothesis, cross-resistance to other targeted and chemotherapies is observed, while importantly, this is lost upon demethylation of the cells. Finally, we clinically validate that DNA methylation changes drive EMT-mediated resistance to sorafenib in patients with advanced hepatocellular carcinoma (HCC). Specifically, we develop a capture-based protocol to interrogate DNA methylation in low amounts of circulating tumor DNA (ctDNA). By interrogating the methylation status in liquid biopsies, longitudinally collected during sorafenib treatment, we assess whether DNA methylation changes also drive EMT and therapy resistance in a clinical setting. Particularly, by monitoring methylation changes in EMT genes, we are able to predict tumor response and acquired resistance to sorafenib.

Conclusions

We propose methylation changes underlying EMT to constitute a common resistance mechanism to cancer therapies. This process can be reversed pharmacologically and monitored non-invasively in ctDNA to predict resistance to treatment.

An Integrated Gene Expression Landscape Profiling Approach to Identify Lung Tumor Endothelial Cell Heterogeneity and Angiogenic Candidates

Heterogeneity of lung tumor endothelial cell (TEC) phenotypes across patients, species (human/mouse), and models (in vivo/in vitro) remains poorly inventoried at the single-cell level. We single-cell RNA (scRNA)-sequenced 56,771 endothelial cells from human/mouse (peri)-tumoral lung and cultured human lung TECs, and detected 17 known and 16 previously unrecognized phenotypes, including TECs putatively regulating immune surveillance. We resolved the canonical tip TECs into a known migratory tip and a putative basement-membrane remodeling breach phenotype. Tip TEC signatures correlated with patient survival, and tip/breach TECs were most sensitive to vascular endothelial growth factor blockade. Only tip TECs were congruent across species/models and shared conserved markers. Integrated analysis of the scRNA-sequenced data with orthogonal multi-omics and meta-analysis data across different human tumors, validated by functional analysis, identified collagen modification as a candidate angiogenic pathway.

Genetic biomarkers in the VEGF pathway predicting response to anti-VEGF therapy in age-related macular degeneration

Objective

Age-related macular degeneration (ARMD) is a leading cause of visual impairment. Intravitreal injections of anti-vascular endothelial growth factor (VEGF) are the standard treatment for wet ARMD. There is however, variability in patient responses, suggesting patient-specific factors influencing drug efficacy. We tested whether single nucleotide polymorphisms (SNPs) in genes encoding VEGF pathway members contribute to therapy response.

Methods and analysis

A retrospective cohort of 281 European wet ARMD patients treated with anti-VEGF was genotyped for 138 tagging SNPs in the VEGF pathway. Per patient, we collected best corrected visual acuity at baseline, after three loading injections and at 12 months. We also registered the injection number and changes in retinal morphology after three loading injections (central foveal thickness (CFT), intraretinal cysts and serous neuroepithelium detachment). Changes in CFT after 3 months were our primary outcome measure. Association of SNPs to response was assessed by binomial logistic regression. Replication was attempted by associating visual acuity changes to genotypes in an independent Japanese cohort.

Results

Association with treatment response was detected for seven SNPs, including in FLT4 (rs55667289: OR=0.746, 95% CI 0.63 to 0.88, p=0.0005) and KDR (rs7691507: OR=1.056, 95% CI 1.02 to 1.10, p=0.005; and rs2305945: OR=0.963, 95% CI 0.93 to 1.00, p=0.0472). Only association with rs55667289 in FLT4 survived multiple testing correction. This SNP was unavailable for testing in the replication cohort. Of six SNPs tested for replication, one was significant although not after multiple testing correction.

Conclusion

Identifying genetic variants that define treatment response can help to develop individualised therapeutic approaches for wet ARMD patients and may point towards new targets in non-responders.

Quiescent Endothelial Cells Upregulate Fatty Acid β-Oxidation for Vasculoprotection via Redox Homeostasis

Little is known about the metabolism of quiescent endothelial cells (QECs). Nonetheless, when dysfunctional, QECs contribute to multiple diseases. Previously, we demonstrated that proliferating endothelial cells (PECs) use fatty acid β-oxidation (FAO) for de novo dNTP synthesis. We report now that QECs are not hypometabolic, but upregulate FAO >3-fold higher than PECs, not to support biomass or energy production but to sustain the tricarboxylic acid cycle for redox homeostasis through NADPH regeneration. Hence, endothelial loss of FAO-controlling CPT1A in CPT1A^ΔEC mice promotes EC dysfunction (leukocyte infiltration, barrier disruption) by increasing endothelial oxidative stress, rendering CPT1A^ΔEC mice more susceptible to LPS and inflammatory bowel disease. Mechanistically, Notch1 orchestrates the use of FAO for redox balance in QECs. Supplementation of acetate (metabolized to acetyl-coenzyme A) restores endothelial quiescence and counters oxidative stress-mediated EC dysfunction in CPT1A^ΔEC mice, offering therapeutic opportunities. Thus, QECs use FAO for vasculoprotection against oxidative stress-prone exposure.

Gene correlation network analysis to identify regulatory factors in idiopathic pulmonary fibrosis

BACKGROUND

Idiopathic pulmonary fibrosis (IPF) is a severe lung disease characterised by extensive pathological changes. The objective for this study was to identify the gene network and regulators underlying disease pathology in IPF and its association with lung function.

METHODS

Lung Tissue Research Consortium dataset with 262 IPF and control subjects (GSE47460) was randomly divided into two non-overlapping groups for cross-validated differential gene expression analysis. Consensus weighted gene coexpression network analysis identified overlapping coexpressed gene modules between both IPF groups. Modules were correlated with lung function (diffusion capacity, DL_CO; forced expiratory volume in 1 s, FEV₁; forced vital capacity, FVC) and enrichment analyses used to identify biological function and transcription factors. Module correlation with miRNA data (GSE72967) identified associated regulators. Clinical relevance in IPF was assessed in a peripheral blood gene expression dataset (GSE93606) to identify modules related to survival.

RESULTS

Correlation network analysis identified 16 modules in IPF. Upregulated modules were associated with cilia, DNA replication and repair, contractile fibres, B-cell and unfolded protein response, and extracellular matrix. Downregulated modules were associated with blood vessels, T-cell and interferon responses, leucocyte activation and degranulation, surfactant metabolism, and cellular metabolic and catabolic processes. Lung function correlated with nine modules (eight with DL_CO, five with FVC). Intermodular network of transcription factors and miRNA showed clustering of fibrosis, immune response and contractile modules. The cilia-associated module was able to predict survival (p=0.0097) in an independent peripheral blood IPF cohort.

CONCLUSIONS

We identified a correlation gene expression network with associated regulators in IPF that provides novel insight into the pathological process of this disease.

Phenotype molding of stromal cells in the lung tumor microenvironment

Cancer cells are embedded in the tumor microenvironment (TME), a complex ecosystem of stromal cells. Here, we present a 52,698-cell catalog of the TME transcriptome in human lung tumors at single-cell resolution, validated in independent samples where 40,250 additional cells were sequenced. By comparing with matching non-malignant lung samples, we reveal a highly complex TME that profoundly molds stromal cells. We identify 52 stromal cell subtypes, including novel subpopulations in cell types hitherto considered to be homogeneous, as well as transcription factors underlying their heterogeneity. For instance, we discover fibroblasts expressing different collagen sets, endothelial cells downregulating immune cell homing and genes coregulated with established immune checkpoint transcripts and correlating with T-cell activity. By assessing marker genes for these cell subtypes in bulk RNA-sequencing data from 1,572 patients, we illustrate how these correlate with survival, while immunohistochemistry for selected markers validates them as separate cellular entities in an independent series of lung tumors. Hence, in providing a comprehensive catalog of stromal cells types and by characterizing their phenotype and co-optive behavior, this resource provides deeper insights into lung cancer biology that will be helpful in advancing lung cancer diagnosis and therapy.

A Single-Cell Transcriptome Atlas of the Aging Drosophila Brain

The diversity of cell types and regulatory states in the brain, and how these change during aging, remains largely unknown. We present a single-cell transcriptome atlas of the entire adult Drosophila melanogaster brain sampled across its lifespan. Cell clustering identified 87 initial cell clusters that are further subclustered and validated by targeted cell-sorting. Our data show high granularity and identify a wide range of cell types. Gene network analyses using SCENIC revealed regulatory heterogeneity linked to energy consumption. During aging, RNA content declines exponentially without affecting neuronal identity in old brains. This single-cell brain atlas covers nearly all cells in the normal brain and provides the tools to study cellular diversity alongside other Drosophila and mammalian single-cell datasets in our unique single-cell analysis platform: SCope (http://scope.aertslab.org). These results, together with SCope, allow comprehensive exploration of all transcriptional states of an entire aging brain.

•

A single-cell atlas of the adult fly brain during aging

•

Network inference reveals regulatory states related to oxidative phosphorylation

•

Cell identity is retained during aging despite exponential decline of gene expression

•

SCope: An online tool to explore and compare single-cell datasets across species

Publisher Correction: Evidence that breast cancer risk at the 2q35 locus is mediated through IGFBP5 regulation

This corrects the article DOI: 10.1038/ncomms5999.

Ischemia-Induced DNA Hypermethylation during Kidney Transplant Predicts Chronic Allograft Injury

Background Ischemia during kidney transplant causes chronic allograft injury and adversely affects outcome, but the underlying mechanisms are incompletely understood. In tumors, oxygen shortage reduces the DNA demethylating activity of the ten-11 translocation (TET) enzymes, yielding hypermethylated genomes that promote tumor progression. We investigated whether ischemia similarly induces DNA hypermethylation in kidney transplants and contributes to chronic injury.Methods We profiled genome-wide DNA methylation in three cohorts of brain-dead donor kidney allograft biopsy specimens: a longitudinal cohort with paired biopsy specimens obtained at allograft procurement (preischemia; n=13), after implantation and reperfusion (postischemia; n=13), and at 3 or 12 months after transplant (n=5 each); a cross-sectional cohort with preimplantation biopsy specimens (n=82); and a cross-sectional cohort with postreperfusion biopsy specimens (n=46).Results Analysis of the paired preischemia and postischemia specimens revealed that methylation increased drastically in all allografts on ischemia. Hypermethylation was caused by loss of 5-hydroxymethylcytosine, the product of TET activity, and it was stable 1 year after transplant. In the preimplantation cohort, CpG hypermethylation directly correlated with ischemia time and for some CpGs, increased 2.6% per additional hour of ischemia. Hypermethylation preferentially affected and reduced the expression of genes involved in suppressing kidney injury and fibrosis. Moreover, CpG hypermethylation in preimplantation specimens predicted chronic injury, particularly fibrosis and glomerulosclerosis, 1 year after transplant. This finding was validated in the independent postreperfusion cohort, in which hypermethylation also predicted reduced allograft function 1 year after transplant, outperforming established clinical variables.Conclusions We highlight a novel epigenetic basis for ischemia-induced chronic allograft injury with biomarker potential.

Genomic and epigenomic analysis of high-risk prostate cancer reveals changes in hydroxymethylation and TET1

The clinical heterogeneity of prostate cancer (PCa) makes it difficult to identify those patients that could benefit from more aggressive treatments. As a contribution to a better understanding of the genomic changes in the primary tumor that are associated with the development of high-risk disease, we performed exome sequencing and copy number determination of a clinically homogeneous cohort of 47 high-risk PCas. We confirmed recurrent mutations in SPOP, PTEN and TP53 among the 850 point mutations we detected. In seven cases, we discovered genomic aberrations in the TET1 (Ten-Eleven Translocation 1) gene which encodes a DNA hydroxymethylase than can modify methylated cytosines in genomic DNA and thus is linked with gene expression changes. TET1 protein levels were reduced in tumor versus non-tumor prostate tissue in 39 of 40 cases. The clinical relevance of changes in TET1 levels was demonstrated in an independent PCa cohort, in which low TET1 mRNA levels were significantly associated with worse metastases-free survival. We also demonstrate a strong reduction in hydroxymethylated DNA in tumor tissue in 27 of 41 cases. Furthermore, we report the first exploratory (h)MeDIP-Seq analyses of eight high-risk PCa samples. This reveals a large heterogeneity in hydroxymethylation changes in tumor versus non-tumor genomes which can be linked with cell polarity.

It's T Time for Normal Blood Vessels

In a new report in Nature, Tian and colleagues (2017) describe that T cell activation in tumors induces blood vessel normalization. Because abnormal vasculature in tumors is known to reduce chemotherapy effectiveness and T lymphocyte infiltration and to increase metastasis, this suggests an exciting novel angle for checkpoint immunotherapy.

Comparative oncogenomics identifies tyrosine kinase FES as a tumor suppressor in melanoma

Identification and functional validation of oncogenic drivers are essential steps toward advancing cancer precision medicine. Here, we have presented a comprehensive analysis of the somatic genomic landscape of the widely used BRAFV600E- and NRASQ61K-driven mouse models of melanoma. By integrating the data with publically available genomic, epigenomic, and transcriptomic information from human clinical samples, we confirmed the importance of several genes and pathways previously implicated in human melanoma, including the tumor-suppressor genes phosphatase and tensin homolog (PTEN), cyclin dependent kinase inhibitor 2A (CDKN2A), LKB1, and others. Importantly, this approach also identified additional putative melanoma drivers with prognostic and therapeutic relevance. Surprisingly, one of these genes encodes the tyrosine kinase FES. Whereas FES is highly expressed in normal human melanocytes, FES expression is strongly decreased in over 30% of human melanomas. This downregulation correlates with poor overall survival. Correspondingly, engineered deletion of Fes accelerated tumor progression in a BRAFV600E-driven mouse model of melanoma. Together, these data implicate FES as a driver of melanoma progression and demonstrate the potential of cross-species oncogenomic approaches combined with mouse modeling to uncover impactful mutations and oncogenic driver alleles with clinical importance in the treatment of human cancer.

Tumors smother their epigenome

DNA hypermethylation is pervasive in tumors, but the factors triggering this modification are largely unknown. We recently demonstrated that the activity of 10-11-translocation methylcytosine dioxygenases, initiators of DNA demethylation, is compromised in hypoxic tumors. The resultant accumulation of methylation inactivates associated genes, linking the tumor microenvironment to epigenetic changes in cancer cells.

Inhibition of the Glycolytic Activator PFKFB3 in Endothelium Induces Tumor Vessel Normalization, Impairs Metastasis, and Improves Chemotherapy

Abnormal tumor vessels promote metastasis and impair chemotherapy. Hence, tumor vessel normalization (TVN) is emerging as an anti-cancer treatment. Here, we show that tumor endothelial cells (ECs) have a hyper-glycolytic metabolism, shunting intermediates to nucleotide synthesis. EC haplo-deficiency or blockade of the glycolytic activator PFKFB3 did not affect tumor growth, but reduced cancer cell invasion, intravasation, and metastasis by normalizing tumor vessels, which improved vessel maturation and perfusion. Mechanistically, PFKFB3 inhibition tightened the vascular barrier by reducing VE-cadherin endocytosis in ECs, and rendering pericytes more quiescent and adhesive (via upregulation of N-cadherin) through glycolysis reduction; it also lowered the expression of cancer cell adhesion molecules in ECs by decreasing NF-κB signaling. PFKFB3-blockade treatment also improved chemotherapy of primary and metastatic tumors.

The H3K9 dimethyltransferases EHMT1/2 protect against pathological cardiac hypertrophy

Cardiac hypertrophic growth in response to pathological cues is associated with reexpression of fetal genes and decreased cardiac function and is often a precursor to heart failure. In contrast, physiologically induced hypertrophy is adaptive, resulting in improved cardiac function. The processes that selectively induce these hypertrophic states are poorly understood. Here, we have profiled 2 repressive epigenetic marks, H3K9me2 and H3K27me3, which are involved in stable cellular differentiation, specifically in cardiomyocytes from physiologically and pathologically hypertrophied rat hearts, and correlated these marks with their associated transcriptomes. This analysis revealed the pervasive loss of euchromatic H3K9me2 as a conserved feature of pathological hypertrophy that was associated with reexpression of fetal genes. In hypertrophy, H3K9me2 was reduced following a miR-217-mediated decrease in expression of the H3K9 dimethyltransferases EHMT1 and EHMT2 (EHMT1/2). miR-217-mediated, genetic, or pharmacological inactivation of EHMT1/2 was sufficient to promote pathological hypertrophy and fetal gene reexpression, while suppression of this pathway protected against pathological hypertrophy both in vitro and in mice. Thus, we have established a conserved mechanism involving a departure of the cardiomyocyte epigenome from its adult cellular identity to a reprogrammed state that is accompanied by reexpression of fetal genes and pathological hypertrophy. These results suggest that targeting miR-217 and EHMT1/2 to prevent H3K9 methylation loss is a viable therapeutic approach for the treatment of heart disease.

Neurogenic Radial Glia-like Cells in Meninges Migrate and Differentiate into Functionally Integrated Neurons in the Neonatal Cortex

Whether new neurons are added in the postnatal cerebral cortex is still debated. Here, we report that the meninges of perinatal mice contain a population of neurogenic progenitors formed during embryonic development that migrate to the caudal cortex and differentiate into Satb2⁺ neurons in cortical layers II-IV. The resulting neurons are electrically functional and integrated into local microcircuits. Single-cell RNA sequencing identified meningeal cells with distinct transcriptome signatures characteristic of (1) neurogenic radial glia-like cells (resembling neural stem cells in the SVZ), (2) neuronal cells, and (3) a cell type with an intermediate phenotype, possibly representing radial glia-like meningeal cells differentiating to neuronal cells. Thus, we have identified a pool of embryonically derived radial glia-like cells present in the meninges that migrate and differentiate into functional neurons in the neonatal cerebral cortex.

Fine-scale mapping of 8q24 locus identifies multiple independent risk variants for breast cancer

Previous genome-wide association studies among women of European ancestry identified two independent breast cancer susceptibility loci represented by single nucleotide polymorphisms (SNPs) rs13281615 and rs11780156 at 8q24. A fine-mapping study across 2.06 Mb (chr8:127,561,724-129,624,067, hg19) in 55,540 breast cancer cases and 51,168 controls within the Breast Cancer Association Consortium was conducted. Three additional independent association signals in women of European ancestry, represented by rs35961416 (OR = 0.95, 95% CI = 0.93-0.97, conditional p = 5.8 × 10(-6) ), rs7815245 (OR = 0.94, 95% CI = 0.91-0.96, conditional p = 1.1 × 10(-6) ) and rs2033101 (OR = 1.05, 95% CI = 1.02-1.07, conditional p = 1.1 × 10(-4) ) were found. Integrative analysis using functional genomic data from the Roadmap Epigenomics, the Encyclopedia of DNA Elements project, the Cancer Genome Atlas and other public resources implied that SNPs rs7815245 in Signal 3, and rs1121948 in Signal 5 (in linkage disequilibrium with rs11780156, r(2)  = 0.77), were putatively functional variants for two of the five independent association signals. The results highlighted multiple 8q24 variants associated with breast cancer susceptibility in women of European ancestry.