Search

Colorectal carcinoma (CRC) is a common cause of mortality 1 , but a comprehensive description of its genomic landscape is lacking 2-9 . Here we perform whole-genome sequencing of 2,023 CRC samples from participants in the UK 100,000 Genomes Project, thereby providing a highly detailed somatic mutational landscape of this cancer. Integrated analyses identify more than 250 putative CRC driver genes, many not previously implicated in CRC or other cancers, including several recurrent changes outside the coding genome. We extend the molecular pathways involved in CRC development, define four new common subgroups of microsatellite-stable CRC based on genomic features and show that these groups have independent prognostic associations. We also characterize several rare molecular CRC subgroups, some with potential clinical relevance, including cancers with both microsatellite and chromosomal instability. We demonstrate a spectrum of mutational profiles across the colorectum, which reflect aetiological differences. These include the role of Escherichia coli pks+ colibactin in rectal cancers 10 and the importance of the SBS93 signature 11-13 , which suggests that diet or smoking is a risk factor. Immune-escape driver mutations 14 are near-ubiquitous in hypermutant tumours and occur in about half of microsatellite-stable CRCs, often in the form of HLA copy number changes. Many driver mutations are actionable, including those associated with rare subgroups (for example, BRCA1 and IDH1), highlighting the role of whole-genome sequencing in optimizing patient care., (© 2024. The Author(s).)

Background: Mutations in the gene MTARC1 (mitochondrial amidoxime-reducing component 1) protect carriers from metabolic dysfunction-associated steatohepatitis (MASH) and cirrhosis. MTARC1 encodes the mARC1 enzyme, which is localized to the mitochondria and has no known MASH-relevant molecular function. Our studies aimed to expand on the published human genetic mARC1 data and to observe the molecular effects of mARC1 modulation in preclinical MASH models., Methods and Results: We identified a novel human structural variant deletion in MTARC1, which is associated with various biomarkers of liver health, including alanine aminotransferase levels. Phenome-wide Mendelian Randomization analyses additionally identified novel putatively causal associations between MTARC1 expression, and esophageal varices and cardiorespiratory traits. We observed that protective MTARC1 variants decreased protein accumulation in in vitro overexpression systems and used genetic tools to study mARC1 depletion in relevant human and mouse systems. Hepatocyte mARC1 knockdown in murine MASH models reduced body weight, liver steatosis, oxidative stress, cell death, and fibrogenesis markers. mARC1 siRNA treatment and overexpression modulated lipid accumulation and cell death consistently in primary human hepatocytes, hepatocyte cell lines, and primary human adipocytes. mARC1 depletion affected the accumulation of distinct lipid species and the expression of inflammatory and mitochondrial pathway genes/proteins in both in vitro and in vivo models., Conclusions: Depleting hepatocyte mARC1 improved metabolic dysfunction-associated steatotic liver disease-related outcomes. Given the functional role of mARC1 in human adipocyte lipid accumulation, systemic targeting of mARC1 should be considered when designing mARC1 therapies. Our data point to plasma lipid biomarkers predictive of mARC1 abundance, such as Ceramide 22:1. We propose future areas of study to describe the precise molecular function of mARC1, including lipid trafficking and subcellular location within or around the mitochondria and endoplasmic reticulum., (Copyright © 2024 The Author(s). Published by Wolters Kluwer Health, Inc. on behalf of the American Association for the Study of Liver Diseases.)

The Cancer Programme of the 100,000 Genomes Project was an initiative to provide whole-genome sequencing (WGS) for patients with cancer, evaluating opportunities for precision cancer care within the UK National Healthcare System (NHS). Genomics England, alongside NHS England, analyzed WGS data from 13,880 solid tumors spanning 33 cancer types, integrating genomic data with real-world treatment and outcome data, within a secure Research Environment. Incidence of somatic mutations in genes recommended for standard-of-care testing varied across cancer types. For instance, in glioblastoma multiforme, small variants were present in 94% of cases and copy number aberrations in at least one gene in 58% of cases, while sarcoma demonstrated the highest occurrence of actionable structural variants (13%). Homologous recombination deficiency was identified in 40% of high-grade serous ovarian cancer cases with 30% linked to pathogenic germline variants, highlighting the value of combined somatic and germline analysis. The linkage of WGS and longitudinal life course clinical data allowed the assessment of treatment outcomes for patients stratified according to pangenomic markers. Our findings demonstrate the utility of linking genomic and real-world clinical data to enable survival analysis to identify cancer genes that affect prognosis and advance our understanding of how cancer genomics impacts patient outcomes., (© 2024. The Author(s).)

