Tumor monocyte content predicts immunochemotherapy outcomes in esophageal adenocarcinoma.

For inoperable esophageal adenocarcinoma (EAC), identifying patients likely to benefit from recently approved immunochemotherapy (ICI+CTX) treatments remains a key challenge. We address this using a uniquely designed window-of-opportunity trial (LUD2015-005), in which 35 inoperable EAC patients received first-line immune checkpoint inhibitors for four weeks (ICI-4W), followed by ICI+CTX. Comprehensive biomarker profiling, including generation of a 65,000-cell single-cell RNA-sequencing atlas of esophageal cancer, as well as multi-timepoint transcriptomic profiling of EAC during ICI-4W, reveals a novel T cell inflammation signature (INCITE) whose upregulation correlates with ICI-induced tumor shrinkage. Deconvolution of pre-treatment gastro-esophageal cancer transcriptomes using our single-cell atlas identifies high tumor monocyte content (TMC) as an unexpected ICI+CTX-specific predictor of greater overall survival (OS) in LUD2015-005 patients and of ICI response in prevalent gastric cancer subtypes from independent cohorts. Tumor mutational burden is an additional independent and additive predictor of LUD2015-005 OS. TMC can improve patient selection for emerging ICI+CTX therapies in gastro-esophageal cancer.

Cell-free DNA TAPS provides multimodal information for early cancer detection.

Multimodal, genome-wide characterization of epigenetic and genetic information in circulating cell-free DNA (cfDNA) could enable more sensitive early cancer detection, but it is technologically challenging. Recently, we developed TET-assisted pyridine borane sequencing (TAPS), which is a mild, bisulfite-free method for base-resolution direct DNA methylation sequencing. Here, we optimized TAPS for cfDNA (cfTAPS) to provide high-quality and high-depth whole-genome cell-free methylomes. We applied cfTAPS to 85 cfDNA samples from patients with hepatocellular carcinoma (HCC) or pancreatic ductal adenocarcinoma (PDAC) and noncancer controls. From only 10 ng of cfDNA (1 to 3 ml of plasma), we generated the most comprehensive cfDNA methylome to date. We demonstrated that cfTAPS provides multimodal information about cfDNA characteristics, including DNA methylation, tissue of origin, and DNA fragmentation. Integrated analysis of these epigenetic and genetic features enables accurate identification of early HCC and PDAC.

Tuning Transcription Factor Availability through Acetylation-Mediated Genomic Redistribution.

It is widely assumed that decreasing transcription factor DNA-binding affinity reduces transcription initiation by diminishing occupancy of sequence-specific regulatory elements. However, in vivo transcription factors find their binding sites while confronted with a large excess of low-affinity degenerate motifs. Here, using the melanoma lineage survival oncogene MITF as a model, we show that low-affinity binding sites act as a competitive reservoir in vivo from which transcription factors are released by mitogen-activated protein kinase (MAPK)-stimulated acetylation to promote increased occupancy of their regulatory elements. Consequently, a low-DNA-binding-affinity acetylation-mimetic MITF mutation supports melanocyte development and drives tumorigenesis, whereas a high-affinity non-acetylatable mutant does not. The results reveal a paradoxical acetylation-mediated molecular clutch that tunes transcription factor availability via genome-wide redistribution and couples BRAF to tumorigenesis. Our results further suggest that p300/CREB-binding protein-mediated transcription factor acetylation may represent a common mechanism to control transcription factor availability.

Author Correction: The effects of mutational processes and selection on driver mutations across cancer types.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Cytosine Methylation Affects the Mutability of Neighboring Nucleotides in Germline and Soma.

