Mutation signatures specific to DNA alkylating agents in yeast and cancers

Direct Measurement of 8OG Syn-Anti Flips in Mutagenic 8OG·A and Long-Range Damage-Dependent Hoogsteen Breathing Dynamics Using 1H CEST NMR

Elucidating how damage impacts DNA dynamics is essential for understanding the mechanisms of damage recognition and repair. Many DNA lesions alter their propensities to form low-populated and short-lived conformational states. However, NMR methods to measure these dynamics require isotopic enrichment, which is difficult for damaged nucleotides. Here, we demonstrate the utility of the 1H chemical exchange saturation transfer (CEST) NMR experiment in measuring the dynamics of oxidatively damaged 8-oxoguanine (8OG) in the mutagenic 8OGsyn·Aanti mismatch. Using 8OG-H7 as an NMR probe of the damaged base, we directly measured 8OG syn-anti flips to form a lowly populated (pop. ∼ 5%) and short-lived (lifetime ∼50 ms) nonmutagenic 8OGanti·Aanti. These exchange parameters were in quantitative agreement with values from 13C off-resonance R1ρ and CEST on the labeled partner adenine. The Watson-Crick-like 8OGsyn·Aanti mismatch also rescued the kinetics of Hoogsteen motions at distant A-T base pairs, which the G·A mismatch had slowed down. The results lend further support for 8OGanti·Aanti as a minor conformational state of 8OG·A, reveal that 8OG damage can impact Hoogsteen dynamics at a distance, and demonstrate the utility of 1H CEST for measuring damage-dependent dynamics in unlabeled DNA.

Direct Measurement of 8OG syn-anti Flips in Mutagenic 8OG•A and Long-Range Damage-Dependent Hoogsteen Breathing Dynamics Using 1H CEST NMR

Elucidating how damage impacts DNA dynamics is essential for understanding the mechanisms of damage recognition and repair. Many DNA lesions alter the propensities to form lowly-populated and short-lived conformational states. However, NMR methods to measure these dynamics require isotopic enrichment, which is difficult for damaged nucleotides. Here, we demonstrate the utility of the 1H chemical exchange saturation transfer (CEST) NMR experiment in measuring the dynamics of oxidatively damaged 8-oxoguanine (8OG) in the mutagenic 8OGsyn•Aanti mismatch. Using 8OG-H7 as an NMR probe of the damaged base, we directly measured 8OG syn-anti flips to form a lowly-populated (pop. ~ 5%) and short-lived (lifetime ~ 50 ms) non-mutagenic 8OGanti•Aanti. These exchange parameters were in quantitative agreement with values from 13C off-resonance R1ρ and CEST on a labeled partner adenine. The Watson-Crick-like 8OGsyn•Aanti mismatch also rescued the kinetics of Hoogsteen motions at distance A-T base pairs, which the G•A mismatch had slowed down. The results lend further support for 8OGanti•Aanti as a minor conformational state of 8OG•A, reveal that 8OG damage can impact Hoogsteen dynamics at a distance, and demonstrate the utility of 1H CEST for measuring damage-dependent dynamics in unlabeled DNA.

Glycidamide-induced hypermutation in yeast single-stranded DNA reveals a ubiquitous clock-like mutational motif in humans

Mutagens often prefer specific nucleotides or oligonucleotide motifs that can be revealed by studying the hypermutation spectra in single-stranded (ss) DNA. We utilized a yeast model to explore mutagenesis by glycidamide, a simple epoxide formed endogenously in humans from the environmental toxicant acrylamide. Glycidamide caused ssDNA hypermutation in yeast predominantly in cytosines and adenines. The most frequent mutations in adenines occurred in the nAt→nGt trinucleotide motif. Base substitutions A→G in this motif relied on Rev1 translesion polymerase activity. Inactivating Rev1 did not alter the nAt trinucleotide preference, suggesting it may be an intrinsic specificity of the chemical reaction between glycidamide and adenine in the ssDNA. We found this mutational motif enriched in published sequencing data from glycidamide-treated mouse cells and ubiquitous in human cancers. In cancers, this motif was positively correlated with the single base substitution (SBS) smoking-associated SBS4 signature, with the clock-like signatures SBS1, SBS5, and was strongly correlated with smoking history and with age of tumor donors. Clock-like feature of the motif was also revealed in cells of human skin and brain. Given its pervasiveness, we propose that this mutational motif reflects mutagenic lesions to adenines in ssDNA from a potentially broad range of endogenous and exogenous agents.