Inherited bone marrow failure associated with heterozygous mutations in GATA2 predisposes toward hematological malignancies, but the mechanisms remain poorly understood. Here, we investigate the mechanistic basis of marrow failure in a zebrafish model of GATA2 deficiency. Single-cell transcriptomics and chromatin accessibility assays reveal that loss of gata2a leads to skewing toward the erythroid lineage at the expense of myeloid cells, associated with loss of cebpa expression and decreased PU.1 and CEBPA transcription factor accessibility in hematopoietic stem and progenitor cells (HSPCs). Furthermore, gata2a mutants show impaired expression of npm1a, the zebrafish NPM1 ortholog. Progressive loss of npm1a in HSPCs is associated with elevated levels of DNA damage in gata2a mutants. Thus, Gata2a maintains myeloid lineage priming through cebpa and protects against genome instability and marrow failure by maintaining expression of npm1a. Our results establish a potential mechanism underlying bone marrow failure in GATA2 deficiency., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 The Author(s). Published by Elsevier Inc. All rights reserved.)

Background: Three-quarters of bladder cancer patients present with early-stage disease (non-muscle-invasive bladder cancer, NMIBC, UICC TNM stages Ta, T1 and Tis); however, most next-generation sequencing studies to date have concentrated on later-stage disease (muscle-invasive BC, stages T2+). We used exome and transcriptome sequencing to comprehensively characterise NMIBCs of all grades and stages to identify prognostic genes and pathways that could facilitate treatment decisions. Tumour grading is based upon microscopy and cellular appearances (grade 1 BCs are less aggressive, and grade 3 BCs are most aggressive), and we chose to also focus on the most clinically complex NMIBC subgroup, those patients with grade 3 pathological stage T1 (G3 pT1) disease., Methods: Whole-exome and RNA sequencing were performed in total on 96 primary NMIBCs including 22 G1 pTa, 14 G3 pTa and 53 G3 pT1s, with both exome and RNA sequencing data generated from 75 of these individual samples. Associations between genomic alterations, expression profiles and progression-free survival (PFS) were investigated., Results: NMIBCs clustered into 3 expression subtypes with different somatic alteration characteristics. Amplifications of ARNT and ERBB2 were significant indicators of worse PFS across all NMIBCs. High APOBEC mutagenesis and high tumour mutation burden were both potential indicators of better PFS in G3pT1 NMIBCs. The expression of individual genes was not prognostic in BCG-treated G3pT1 NMIBCs; however, downregulated interferon-alpha and gamma response pathways were significantly associated with worse PFS (adjusted p-value < 0.005)., Conclusions: Multi-omic data may facilitate better prognostication and selection of therapeutic interventions in patients with G3pT1 NMIBC. These findings demonstrate the potential for improving the management of high-risk NMIBC patients and warrant further prospective validation., (© 2022. The Author(s).)

We report an autosomal recessive, multi-organ tumor predisposition syndrome, caused by bi-allelic loss-of-function germline variants in the base excision repair (BER) gene MBD4. We identified five individuals with bi-allelic MBD4 variants within four families and these individuals had a personal and/or family history of adenomatous colorectal polyposis, acute myeloid leukemia, and uveal melanoma. MBD4 encodes a glycosylase involved in repair of G:T mismatches resulting from deamination of 5'-methylcytosine. The colorectal adenomas from MBD4-deficient individuals showed a mutator phenotype attributable to mutational signature SBS1, consistent with the function of MBD4. MBD4-deficient polyps harbored somatic mutations in similar driver genes to sporadic colorectal tumors, although AMER1 mutations were more common and KRAS mutations less frequent. Our findings expand the role of BER deficiencies in tumor predisposition. Inclusion of MBD4 in genetic testing for polyposis and multi-tumor phenotypes is warranted to improve disease management., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 The Author(s). Published by Elsevier Inc. All rights reserved.)

