DNA replication errors that persist as mismatch mutations make up the molecular fingerprint of mismatch repair (MMR)-deficient tumors and convey them with resistance to standard therapy. Using whole-genome and whole-exome sequencing, we here confirm an MMR-deficient mutation signature that is distinct from other tumor genomes, but surprisingly similar to germ-line DNA, indicating that a substantial fraction of human genetic variation arises through mutations escaping MMR. Moreover, we identify a large set of recurrent indels that may serve to detect microsatellite instability (MSI). Indeed, using endometrial tumors with immunohistochemically proven MMR deficiency, we optimize a novel marker set capable of detecting MSI and show it to have greater specificity and selectivity than standard MSI tests. Additionally, we show that recurrent indels are enriched for the 'DNA double-strand break repair by homologous recombination' pathway. Consequently, DSB repair is reduced in MMR-deficient tumors, triggering a dose-dependent sensitivity of MMR-deficient tumor cultures to DSB inducers.

Hypermethylation of the promoters of tumour suppressor genes represses transcription of these genes, conferring growth advantages to cancer cells. How these changes arise is poorly understood. Here we show that the activity of oxygen-dependent ten-eleven translocation (TET) enzymes is reduced by tumour hypoxia in human and mouse cells. TET enzymes catalyse DNA demethylation through 5-methylcytosine oxidation. This reduction in activity occurs independently of hypoxia-associated alterations in TET expression, proliferation, metabolism, hypoxia-inducible factor activity or reactive oxygen species, and depends directly on oxygen shortage. Hypoxia-induced loss of TET activity increases hypermethylation at gene promoters in vitro. In patients, tumour suppressor gene promoters are markedly more methylated in hypoxic tumour tissue, independent of proliferation, stromal cell infiltration and tumour characteristics. Our data suggest that up to half of hypermethylation events are due to hypoxia, with these events conferring a selective advantage. Accordingly, increased hypoxia in mouse breast tumours increases hypermethylation, while restoration of tumour oxygenation abrogates this effect. Tumour hypoxia therefore acts as a novel regulator of DNA methylation.

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.

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.

The 10q26 locus in the second intron of FGFR2 is the locus most strongly associated with estrogen-receptor-positive breast cancer in genome-wide association studies. We conducted fine-scale mapping in case-control studies genotyped with a custom chip (iCOGS), comprising 41 studies (n = 89,050) of European ancestry, 9 Asian ancestry studies (n = 13,983), and 2 African ancestry studies (n = 2,028) from the Breast Cancer Association Consortium. We identified three statistically independent risk signals within the locus. Within risk signals 1 and 3, genetic analysis identified five and two variants, respectively, highly correlated with the most strongly associated SNPs. By using a combination of genetic fine mapping, data on DNase hypersensitivity, and electrophoretic mobility shift assays to study protein-DNA binding, we identified rs35054928, rs2981578, and rs45631563 as putative functional SNPs. Chromatin immunoprecipitation showed that FOXA1 preferentially bound to the risk-associated allele (C) of rs2981578 and was able to recruit ERα to this site in an allele-specific manner, whereas E2F1 preferentially bound the risk variant of rs35054928. The risk alleles were preferentially found in open chromatin and bound by Ser5 phosphorylated RNA polymerase II, suggesting that the risk alleles are associated with changes in transcription. Chromatin conformation capture demonstrated that the risk region was able to interact with the promoter of FGFR2, the likely target gene of this risk region. A role for FOXA1 in mediating breast cancer susceptibility at this locus is consistent with the finding that the FGFR2 risk locus primarily predisposes to estrogen-receptor-positive disease.

Congenital heart defects (CHDs) are the most common birth defect worldwide and are a leading cause of neonatal mortality. Nonsyndromic atrioventricular septal defects (AVSDs) are an important subtype of CHDs for which the genetic architecture is poorly understood. We performed exome sequencing in 13 parent-offspring trios and 112 unrelated individuals with nonsyndromic AVSDs and identified five rare missense variants (two of which arose de novo) in the highly conserved gene NR2F2, a very significant enrichment (p = 7.7 × 10(-7)) compared to 5,194 control subjects. We identified three additional CHD-affected families with other variants in NR2F2 including a de novo balanced chromosomal translocation, a de novo substitution disrupting a splice donor site, and a 3 bp duplication that cosegregated in a multiplex family. NR2F2 encodes a pleiotropic developmental transcription factor, and decreased dosage of NR2F2 in mice has been shown to result in abnormal development of atrioventricular septa. Via luciferase assays, we showed that all six coding sequence variants observed in individuals significantly alter the activity of NR2F2 on target promoters.

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.