Methylated cytosines deaminate at higher rates than unmethylated cytosines, and the lesions they produce are repaired less efficiently. As a result, methylated cytosines are mutational hotspots. Here, combining rare polymorphism and base-resolution methylation data in humans, Arabidopsis thaliana, and rice (Oryza sativa), we present evidence that methylation state affects mutation dynamics not only at the focal cytosine but also at neighboring nucleotides. In humans, contrary to prior suggestions, we find that nucleotides in the close vicinity (±3 bp) of methylated cytosines mutate less frequently. Reduced mutability around methylated CpGs is also observed in cancer genomes, considering single nucleotide variants alongside tissue-of-origin-matched methylation data. In contrast, methylation is associated with increased neighborhood mutation risk in A. thaliana and rice. The difference in neighborhood mutation risk is less pronounced further away from the focal CpG and modulated by regional GC content. Our results are consistent with a model where altered risk at neighboring bases is linked to lesion formation at the focal CpG and subsequent long-patch repair. Our findings indicate that cytosine methylation has a broader mutational footprint than is commonly assumed.

Bisulfite-free direct detection of 5-methylcytosine and 5-hydroxymethylcytosine at base resolution.

Bisulfite sequencing has been the gold standard for mapping DNA modifications including 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC) for decades1-4. However, this harsh chemical treatment degrades the majority of the DNA and generates sequencing libraries with low complexity2,5,6. Here, we present a bisulfite-free and base-level-resolution sequencing method, TET-assisted pyridine borane sequencing (TAPS), for detection of 5mC and 5hmC. TAPS combines ten-eleven translocation (TET) oxidation of 5mC and 5hmC to 5-carboxylcytosine (5caC) with pyridine borane reduction of 5caC to dihydrouracil (DHU). Subsequent PCR converts DHU to thymine, enabling a C-to-T transition of 5mC and 5hmC. TAPS detects modifications directly with high sensitivity and specificity, without affecting unmodified cytosines. This method is nondestructive, preserving DNA fragments over 10 kilobases long. We applied TAPS to the whole-genome mapping of 5mC and 5hmC in mouse embryonic stem cells and show that, compared with bisulfite sequencing, TAPS results in higher mapping rates, more even coverage and lower sequencing costs, thus enabling higher quality, more comprehensive and cheaper methylome analyses.

Single cell RNA-seq reveals profound transcriptional similarity between Barrett's oesophagus and oesophageal submucosal glands.

Barrett's oesophagus is a precursor of oesophageal adenocarcinoma. In this common condition, squamous epithelium in the oesophagus is replaced by columnar epithelium in response to acid reflux. Barrett's oesophagus is highly heterogeneous and its relationships to normal tissues are unclear. Here we investigate the cellular complexity of Barrett's oesophagus and the upper gastrointestinal tract using RNA-sequencing of single cells from multiple biopsies from six patients with Barrett's oesophagus and two patients without oesophageal pathology. We find that cell populations in Barrett's oesophagus, marked by LEFTY1 and OLFM4, exhibit a profound transcriptional overlap with oesophageal submucosal gland cells, but not with gastric or duodenal cells. Additionally, SPINK4 and ITLN1 mark cells that precede morphologically identifiable goblet cells in colon and Barrett's oesophagus, potentially aiding the identification of metaplasia. Our findings reveal striking transcriptional relationships between normal tissue populations and cells in a premalignant condition, with implications for clinical practice.

Mutational signature distribution varies with DNA replication timing and strand asymmetry.

BACKGROUND: DNA replication plays an important role in mutagenesis, yet little is known about how it interacts with other mutagenic processes. Here, we use somatic mutation signatures-each representing a mutagenic process-derived from 3056 patients spanning 19 cancer types to quantify the strand asymmetry of mutational signatures around replication origins and between early and late replicating regions. RESULTS: We observe that most of the detected mutational signatures are significantly correlated with the timing or direction of DNA replication. The properties of these associations are distinct for different signatures and shed new light on several mutagenic processes. For example, our results suggest that oxidative damage to the nucleotide pool substantially contributes to the mutational landscape of esophageal adenocarcinoma. CONCLUSIONS: Together, our results indicate an interaction between DNA replication, the associated damage repair, and most mutagenic processes.

DNA Modifications: Naturally More Error Prone?