The ubiquitous host protein, CCCTC-binding factor (CTCF), is an essential regulator of cellular transcription and functions to maintain epigenetic boundaries, stabilise chromatin loops and regulate splicing of alternative exons. We have previously demonstrated that CTCF binds to the E2 open reading frame (ORF) of human papillomavirus (HPV) 18 and functions to repress viral oncogene expression in undifferentiated keratinocytes by co-ordinating an epigenetically repressed chromatin loop within HPV episomes. Keratinocyte differentiation disrupts CTCF-dependent chromatin looping of HPV18 episomes promoting induction of enhanced viral oncogene expression. To further characterise CTCF function in HPV transcription control we utilised direct, long-read Nanopore RNA-sequencing which provides information on the structure and abundance of full-length transcripts. Nanopore analysis of primary human keratinocytes containing HPV18 episomes before and after synchronous differentiation allowed quantification of viral transcript species, including the identification of low abundance novel transcripts. Comparison of transcripts produced in wild type HPV18 genome-containing cells to those identified in CTCF-binding deficient genome-containing cells identifies CTCF as a key regulator of differentiation-dependent late promoter activation, required for efficient E1^E4 and L1 protein expression. Furthermore, our data show that CTCF binding at the E2 ORF promotes usage of the downstream weak splice donor (SD) sites SD3165 and SD3284, to the dominant E4 splice acceptor site at nucleotide 3434. These findings demonstrate that in the HPV life cycle both early and late virus transcription programmes are facilitated by recruitment of CTCF to the E2 ORF., Competing Interests: The authors have declared that no competing interests exist.

Introduction: Research demonstrates strong evidence that circulating tumour cells (CTCs) can provide diagnostic and/or prognostic biomarkers in head and neck squamous cell carcinoma (HNSCC) and a potential tool for therapeutic stratification. However, the question still remains as to the optimum method of CTC enrichment and how this can be translated into clinical practice. We aimed to evaluate the Parsortix microfluidic device for CTC enrichment and characterisation in HNSCC, seeking to optimise a sample collection and processing protocol that preserves CTC integrity and phenotype., Method: Spiking experiments of the FaDu and SCC040 HNSCC cell lines were used to determine the Parsortix capture rate of rare "CTC-like" cells. Capture rates of cancer cells spiked into EDTA blood collections tubes (BCTs) were compared to the Transfix fixative BCT and Cytodelics whole blood freezing protocol. The Lexogen Quantseq library preparation was used to profile gene expression of unfixed cells before and after microfluidic enrichment and enriched cell line spiked Transfix blood samples. An antibody panel was optimised to enable immunofluorescence microscopy CTC detection in HNSCC patient Transfix blood samples, using epithelial (EpCAM) and mesenchymal (N-cadherin) CTC markers., Results: Across a spiked cell concentration range of 9-129 cells/mL, Parsortix demonstrated a mean cell capture rate of 53.5% for unfixed cells, with no significant relationship between spiked cell concentration and capture rate. Samples preserved in Transfix BCTs demonstrated significantly increased capture rates at 0 h (time to processing) compared to EDTA BCTs (65.3% vs. 51.0%). Capture rates in Transfix BCTs were maintained at 24 h and 72 h timepoints, but dropped significantly in EDTA BCTs. Gene expression profiling revealed that microfluidic enrichment of unfixed cell lines caused downregulation of RNA processing/binding gene pathways and upregulation of genes involved in cell injury, apoptosis and oxidative stress. RNA was successfully extracted and sequenced from Transfix preserved cells enriched using Parsortix, demonstrating epithelial specific transcripts from spiked cells. In a proof-of-concept cohort of four patients with advanced HNSCC, CTCs were successfully identified and visualised with epithelial and epithelial-mesenchymal phenotypes., Conclusion: We have optimised a protocol for detection of CTCs in HNSCC with the Parsortix microfluidic device, using Transfix BCTs. We report a significant benefit, both in terms of cell capture rates and preserving cell phenotype, for using a fixative BCT- particularly if samples are stored before processing. In the design of large cohort multi-site clinical trials, such data are of paramount importance.

