Molecular landscape of proliferative verrucous leukoplakia: a systematic review.

Proliferative verrucous leukoplakia (PVL) is a rare oral potentially malignant disorder characterised by multifocal origin and unpredictable long-term evolution to oral squamous cell carcinoma (OSCC) or oral verrucous carcinoma (OVC). Currently no predictive biomarkers are in clinical use. We aimed to explore the genomic profile of PVL. A total of 685 cases in 26 studies were included in this review. Genomic data were presented in 15% of studies and biomarker analysis was reported in 85% of studies. At first clinical presentation, PVL is characterised by a high loss of heterozygosity (LOH), similar to OSCC, and low copy number alterations (CNA). As these progress, more CNAs and mutations in CDKN2A and alterations to ELAVL1 expression are noted, but no TP53 mutations are identified. There is significantly lower LOH at 17p in early PVL compared with OSCC (p = 0.037). Deletions in chromosomal loci 17q12, 5q31.1 and amplifications in 7q11.2, 7q22 are shared between early lesions and OVC. PVL shows CNAs at 11q31. WNT signalling pathway genes (SUZ12, CTTN and FOLR3) are enriched in CN-altered regions. PVL stroma shows significantly lower α-SMA and higher CD34 expression than OVC and OSCC. The exact genomic landscape is currently unclear, and further studies are necessary to unravel this mystery.

Single-cell analysis reveals prognostic fibroblast subpopulations linked to molecular and immunological subtypes of lung cancer.

Fibroblasts are poorly characterised cells that variably impact tumour progression. Here, we use single cell RNA-sequencing, multiplexed immunohistochemistry and digital cytometry (CIBERSORTx) to identify and characterise three major fibroblast subpopulations in human non-small cell lung cancer: adventitial, alveolar and myofibroblasts. Alveolar and adventitial fibroblasts (enriched in control tissue samples) localise to discrete spatial niches in histologically normal lung tissue and indicate improved overall survival rates when present in lung adenocarcinomas (LUAD). Trajectory inference identifies three phases of control tissue fibroblast activation, leading to myofibroblast enrichment in tumour samples: initial upregulation of inflammatory cytokines, followed by stress-response signalling and ultimately increased expression of fibrillar collagens. Myofibroblasts correlate with poor overall survival rates in LUAD, associated with loss of epithelial differentiation, TP53 mutations, proximal molecular subtypes and myeloid cell recruitment. In squamous carcinomas myofibroblasts were not prognostic despite being transcriptomically equivalent. These findings have important implications for developing fibroblast-targeting strategies for cancer therapy.

The spatiotemporal dynamics of microglia across the human lifespan.

Microglia, the brain's resident macrophages, shape neural development and are key neuroimmune hubs in the pathological signatures of neurodevelopmental disorders. Despite the importance of microglia, their development has not been carefully examined in the human brain, and most of our knowledge derives from rodents. We aimed to address this gap in knowledge by establishing an extensive collection of 97 post-mortem tissues in order to enable quantitative, sex-matched, detailed analysis of microglia across the human lifespan. We identify the dynamics of these cells in the human telencephalon, describing waves in microglial density across gestation, infancy, and childhood, controlled by a balance of proliferation and apoptosis, which track key neurodevelopmental milestones. These profound changes in microglia are also observed in bulk RNA-seq and single-cell RNA-seq datasets. This study provides a detailed insight into the spatiotemporal dynamics of microglia across the human lifespan and serves as a foundation for elucidating how microglia contribute to shaping neurodevelopment in humans.

Genomic Analysis of Response to Neoadjuvant Chemotherapy in Esophageal Adenocarcinoma.