Changes in metabolic landscapes shape divergent but distinct mutational signatures and cytotoxic consequences of redox stress

Mutational signatures discerned in cancer genomes, in aging tissues and in cells exposed to toxic agents, reflect complex processes underlying transformation of cells from normal to dysfunctional. Due to its ubiquitous and chronic nature, redox stress contributions to cellular makeover remain equivocal. The deciphering of a new mutational signature of an environmentally-relevant oxidizing agent, potassium bromate, in yeast single strand DNA uncovered a surprising heterogeneity in the mutational signatures of oxidizing agents. NMR-based analysis of molecular outcomes of redox stress revealed profound dissimilarities in metabolic landscapes following exposure to hydrogen peroxide versus potassium bromate. The predominance of G to T substitutions in the mutational spectra distinguished potassium bromate from hydrogen peroxide and paraquat and mirrored the observed metabolic changes. We attributed these changes to the generation of uncommon oxidizing species in a reaction with thiol-containing antioxidants; a nearly total depletion of intracellular glutathione and a paradoxical augmentation of potassium bromate mutagenicity and toxicity by antioxidants. Our study provides the framework for understanding multidimensional processes triggered by agents collectively known as oxidants. Detection of increased mutational loads associated with potassium bromate-related mutational motifs in human tumors may be clinically relevant as a biomarker of this distinct type of redox stress.

The interplay between noncoding RNAs and drug resistance in hepatocellular carcinoma: the big impact of little things

Hepatocellular carcinoma (HCC) is the leading cause of cancer-related death in people, and a common primary liver cancer. Lacking early diagnosis and a high recurrence rate after surgical resection, systemic treatment is still an important treatment method for advanced HCC. Different drugs have distinct curative effects, side effects and drug resistance due to different properties. At present, conventional molecular drugs for HCC have displayed some limitations, such as adverse drug reactions, insensitivity to some medicines, and drug resistance. Noncoding RNAs (ncRNAs), including microRNAs (miRNAs), long noncoding RNAs (lncRNAs) and circular RNAs (circRNAs), have been well documented to be involved in the occurrence and progression of cancer. Novel biomarkers and therapeutic targets, as well as research into the molecular basis of drug resistance, are urgently needed for the management of HCC. We review current research on ncRNAs and consolidate the known roles regulating drug resistance in HCC and examine the potential clinical applications of ncRNAs in overcoming drug resistance barriers in HCC based on targeted therapy, cell cycle non-specific chemotherapy and cell cycle specific chemotherapy.

Identification of Some Glutamic Acid Derivatives with Biological Potential by Computational Methods

Glutamic acid is a non-essential amino acid involved in multiple metabolic pathways. Of high importance is its relationship with glutamine, an essential fuel for cancer cell development. Compounds that can modify glutamine or glutamic acid behaviour in cancer cells have resulted in attractive anticancer therapeutic alternatives. Based on this idea, we theoretically formulated 123 glutamic acid derivatives using Biovia Draw. Suitable candidates for our research were selected among them. For this, online platforms and programs were used to describe specific properties and their behaviour in the human organism. Nine compounds proved to have suitable or easy to optimise properties. The selected compounds showed cytotoxicity against breast adenocarcinoma, lung cancer cell lines, colon carcinoma, and T cells from acute leukaemia. Compound 2Ba5 exhibited the lowest toxicity, and derivative 4Db6 exhibited the most intense bioactivity. Molecular docking studies were also performed. The binding site of the 4Db6 compound in the glutamine synthetase structure was determined, with the D subunit and cluster 1 being the most promising. In conclusion, glutamic acid is an amino acid that can be manipulated very easily. Therefore, molecules derived from its structure have great potential to become innovative drugs, and further research on these will be conducted.