Background: Targeting of MDSCs is a major clinical challenge in the era of immunotherapy. Antibodies which deplete MDSCs in murine models can reactivate T cell responses. In humans such approaches have not developed due to difficulties in identifying targets amenable to clinical translation., Methods: RNA-sequencing of M-MDSCs and G-MDSCs from cancer patients was undertaken. Flow cytometry and immunohistochemistry of blood and tumours determined MDSC CD33 expression. MDSCs were treated with Gemtuzumab ozogamicin and internalisation kinetics, and cell death mechanisms determined by flow cytometry, confocal microscopy and electron microscopy. Effects on T cell proliferation and CAR-T cell anti-tumour cytotoxicity were identified in the presence of Gemtuzumab ozogamicin., Findings: RNA-sequencing of human M-MDSCs and G-MDSCs identified transcriptomic differences, but that CD33 is a common surface marker. Flow cytometry indicated CD33 expression is higher on M-MDSCs, and CD33+ MDSCs are found in the blood and tumours regardless of cancer subtype. Treatment of human MDSCs leads to Gemtuzumab ozogamicin internalisation, increased p-ATM, and cell death; restoring T cell proliferation. Anti-GD2-/mesothelin-/EGFRvIII-CAR-T cell activity is enhanced in combination with the anti-MDSC effects of Gemtuzumab ozogamicin., Interpretation: The study identifies that M-MDSCs and G-MDSCs are transcriptomically different but CD33 is a therapeutic target on peripheral and infiltrating MDSCs across cancer subtypes. The immunotoxin Gemtuzumab ozogamicin can deplete MDSCs providing a translational approach to reactivate T cell and CAR-T cell responses against multiple cancers. In the rare conditions of HLH/MAS gemtuzumab ozogamicin provides a novel anti-myeloid strategy. FUND: This work was supported by Cancer Research UK, CCLG, Treating Children with Cancer, and the alumni and donors to the University of Birmingham., (Copyright © 2019 The Authors. Published by Elsevier B.V. All rights reserved.)

Severe Malarial Anemia (SMA), a life-threatening childhood Plasmodium falciparum malaria syndrome requiring urgent blood transfusion, exhibits inflammatory and hemolytic pathology. Differentiating between hypo-haptoglobinemia due to hemolysis or that of genetic origin is key to understand SMA pathogenesis. We hypothesized that while malaria-induced hypo-haptoglobinemia should reverse at recovery, that of genetic etiology should not. We carried-out a case-control study of children living under hyper-endemic holoendemic malaria burden in the sub-Saharan metropolis of Ibadan, Nigeria. We show that hypo-haptoglobinemia is a risk factor for childhood SMA and not solely due to intravascular hemolysis from underlying schizogony. In children presenting with SMA, hypo-haptoglobinemia remains through convalescence to recovery suggesting a genetic cause. We identified a haptoglobin gene variant, rs12162087 (g.-1203G > A, frequency = 0.67), to be associated with plasma haptoglobin levels (p = 8.5 × 10 -6 ). The Homo-Var:(AA) is associated with high plasma haptoglobin while the reference Homo-Ref:(GG) is associated with hypo-haptoglobinemia (p = 2.3 × 10 -6 ). The variant is associated with SMA, with the most support for a risk effect for Homo-Ref genotype. Our insights on regulatory haptoglobin genotypes and hypo-haptoglobinemia suggest that haptoglobin screening could be part of risk-assessment algorithms to prevent rapid disease progression towards SMA in regions with no-access to urgent blood transfusion where SMA accounts for high childhood mortality rates.

The complex life cycle of oncogenic human papillomavirus (HPV) initiates in undifferentiated basal epithelial keratinocytes where expression of the E6 and E7 oncogenes is restricted. Upon epithelial differentiation, E6/E7 transcription is increased through unknown mechanisms to drive cellular proliferation required to support virus replication. We report that the chromatin-organising CCCTC-binding factor (CTCF) promotes the formation of a chromatin loop in the HPV genome that epigenetically represses viral enhancer activity controlling E6/E7 expression. CTCF-dependent looping is dependent on the expression of the CTCF-associated Yin Yang 1 (YY1) transcription factor and polycomb repressor complex (PRC) recruitment, resulting in trimethylation of histone H3 at lysine 27. We show that viral oncogene up-regulation during cellular differentiation results from YY1 down-regulation, disruption of viral genome looping, and a loss of epigenetic repression of viral enhancer activity. Our data therefore reveal a key role for CTCF-YY1-dependent looping in the HPV life cycle and identify a regulatory mechanism that could be disrupted in HPV carcinogenesis., Competing Interests: The authors have declared that no competing interests exist.

Arginine is a semi-essential amino acid that plays a key role in cell survival and proliferation in normal and malignant cells. BCT-100, a pegylated (PEG) recombinant human arginase, can deplete arginine and starve malignant cells of the amino acid. Acute lymphoblastic leukemia (ALL) is the most common cancer of childhood, yet for patients with high risk or relapsed disease prognosis remains poor. We show that BCT-100 is cytotoxic to ALL blasts from patients in vitro by necrosis, and is synergistic in combination with dexamethasone. Against ALL xenografts, BCT-100 leads to a reduction in ALL engraftment and a prolongation of survival. ALL blasts express the arginine transporter CAT-1, yet the majority of blasts are arginine auxotrophic due to deficiency in either argininosuccinate synthase (ASS) or ornithine transcarbamylase (OTC). Although endogenous upregulation or retroviral transduced increases in ASS or OTC may promote ALL survival under moderately low arginine conditions, expression of these enzymes cannot prevent BCT-100 cytotoxicity at arginine depleting doses. RNA-sequencing of ALL blasts and supporting stromal cells treated with BCT-100 identifies a number of candidate pathways which are altered in the presence of arginine depletion. Therefore, BCT-100 provides a new clinically relevant therapeutic approach to target arginine metabolism in ALL., (© 2017 The Authors International Journal of Cancer published by John Wiley & Sons Ltd on behalf of UICC.)