Genetic studies (in particular linkage and association studies) identify chromosomal regions involved in a disease or phenotype of interest, but those regions often contain many candidate genes, only a few of which can be followed-up for biological validation. Recently, computational methods to identify (prioritize) the most promising candidates within a region have been proposed, but they are usually not applicable to cases where little is known about the phenotype (no or few confirmed disease genes, fragmentary understanding of the biological cascades involved). We seek to overcome this limitation by replacing knowledge about the biological process by experimental data on differential gene expression between affected and healthy individuals. Considering the problem from the perspective of a gene/protein network, we assess a candidate gene by considering the level of differential expression in its neighborhood under the assumption that strong candidates will tend to be surrounded by differentially expressed neighbors. We define a notion of soft neighborhood where each gene is given a contributing weight, which decreases with the distance from the candidate gene on the protein network. To account for multiple paths between genes, we define the distance using the Laplacian exponential diffusion kernel. We score candidates by aggregating the differential expression of neighbors weighted as a function of distance. Through a randomization procedure, we rank candidates by p-values. We illustrate our approach on four monogenic diseases and successfully prioritize the known disease causing genes.

Genomic imbalances are a major cause of constitutional and acquired disorders. Therefore, aneuploidy screening has become the cornerstone of preimplantation, prenatal and postnatal genetic diagnosis, as well as a routine aspect of the diagnostic workup of many acquired disorders. Recently, array comparative genomic hybridization (array CGH) has been introduced as a rapid and high-resolution method for the detection of both benign and disease-causing genomic copy-number variations. Until now, array CGH has been performed using a significant quantity of DNA derived from a pool of cells. Here, we present an array CGH method that accurately detects chromosomal imbalances from a single lymphoblast, fibroblast and blastomere within a single day. Trisomy 13, 18, 21 and monosomy X, as well as normal ploidy levels of all other chromosomes, were accurately determined from single fibroblasts. Moreover, we showed that a segmental deletion as small as 34 Mb could be detected. Finally, we demonstrated the possibility to detect aneuploidies in single blastomeres derived from preimplantation embryos. This technique offers new possibilities for genetic analysis of single cells in general and opens the route towards aneuploidy screening and detection of unbalanced translocations in preimplantation embryos in particular.

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.

Congenital heart disease (CHD) occurs in almost 1% of newborn children and is considered a multifactorial disorder. CHD may segregate in families due to significant contribution of genetic factors in the disease etiology. The aim of the study was to identify pathophysiological mechanisms in families segregating CHD.

Congenital heart defects (CHDs) have a neonatal incidence of 0.8-1% (refs. 1,2). Despite abundant examples of monogenic CHD in humans and mice, CHD has a low absolute sibling recurrence risk (∼2.7%), suggesting a considerable role for de novo mutations (DNMs) and/or incomplete penetrance. De novo protein-truncating variants (PTVs) have been shown to be enriched among the 10% of 'syndromic' patients with extra-cardiac manifestations. We exome sequenced 1,891 probands, including both syndromic CHD (S-CHD, n = 610) and nonsyndromic CHD (NS-CHD, n = 1,281). In S-CHD, we confirmed a significant enrichment of de novo PTVs but not inherited PTVs in known CHD-associated genes, consistent with recent findings. Conversely, in NS-CHD we observed significant enrichment of PTVs inherited from unaffected parents in CHD-associated genes. We identified three genome-wide significant S-CHD disorders caused by DNMs in CHD4, CDK13 and PRKD1. Our study finds evidence for distinct genetic architectures underlying the low sibling recurrence risk in S-CHD and NS-CHD.

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.

