Corrigendum: Cholangiocytes act as facultative liver stem cells during impaired hepatocyte regeneration.

This corrects the article DOI: 10.1038/nature23015.

The Scottish COVID Cancer Immunity Prevalence Study: A Longitudinal Study of SARS-CoV-2 Immune Response in Patients Receiving Anti-Cancer Treatment.

BACKGROUND: Cancer and anti-cancer treatment (ACT) may be risk factors for severe SARS-CoV-2 infection and limited vaccine efficacy. Long-term longitudinal studies are needed to evaluate these risks. The Scottish COVID cancer immunity prevalence (SCCAMP) study characterizes the incidence and outcomes of SARS-CoV-2 infection and vaccination in patients with solid tumors undergoing ACT. This preliminary analysis includes 766 patients recruited since May 2020. METHODS: Patients with solid-organ cancers attending secondary care for active ACT consented to the collection of routine electronic health record data and serial blood samples over 12 months. Blood samples were tested for total SARS-CoV-2 antibody. RESULTS: A total of 766 participants were recruited between May 28, 2020 and October 31, 2021. Most received cytotoxic chemotherapy (79%). Among the participants, 48 (6.3%) were tested positive for SARS-CoV-2 by PCR. Infection rates were unaffected by ACT, largely aligning with the local population. Mortality proportion was not higher with a recent positive SARS-CoV-2 PCR (10.4% vs 10.6%). Multivariate analysis revealed lower infection rates in vaccinated patients regardless of chemotherapy (HR 0.307 [95% CI, 0.144-0.6548]) or immunotherapy (HR 0.314 [95% CI, 0.041-2.367]) treatment. A total of 96.3% of patients successfully raised SARS-CoV-2 antibodies after >2 vaccines. This was independent of the treatment type. CONCLUSION: This is the largest on-going longitudinal real-world dataset of patients undergoing ACT during the early stages of the COVID-19 pandemic. This preliminary analysis demonstrates that patients with solid tumors undergoing ACT have high protection from SARS-CoV-2 infection following COVID-19 vaccination. The SCCAMP study will evaluate long-term COVID-19 antibody trends, focusing on specific ACTs and patient subgroups.

Whole exome sequencing of low grade serous ovarian carcinoma identifies genomic events associated with clinical outcome.

OBJECTIVES: Low-grade serous ovarian carcinoma (LGSOC) is a distinct, rare, ovarian cancer type characterised by younger patient age and intrinsic chemoresistance. Understanding the molecular landscape is crucial for optimising targeted therapy. METHODS: Genomic data from whole exome sequencing of tumour tissue was analysed in a LGSOC cohort with detailed clinical annotation. RESULTS: 63 cases were analysed and three subgroups identified based on single nucleotide variants: canonical MAPK mutant (cMAPKm: 52%, KRAS/BRAF/NRAS), MAPK-associated gene mutation (MAPK-assoc: 27%) and MAPK wild-type (MAPKwt: 21%). NOTCH pathway disruption occurred across all subgroups. Tumour mutational burden (TMB), mutational signatures and recurrent copy number (CN) changes varied across the cohort with co-occurrence of chromosome 1p loss and 1q gain (CN Chr1pq) a recurrent feature. Low TMB and CN Chr1pq were associated with inferior disease-specific survival (HR 6.43; p < 0.001 and HR 3.29, p = 0.011 respectively). Stepwise genomic classification in relation to outcome resulted in four groups (TMB low; CN Chr1pq; MAPKwt/MAPKassoc; cMAPKm). 5 year disease-specific survival was 46%, 55%, 79% and 100% respectively for these groups. The two most favourable genomic subgroups were enriched for the SBS10b mutational signature, particularly the cMAPKm subgroup. CONCLUSIONS: LGSOC comprises multiple genomic subgroups with distinct clinical and molecular features. Chr1pq CN arm disruption and TMB represent promising methods to identify individuals with poorer prognosis. Further investigation of the molecular basis for these observations is required. MAPKwt cases represent around a fifth of patients. NOTCH inhibitors represent a candidate therapeutic strategy worthy of exploration across these cases.

Cholangiocytes act as facultative liver stem cells during impaired hepatocyte regeneration.

