Increase In Mutation Frequency Research Articles

Identification of a Macrophage marker gene signature to evaluate immune infiltration and therapeutic response in hepatocellular carcinoma

BackgroundOnly a minority of hepatocellular carcinoma (HCC) patients can benefit from systemic regimens. Macrophages, which abundantly infiltrate in HCC, could mediate tumour microenvironment remodelling and immune escape, proving to be powerful weapons in combating HCC. Thus, a deeper understanding of macrophages is necessary for improving existing antitumour treatments. MethodsWith a series of bioinformatic approaches, we comprehensively explored the role of macrophage-related genes in human HCCs from multiple single-cell and bulk RNA sequencing datasets. Unsupervised clustering was performed to cluster the macrophage marker genes (MMGs). GSVA and functional enrichment analysis were used to elucidate the functional differences among the MMG-associated clusters. Subsequently, a component analysis algorithm was used to construct a Macrosig score, and the prognosis, biological characteristics, mutation profile, TME cell infiltration status and drug response of patients with different Macrosig scores were further analysed. ResultsWe identified 13 MMGs in 574 HCC samples, based on which three MMG-associated clusters were defined. Overall survival time, clinicopathological features and immune infiltration scores differed among the different clusters. On this basis, 12 hub genes were identified among these clusters; subsequently, a scoring system was constructed to determine the Macrosig score. Importantly, patients with low-Macrosig scores, characterized by increased immune infiltration, increased mutation frequency and increased immune checkpoint expression, including CTLA-4, LAG3, PDCD1 and TIGIT, exhibited enhanced efficacy of immunotherapy when validated in an external database. Moreover, a low-Macrosig score indicates increased sensitivity to AZD.2281, A.443654, ABT.263, ABT.888, AG.014699 and ATRA, while a high Macrosig score indicates increased sensitivity to AZD6482, AKT inhibitor VIII, AS601245, AZ628, AZD.0530 and AZD6244. ConclusionsA novel scoring system was constructed to guide more effective prognostic evaluation and tailoring therapeutic regimens for HCC patients.

Phage resistance in Klebsiella pneumoniae and bidirectional effects impacting antibiotic susceptibility

ObjectivesBacteriophage (phage) therapy is a promising anti-infective option to combat antimicrobial resistance. However, the clinical utilization of phage therapy has been severely compromised by the potential emergence of phage resistance. Although certain phage resistance mechanisms can restore bacterial susceptibility to certain antibiotics, a lack of knowledge of phage resistance mechanisms hinders optimal use of phages and their combination with antibiotics. MethodsGenome-wide transposon screening was performed with a mutant library of Klebsiella pneumoniae MKP103 to identify phage pKMKP103_1-resistant mutants. Phage-resistant phenotypes were evaluated by time-kill kinetics and efficiency of plating assays. Phage resistance mechanisms were investigated with adsorption, one-step growth, and mutation frequency assays. Antibiotic susceptibility was determined with broth microdilution and population analysis profiles. ResultsWe observed a repertoire of phage resistance mechanisms in K pneumoniae, such as disruption of phage binding (fhuA::Tn and tonB::Tn), extension of the phage latent period (mnmE::Tn and rpoN::Tn), and increased mutation frequency (mutS::Tn and mutL::Tn). Notably, in contrast to the prevailing view that phage resistance re-sensitizes antibiotic-resistant bacteria, we observed a bidirectional steering effect on bacterial antibiotic susceptibility. Specifically, rpoN::Tn increased susceptibility to colistin while mutS::Tn and mutL::Tn increased resistance to rifampicin and colistin. DiscussionOur findings demonstrate that K pneumoniae employs multiple strategies to overcome phage infection, which may result in enhanced or reduced antibiotic susceptibility. Mechanism-guided phage steering should be incorporated into phage therapy to better inform clinical decisions on phage-antibiotic combinations.

Possible contribution of 8-hydroxydeoxyguanosine to gene mutations in the kidney DNA of gpt delta rats following potassium bromate treatment