Oligonucleotide synthesis under mild deprotection conditions

Over a hundred non-canonical nucleotides have been found in DNA and RNA. Many of them are sensitive toward nucleophiles. Because known oligonucleotide synthesis technologies require nucleophilic conditions for deprotection, currently there is no suitable technology for their synthesis. The recently disclosed method regarding the use of 1,3-dithian-2-yl-methyl (Dim) for phosphate protection and 1,3-dithian-2-yl-methoxycarbonyl (Dmoc) for amino protection can solve the problem. With Dim-Dmoc protection, oligodeoxynucleotide (ODN) deprotection can be achieved with NaIO4 followed by aniline. Some sensitive groups have been determined to be stable under these conditions. Besides serving as a base, aniline also serves as a nucleophilic scavenger, which prevents deprotection side products from reacting with ODN. For this reason, excess aniline is needed. Here, we report the use of alkyl Dim (aDim) and alkyl Dmoc (aDmoc) for ODN synthesis. With aDim-aDmoc protection, deprotection is achieved with NaIO4 followed by K2CO3. No nucleophilic scavenger such as aniline is needed. Over 10 ODNs including one containing the highly sensitive N4-acetylcytidine were synthesized. Work on extending the method for sensitive RNA synthesis is in progress.

Widely spaced and divergent inverted repeats become a potent source of chromosomal rearrangements in long single-stranded DNA regions

DNA inverted repeats (IRs) are widespread across many eukaryotic genomes. Their ability to form stable hairpin/cruciform secondary structures is causative in triggering chromosome instability leading to several human diseases. Distance and sequence divergence between IRs are inversely correlated with their ability to induce gross chromosomal rearrangements (GCRs) because of a lesser probability of secondary structure formation and chromosomal breakage. In this study, we demonstrate that structural parameters that normally constrain the instability of IRs are overcome when the repeats interact in single-stranded DNA (ssDNA). We established a system in budding yeast whereby >73 kb of ssDNA can be formed in cdc13-707fs mutants. We found that in ssDNA, 12 bp or 30 kb spaced Alu-IRs show similarly high levels of GCRs, while heterology only beyond 25% suppresses IR-induced instability. Mechanistically, rearrangements arise after cis-interaction of IRs leading to a DNA fold-back and the formation of a dicentric chromosome, which requires Rad52/Rad59 for IR annealing as well as Rad1-Rad10, Slx4, Msh2/Msh3 and Saw1 proteins for nonhomologous tail removal. Importantly, using structural characteristics rendering IRs permissive to DNA fold-back in yeast, we found that ssDNA regions mapped in cancer genomes contain a substantial number of potentially interacting and unstable IRs.

Overcoming Cancer Multi-drug Resistance (MDR): Reasons, mechanisms, nanotherapeutic solutions, and challenges

Multi-drug resistance (MDR) in cancer cells, either intrinsic or acquired through various mechanisms, significantly hinders the therapeutic efficacy of drugs. Typically, the reduced therapeutic performance of various drugs is predominantly due to the inherent over expression of ATP-binding cassette (ABC) transporter proteins on the cell membrane, resulting in the deprived uptake of drugs, augmenting drug detoxification, and DNA repair. In addition to various physiological abnormalities and extensive blood flow, MDR cancer phenotypes exhibit improved apoptotic threshold and drug efflux efficiency. These severe consequences have substantially directed researchers in the fabrication of various advanced therapeutic strategies, such as co-delivery of drugs along with various generations of MDR inhibitors, augmented dosage regimens and frequency of administration, as well as combinatorial treatment options, among others. In this review, we emphasize different reasons and mechanisms responsible for MDR in cancer, including but not limited to the known drug efflux mechanisms mediated by permeability glycoprotein (P-gp) and other pumps, reduced drug uptake, altered DNA repair, and drug targets, among others. Further, an emphasis on specific cancers that share pathogenesis in executing MDR and effluxed drugs in common is provided. Then, the aspects related to various nanomaterials-based supramolecular programmable designs (organic- and inorganic-based materials), as well as physical approaches (light- and ultrasound-based therapies), are discussed, highlighting the unsolved issues and future advancements. Finally, we summarize the review with interesting perspectives and future trends, exploring further opportunities to overcome MDR.