Epigenetic DNA modifications are essential for normal cell function in vertebrates, but they can also be hotspots of mutagenesis. Methylcytosine in particular has long been known to be less stable than other nucleotides and spontaneously deaminates to thymine. Beyond this well-established phenomenon, however, the influence of epigenetic marks on mutagenesis has recently become an active field of investigation. In this review, we summarize current knowledge of the interactions between different DNA modifications and other mutagenic processes. External mutagens, such as UV light or smoking carcinogens, affect modified cytosines differently from unmodified ones, and modified cytosine can in some cases be protective rather than mutagenic. Notably, cell-intrinsic processes, such as DNA replication, also appear to influence the mutagenesis of modified cytosines. Altogether, evidence is accumulating to show that epigenetic changes have a profound influence on tissue-specific mutation accumulation.

The effects of mutational processes and selection on driver mutations across cancer types.

Epidemiological evidence has long associated environmental mutagens with increased cancer risk. However, links between specific mutation-causing processes and the acquisition of individual driver mutations have remained obscure. Here we have used public cancer sequencing data from 11,336 cancers of various types to infer the independent effects of mutation and selection on the set of driver mutations in a cancer type. First, we detect associations between a range of mutational processes, including those linked to smoking, ageing, APOBEC and DNA mismatch repair (MMR) and the presence of key driver mutations across cancer types. Second, we quantify differential selection between well-known alternative driver mutations, including differences in selection between distinct mutant residues in the same gene. These results show that while mutational processes have a large role in determining which driver mutations are present in a cancer, the role of selection frequently dominates.

DNA Replication and associated repair pathways are involved in the mutagenesis of methylated cytosine.

Transitions of cytosine to thymine in CpG dinucleotides are the most frequent type of mutations observed in cancer. This increased mutability is commonly explained by the presence of 5-methylcytosine (5mC) and its spontaneous hydrolytic deamination into thymine. Here, we describe observations that question whether spontaneous deamination alone causes the elevated mutagenicity of 5mC. Tumours with somatic mutations in DNA mismatch-repair genes or in the proofreading domain of DNA polymerase ε (Pol ε) exhibit more 5mC to T transitions than would be expected, given the kinetics of hydrolytic deamination. This enrichment is asymmetrical around replication origins with a preference for the leading strand template, in particular in methylated cytosines flanked by guanines (GCG). Notably, GCG to GTG mutations also exhibit strand asymmetry in mismatch-repair and Pol ε wild-type tumours. Together, these findings suggest that mis-incorporation of A opposite 5mC during replication of the leading strand might be a contributing factor in the mutagenesis of methylated cytosine.

The importance of p53 pathway genetics in inherited and somatic cancer genomes.

Decades of research have shown that mutations in the p53 stress response pathway affect the incidence of diverse cancers more than mutations in other pathways. However, most evidence is limited to somatic mutations and rare inherited mutations. Using newly abundant genomic data, we demonstrate that commonly inherited genetic variants in the p53 pathway also affect the incidence of a broad range of cancers more than variants in other pathways. The cancer-associated single nucleotide polymorphisms (SNPs) of the p53 pathway have strikingly similar genetic characteristics to well-studied p53 pathway cancer-causing somatic mutations. Our results enable insights into p53-mediated tumour suppression in humans and into p53 pathway-based cancer surveillance and treatment strategies.

Chromatin organization is a major influence on regional mutation rates in human cancer cells.

Cancer genome sequencing provides the first direct information on how mutation rates vary across the human genome in somatic cells. Testing diverse genetic and epigenetic features, here we show that mutation rates in cancer genomes are strikingly related to chromatin organization. Indeed, at the megabase scale, a single feature­levels of the heterochromatin-associated histone modification H3K9me3­can account for more than 40% of mutation-rate variation, and a combination of features can account for more than 55%. The strong association between mutation rates and chromatin organization is upheld in samples from different tissues and for different mutation types. This suggests that the arrangement of the genome into heterochromatin- and euchromatin-like domains is a dominant influence on regional mutation-rate variation in human somatic cells.

Dosage sensitivity shapes the evolution of copy-number varied regions.