8-Hydroxydeoxyguanosine (8-OHdG) is well known not only as an effective biomarker of oxidative stress but also as a mutagenic DNA modification. Incorporation of dAMP at the opposite site of 8-OHdG induces G>T or A>C transversions. However, in vivo analyses of gene mutations caused by potassium bromate (KBrO3), which can induce 8-OHdG at carcinogenic target sites, showed that G>T was prominent in the small intestines of mice, but not in the kidneys of rats. Because KBrO3 was a much clearer carcinogen in the kidneys of rats, detailed analyses of gene mutations in the kidney DNA of rats treated with KBrO3 could improve our understanding of oxidative stress-mediated carcinogenesis. In the current study, site-specific reporter gene mutation assays were performed in the kidneys of gpt delta rats treated with KBrO3. Groups of 5 gpt delta rats were treated with KBrO3 at concentrations of 0, 125, 250, or 500 ppm in the drinking water for 9 weeks. At necropsy, the kidneys were macroscopically divided into the cortex and medulla. 8-OHdG levels in DNA extracted from the cortex were dramatically elevated at concentrations of 250 ppm and higher compared with those from the medulla. Cortex-specific increases in mutant frequencies in gpt and red/gam genes were found at 500 ppm. Mutation spectrum and sequence analyses of their mutants demonstrated significant elevations in A>T transversions in the gpt gene and single base deletions at guanine or adenine in the gpt or red/gam genes. While A>T transversions and single base deletions of adenine may result from the oxidized modification of adenine, the contribution of 8-OHdG to gene mutations was limited despite possible participation of the 8-OHdG repair process in guanine deletion.

Using Selective Enzymes to Measure Noncanonical DNA Building Blocks: dUTP, 5-Methyl-dCTP, and 5-Hydroxymethyl-dCTP.

Cells maintain a fine-tuned balance of deoxyribonucleoside 5'-triphosphates (dNTPs), a crucial factor in preserving genomic integrity. Any alterations in the nucleotide pool's composition or chemical modifications to nucleotides before their incorporation into DNA can lead to increased mutation frequency and DNA damage. In addition to the chemical modification of canonical dNTPs, the cellular de novo dNTP metabolism pathways also produce noncanonical dNTPs. To keep their levels low and prevent them from incorporating into the DNA, these noncanonical dNTPs are removed from the dNTP pool by sanitizing enzymes. In this study, we introduce innovative protocols for the high-throughput fluorescence-based quantification of dUTP, 5-methyl-dCTP, and 5-hydroxymethyl-dCTP. To distinguish between noncanonical dNTPs and their canonical counterparts, specific enzymes capable of hydrolyzing either the canonical or noncanonical dNTP analogs are employed. This approach provides a more precise understanding of the composition and noncanonical constituents of dNTP pools, facilitating a deeper comprehension of DNA metabolism and repair. It is also crucial for accurately interpreting mutational patterns generated through the next-generation sequencing of biological samples.

SETD2 mutations do not contribute to clonal fitness in response to chemotherapy in childhood B cell acute lymphoblastic leukemia

Mutations in genes encoding epigenetic regulators are commonly observed at relapse in B cell acute lymphoblastic leukemia (B-ALL). Loss-of-function mutations in SETD2, an H3K36 methyltransferase, have been observed in B-ALL and other cancers. Previous studies on mutated SETD2 in solid tumors and acute myelogenous leukemia support a role in promoting resistance to DNA damaging agents. We did not observe chemoresistance, an impaired DNA damage response, nor increased mutation frequency in response to thiopurines using CRISPR-mediated knockout in wild-type B-ALL cell lines. Likewise, restoration of SETD2 in cell lines with hemizygous mutations did not increase sensitivity. SETD2 mutations affected the chromatin landscape and transcriptional output that was unique to each cell line. Collectively our data does not support a role for SETD2 mutations in driving clonal evolution and relapse in B-ALL, which is consistent with the lack of enrichment of SETD2 mutations at relapse in most studies.

Assessment of the in vivo genotoxic potential of three smoke flavoring primary product mixtures.

Smoke flavorings are mixtures generated from wood pyrolysis that are filtered to remove tar and are often considered healthier alternatives to conventional smoking processes. While the latter is mostly unregulated, smoke-flavoring primary products (SFPPs) are undergoing the 10-year required re-evaluation in the European Union (EU). To comply with recent smoke flavor guidance, in vivo micronucleus studies in rats and transgenic rodent (TGR) mutation assays in Muta™Mice were conducted on three SFPPs. For most studies, typical limit doses were exceeded to comply with regulatory requests. Exposure to SFPPs by oral gavage did not result in significant increases in bone marrow micronucleus formation. Except for one group, exposure to SFPPs via feed for 28 days did not result in significant increases in mutant frequency (MF) in the glandular stomach or liver. One group exposed to a maximal feasible dietary dose of 50,000 ppm (>10,000 mg/kg bodyweight per day) exhibited a statistically significant increase in liver MF; however, the MF in all mice in this group were within the historical vehicle control 95% quantile confidence intervals and therefore not considered biologically relevant. Based on estimates of human dietary exposure to each SFPP, the margin of exposure (MOE) values in the TGR assays exceed 10,000. The MOE for one unintentionally present constituent, 2,5(H)-furanone, also exceeds 10,000. Collectively, these data indicate that these SFPPs pose no genotoxic risk and are safe alternatives to conventional smoking.