Mutagenesis induced by protonation of single-stranded DNA is linked to glycolytic sugar metabolism

Mutagenesis can be thought of as random, in the sense that the occurrence of each mutational event cannot be predicted with precision in space or time. However, when sufficiently large numbers of mutations are analyzed, recurrent patterns of base changes called mutational signatures can be identified. To date, some 60 single base substitution or SBS signatures have been derived from analysis of cancer genomics data. We recently reported that the ubiquitous signature SBS5 matches the pattern of single nucleotide polymorphisms (SNPs) in humans and has analogs in many species. Using a temperature-sensitive single-stranded DNA (ssDNA) mutation reporter system, we also showed that a similar mutational pattern in yeast is dependent on error-prone translesion DNA synthesis (TLS) and glycolytic sugar metabolism. Here, we further investigated mechanisms that are responsible for this form of mutagenesis in yeast. We first confirmed that excess sugar metabolism leads to increased mutation rate, which was detectable by fluctuation assay. Since glycolysis is known to produce excess protons, we then investigated the effects of experimental manipulations on pH and mutagenesis. We hypothesized that yeast metabolizing 8% glucose would produce more excess protons than cells metabolizing 2% glucose. Consistent with this, cells metabolizing 8% glucose had lower intracellular and extracellular pH values. Similarly, deletion of vma3 (encoding a vacuolar H⁺-ATPase subunit) increased mutagenesis. We also found that treating cells with edelfosine (which renders membranes more permeable, including to protons) or culturing in low pH media increased mutagenesis. Analysis of the mutational pattern attributable to 20 µM edelfosine treatment revealed similarity to the SBS5-like TLS- and glycolysis-dependant mutational patterns previously observed in ssDNA. Altogether, our results agree with multiple biochemical studies showing that protonation of nitrogenous bases can alter base pairing so as to stabilize some mispairs, and shed new light on a common form of intrinsic mutagenesis.

Analyses of mutational patterns induced by formaldehyde and acetaldehyde reveal similarity to a common mutational signature

Formaldehyde and acetaldehyde are reactive small molecules produced endogenously in cells as well as being environmental contaminants. Both of these small aldehydes are classified as human carcinogens, since they are known to damage DNA and exposure is linked to cancer incidence. However, the mutagenic properties of formaldehyde and acetaldehyde remain incompletely understood, at least in part because they are relatively weak mutagens. Here, we use a highly sensitive yeast genetic reporter system featuring controlled generation of long single-stranded DNA regions to show that both small aldehydes induced mutational patterns characterized by predominantly C/G → A/T, C/G → T/A, and T/A → C/G substitutions, each in similar proportions. We observed an excess of C/G → A/T transversions when compared to mock-treated controls. Many of these C/G → A/T transversions occurred at TC/GA motifs. Interestingly, the formaldehyde mutational pattern resembles single base substitution signature 40 from the Catalog of Somatic Mutations in Cancer. Single base substitution signature 40 is a mutational signature of unknown etiology. We also noted that acetaldehyde treatment caused an excess of deletion events longer than 4 bases while formaldehyde did not. This latter result could be another distinguishing feature between the mutational patterns of these simple aldehydes. These findings shed new light on the characteristics of 2 important, commonly occurring mutagens.

Acetaldehyde makes a distinct mutation signature in single-stranded DNA

Acetaldehyde (AA), a by-product of ethanol metabolism, is acutely toxic due to its ability to react with various biological molecules including DNA and proteins, which can greatly impede key processes such as replication and transcription and lead to DNA damage. As such AA is classified as a group 1 carcinogen by the International Agency for Research on Cancer (IARC). Previous in vitro studies have shown that AA generates bulky adducts on DNA, with signature guanine-centered (GG→TT) mutations. However, due to its weak mutagenicity, short chemical half-life, and the absence of powerful genetic assays, there is considerable variability in reporting the mutagenic effects of AA in vivo. Here, we used an established yeast genetic reporter system and demonstrate that AA treatment is highly mutagenic to cells and leads to strand-biased mutations on guanines (G→T) at a high frequency on single stranded DNA (ssDNA). We further demonstrate that AA-derived mutations occur through lesion bypass on ssDNA by the translesion polymerase Polζ. Finally, we describe a unique mutation signature for AA, which we then identify in several whole-genome and -exome sequenced cancers, particularly those associated with alcohol consumption. Our study proposes a key mechanism underlying carcinogenesis by acetaldehyde—mutagenesis of single-stranded DNA.