Neoadjuvant therapy followed by surgery is the standard of care for locally advanced esophageal adenocarcinoma (EAC). Unfortunately, response to neoadjuvant chemotherapy (NAC) is poor (20-37%), as is the overall survival benefit at five years (9%). The EAC genome is complex and heterogeneous between patients, and it is not yet understood whether specific mutational patterns may result in chemotherapy sensitivity or resistance. To identify associations between genomic events and response to NAC in EAC, a comparative genomic analysis was performed in 65 patients with extensive clinical and pathological annotation using whole-genome sequencing (WGS). We defined response using Mandard Tumor Regression Grade (TRG), with responders classified as TRG1-2 (n = 27) and non-responders classified as TRG4-5 (n =38). We report a higher non-synonymous mutation burden in responders (median 2.08/Mb vs. 1.70/Mb, p = 0.036) and elevated copy number variation in non-responders (282 vs. 136/patient, p < 0.001). We identified copy number variants unique to each group in our cohort, with cell cycle (CDKN2A, CCND1), c-Myc (MYC), RTK/PIK3 (KRAS, EGFR) and gastrointestinal differentiation (GATA6) pathway genes being specifically altered in non-responders. Of note, NAV3 mutations were exclusively present in the non-responder group with a frequency of 22%. Thus, lower mutation burden, higher chromosomal instability and specific copy number alterations are associated with resistance to NAC.

Single-nucleotide Fcγ receptor polymorphisms do not impact obinutuzumab/rituximab outcome in patients with lymphoma.