MiRNA-767 and its binding site polymorphism in the mTOR gene act as potential biomarkers for female reproductive cancers.

Aims: The present study aimed to understand the relationship between the mTOR gene SNP (rs2536) and reproductive cancer risk. The expression level of miRNA-767 was also assessed. Methods: 700 tumor samples (300 breast, 200 ovarian and 200 cervical cancers), along with adjacent uninvolved control tissue, were used. rs2536 was screened using Tetra-ARMS PCR and expression level of miRNA-767 was assessed using quantitative PCR. Results: The frequency of the homozygous mutant genotype of rs2536 was observed significantly higher in breast (p<0.04), ovarian (p<0.005) and cervical (p<0.003) cancers. Significant downregulation of miRNA-767 was observed in tumors compared with controls. Conclusion: The present study demonstrates that increased mutant frequency of rs2536 and deregulation of miRNA-767 are associated with increased reproductive cancer risk.

Advances in using adaptive laboratory evolution technology for engineering of photosynthetic cyanobacteria

Cyanobacteria are the only prokaryotes capable of oxygenic photosynthesis, which have potential to serve as "autotrophic cell factories". However, the synthesis of biofuels and chemicals using cyanobacteria as chassis are suffered from poor stress tolerance and low yield, resulting in low economic feasibility for industrial production. Thus, it's urgent to construct new cyanobacterial chassis by means of synthetic biology. In recent years, adaptive laboratory evolution (ALE) has made great achievements in chassis engineering, including optimizing growth rate, increasing tolerance, enhancing substrate utilization and increasing product yield. ALE has also made some progress in improving the tolerance of cyanobacteria to high light intensity, heavy metal ions, high concentrations of salt and organic solvents. However, the engineering efficiency of ALE strategy in cyanobacteria is generally low, and the molecular mechanisms underpinning the tolerance to various stresses have not been fully elucidated. To this end, this review summarizes the ALE-associated technical strategies and their applications in cyanobacteria chassis engineering, following by discussing how to construct larger ALE mutation library, increase mutation frequency of strains and shorten evolution time. Moreover, exploration of the construction principles and strategies for constructing multi-stress tolerant cyanobacteria, and efficient analysis the mutant libraries of evolved strains as well as construction of strains with high yield and strong robustness are discussed, with the aim to facilitate the engineering of cyanobacteria chassis and the application of engineered cyanobacteria in the future.

Comprehensive characterization of pyroptosis phenotypes with distinct tumor immune profiles in gastric cancer to aid immunotherapy.

Pyroptosis is a form of programmed cell death that is essential for immunity. Herein, this study was conducted to uncover the implication of pyroptosis in immunomodulation and tumor microenvironment (TME) in gastric cancer. Prognostic pyroptosis-related genes were extracted to identify different pyroptosis phenotypes and pyroptosis genomic phenotypes via unsupervised clustering analysis in the gastric cancer meta-cohort cohort (GSE15459, GSE62254, GSE84437, GSE26253 and TCGA-STAD). The activation of hallmark gene sets was quantified by GSVA and immune cell infiltration was estimated via ssGSEA and CIBERSORT. Through PCA algorithm, pyroptosis score was conducted. The predictors of immune response (TMB and IPS) and genetic mutations were evaluated. The efficacy of pyroptosis score in predicting immune response was verified in two anti-PD-1 therapy cohorts. Three different pyroptosis phenotypes with different prognosis, biological pathways and tumor immune microenvironment were established among 1275 gastric cancer patients, corresponding to three immune phenotypes: immune-inflamed, immune-desert, and immune-excluded. According to the pyroptosis score, patients were separated into high and low pyroptosis score groups. Low pyroptosis score indicated favorable survival outcomes, enhanced immune responses, and increased mutation frequency. Moreover, low pyroptosis score patients displayed more clinical benefits from anti-PD-1 and prolonged survival time. Our findings uncovered a nonnegligible role of pyroptosis in immunomodulation and TME multiformity and complicacy in gastric cancer. Quantifying the pyroptosis score in individual tumors may tailor more effective immunotherapeutic strategies.