Cellular heterogeneity in DNA alkylation repair increases population genetic plasticity

DNA repair mechanisms fulfil a dual role, as they are essential for cell survival and genome maintenance. Here, we studied how cells regulate the interplay between DNA repair and mutation. We focused on the adaptive response that increases the resistance of Escherichia coli cells to DNA alkylation damage. Combination of single-molecule imaging and microfluidic-based single-cell microscopy showed that noise in the gene activation timing of the master regulator Ada is accurately propagated to generate a distinct subpopulation of cells in which all proteins of the adaptive response are essentially absent. Whereas genetic deletion of these proteins causes extreme sensitivity to alkylation stress, a temporary lack of expression is tolerated and increases genetic plasticity of the whole population. We demonstrated this by monitoring the dynamics of nascent DNA mismatches during alkylation stress as well as the frequency of fixed mutations that are generated by the distinct subpopulations of the adaptive response. We propose that stochastic modulation of DNA repair capacity by the adaptive response creates a viable hypermutable subpopulation of cells that acts as a source of genetic diversity in a clonal population.

Identification of copy number variation-driven molecular subtypes informative for prognosis and treatment in pancreatic adenocarcinoma of a Chinese cohort

BACKGROUND

Pancreatic adenocarcinoma (PAAD) is one of the most lethal carcinomas, and the current histopathological classifications are of limited use in clinical decision-making. There is an unmet need to identify new biomarkers for prognosis-informative molecular subtyping and ultimately for precision medicine.

METHODS

We profiled genomic alterations for 608 PAAD patients in a Chinese cohort, including somatic mutations, pathogenic germline variants and copy number variations (CNV). Using the CNV information, we performed unsupervised consensus clustering of these patients, differential CNV analysis and functional/pathway enrichment analysis. Cox regression was conducted for progression-free survival analysis, the elastic net algorithm used for prognostic model construction, and rank-based gene set enrichment analysis for exploring tumor microenvironments.

FINDINGS

Our data did not support prognostic value of point mutations in either highly mutated genes (such as KRAS, TP53, CDKN2A and SMAD4) or homologous recombination repair genes. Instead, associated with worse prognosis were amplified genes involved in DNA repair and receptor tyrosine kinase (RTK) related signalings. Motivated by this observation, we categorized patients into four molecular subtypes (namely repair-deficient, proliferation-active, repair-proficient and repair-enhanced) that differed in prognosis, and also constructed a prognostic model that can stratify patients with low or high risk of relapse. Finally, we analyzed publicly available datasets, not only reinforcing the prognostic value of our identified genes in DNA repair and RTK related signalings, but also identifying tumor microenvironment correlates with prognostic risks.

INTERPRETATION

Together with the evidence from genomic footprint analysis, we suggest that repair-deficient and proliferation-active subtypes are better suited for DNA damage therapies, while immunotherapy is highly recommended for repair-proficient and repair-enhanced subtypes. Our results represent a significant step in molecular subtyping, diagnosis and management for PAAD patients.

FUNDING

This work was supported by the National Natural Science Foundation of China (grant numbers 81470894, 81502695, 81672325, 81871906, 82073326, 82103482 and 32170663), the Shanghai Sailing Program (grant number 20YF1426900), and the Program for Professor of Special Appointment (Eastern Scholar) at Shanghai Institutions of Higher Learning (awarded to H.F.).

The Shu complex prevents mutagenesis and cytotoxicity of single-strand specific alkylation lesions