After liver injury, regeneration occurs through self-replication of hepatocytes. In severe liver injury, hepatocyte proliferation is impaired-a feature of human chronic liver disease. It is unclear whether other liver cell types can regenerate hepatocytes. Here we use two independent systems to impair hepatocyte proliferation during liver injury to evaluate the contribution of non-hepatocytes to parenchymal regeneration. First, loss of ß1-integrin in hepatocytes with liver injury triggered a ductular reaction of cholangiocyte origin, with approximately 25% of hepatocytes being derived from a non-hepatocyte origin. Second, cholangiocytes were lineage traced with concurrent inhibition of hepatocyte proliferation by ß1-integrin knockdown or p21 overexpression, resulting in the significant emergence of cholangiocyte-derived hepatocytes. We describe a model of combined liver injury and inhibition of hepatocyte proliferation that causes physiologically significant levels of regeneration of functional hepatocytes from biliary cells.

Activation of transcription factor circuity in 2i-induced ground state pluripotency is independent of repressive global epigenetic landscapes.

Mouse embryonic stem cells (mESCs) cultured with MEK/ERK and GSK3ß (2i) inhibitors transition to ground state pluripotency. Gene expression changes, redistribution of histone H3K27me3 profiles and global DNA hypomethylation are hallmarks of 2i exposure, but it is unclear whether epigenetic alterations are required to achieve and maintain ground state or occur as an outcome of 2i signal induced changes. Here we show that ESCs with three epitypes, WT, constitutively methylated, or hypomethylated, all undergo comparable morphological, protein expression and transcriptome changes independently of global alterations of DNA methylation levels or changes in H3K27me3 profiles. Dazl and Fkbp6 expression are induced by 2i in all three epitypes, despite exhibiting hypermethylated promoters in constitutively methylated ESCs. We identify a number of activated gene promoters that undergo 2i dependent loss of H3K27me3 in all three epitypes, however genetic and pharmaceutical inhibition experiments show that H3K27me3 is not required for their silencing in non-2i conditions. By separating and defining their contributions, our data suggest that repressive epigenetic systems play minor roles in mESC self-renewal and naïve ground state establishment by core sets of dominant pluripotency associated transcription factor networks, which operate independently from these epigenetic processes.

CpG islands influence chromatin structure via the CpG-binding protein Cfp1.

CpG islands (CGIs) are prominent in the mammalian genome owing to their GC-rich base composition and high density of CpG dinucleotides. Most human gene promoters are embedded within CGIs that lack DNA methylation and coincide with sites of histone H3 lysine 4 trimethylation (H3K4me3), irrespective of transcriptional activity. In spite of these intriguing correlations, the functional significance of non-methylated CGI sequences with respect to chromatin structure and transcription is unknown. By performing a search for proteins that are common to all CGIs, here we show high enrichment for Cfp1, which selectively binds to non-methylated CpGs in vitro. Chromatin immunoprecipitation of a mono-allelically methylated CGI confirmed that Cfp1 specifically associates with non-methylated CpG sites in vivo. High throughput sequencing of Cfp1-bound chromatin identified a notable concordance with non-methylated CGIs and sites of H3K4me3 in the mouse brain. Levels of H3K4me3 at CGIs were markedly reduced in Cfp1-depleted cells, consistent with the finding that Cfp1 associates with the H3K4 methyltransferase Setd1 (refs 7, 8). To test whether non-methylated CpG-dense sequences are sufficient to establish domains of H3K4me3, we analysed artificial CpG clusters that were integrated into the mouse genome. Despite the absence of promoters, the insertions recruited Cfp1 and created new peaks of H3K4me3. The data indicate that a primary function of non-methylated CGIs is to genetically influence the local chromatin modification state by interaction with Cfp1 and perhaps other CpG-binding proteins.

Dynamic changes in 5-hydroxymethylation signatures underpin early and late events in drug exposed liver.