Unbiased whole genome detection of ultrarare off-target mutations in genome-edited cell populations by HiFi sequencing.

DNA base editors (BEs) composed of a nuclease-deficient Cas9 fused to a DNA-modifying enzyme can achieve on-target mutagenesis without creating double-strand DNA breaks (DSBs). As a result, BEs generate far less DNA damage than traditional nuclease-proficient Cas9 systems, which do rely on the creation of DSBs to achieve on-target mutagenesis. The inability of BEs to create DSBs makes the detection of their undesired off-target effects very difficult. PacBio HiFi sequencing can efficiently detect ultrarare mutations resulting from chemical mutagenesis in whole genomes with a sensitivity ~1 × 10-8 mutations per base pair. In this proof-of-principle study, we evaluated whether this technique could also detect the on- and off-target mutations generated by a cytosine-to-thymine (C>T) BE targeting the LacZ gene in Escherichia coli (E. coli). HiFi sequencing detected on-target mutant allele fractions ranging from ~7% to ~63%, depending on the single-guide RNA (sgRNA) used, while no on-target mutations were detected in controls lacking the BE. The presence of the BE resulted in a ~3-fold increase in mutation frequencies compared to controls lacking the BE, irrespective of the sgRNA used. These increases were mostly composed of C:G>T:A substitutions distributed throughout the genome. Our results demonstrate that HiFi sequencing can efficiently identify on- and off-target mutations in cell populations that have undergone genome editing.

Duplex sequencing provides detailed characterization of mutation frequencies and spectra in the bone marrow of MutaMouse males exposed to procarbazine hydrochloride

Mutagenicity testing is an essential component of health safety assessment. Duplex Sequencing (DS), an emerging high-accuracy DNA sequencing technology, may provide substantial advantages over conventional mutagenicity assays. DS could be used to eliminate reliance on standalone reporter assays and provide mechanistic information alongside mutation frequency (MF) data. However, the performance of DS must be thoroughly assessed before it can be routinely implemented for standard testing. We used DS to study spontaneous and procarbazine (PRC)-induced mutations in the bone marrow (BM) of MutaMouse males across a panel of 20 diverse genomic targets. Mice were exposed to 0, 6.25, 12.5, or 25 mg/kg-bw/day for 28 days by oral gavage and BM sampled 42 days post-exposure. Results were compared with those obtained using the conventional lacZ viral plaque assay on the same samples. DS detected significant increases in mutation frequencies and changes to mutation spectra at all PRC doses. Low intra-group variability within DS samples allowed for detection of increases at lower doses than the lacZ assay. While the lacZ assay initially yielded a higher fold-change in mutant frequency than DS, inclusion of clonal mutations in DS mutation frequencies reduced this discrepancy. Power analyses suggested that three animals per dose group and 500 million duplex base pairs per sample is sufficient to detect a 1.5-fold increase in mutations with > 80% power. Overall, we demonstrate several advantages of DS over classical mutagenicity assays and provide data to support efforts to identify optimal study designs for the application of DS as a regulatory test.

High frequency of Voltage-gated sodium channel (VGSC) gene mutations in Aedes albopictus (Diptera: Culicidae) suggest rapid insecticide resistance evolution in Shanghai, China.

Dengue fever is an infectious disease that is imported into Shanghai, China and requires prevention and control measures. Controlling the vector Aedes albopictus through insecticide use is a key approach to dengue control. However, the rapid evolution of insecticide resistance in Ae. albopictus has raised concerns about the failure of dengue control efforts. Knockdown resistance (kdr) caused by point mutations in the voltage-gated sodium channel (VGSC) gene is a primary mechanism of pyrethroid resistance. In this study, we investigated the kdr mutations of Ae. albopictus in Shanghai and evaluated the trend in its evolution. We collected 17 populations of Ae. albopictus from 15 districts in Shanghai in 2020, extracted genomic DNA from individual mosquitoes, and amplified Domain II, III, and IV in VGSC using PCR. Following sequencing, we obtained 658 VGSC sequences. We detected the nonsynonymous mutations V1016G, I1532T, and F1534S/C/I, among which V1016G and F1534C/I were reported in Shanghai for the first time and F1534I was a novel mutant allele in Ae. albopictus. The overall mutation frequency was 84.65%, with individual mutation frequencies ranging from 46.81% to 100%, excluding the Fengxian District population, which had a frequency of 0%. The V1016G and I1532T mutation types accounted for 7.14% and 3.42%, respectively. The mutant allele at codon 1534 accounted for 63.98% of all mutations, including TCC/S (62.77%), TGC/C (1.06%), and ATC/I (0.15%). We identified and classified five intron types in Domain III by length, including A (83 bp, 12.07%), B (68 bp, 87.30%), C (80 bp, 0.16%), D (72 bp, 0.16%), and E (70 bp, 0.31%). Individuals with intron B had a significant mutation tendency at codon 1534 relative to intron A (chi-square test, p < 0.0001). We found no correlation between mutation frequency and the amount of pyrethroid used (Pearson correlation, p = 0.4755). In recent years, kdr mutations in the Ae. albopictus population in Shanghai have rapidly evolved, as evidenced by an increase in mutation types and significantly increased mutation frequency. The F1534I/ATC mutant allele was found to be a novel mutation, F1534C/TGC was reported for the first time in Shanghai, and intron B in Domain III was significantly associated with mutation frequency at codon 1534. Continuous monitoring of resistance changes and strict regulation of insecticide use are required.