Three-methyl cytosine (3meC) are toxic DNA lesions, blocking base pairing. Bacteria and humans express members of the AlkB enzymes family, which directly remove 3meC. However, other organisms, including budding yeast, lack this class of enzymes. It remains an unanswered evolutionary question as to how yeast repairs 3meC, particularly in single-stranded DNA. The yeast Shu complex, a conserved homologous recombination factor, aids in preventing replication-associated mutagenesis from DNA base damaging agents such as methyl methanesulfonate (MMS). We found that MMS-treated Shu complex-deficient cells exhibit a genome-wide increase in A:T and G:C substitutions mutations. The G:C substitutions displayed transcriptional and replicational asymmetries consistent with mutations resulting from 3meC. Ectopic expression of a human AlkB homolog in Shu-deficient yeast rescues MMS-induced growth defects and increased mutagenesis. Thus, our work identifies a novel homologous recombination-based mechanism mediated by the Shu complex for coping with alkylation adducts.

Genomic and evolutionary classification of lung cancer in never smokers

Lung cancer in never smokers (LCINS) is a common cause of cancer mortality but its genomic landscape is poorly characterized. Here high-coverage whole-genome sequencing of 232 LCINS showed 3 subtypes defined by copy number aberrations. The dominant subtype (piano), which is rare in lung cancer in smokers, features somatic UBA1 mutations, germline AR variants and stem cell-like properties, including low mutational burden, high intratumor heterogeneity, long telomeres, frequent KRAS mutations and slow growth, as suggested by the occurrence of cancer drivers' progenitor cells many years before tumor diagnosis. The other subtypes are characterized by specific amplifications and EGFR mutations (mezzo-forte) and whole-genome doubling (forte). No strong tobacco smoking signatures were detected, even in cases with exposure to secondhand tobacco smoke. Genes within the receptor tyrosine kinase-Ras pathway had distinct impacts on survival; five genomic alterations independently doubled mortality. These findings create avenues for personalized treatment in LCINS.

From RNA World to SARS-CoV-2: The Edited Story of RNA Viral Evolution

The current SARS-CoV-2 pandemic underscores the importance of understanding the evolution of RNA genomes. While RNA is subject to the formation of similar lesions as DNA, the evolutionary and physiological impacts RNA lesions have on viral genomes are yet to be characterized. Lesions that may drive the evolution of RNA genomes can induce breaks that are repaired by recombination or can cause base substitution mutagenesis, also known as base editing. Over the past decade or so, base editing mutagenesis of DNA genomes has been subject to many studies, revealing that exposure of ssDNA is subject to hypermutation that is involved in the etiology of cancer. However, base editing of RNA genomes has not been studied to the same extent. Recently hypermutation of single-stranded RNA viral genomes have also been documented though its role in evolution and population dynamics. Here, we will summarize the current knowledge of key mechanisms and causes of RNA genome instability covering areas from the RNA world theory to the SARS-CoV-2 pandemic of today. We will also highlight the key questions that remain as it pertains to RNA genome instability, mutations accumulation, and experimental strategies for addressing these questions.

A selection-based next generation sequencing approach to develop robust, genotype-specific mutation profiles in Saccharomyces cerevisiae

Distinct mutation signatures arise from environmental exposures and/or from defects in metabolic pathways that promote genome stability. The presence of a particular mutation signature can therefore predict the underlying mechanism of mutagenesis. These insults to the genome often alter dNTP pools, which itself impacts replication fidelity. Therefore, the impact of altered dNTP pools should be considered when making mechanistic predictions based on mutation signatures. We developed a targeted deep-sequencing approach on the CAN1 gene in Saccharomyces cerevisiae to define information-rich mutational profiles associated with distinct rnr1 backgrounds. Mutations in the activity and selectivity sites of rnr1 lead to elevated and/or unbalanced dNTP levels, which compromises replication fidelity and increases mutation rates. The mutation spectra of rnr1Y285F and rnr1Y285A alleles were characterized previously; our analysis was consistent with this prior work but the sequencing depth achieved in our study allowed a significantly more robust and nuanced computational analysis of the variants observed, generating profiles that integrated information about mutation spectra, position effects, and sequence context. This approach revealed previously unidentified, genotype-specific mutation profiles in the presence of even modest changes in dNTP pools. Furthermore, we identified broader sequence contexts and nucleotide motifs that influenced variant profiles in different rnr1 backgrounds, which allowed specific mechanistic predictions about the impact of altered dNTP pools on replication fidelity.