Dosage sensitivity is an important evolutionary force which impacts on gene dispensability and duplicability. The newly available data on human copy-number variation (CNV) allow an analysis of the most recent and ongoing evolution. Provided that heterozygous gene deletions and duplications actually change gene dosage, we expect to observe negative selection against CNVs encompassing dosage sensitive genes. In this study, we make use of several sources of population genetic data to identify selection on structural variations of dosage sensitive genes. We show that CNVs can directly affect expression levels of contained genes. We find that genes encoding members of protein complexes exhibit limited expression variation and overlap significantly with a manually derived set of dosage sensitive genes. We show that complexes and other dosage sensitive genes are underrepresented in CNV regions, with a particular bias against frequent variations and duplications. These results suggest that dosage sensitivity is a significant force of negative selection on regions of copy-number variation.

Protein interactions in human genetic diseases.

We present a novel method that combines protein structure information with protein interaction data to identify residues that form part of an interaction interface. Our prediction method can retrieve interaction hotspots with an accuracy of 60% (at a 20% false positive rate). The method was applied to all mutations in the Online Mendelian Inheritance in Man (OMIM) database, predicting 1,428 mutations to be related to an interaction defect. Combining predicted and hand-curated sets, we discuss how mutations affect protein interactions in general.

Reuse of structural domain-domain interactions in protein networks.

BACKGROUND: Protein interactions are thought to be largely mediated by interactions between structural domains. Databases such as iPfam relate interactions in protein structures to known domain families. Here, we investigate how the domain interactions from the iPfam database are distributed in protein interactions taken from the HPRD, MPact, BioGRID, DIP and IntAct databases. RESULTS: We find that known structural domain interactions can only explain a subset of 4-19% of the available protein interactions, nevertheless this fraction is still significantly bigger than expected by chance. There is a correlation between the frequency of a domain interaction and the connectivity of the proteins it occurs in. Furthermore, a large proportion of protein interactions can be attributed to a small number of domain interactions. We conclude that many, but not all, domain interactions constitute reusable modules of molecular recognition. A substantial proportion of domain interactions are conserved between E. coli, S. cerevisiae and H. sapiens. These domains are related to essential cellular functions, suggesting that many domain interactions were already present in the last universal common ancestor. CONCLUSION: Our results support the concept of domain interactions as reusable, conserved building blocks of protein interactions, but also highlight the limitations currently imposed by the small number of available protein structures.

An introduction to hidden Markov models.

This unit introduces the concept of hidden Markov models in computational biology. It describes them using simple biological examples, requiring as little mathematical knowledge as possible. The unit also presents a brief history of hidden Markov models and an overview of their current applications before concluding with a discussion of their limitations.

Pfam: clans, web tools and services.

Pfam is a database of protein families that currently contains 7973 entries (release 18.0). A recent development in Pfam has enabled the grouping of related families into clans. Pfam clans are described in detail, together with the new associated web pages. Improvements to the range of Pfam web tools and the first set of Pfam web services that allow programmatic access to the database and associated tools are also presented. Pfam is available on the web in the UK (http://www.sanger.ac.uk/Software/Pfam/), the USA (http://pfam.wustl.edu/), France (http://pfam.jouy.inra.fr/) and Sweden (http://pfam.cgb.ki.se/).

Visualizing profile-profile alignment: pairwise HMM logos.

UNLABELLED: The availability of advanced profile-profile comparison tools, such as PRC or HHsearch demands sophisticated visualization tools not presently available. We introduce an approach built upon the concept of HMM logos. The method illustrates the similarities of pairs of protein family profiles in an intuitive way. Two HMM logos, one for each profile, are drawn one upon the other. The aligned states are then highlighted and connected. AVAILABILITY: A web interface offering online creation of pairwise HMM logos is available at http://www.sanger.ac.uk/Software/analysis/logomat-p. Furthermore, software developers may download a Perl package that includes methods for creation of pairwise HMM logos locally. CONTACT: bsb@sanger.ac.uk.