Characterization of POLE c.1373A > T p.(Tyr458Phe), causing high cancer risk

The cancer syndrome polymerase proofreading-associated polyposis results from germline mutations in the POLE and POLD1 genes. Mutations in the exonuclease domain of these genes are associated with hyper- and ultra-mutated tumors with a predominance of base substitutions resulting from faulty proofreading during DNA replication. When a new variant is identified by gene testing of POLE and POLD1, it is important to verify whether the variant is associated with PPAP or not, to guide genetic counseling of mutation carriers. In 2015, we reported the likely pathogenic (class 4) germline POLE c.1373A > T p.(Tyr458Phe) variant and we have now characterized this variant to verify that it is a class 5 pathogenic variant. For this purpose, we investigated (1) mutator phenotype in tumors from two carriers, (2) mutation frequency in cell-based mutagenesis assays, and (3) structural consequences based on protein modeling. Whole-exome sequencing of two tumors identified an ultra-mutator phenotype with a predominance of base substitutions, the majority of which are C > T. A SupF mutagenesis assay revealed increased mutation frequency in cells overexpressing the variant of interest as well as in isogenic cells encoding the variant. Moreover, exonuclease repair yeast-based assay supported defect in proofreading activity. Lastly, we present a homology model of human POLE to demonstrate structural consequences leading to pathogenic impact of the p.(Tyr458Phe) mutation. The three lines of evidence, taken together with updated co-segregation and previously published data, allow the germline variant POLE c.1373A > T p.(Tyr458Phe) to be reclassified as a class 5 variant. That means the variant is associated with PPAP.

Abstract C016: Genomic analysis of patients across continental ancestry groups

Abstract People with different ancestries inherit different risks and encounter different environmental exposures resulting in different somatic profiles. A lack of knowledge about ancestry-specific alterations is a major barrier to implementing precision medicine, leading to inequities in genetic testing, targeted treatment and clinical trial design for cancer patients from the underserved populations. However, a limited number of cancer cases from non-European populations have been sequenced in research setting, and paired normal samples are often not collected in routine clinical care. Here, we inferred African (AFR), South Asian (SAS), East Asian (EAS), American (AMR) and European (EUR) ancestry from comprehensive genomic profiling (CGP) of over 200,000 tumors at Foundation Medicine and identified 18717, 8588, 6594 patients with significant proportions of AFR, AMR, EAS ancestry, respectively. Using logistic regressions with age of diagnosis, gender and tumor mutation burden (TMB) as covariates, we identified 165 ancestry-associated genes across 14 common cancer types. We found that in colorectal cancer, KRAS (OR=1.5, FDR corrected p=4x10-44) and APC (OR=1.4, FDR corrected p=1x10-16) mutations were enriched in patients with AFR ancestry, whereas BRAF mutations were depleted in both AFR (OR=0.5, FDR corrected p=5x10-24) and EAS (OR=0.6, FDR corrected p=6x10-5) patients. CDK12 mutations were enriched in both AFR (OR=1.6, FDR corrected p=2x10-4) and EAS (OR=3.1, FDR corrected p=2x10-7) prostate cancer patients. For patients with AMR ancestry, increased mutation frequency was observed in SETD2 in pancreatic cancer (OR=3.4, FDR corrected p=1x10-4), and VHL in renal clear cell carcinoma (OR=1.9, FDR corrected p=2x10-4). For renal cell carcinoma patients with AFR ancestry, mutations in NF2 were enriched (OR=2.5, FDR =8x10-6). In GBM, TERT alterations were associated with EUR ancestry (OR=1.3, FDR corrected p=0.003). Interestingly, TMB-high status was associated with EUR ancestry in melanoma (OR=5.9, p=6x10-36) and anti-correlated with AFR ancestry in endometrial carcinoma (OR=0.7, p=1x10-6). Our large-scale ancestry analysis characterized the somatic landscape of non-European cancer patients. Future directions include quantifying the germline heritability of identified ancestry-associated mutations using a local ancestry risk score (LRS). Citation Format: Jian Carrot-Zhang, Justin New, Amy Xie, Garrett Frampton. Genomic analysis of patients across continental ancestry groups [abstract]. In: Proceedings of the 15th AACR Conference on the Science of Cancer Health Disparities in Racial/Ethnic Minorities and the Medically Underserved; 2022 Sep 16-19; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2022;31(1 Suppl):Abstract nr C016.