Tracing Lung Cancer Risk Factors Through Mutational Signatures in Never-Smokers

Epidemiologic studies often rely on questionnaire data, exposure measurement tools, and/or biomarkers to identify risk factors and the underlying carcinogenic processes. An emerging and promising complementary approach to investigate cancer etiology is the study of somatic "mutational signatures" that endogenous and exogenous processes imprint on the cellular genome. These signatures can be identified from a complex web of somatic mutations thanks to advances in DNA sequencing technology and analytical algorithms. This approach is at the core of the Sherlock-Lung study (2018-ongoing), a retrospective case-only study of over 2,000 lung cancers in never-smokers (LCINS), using different patterns of mutations observed within LCINS tumors to trace back possible exposures or endogenous processes. Whole genome and transcriptome sequencing, genome-wide methylation, microbiome, and other analyses are integrated with data from histological and radiological imaging, lifestyle, demographic characteristics, environmental and occupational exposures, and medical records to classify LCINS into subtypes that could reveal distinct risk factors. To date, we have received samples and data from 1,370 LCINS cases from 17 study sites worldwide and whole-genome sequencing has been completed on 1,257 samples. Here, we present the Sherlock-Lung study design and analytical strategy, also illustrating some empirical challenges and the potential for this approach in future epidemiologic studies.

Genetic toxicity testing using human in vitro organotypic airway cultures: Assessing DNA damage with the CometChip and mutagenesis by Duplex Sequencing

The organotypic human air-liquid-interface (ALI) airway tissue model has been used as an in vitro cell culture system for evaluating the toxicity of inhaled substances. ALI airway cultures are highly differentiated, which has made it challenging to evaluate genetic toxicology endpoints. In the current study, we assayed DNA damage with the high-throughput CometChip assay and quantified mutagenesis with Duplex Sequencing, an error-corrected next-generation sequencing method capable of detecting a single mutation per 107 base pairs. Fully differentiated human ALI airway cultures were treated from the basolateral side with 6.25 to 100 μg/mL ethyl methanesulfonate (EMS) over a period of 28 days. CometChip assays were conducted after 3 and 28 days of treatment, and Duplex Sequencing after 28 days of treatment. Treating the airway cultures with EMS resulted in time- and concentration-dependent increases in DNA damage and a concentration-dependent increase in mutant frequency. The mutations observed in the EMS-treated cultures were predominantly C → T transitions and exhibited a unique trinucleotide signature relative to the negative control. Measurement of physiological endpoints indicated that the EMS treatments had no effect on anti-p63-positive basal cell frequency, but produced concentration-responsive increases in cytotoxicity and perturbations in cell morphology, along with concentration-responsive decreases in culture viability, goblet cell and anti-Ki67-positive proliferating cell frequency, cilia beating frequency, and mucin secretion. The results indicate that a unified 28-day study can be used to measure several important safety endpoints in physiologically relevant human in vitro ALI airway cultures, including DNA damage, mutagenicity, and tissue-specific general toxicity.

Favorable response to immunotherapy in a pancreatic neuroendocrine tumor with temozolomide-induced high tumor mutational burden

Neuroendocrine neoplasm of the pancreas is a rare tumor with limited treatment options. Among such tumors, treatment for pancreatic neuroendocrine tumor (PanNET) G3 is the most difficult. Temozolomide (TMZ) is commonly used to treat PanNET. However, TMZ may cause tumor gene alkylation, which induces drug resistance and rapid disease progression. Herein, we present a case of a female who was diagnosed with PanNET G3 and achieved a partial response to toripalimab, an anti‐programmed cell death‐ligand 1 (anti‐PD‐L1) monoclonal antibody, after multiple cycles of TMZ treatment. Genomic profiling revealed that compared with the patient's samples collected at baseline, the post‐TMZ‐treatment samples had markedly higher levels of tumor mutational burden (TMB) associated with characteristic alkylating mutational signature representing a positive correlation with favorable response to anti‐PD‐1 treatment. In addition, we observed a germline truncating mutation of MUTYH (W156*) that was considered to be pathogenic and potentially conferred to genomic instability. This case suggests that anti‐PD‐1 therapy could be a treatment option for PanNET patients with increased TMB after TMZ‐based treatment.