During ageing, kidney function decreases due to renal tubular atrophy, interstitial fibrosis, glomerulosclerosis and arteriosclerosis. Recently, changes in DNA methylation were shown to contribute to various ageing processes. However, it is unknown whether such changes also contribute to age-related kidney dysfunction. To assess this, we profiled genome-wide changes in DNA methylation (over 800 000 CpG sites) in 95 renal biopsies obtained prior to kidney transplantation from donors aged 16 to 73 years. Donor age significantly associated with the methylation of 92 778 CpGs (false discovery rate under 0.05), corresponding to 10 285 differentially methylated regions. These regions were most frequently located in genes involved in the Wnt/beta-catenin signaling pathway. Using an independent cohort of 67 biopsies, we autonomously validated these findings. Interestingly, the methylation status of these 92 778 age-related CpGs was associated with glomerulosclerosis (34.4% of CpGs at a false discovery rate under 0.05) and interstitial fibrosis (0.9%) and graft function at one year after transplantation, but not with tubular atrophy and arteriosclerosis. No association was observed with any of these pathologies at the time of transplantation (0% at a false discovery rate under 0.05). Thus, age-associated changes in DNA methylation at the time of transplantation predict future injury of transplanted kidneys. Specifically, our epigenome-wide association study demonstrates that epigenetic renal ageing is implicated in progressive fibrosis in both the glomerulus and the interstitium.

Aflibercept (ziv-aflibercept) is an anti-angiogenic agent recently approved in combination with FOLFIRI for the treatment of metastatic colorectal cancer (mCRC) patients previously treated with oxaliplatin. Despite heterogeneity in response to aflibercept, no biomarkers for efficacy or adverse effects have been identified. Here we present biomarker data from the randomised phase II AFFIRM trial assessing aflibercept in combination with mFOLFOX6 first line in mCRC.

Genome-wide DNA methylation reprogramming occurs in mouse primordial germ cells (PGCs) and preimplantation embryos, but the precise dynamics and biological outcomes are largely unknown. We have carried out whole-genome bisulfite sequencing (BS-Seq) and RNA-Seq across key stages from E6.5 epiblast to E16.5 PGCs. Global loss of methylation takes place during PGC expansion and migration with evidence for passive demethylation, but sequences that carry long-term epigenetic memory (imprints, CpG islands on the X chromosome, germline-specific genes) only become demethylated upon entry of PGCs into the gonads. The transcriptional profile of PGCs is tightly controlled despite global hypomethylation, with transient expression of the pluripotency network, suggesting that reprogramming and pluripotency are inextricably linked. Our results provide a framework for the understanding of the epigenetic ground state of pluripotency in the germline.

Preservation of bone mass is crucial for healthy ageing and largely depends on adequate responses of matrix-embedded osteocytes. These cells control bone formation and resorption concurrently by secreting the WNT/β-catenin antagonist sclerostin (SOST). Osteocytes reside within a low oxygen microenvironment, but whether and how oxygen sensing regulates their function remains elusive. Here, we show that conditional deletion of the oxygen sensor prolyl hydroxylase (PHD) 2 in osteocytes results in a high bone mass phenotype, which is caused by increased bone formation and decreased resorption. Mechanistically, enhanced HIF-1α signalling increases Sirtuin 1-dependent deacetylation of the Sost promoter, resulting in decreased sclerostin expression and enhanced WNT/β-catenin signalling. Additionally, genetic ablation of PHD2 in osteocytes blunts osteoporotic bone loss induced by oestrogen deficiency or mechanical unloading. Thus, oxygen sensing by PHD2 in osteocytes negatively regulates bone mass through epigenetic regulation of sclerostin and targeting PHD2 elicits an osteo-anabolic response in osteoporotic models.

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.

Congenital heart defects (CHDs) are the most common major developmental anomalies and the most frequent cause for perinatal mortality, but their etiology remains often obscure. We identified a locus for CHDs on 6q24-q25. Genotype-phenotype correlations in 12 patients carrying a chromosomal deletion on 6q delineated a critical 850 kb region on 6q25.1 harboring five genes. Bioinformatics prioritization of candidate genes in this locus for a role in CHDs identified the TGF-beta-activated kinase 1/MAP3K7 binding protein 2 gene (TAB2) as the top-ranking candidate gene. A role for this candidate gene in cardiac development was further supported by its conserved expression in the developing human and zebrafish heart. Moreover, a critical, dosage-sensitive role during development was demonstrated by the cardiac defects observed upon titrated knockdown of tab2 expression in zebrafish embryos. To definitively confirm the role of this candidate gene in CHDs, we performed mutation analysis of TAB2 in 402 patients with a CHD, which revealed two evolutionarily conserved missense mutations. Finally, a balanced translocation was identified, cosegregating with familial CHD. Mapping of the breakpoints demonstrated that this translocation disrupts TAB2. Taken together, these data clearly demonstrate a role for TAB2 in human cardiac development.

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.