The role of migration in mutant dynamics in fragmented populations.

Mutant dynamics in fragmented populations have been studied extensively in evolutionary biology. Yet, open questions remain, both experimentally and theoretically. Some of the fundamental properties predicted by models still need to be addressed experimentally. We contribute to this by using a combination of experiments and theory to investigate the role of migration in mutant distribution. In the case of neutral mutants, while the mean frequency of mutants is not influenced by migration, the probability distribution is. To address this empirically, we performed in vitro experiments, where mixtures of GFP-labelled ("mutant") and non-labelled ("wid-type") murine cells were grown in wells (demes), and migration was mimicked via cell transfer from well to well. In the presence of migration, we observed a change in the skewedness of the distribution of the mutant frequencies in the wells, consistent with previous and our own model predictions. In the presence of de novo mutant production, we used modelling to investigate the level at which disadvantageous mutants are predicted to exist, which has implications for the adaptive potential of the population in case of an environmental change. In panmictic populations, disadvantageous mutants can persist around a steady state, determined by the rate of mutant production and the selective disadvantage (selection-mutation balance). In a fragmented system that consists of demes connected by migration, a steady-state persistence of disadvantageous mutants is also observed, which, however, is fundamentally different from the mutation-selection balance and characterized by higher mutant levels. The increase in mutant frequencies above the selection-mutation balance can be maintained in small ( ) demes as long as the migration rate is sufficiently small. The migration rate above which the mutants approach the selection-mutation balance decays exponentially with . The observed increase in the mutant numbers is not explained by the change in the effective population size. Implications for evolutionary processes in diseases are discussed, where the pre-existence of disadvantageous drug-resistant mutant cells or pathogens drives the response of the disease to treatments.

In vitro and in vivo assessments of the genotoxic potential of 3-chloroallyl alcohol.

3-Chloroallyl alcohol (3-CAA) can be found in the environment following the application of plant protection products. 3-CAA is formed in groundwater following the injection of 1,3-dichloropropene, a fumigant used to control nematodes. 3-CAA is also formed, in leafy crops, as a glycoside conjugate following application of the herbicide, clethodim. Human exposure may occur from groundwater used as drinking water or through dietary consumption. To characterize 3-CAA's potential to cause genotoxicity in mammals, in vitro and in vivo studies were conducted. 3-CAA was negative in an Ames test and positive in a mouse lymphoma forward mutation assay. 3-CAA was negative in an acute in vivo CD-1 mouse bone marrow micronucleus assay when administered up to a dose level of 125 mg/kg/day for two consecutive days. In a combined gene mutation assay and erythrocyte micronucleus assay, using transgenic Big Blue® Fischer 344 rats, 3-CAA was administered via drinking water at targeted dose levels of 0, 10, 30, and 100 mg/kg/day for 29 days. Peripheral blood samples, collected at the end of treatment, were analyzed for micronucleus induction in reticulocytes using flow cytometry. Liver and bone marrow samples, collected 2 days after the termination of the treatment, were analyzed for the induction of mutations at the cII locus. 3-CAA did not induce an increase in mutant frequency or micronuclei under the experimental conditions. In conclusion, the mutagenic response observed in the in vitro mouse lymphoma assay is not confirmed in the whole animal. 3-CAA is not considered to pose a mutagenic risk.

Multiple Myeloma Evolution Is Characterized By Dynamic Epigenetic Landscapes