Aberrant DNA methylation is a common feature of neoplastic lesions, and early detection of such changes may provide powerful mechanistic insights and biomarkers for carcinogenesis. Here, we investigate dynamic changes in the mouse liver DNA methylome associated with short (1 day) and prolonged (7, 28 and 91 days) exposure to the rodent liver non-genotoxic carcinogen, phenobarbital (PB). We find that the distribution of 5mC/5hmC is highly consistent between untreated individuals of a similar age; yet, changes during liver maturation in a transcriptionally dependent manner. Following drug treatment, we identify and validate a series of differentially methylated or hydroxymethylated regions: exposure results in staged transcriptional responses with distinct kinetic profiles that strongly correlate with promoter proximal region 5hmC levels. Furthermore, reciprocal changes for both 5mC and 5hmC in response to PB suggest that active demethylation may be taking place at each set of these loci via a 5hmC intermediate. Finally, we identify potential early biomarkers for non-genotoxic carcinogenesis, including several genes aberrantly expressed in liver cancer. Our work suggests that 5hmC profiling can be used as an indicator of cell states during organ maturation and drug-induced responses and provides novel epigenetic signatures for non-genotoxic carcinogen exposure.

Comparative analysis of affinity-based 5-hydroxymethylation enrichment techniques.

The epigenetic modification of 5-hydroxymethylcytosine (5hmC) is receiving great attention due to its potential role in DNA methylation reprogramming and as a cell state identifier. Given this interest, it is important to identify reliable and cost-effective methods for the enrichment of 5hmC marked DNA for downstream analysis. We tested three commonly used affinity-based enrichment techniques; (i) antibody, (ii) chemical capture and (iii) protein affinity enrichment and assessed their ability to accurately and reproducibly report 5hmC profiles in mouse tissues containing high (brain) and lower (liver) levels of 5hmC. The protein-affinity technique is a poor reporter of 5hmC profiles, delivering 5hmC patterns that are incompatible with other methods. Both antibody and chemical capture-based techniques generate highly similar genome-wide patterns for 5hmC, which are independently validated by standard quantitative PCR (qPCR) and glucosyl-sensitive restriction enzyme digestion (gRES-qPCR). Both antibody and chemical capture generated profiles reproducibly link to unique chromatin modification profiles associated with 5hmC. However, there appears to be a slight bias of the antibody to bind to regions of DNA rich in simple repeats. Ultimately, the increased specificity observed with chemical capture-based approaches makes this an attractive method for the analysis of locus-specific or genome-wide patterns of 5hmC.

Whole exome sequencing reveals diverse genomic relatedness between paired concurrent endometrial and ovarian carcinomas.

INTRODUCTION: Concurrent non-serous endometrial and ovarian tumours are often managed clinically as two separate primary tumours, but almost all exhibit evidence of a genomic relationship. METHODOLOGY: This study investigates the extent of relatedness using whole exome sequencing, which was performed on paired non-serous endometrial and ovarian carcinomas from 27 patients with concurrent tumours in a cohort with detailed clinicopathological annotation. Four whole exome sequencing-derived parameters (mutation, mutational burden, mutational signatures and mutant allele tumour heterogeneity scores) were used to develop a novel unsupervised model for the assessment of genomic similarity to infer genomic relatedness of paired tumours. RESULTS: This novel model demonstrated genomic relatedness across all four parameters in all tumour pairs. Mutations in PTEN, ARID1A, CTNNB1, KMT2D and PIK3CA occurred most frequently and 24 of 27 (89 %) tumour pairs shared identical mutations in at least one of these genes, with all pairs sharing mutations in a number of other genes. Ovarian endometriosis, CTNNB1 exon 3 mutation, and progression and death from disease were present across the similarity ranking. Mismatch repair deficiency was associated with less genomically similar pairs. DISCUSSION: Although there was diversity across the cohort, the presence of genomic similarity in all paired tumours supports the hypothesis that concurrent non-serous endometrial and ovarian carcinomas are genomically related and may have arisen from a common origin.

The role of liquid biopsy in management of the neck with indeterminate response on post-treatment imaging following non-surgical management of oropharyngeal cancer.