Single-nucleotide polymorphisms (SNPs) have been shown to influence Fcγ receptor (FcγR) affinity and activity, but their effect on treatment response is unclear. We assessed their importance in the efficacy of obinutuzumab or rituximab combined with chemotherapy in untreated advanced follicular lymphoma (FL) and diffuse large B-cell lymphoma (DLBCL) in the GALLIUM (www.clinicaltrials.gov #NCT01332968) and GOYA (#NCT01287741) trials, respectively. Genomic DNA was extracted from patients enrolled in GALLIUM (n = 1202) and GOYA (n = 1418). Key germline SNPs, FCGR2A R131H (rs1801274), FCGR3A F158V (rs396991), and FCGR2B I232T (rs1050501), were genotyped and assessed for their impact on investigator-assessed progression-free survival (PFS). In both cohorts there was no prognostic effect of FCGR2A or FCGR3A. In FL, FCGR2B was associated with favorable PFS in univariate and multivariate analyses comparing I232T with I232I, with a more modest association for rituximab-treated (univariate: hazard ratio [HR], 0.78; 95% confidence interval [CI], 0.54-1.14; P = .21) vs obinutuzumab-treated patients (HR, 0.56; 95% CI, 0.34-0.91; P = .02). Comparing T232T with I232I, an association was found for obinutuzumab (univariate: HR, 2.76; 95% CI, 1.02-7.5; P = .0459). Neither observation retained significance after multiple-test adjustment. FCGR2B was associated with poorer PFS in multivariate analyses comparing T232T with I232I in rituximab- but not obinutuzumab-treated patients with DLBCL (HR, 4.40; 95% CI, 1.71-11.32; P = .002; multiple-test-adjusted P = .03); however, this genotype was rare (n = 13). This study shows that FcγR genotype is not associated with response to rituximab/obinutuzumab plus chemotherapy in treatment-naive patients with advanced FL or DLBCL.

p53 is regulated by aerobic glycolysis in cancer cells by the CtBP family of NADH-dependent transcriptional regulators.

High rates of glycolysis in cancer cells are a well-established characteristic of many human tumors, providing rapidly proliferating cancer cells with metabolites that can be used as precursors for anabolic pathways. Maintenance of high glycolytic rates depends on the lactate dehydrogenase-catalyzed regeneration of NAD+ from GAPDH-generated NADH because an increased NADH:NAD+ ratio inhibits GAPDH. Here, using human breast cancer cell models, we identified a pathway in which changes in the extramitochondrial-free NADH:NAD+ ratio signaled through the CtBP family of NADH-sensitive transcriptional regulators to control the abundance and activity of p53. NADH-free forms of CtBPs cooperated with the p53-binding partner HDM2 to suppress p53 function, and loss of these forms in highly glycolytic cells resulted in p53 accumulation. We propose that this pathway represents a "glycolytic stress response" in which the initiation of a protective p53 response by an increased NADH:NAD+ ratio enables cells to avoid cellular damage caused by mismatches between metabolic supply and demand.

Clinical significance of TP53, BIRC3, ATM and MAPK-ERK genes in chronic lymphocytic leukaemia: data from the randomised UK LRF CLL4 trial.

Despite advances in chronic lymphocytic leukaemia (CLL) treatment, globally chemotherapy remains a central treatment modality, with chemotherapy trials representing an invaluable resource to explore disease-related/genetic features contributing to long-term outcomes. In 499 LRF CLL4 cases, a trial with >12 years follow-up, we employed targeted resequencing of 22 genes, identifying 623 mutations. After background mutation rate correction, 11/22 genes were recurrently mutated at frequencies between 3.6% (NFKBIE) and 24% (SF3B1). Mutations beyond Sanger resolution (<12% VAF) were observed in all genes, with KRAS mutations principally composed of these low VAF variants. Firstly, employing orthogonal approaches to confirm <12% VAF TP53 mutations, we assessed the clinical impact of TP53 clonal architecture. Whilst ≥ 12% VAF TP53mut cases were associated with reduced PFS and OS, we could not demonstrate a difference between <12% VAF TP53 mutations and either wild type or ≥12% VAF TP53mut cases. Secondly, we identified biallelic BIRC3 lesions (mutation and deletion) as an independent marker of inferior PFS and OS. Finally, we observed that mutated MAPK-ERK genes were independent markers of poor OS in multivariate survival analysis. In conclusion, our study supports using targeted resequencing of expanded gene panels to elucidate the prognostic impact of gene mutations.

Single-cell analysis based dissection of clonality in myelofibrosis.

Cancer development is an evolutionary genomic process with parallels to Darwinian selection. It requires acquisition of multiple somatic mutations that collectively cause a malignant phenotype and continuous clonal evolution is often linked to tumor progression. Here, we show the clonal evolution structure in 15 myelofibrosis (MF) patients while receiving treatment with JAK inhibitors (mean follow-up 3.9 years). Whole-exome sequencing at multiple time points reveal acquisition of somatic mutations and copy number aberrations over time. While JAK inhibition therapy does not seem to create a clear evolutionary bottleneck, we observe a more complex clonal architecture over time, and appearance of unrelated clones. Disease progression associates with increased genetic heterogeneity and gain of RAS/RTK pathway mutations. Clonal diversity results in clone-specific expansion within different myeloid cell lineages. Single-cell genotyping of circulating CD34 + progenitor cells allows the reconstruction of MF phylogeny demonstrating loss of heterozygosity and parallel evolution as recurrent events.

Clinical significance of DNA methylation in chronic lymphocytic leukemia patients: results from 3 UK clinical trials.

Chronic lymphocytic leukemia patients with mutated immunoglobulin heavy-chain genes (IGHV-M), particularly those lacking poor-risk genomic lesions, often respond well to chemoimmunotherapy (CIT). DNA methylation profiling can subdivide early-stage patients into naive B-cell-like CLL (n-CLL), memory B-cell-like CLL (m-CLL), and intermediate CLL (i-CLL), with differing times to first treatment and overall survival. However, whether DNA methylation can identify patients destined to respond favorably to CIT has not been ascertained. We classified treatment-naive patients (n = 605) from 3 UK chemo and CIT clinical trials into the 3 epigenetic subgroups, using pyrosequencing and microarray analysis, and performed expansive survival analysis. The n-CLL, i-CLL, and m-CLL signatures were found in 80% (n = 245/305), 17% (53/305), and 2% (7/305) of IGHV-unmutated (IGHV-U) cases, respectively, and in 9%, (19/216), 50% (108/216), and 41% (89/216) of IGHV-M cases, respectively. Multivariate Cox proportional analysis identified m-CLL as an independent prognostic factor for overall survival (hazard ratio [HR], 0.46; 95% confidence interval [CI], 0.24-0.87; P = .018) in CLL4, and for progression-free survival (HR, 0.25; 95% CI, 0.10-0.57; P = .002) in ARCTIC and ADMIRE patients. The analysis of epigenetic subgroups in patients entered into 3 first-line UK CLL trials identifies m-CLL as an independent marker of prolonged survival and may aid in the identification of patients destined to demonstrate prolonged survival after CIT.

An optimised tissue disaggregation and data processing pipeline for characterising fibroblast phenotypes using single-cell RNA sequencing.

Single-cell RNA sequencing (scRNA-Seq) provides a valuable platform for characterising multicellular ecosystems. Fibroblasts are a heterogeneous cell type involved in many physiological and pathological processes, but remain poorly-characterised. Analysis of fibroblasts is challenging: these cells are difficult to isolate from tissues, and are therefore commonly under-represented in scRNA-seq datasets. Here, we describe an optimised approach for fibroblast isolation from human lung tissues. We demonstrate the potential for this procedure in characterising stromal cell phenotypes using scRNA-Seq, analyse the effect of tissue disaggregation on gene expression, and optimise data processing to improve clustering quality. We also assess the impact of in vitro culture conditions on stromal cell gene expression and proliferation, showing that altering these conditions can skew phenotypes.

Stem cell-like breast cancer cells with acquired resistance to metformin are sensitive to inhibitors of NADH-dependent CtBP dimerization.

Altered flux through major metabolic pathways is a hallmark of cancer cells and provides opportunities for therapy. Stem cell-like cancer (SCLC) cells can cause metastasis and therapy resistance. They possess metabolic plasticity, theoretically enabling resistance to therapies targeting a specific metabolic state. The C-terminal binding protein (CtBP) transcriptional regulators are potential therapeutic targets in highly glycolytic cancer cells, as they are activated by the glycolytic coenzyme nicotinamide adenine dinucleotide (NADH). However, SCLC cells commonly exist in an oxidative state with low rates of glycolysis. Metformin inhibits complex I of the mitochondrial electron transport chain; it can kill oxidative SCLC cells and has anti-cancer activity in patients. SCLC cells can acquire resistance to metformin through increased glycolysis. Given the potential for long-term metformin therapy, we have studied acquired metformin resistance in cells from the claudin-low subtype of breast cancer. Cells cultured for 8 weeks in sub-IC50 metformin concentration proliferated comparably to untreated cells and exhibited higher rates of glucose uptake. SCLC cells were enriched in metformin-adapted cultures. These SCLC cells acquired sensitivity to multiple methods of inhibition of CtBP function, including a cyclic peptide inhibitor of NADH-induced CtBP dimerization. Single-cell mRNA sequencing identified a reprogramming of epithelial-mesenchymal and stem cell gene expression in the metformin-adapted SCLC cells. These SCLC cells demonstrated an acquired dependency on one of these genes, Tenascin C. Thus, in addition to acquisition of sensitivity to glycolysis-targeting therapeutic strategies, the reprograming of gene expression in the metformin-adapted SCLC cells renders them sensitive to potential therapeutic approaches not directly linked to cell metabolism.

The effects of restricted glycolysis on stem-cell like characteristics of breast cancer cells.

Altered glycolysis is a characteristic of many cancers, and can also be associated with changes in stem cell-like cancer (SCLC) cell populations. We therefore set out to directly examine the effect of glycolysis on SCLC cell phenotype, using a model where glycolysis is stably reduced by adapting the cells to a sugar source other than glucose. Restricting glycolysis using this approach consistently resulted in cells with increased oncogenic potential; including an increase in SCLC cells, proliferation in 3D matrigel, invasiveness, chemoresistance, and altered global gene expression. Tumorigenicity in vivo was also markedly increased. SCLC cells exhibited increased dependence upon alternate metabolic pathways. They also became c-KIT dependent, indicating that their apparent state of maturation is regulated by glycolysis. Single-cell mRNA sequencing identified altered networks of metabolic-, stem- and signaling- gene expression within SCLC-enriched populations in response to glycolytic restriction. Therefore, reduced glycolysis, which may occur in niches within tumors where glucose availability is limiting, can promote tumor aggressiveness by increasing SCLC cell populations, but can also introduce novel, potentially exploitable, vulnerabilities in SCLC cells.

Subclonal Evolution of Cancer-Related Gene Mutations in p53 Immunopositive Patches in Human Skin.

Normal sun-exposed skin contains numerous epidermal patches that stain positive for p53 protein (p53 immunopositive patches, PIPs), which are considered potential early precursors of skin cancer. Although the TP53 gene is mutated in many PIPs, it is unclear whether PIPs contain any other cancer-related mutations. Here we report that PIPs, predominantly <3,000 p53 immunopositive cells in size, within normal chronically exposed skin contain mutations in multiple genes that are mutated in cutaneous squamous cell cancers. These mutations in the PIPs were not detected within the non-PIP epidermis of corresponding normal chronically exposed skin. Although some of these genetic alterations are clonal in the PIPs, many of the mutations are subclonal within these lesions. Similar mutations are seen in later precancers (actinic keratoses and Bowen's disease). Our results demonstrate that PIPs in chronically exposed skin contain multiple mutations in cancer-related genes. In addition, the results indicate that the clonal evolution of mutations that are seen within later precancerous lesions and in established malignancy can also occur in PIPs within normal human skin.

Antibody Tumor Targeting Is Enhanced by CD27 Agonists through Myeloid Recruitment.

Monoclonal antibodies (mAbs) can destroy tumors by recruiting effectors such as myeloid cells, or targeting immunomodulatory receptors to promote cytotoxic T cell responses. Here, we examined the therapeutic potential of combining a direct tumor-targeting mAb, anti-CD20, with an extended panel of immunomodulatory mAbs. Only the anti-CD27/CD20 combination provided cures. This was apparent in multiple lymphoma models, including huCD27 transgenic mice using the anti-huCD27, varlilumab. Detailed mechanistic analysis using single-cell RNA sequencing demonstrated that anti-CD27 stimulated CD8+ T and natural killer cells to release myeloid chemo-attractants and interferon gamma, to elicit myeloid infiltration and macrophage activation. This study demonstrates the therapeutic advantage of using an immunomodulatory mAb to regulate lymphoid cells, which then recruit and activate myeloid cells for enhanced killing of mAb-opsonized tumors.

Evaluation of High-Throughput Genomic Assays for the Fc Gamma Receptor Locus.

Cancer immunotherapy has been revolutionised by the use monoclonal antibodies (mAb) that function through their interaction with Fc gamma receptors (FcγRs). The low-affinity FcγR genes are highly homologous, map to a complex locus at 1p23 and harbour single nucleotide polymorphisms (SNPs) and copy number variation (CNV) that can impact on receptor function and response to therapeutic mAbs. This complexity can hinder accurate characterisation of the locus. We therefore evaluated and optimised a suite of assays for the genomic analysis of the FcγR locus amenable to peripheral blood mononuclear cells and formalin-fixed paraffin-embedded (FFPE) material that can be employed in a high-throughput manner. Assessment of TaqMan genotyping for FCGR2A-131H/R, FCGR3A-158F/V and FCGR2B-232I/T SNPs demonstrated the need for additional methods to discriminate genotypes for the FCGR3A-158F/V and FCGR2B-232I/T SNPs due to sequence homology and CNV in the region. A multiplex ligation-dependent probe amplification assay provided high quality SNP and CNV data in PBMC cases, but there was greater data variability in FFPE material in a manner that was predicted by the BIOMED-2 multiplex PCR protocol. In conclusion, we have evaluated a suite of assays for the genomic analysis of the FcγR locus that are scalable for application in large clinical trials of mAb therapy. These assays will ultimately help establish the importance of FcγR genetics in predicting response to antibody therapeutics.

Fcγ receptors: genetic variation, function, and disease.

Fcγ receptors (FcγRs) are key immune receptors responsible for the effective control of both humoral and innate immunity and are central to maintaining the balance between generating appropriate responses to infection and preventing autoimmunity. When this balance is lost, pathology results in increased susceptibility to cancer, autoimmunity, and infection. In contrast, optimal FcγR engagement facilitates effective disease resolution and response to monoclonal antibody immunotherapy. The underlying genetics of the FcγR gene family are a central component of this careful balance. Complex in humans and generated through ancestral duplication events, here we review the evolution of the gene family in mammals, the potential importance of copy number, and functionally relevant single nucleotide polymorphisms, as well as discussing current approaches and limitations when exploring genetic variation in this region.

Exome sequence read depth methods for identifying copy number changes.

Copy number variants (CNVs) play important roles in a number of human diseases and in pharmacogenetics. Powerful methods exist for CNV detection in whole genome sequencing (WGS) data, but such data are costly to obtain. Many disease causal CNVs span or are found in genome coding regions (exons), which makes CNV detection using whole exome sequencing (WES) data attractive. If reliably validated against WGS-based CNVs, exome-derived CNVs have potential applications in a clinical setting. Several algorithms have been developed to exploit exome data for CNV detection and comparisons made to find the most suitable methods for particular data samples. The results are not consistent across studies. Here, we review some of the exome CNV detection methods based on depth of coverage profiles and examine their performance to identify problems contributing to discrepancies in published results. We also present a streamlined strategy that uses a single metric, the likelihood ratio, to compare exome methods, and we demonstrated its utility using the VarScan 2 and eXome Hidden Markov Model (XHMM) programs using paired normal and tumour exome data from chronic lymphocytic leukaemia patients. We use array-based somatic CNV (SCNV) calls as a reference standard to compute prevalence-independent statistics, such as sensitivity, specificity and likelihood ratio, for validation of the exome-derived SCNVs. We also account for factors known to influence the performance of exome read depth methods, such as CNV size and frequency, while comparing our findings with published results.

ATM mutation rather than BIRC3 deletion and/or mutation predicts reduced survival in 11q-deleted chronic lymphocytic leukemia: data from the UK LRF CLL4 trial.

ATM mutation and BIRC3 deletion and/or mutation have independently been shown to have prognostic significance in chronic lymphocytic leukemia. However, the relative clinical importance of these abnormalities in patients with a deletion of 11q encompassing the ATM gene has not been established. We screened a cohort of 166 patients enriched for 11q-deletions for ATM mutations and BIRC3 deletion and mutation and determined the overall and progression-free survival among the 133 of these cases treated within the UK LRF CLL4 trial. SNP6.0 profiling demonstrated that BIRC3 deletion occurred in 83% of 11q-deleted cases and always co-existed with ATM deletion. For the first time we have demonstrated that 40% of BIRC3-deleted cases have concomitant deletion and mutation of ATM. While BIRC3 mutations were rare, they exclusively occurred with BIRC3 deletion and a wild-type residual ATM allele. In 11q-deleted cases, we confirmed that ATM mutation was associated with a reduced overall and progression-free survival comparable to that seen with TP53 abnormalities, whereas BIRC3 deletion and/or mutation had no impact on overall and progression-free survival. In conclusion, in 11q-deleted patients treated with first-line chemotherapy, ATM mutation rather than BIRC3 deletion and/or mutation identifies a subgroup with a poorer outcome.