The biological mechanisms driving MM progression and multi-drug resistance remain unclear. To better characterize the molecular mechanisms driving malignant transformation and disease progression, we generated a new database comprising 1,016 bone marrow biopsies from patients across the whole MM spectrum. These samples were clinically, genetically (844 RNA-Seq, 870 WES, 7 scRNA-Seq), and epigenetically (11 scATAC-Seq, 8 CUT&Tag) characterized. While we have observed increased mutation frequency in MM driver genes in later disease stages, no specific mutations statistically correlated with malignant transformation, nor with the development of refractory disease. Enrichment analysis of transcriptional data identified biological pathways and underlying molecular events associated with MM progression and therapy resistance: cell adhesion, inflammatory cytokines, TME-dependence and hematopoietic cell identity were under-expressed in NDMM vs. premalignant disease; whereas cell cycle, energy metabolism, DNA repair, and protein/RNA synthesis/degradation were overexpressed in LRMM compared to NDMM. Importantly, this transcriptional dysregulation was anecdotally confirmed in patients with multiple sequential samples. To identify the main biological pathways responsible for the inter-sample transcriptomic variability of our cohort, we performed a principal component analysis (PCA) using single sample gene set enrichment analysis (ssGSEA) NES of Cancer Hallmarks. The first component (PC1) was associated with aforementioned hallmarks differentially expressed in NDMM and LRMM samples. This suggested PC1 as a possible prognostic biomarker, which was confirmed by Cox Proportional Hazard models, associated with shorter TTP in SMOL, PFS in NDMM, and OS in active MM. To elucidate the mechanisms driving such aberrant transcriptional programs, we applied a robust dimensionality reduction technique (t-SNE) and a clustering method (fuzzy c-means) to RNA-Seq data to identify clusters of co-expressing genes, thus characterizing MM transcriptomic topology. In total, 16,738 genes were grouped in 500 gene clusters putatively controlled by shared transcriptional programs. Enrichment analysis showed that gene clusters under-expressed during malignant transformation were enriched for H3K27me3, whereas gene clusters overexpressed in refractory disease were enriched for H3K27ac. We validated these findings using scATAC-Seq in 11 MM samples (1 MGUS, 3 SMOL, 4 NDMM, 1 ERMM and 2 LRMM), which confirmed that a greater decrease in chromatin accessibility occurred among the most under-expressed gene clusters during malignant transformation. scATAC-Seq data also revealed the differential accessibility of B cell master regulators' binding motifs across disease progression (e.g., SPI1 in MGUS, FOS in SMOL, the IRF and POU family in NDMM, RUNX2 and RFX5 in LRMM) and suggested MM pioneer TFs by correlating the accessibility of TF's genes to their targets' expression (e.g., SPI1, BACH2, RELA, IRF4, among others). CUT&Tag data from 4 NDMM and 4 LRMM samples confirmed that the strongest H3K27ac signal was found in the topological region overexpressed in refractory disease, with a positive correlation between H3K27ac signal and gene expression. Super enhancer (SE) analysis demonstrated that, while genes associated with active SEs showed higher level of expression than genes subjected to regular enhancers (RE) in all samples, LRMM presented increased expression of genes under both SE/RE control compared to their counterparts in NDMM. These observations, in conjunction with an increased expression of transcription regulation-associated DNA-binding proteins (DBPs; e.g., CDK9, CTCF, MED1, YY1), suggests that the aberrant gene expression in LRMM may be due to differential expression of these DBPs. In conclusion, MM progression, characterized by transcriptional dysregulation in response to TME and therapeutic stress, is epigenetically driven by DBPs whose activity is altered by genetic, cytogenetic, and epigenetic events. Future studies will determine the role of such masterregulators in MM biology and drug resistance, as well as allow us to propose patient-specific strategies for therapy re-sensitization.

Inhibition of ABL1 by tyrosine kinase inhibitors leads to a downregulation of MLH1 by Hsp70-mediated lysosomal protein degradation.