Similarity between mutation spectra in hypermutated genomes of rubella virus and in SARS-CoV-2 genomes accumulated during the COVID-19 pandemic

Genomes of tens of thousands of SARS-CoV2 isolates have been sequenced across the world and the total number of changes (predominantly single base substitutions) in these isolates exceeds ten thousand. We compared the mutational spectrum in the new SARS-CoV-2 mutation dataset with the previously published mutation spectrum in hypermutated genomes of rubella-another positive single stranded (ss) RNA virus. Each of the rubella virus isolates arose by accumulation of hundreds of mutations during propagation in a single subject, while SARS-CoV-2 mutation spectrum represents a collection events in multiple virus isolates from individuals across the world. We found a clear similarity between the spectra of single base substitutions in rubella and in SARS-CoV-2, with C to U as well as A to G and U to C being the most prominent in plus strand genomic RNA of each virus. Of those, U to C changes universally showed preference for loops versus stems in predicted RNA secondary structure. Similarly, to what was previously reported for rubella virus, C to U changes showed enrichment in the uCn motif, which suggested a subclass of APOBEC cytidine deaminase being a source of these substitutions. We also found enrichment of several other trinucleotide-centered mutation motifs only in SARS-CoV-2-likely indicative of a mutation process characteristic to this virus. Altogether, the results of this analysis suggest that the mutation mechanisms that lead to hypermutation of the rubella vaccine virus in a rare pathological condition may also operate in the background of the SARS-CoV-2 viruses currently propagating in the human population.

Similarity between mutation spectra in hypermutated genomes of rubella virus and in SARS-CoV-2 genomes accumulated during the COVID-19 pandemic

Genomes of tens of thousands of SARS-CoV2 isolates have been sequenced across the world and the total number of changes (predominantly single base substitutions) in these isolates exceeds ten thousand. We compared the mutational spectrum in the new SARS-CoV-2 mutation dataset with the previously published mutation spectrum in hypermutated genomes of rubella - another positive single stranded (ss) RNA virus. Each of the rubella isolates arose by accumulation of hundreds of mutations during propagation in a single subject, while SARS-CoV-2 mutation spectrum represents a collection events in multiple virus isolates from individuals across the world. We found a clear similarity between the spectra of single base substitutions in rubella and in SARS-CoV-2, with C to U as well as A to G and U to C being the most prominent in plus strand genomic RNA of each virus. Of those, U to C changes universally showed preference for loops versus stems in predicted RNA secondary structure. Similarly, to what was previously reported for rubella, C to U changes showed enrichment in the uCn motif, which suggested a subclass of APOBEC cytidine deaminase being a source of these substitutions. We also found enrichment of several other trinucleotide-centered mutation motifs only in SARS-CoV-2 - likely indicative of a mutation process characteristic to this virus. Altogether, the results of this analysis suggest that the mutation mechanisms that lead to hypermutation of the rubella vaccine virus in a rare pathological condition may also operate in the background of the SARS-CoV-2 viruses currently propagating in the human population.

Hypermutation in single-stranded DNA

Regions of genomic DNA can become single-stranded in the course of normal replication and transcription as well as during DNA repair. Abnormal repair and replication intermediates can contain large stretches of persistent single-stranded DNA, which is extremely vulnerable to DNA damaging agents and hypermutation. Since such single-stranded DNA spans only a fraction of the genome at a given instance, hypermutation in these regions leads to tightly-spaced mutation clusters. This phenomenon of hypermutation in single-stranded DNA has been documented in several experimental models as well as in cancer genomes. Recently, hypermutated single-stranded RNA viral genomes also have been documented. Moreover, indications of hypermutation in single-stranded DNA may also be found in the human germline. This review will summarize key current knowledge and the recent developments in understanding the diverse mechanisms and sources of ssDNA hypermutation.