INTRODUCTION: This study aimed to determine if post-treatment HPV cell-free DNA (cfDNA) can assist in the decision-making process for salvage neck dissection in patients following non-surgical treatment of oropharyngeal squamous cell carcinoma (OPSCC) with a partial response in the neck on imaging at 12 weeks post-treatment. METHODS: 86 patients who completed treatment were prospectively recruited through the regional multidisciplinary team (MDT). Treatment response was categorised as complete response (CR), partial response (PR) or progressive disease on 12-week post-treatment imaging. Pre- and post-treatment blood samples were assessed for HPV cfDNA through droplet digital PCR (ddPCR). RESULTS: Eight patients had an isolated partial response in the neck. One (12.5%) had detectable HPV cfDNA (22.96 copies/ml) at ∼12 weeks post-treatment with positive disease on subsequent neck dissection (positive predictive value; PPV = 100%). Of the seven patients with undetectable HPV cfDNA, two patients had evidence of regional disease recurrence at 23.9 and 27.4 months respectively (negative predictive value; NPV = 71%). CONCLUSION: The detection of HPV cfDNA may help target salvage therapy in patients with a partial response in the neck. Follow-up studies in larger cohorts would be required to further validate the use of post-treatment HPV cfDNA in the management of OPSCC.

Distinct histopathological features are associated with molecular subtypes and outcome in low grade serous ovarian carcinoma.

Low grade serous ovarian carcinoma (LGSOC) demonstrates unique clinical and molecular features compared to other ovarian cancer types. The relationship between common histological features of LGSOC and molecular events, such as hormone receptor expression patterns and MAPK gene mutation status, remains poorly understood. Recent data suggest some of these molecular features may be biomarkers of response to recently introduced biologically-targeted therapies, namely endocrine therapy and MEK inhibitors. We utilize a cohort of 63 pathologically-confirmed LGSOC cases with whole exome sequencing and hormone receptor expression data to investigate these relationships. LGSOC cases demonstrated uniformly high oestrogen receptor (ER) expression, but variable progesterone receptor (PR) expression intensity. 60% and 37% of cases demonstrated micropapillary and macropapillary patterns of stromal invasion, respectively. 63% of cases demonstrated desmoplasia, which was significantly associated with advanced disease stage and visible residual disease after cytoreductive surgery. MAPK-mutant cases (KRAS, BRAF, NRAS) more frequently demonstrated macropapillary stromal invasion, while Chr1p loss was associated with desmoplasia and low PR expression. Presence of micropapillary stromal invasion and low PR expression were associated with significantly poorer survival after accounting for stage and residual disease status. Together, these data identify novel relationships between histopathological features and molecularly-defined subgroups in LGSOC.

Longitudinal measurement of HPV copy number in cell-free DNA is associated with patient outcomes in HPV-positive oropharyngeal cancer.

INTRODUCTION: Oropharyngeal squamous cell carcinoma (OPSCC) is increasing in global prevalence and is divided into two types dependent on association with human papillomavirus (HPV). Assay of HPV copy number in plasma cell-free DNA (cfDNA) provides a minimally invasive method for detecting and monitoring tumour-derived HPV, with potential for enhancing clinical care. MATERIALS AND METHODS: In a prospectively recruited cohort of 104 OPSCC patients, we evaluate the utility of cfDNA droplet digital PCR (ddPCR) as a method for characterisation and longitudinal monitoring of patients with OPSCC. RESULTS: ddPCR assay of pre-treatment plasma cfDNA for five HPV types showed overall 95% concordance with p16 immunohistochemistry and PCR analysis of tumour tissue. Longitudinal sampling in 48 HPV+ve patients, with median follow-up of 20 months, was strongly associated with patient outcomes. Persistently elevated cfDNA-HPV post-treatment was associated with treatment failure (2/2 patients) and an increase of cfDNA-HPV in patients whose HPV levels were initially undetectable post-treatment was associated with disease recurrence (5/6 patients). No recurrence was observed in patients in whom cfDNA-HPV was undetectable in all post-treatment samples. In two patients, sequential HPV measurement could have avoided surgical intervention which did not confirm recurrence. CONCLUSION: The high concordance of pre-treatment plasma cfDNA-HPV analysis with tissue-based assays, together with the clinical associations of sequentially measured post-treatment cfDNA-HPV copy number add to a growing body of evidence that suggest utility of cfDNA-HPV ddPCR in management of OPSCC. Standardised clinical trials based on these data are now needed to assess the impact of such testing on overall patient outcomes.

Human biliary epithelial cells from discarded donor livers rescue bile duct structure and function in a mouse model of biliary disease.