The DNA mismatch repair (MMR) pathway and its regulation are critical for genomic stability. Mismatch repair (MMR) follows replication and repairs misincorporated bases and small insertions or deletions that are not recognized and removed by the proofreading polymerase. Cells deficient in MMR exhibit an increased overall mutation rate and increased expansion and contraction of short repeat sequences in the genome termed microsatellite instability (MSI). MSI is often a clinical measure of genome stability in tumors and is used to determine the course of treatment. MMR is also critical for inducing apoptosis after alkylation damage from environmental agents or DNA-damaging chemotherapy. MLH1 is essential for MMR, and loss or mutation of MLH1 leads to defective MMR, increased mutation frequency, and MSI. In this study, we report that tyrosine kinase inhibitors, imatinib and nilotinib, lead to decreased MLH1 protein expression but not decreased MLH1 mRNA levels. Of the seven cellular targets of Imatinib and nilotinib, we show that silencing of ABL1 also reduces MLH1 protein expression. Treatment with tyrosine kinase inhibitors or silencing of ABL1 results in decreased apoptosis after treatment with alkylating agents, suggesting the level of MLH1 reduction is sufficient to disrupt MMR function. We also report MLH1 is tyrosine phosphorylated by ABL1. We demonstrate that MLH1 downregulation by ABL1 knockdown or inhibition requires chaperone protein Hsp70 and that MLH1 degradation can be abolished with the lysosomal inhibitor bafilomycin. Taken together, we propose that ABL1 prevents MLH1 from being targeted for degradation by the chaperone Hsp70 and that in the absence of ABL1 activity at least a portion of MLH1 is degraded through the lysosome. This study represents an advance in understanding MMR pathway regulation and has important clinical implications as MMR status is used in the clinic to inform patient treatment, including the use of immunotherapy.

Mutational profiles of head and neck squamous cell carcinomas based upon human papillomavirus status in the Veterans Affairs National Precision Oncology Program.

Patients with advanced head and neck squamous cell carcinoma (HNSCC) associated with human papillomavirus (HPV) demonstrate favorable clinical outcomes compared to patients bearing HPV-negative HNSCC. We sought to characterize the association between HPV status and mutational profiles among patients served by the Veterans Health Administration (VHA). We performed a retrospective analysis of all Veterans with primary HNSCC tumors who underwent next-generation sequencing (NGS) through the VHA's National Precision Oncology Program between July 2016 and February 2019. HPV status was determined by clinical pathology reports of p16 immunohistochemical staining; gene variant pathogenicity was classified using OncoKB, an online precision oncology knowledge database, and mutation frequencies were compared using Fisher's exact test. A total of 79 patients met inclusion criteria, of which 48 (60.8%) had p16-positive tumors. Patients with p16-negative HNSCC were more likely to have mutations in TP53 (p < 0.0001), and a trend towards increased mutation frequency was observed within NOTCH1 (p = 0.032) and within the composite CDK/Rb pathway (p = 0.065). Mutations in KRAS, NRAS, HRAS, and FBXW7 were exclusively identified within p16-positive tumors, and a trend towards increased frequency was observed within the PI3K pathway (p = 0.051). No difference in overall mutational burden was observed between the two groups. In accordance with the previous studies, no clear molecular basis for improved prognosis among patients harboring HPV-positive disease has been elucidated. Though no targeted therapies are approved based upon HPV-status, current efforts to trial PI3K inhibitors and mTOR inhibitors across patients with HPV-positive disease bear genomic rationale based upon the current findings.

Loss of ALK hotspot mutations in relapsed neuroblastoma.

ALK is the most commonly mutated oncogene in neuroblastoma with increased mutation frequency reported at relapse. Here we report the loss of an ALK mutation in two patients at relapse and a paired neuroblastoma cell line at relapse. ALK detection methods including Sanger sequencing, targeted next-generation sequencing and a new ALK Agena MassARRAY technique were used to detect common hotspot ALK variants in tumors at diagnosis and relapse from two high-risk neuroblastoma patients. Copy number analysis including single nucleotide polymorphism array and array comparative genomic hybridization confirmed adequate tumor cell content in DNA used for mutation testing. Case 1 presented with an ALK F1174L mutation at diagnosis with a variant allele frequency (VAF) ranging between 23.5% and 28.5%, but the mutation was undetectable at relapse. Case 2 presented with an ALK R1257Q mutation at diagnosis (VAF=39%-47.4%) which decreased to <0.01% at relapse. Segmental chromosomal aberrations were maintained between diagnosis and relapse confirming sufficient tumor cell content for mutation detection. The diagnostic SKNBE1n cell line harbors an ALK F1174S mutation, which was lost in the relapsed SKNBE2c cell line. To our knowledge, these are the first reported cases of loss of ALK mutations at relapse in neuroblastoma in the absence of ALK inhibitor therapy, reflecting intra-tumoral spatial and temporal heterogeneity. As ALK inhibitors are increasingly used in the treatment of refractory/relapsed neuroblastoma, our study highlights the importance of confirming whether an ALK mutation detected at diagnosis is still present in clones leading to relapse.