Chronic obstructive pulmonary disease (COPD) is characterized by reduced lung function and is the third leading cause of death globally. Through genome-wide association discovery in 48,943 individuals, selected from extremes of the lung function distribution in UK Biobank, and follow-up in 95,375 individuals, we increased the yield of independent signals for lung function from 54 to 97. A genetic risk score was associated with COPD susceptibility (odds ratio per 1 s.d. of the risk score (∼6 alleles) (95% confidence interval) = 1.24 (1.20-1.27), P = 5.05 × 10 -49 ), and we observed a 3.7-fold difference in COPD risk between individuals in the highest and lowest genetic risk score deciles in UK Biobank. The 97 signals show enrichment in genes for development, elastic fibers and epigenetic regulation pathways. We highlight targets for drugs and compounds in development for COPD and asthma (genes in the inositol phosphate metabolism pathway and CHRM3) and describe targets for potential drug repositioning from other clinical indications.

Background: Association studies have identified a number of loci that contribute to an increased body mass index (BMI), the strongest of which is in the first intron of the FTO gene on human chromosome 16q12.2. However, this region is both non-coding and under strong linkage disequilibrium, making it recalcitrant to functional interpretation. Furthermore, the FTO gene is located within a complex cis-regulatory landscape defined by a topologically associated domain that includes the IRXB gene cluster, a trio of developmental regulators. Consequently, at least three genes in this interval have been implicated in the aetiology of obesity., Methods: Here, we sequence a 2 Mb region encompassing the FTO, RPGRIP1L and IRXB cluster genes in 284 individuals from a well-characterised study group of Danish men containing extremely overweight young adults and controls. We further replicate our findings both in an expanded male cohort and an independent female study group. Finally, we compare our variant data with a previous study describing IRX3 and FTO interactions in this region., Results: We obtain deep coverage across the entire region, allowing accurate and unequivocal determination of almost every single nucleotide polymorphism and short insertion/deletion. As well as confirming previous findings across the interval, we identify a further novel age-dependent association upstream of IRX5 that imposes a similar burden on BMI to the FTO locus., Conclusions: Our findings are consistent with the hypothesis that chromatin architectures play a role in regulating gene expression levels across topological domains while our targeted sequence approach represents a widely applicable methodology for high-resolution analysis of regional variation across candidate genomic loci.

Background: Determining the function of regulatory elements is fundamental for our understanding of development, disease and evolution. However, the sequence features that mediate these functions are often unclear and the prediction of tissue-specific expression patterns from sequence alone is non-trivial. Previous functional studies have demonstrated a link between PBX-HOX and MEIS/PREP binding interactions and hindbrain enhancer activity, but the defining grammar of these sites, if any exists, has remained elusive., Results: Here, we identify a shared sequence signature (syntax) within a heterogeneous set of conserved vertebrate hindbrain enhancers composed of spatially co-occurring PBX-HOX and MEIS/PREP transcription factor binding motifs. We use this syntax to accurately predict hindbrain enhancers in 89% of cases (67/75 predicted elements) from a set of conserved non-coding elements (CNEs). Furthermore, mutagenesis of the sites abolishes activity or generates ectopic expression, demonstrating their requirement for segmentally restricted enhancer activity in the hindbrain. We refine and use our syntax to predict over 3,000 hindbrain enhancers across the human genome. These sequences tend to be located near developmental transcription factors and are enriched in known hindbrain activating elements, demonstrating the predictive power of this simple model., Conclusion: Our findings support the theory that hundreds of CNEs, and perhaps thousands of regions across the human genome, function to coordinate gene expression in the developing hindbrain. We speculate that deeply conserved sequences of this kind contributed to the co-option of new genes into the hindbrain gene regulatory network during early vertebrate evolution by linking patterns of hox expression to downstream genes involved in segmentation and patterning, and evolutionarily newer instances may have continued to contribute to lineage-specific elaboration of the hindbrain.