Biliary diseases can cause inflammation, fibrosis, bile duct destruction, and eventually liver failure. There are no curative treatments for biliary disease except for liver transplantation. New therapies are urgently required. We have therefore purified human biliary epithelial cells (hBECs) from human livers that were not used for liver transplantation. hBECs were tested as a cell therapy in a mouse model of biliary disease in which the conditional deletion of Mdm2 in cholangiocytes causes senescence, biliary strictures, and fibrosis. hBECs are expandable and phenotypically stable and help restore biliary structure and function, highlighting their regenerative capacity and a potential alternative to liver transplantation for biliary disease.

Integrated molecular characterisation of endometrioid ovarian carcinoma identifies opportunities for stratification.

Endometrioid ovarian carcinoma (EnOC) is an under-investigated ovarian cancer type. Recent studies have described disease subtypes defined by genomics and hormone receptor expression patterns; here, we determine the relationship between these subtyping layers to define the molecular landscape of EnOC with high granularity and identify therapeutic vulnerabilities in high-risk cases. Whole exome sequencing data were integrated with progesterone and oestrogen receptor (PR and ER) expression-defined subtypes in 90 EnOC cases following robust pathological assessment, revealing dominant clinical and molecular features in the resulting integrated subtypes. We demonstrate significant correlation between subtyping approaches: PR-high (PR + /ER + , PR + /ER-) cases were predominantly CTNNB1-mutant (73.2% vs 18.4%, P < 0.001), while PR-low (PR-/ER + , PR-/ER-) cases displayed higher TP53 mutation frequency (38.8% vs 7.3%, P = 0.001), greater genomic complexity (P = 0.007) and more frequent copy number alterations (P = 0.001). PR-high EnOC patients experience favourable disease-specific survival independent of clinicopathological and genomic features (HR = 0.16, 95% CI 0.04-0.71). TP53 mutation further delineates the outcome of patients with PR-low tumours (HR = 2.56, 95% CI 1.14-5.75). A simple, routinely applicable, classification algorithm utilising immunohistochemistry for PR and p53 recapitulated these subtypes and their survival profiles. The genomic profile of high-risk EnOC subtypes suggests that inhibitors of the MAPK and PI3K-AKT pathways, alongside PARP inhibitors, represent promising candidate agents for improving patient survival. Patients with PR-low TP53-mutant EnOC have the greatest unmet clinical need, while PR-high tumours-which are typically CTNNB1-mutant and TP53 wild-type-experience excellent survival and may represent candidates for trials investigating de-escalation of adjuvant chemotherapy to agents such as endocrine therapy.

A human pluripotent stem cell model for the analysis of metabolic dysfunction in hepatic steatosis.

Nonalcoholic fatty liver disease (NAFLD) is currently the most prevalent form of liver disease worldwide. This term encompasses a spectrum of pathologies, from benign hepatic steatosis to non-alcoholic steatohepatitis, which have, to date, been challenging to model in the laboratory setting. Here, we present a human pluripotent stem cell (hPSC)-derived model of hepatic steatosis, which overcomes inherent challenges of current models and provides insights into the metabolic rewiring associated with steatosis. Following induction of macrovesicular steatosis in hepatocyte-like cells using lactate, pyruvate, and octanoate (LPO), respirometry and transcriptomic analyses revealed compromised electron transport chain activity. 13C isotopic tracing studies revealed enhanced TCA cycle anaplerosis, with concomitant development of a compensatory purine nucleotide cycle shunt leading to excess generation of fumarate. This model of hepatic steatosis is reproducible, scalable, and overcomes the challenges of studying mitochondrial metabolism in currently available models.

Neutrophils Fuel Effective Immune Responses through Gluconeogenesis and Glycogenesis.

Neutrophils can function and survive in injured and infected tissues, where oxygen and metabolic substrates are limited. Using radioactive flux assays and LC-MS tracing with U-13C glucose, glutamine, and pyruvate, we observe that neutrophils require the generation of intracellular glycogen stores by gluconeogenesis and glycogenesis for effective survival and bacterial killing. These metabolic adaptations are dynamic, with net increases in glycogen stores observed following LPS challenge or altitude-induced hypoxia. Neutrophils from patients with chronic obstructive pulmonary disease have reduced glycogen cycling, resulting in impaired function. Metabolic specialization of neutrophils may therefore underpin disease pathology and allow selective therapeutic targeting.

Non-alcoholic fatty liver disease (NAFLD) is associated with dynamic changes in DNA hydroxymethylation.

Non-alcoholic fatty liver disease (NAFLD) is now the commonest cause of liver disease in developed countries affecting 25-33% of the general population and up to 75% of those with obesity. Recent data suggest that alterations in DNA methylation may be related to NAFLD pathogenesis and progression and we have previously shown that dynamic changes in the cell lineage identifier 5-hydroxymethylcytosine (5hmC) may be important in the pathogenesis of liver disease. We used a model of diet-induced obesity, maintaining male mice on a high-fat diet (HFD) to generate hepatic steatosis. We profiled hepatic gene expression, global and locus-specific 5hmC and additionally investigated the effects of weight loss on the phenotype. HFD led to increased weight gain, fasting hyperglycaemia, glucose intolerance, insulin resistance and hepatic periportal macrovesicular steatosis. Diet-induced hepatic steatosis associated with reversible 5hmC changes at a discrete number of functionally important genes. We propose that 5hmC profiles are a useful signature of gene transcription and a marker of cell state in NAFLD and suggest that 5hmC profiles hold potential as a biomarker of abnormal liver physiology.

Molecular stratification of endometrioid ovarian carcinoma predicts clinical outcome.

Endometrioid ovarian carcinoma (EnOC) demonstrates substantial clinical and molecular heterogeneity. Here, we report whole exome sequencing of 112 EnOC cases following rigorous pathological assessment. We detect a high frequency of mutation in CTNNB1 (43%), PIK3CA (43%), ARID1A (36%), PTEN (29%), KRAS (26%), TP53 (26%) and SOX8 (19%), a recurrently-mutated gene previously unreported in EnOC. POLE and mismatch repair protein-encoding genes were mutated at lower frequency (6%, 18%) with significant co-occurrence. A molecular taxonomy is constructed, identifying clinically distinct EnOC subtypes: cases with TP53 mutation demonstrate greater genomic complexity, are commonly FIGO stage III/IV at diagnosis (48%), are frequently incompletely debulked (44%) and demonstrate inferior survival; conversely, cases with CTNNB1 mutation, which is mutually exclusive with TP53 mutation, demonstrate low genomic complexity and excellent clinical outcome, and are predominantly stage I/II at diagnosis (89%) and completely resected (87%). Moreover, we identify the WNT, MAPK/RAS and PI3K pathways as good candidate targets for molecular therapeutics in EnOC.

DNA Methylation Directs Polycomb-Dependent 3D Genome Re-organization in Naive Pluripotency.

The DNA hypomethylation that occurs when embryonic stem cells (ESCs) are directed to the ground state of naive pluripotency by culturing in two small molecule inhibitors (2i) results in redistribution of polycomb (H3K27me3) away from its target loci. Here, we demonstrate that 3D genome organization is also altered in 2i, with chromatin decompaction at polycomb target loci and a loss of long-range polycomb interactions. By preventing DNA hypomethylation during the transition to the ground state, we are able to restore to ESC in 2i the H3K27me3 distribution, as well as polycomb-mediated 3D genome organization that is characteristic of primed ESCs grown in serum. However, these cells retain the functional characteristics of 2i ground-state ESCs. Our findings demonstrate the central role of DNA methylation in shaping major aspects of 3D genome organization but caution against assuming causal roles for the epigenome and 3D genome in gene regulation and function in ESCs.

Affinity-Based Enrichment Techniques for the Genome-Wide Analysis of 5-Hydroxymethylcytosine.

Since its initial characterization in 2009 there has been a great degree of interest in comparative profiling of 5-hydroxymethylcytosine (5hmC) nucleotides in vertebrate DNA. Through a host of genome-wide studies the distribution of 5hmC has been mapped in a range of cell lines, tissue types and organisms; the majority of which have been generated through affinity-based methods for 5hmC enrichment. Although recent advances in the field have resulted in the ability to investigate the levels of both methylated and hydroxymethylated cytosines at single base resolution, such studies are still relatively cost-prohibitive as well as technically challenging. As such affinity-based methods for the enrichment of 5hmC-modified DNA fragments represent a cost-effective and highly informative method for profiling 5hmC residency in genomic DNA. Here we will outline protocols for two independent affinity based methods to generate 5hmC enriched fractions for subsequent locus specific and genome-wide analysis; immunoprecipitation using highly specific 5hmC antibodies as well as a chemical capture based method.