Frequencies and spectra of aflatoxin B1-induced mutations in liver genomes of NEIL1-deficient mice as revealed by duplex sequencing.

Increased risk for the development of hepatocellular carcinoma (HCC) is driven by a number of etiological factors including hepatitis viral infection and dietary exposures to foods contaminated with aflatoxin-producing molds. Intracellular metabolic activation of aflatoxin B1 (AFB1) to a reactive epoxide generates highly mutagenic AFB1-Fapy-dG adducts. Previously, we demonstrated that repair of AFB1-Fapy-dG adducts can be initiated by the DNA glycosylase NEIL1 and that male Neil1-/- mice were significantly more susceptible to AFB1-induced HCC relative to wild-type mice. To investigate the mechanisms underlying this enhanced carcinogenesis, WT and Neil1-/- mice were challenged with a single, 4 mg/kg dose of AFB1 and frequencies and spectra of mutations were analyzed in liver DNAs 2.5 months post-injection using duplex sequencing. The analyses of DNAs from AFB1-challenged mice revealed highly elevated mutation frequencies in the nuclear genomes of both males and females, but not the mitochondrial genomes. In both WT and Neil1-/- mice, mutation spectra were highly similar to the AFB1-specific COSMIC signature SBS24. Relative to wild-type, the NEIL1 deficiency increased AFB1-induced mutagenesis with concomitant elevated HCCs in male Neil1-/- mice. Our data establish a critical role of NEIL1 in limiting AFB1-induced mutagenesis and ultimately carcinogenesis.

Targeted accurate RNA consensus sequencing (tARC-seq) reveals mechanisms of replication error affecting SARS-CoV-2 divergence.

RNA viruses, like SARS-CoV-2, depend on their RNA-dependent RNA polymerases (RdRp) for replication, which is error prone. Monitoring replication errors is crucial for understanding the virus's evolution. Current methods lack the precision to detect rare de novo RNA mutations, particularly in low-input samples such as those from patients. Here we introduce a targeted accurate RNA consensus sequencing method (tARC-seq) to accurately determine the mutation frequency and types in SARS-CoV-2, both in cell culture and clinical samples. Our findings show an average of 2.68 × 10-5 de novo errors per cycle with a C > T bias that cannot be solely attributed to APOBEC editing. We identified hotspots and cold spots throughout the genome, correlating with high or low GC content, and pinpointed transcription regulatory sites as regions more susceptible to errors. tARC-seq captured template switching events including insertions, deletions and complex mutations. These insights shed light on the genetic diversity generation and evolutionary dynamics of SARS-CoV-2.

Access to COVID-19 testing by individuals with housing insecurity during the early days of the COVID-19 pandemic in the United States: a scoping review.

Introduction: The COVID-19 pandemic focused attention on healthcare disparities and inequities faced by individuals within marginalized and structurally disadvantaged groups in the United States. These individuals bore the heaviest burden across this pandemic as they faced increased risk of infection and difficulty in accessing testing and medical care. Individuals experiencing housing insecurity are a particularly vulnerable population given the additional barriers they face. In this scoping review, we identify some of the barriers this high-risk group experienced during the early days of the pandemic and assess novel solutions to overcome these barriers. Methods: A scoping review was performed following PRISMA-Sc guidelines looking for studies focusing on COVID-19 testing among individuals experiencing housing insecurity. Barriers as well as solutions to barriers were identified as applicable and summarized using qualitative methods, highlighting particular ways that proved effective in facilitating access to testing access and delivery. Results: Ultimately, 42 studies were included in the scoping review, with 143 barriers grouped into four categories: lack of cultural understanding, systemic racism, and stigma; medical care cost, insurance, and logistics; immigration policies, language, and fear of deportation; and other. Out of these 42 studies, 30 of these studies also suggested solutions to address them. Conclusion: A paucity of studies have analyzed COVID-19 testing barriers among those experiencing housing insecurity, and this is even more pronounced in terms of solutions to address those barriers. Expanding resources and supporting investigators within this space is necessary to ensure equitable healthcare delivery.

RADx-UP Testing Core: Access to COVID-19 Diagnostics in Community-Engaged Research with Underserved Populations.

Research on the COVID-19 pandemic revealed a disproportionate burden of COVID-19 infection and death among underserved populations and exposed low rates of SARS-CoV-2 testing in these communities. A landmark National Institutes of Health (NIH) funding initiative, the Rapid Acceleration of Diagnostics-Underserved Populations (RADx-UP) program, was developed to address the research gap in understanding the adoption of COVID-19 testing in underserved populations. This program is the single largest investment in health disparities and community-engaged research in the history of the NIH. The RADx-UP Testing Core (TC) provides community-based investigators with essential scientific expertise and guidance on COVID-19 diagnostics. This commentary describes the first 2 years of the TC's experience, highlighting the challenges faced and insights gained to safely and effectively deploy large-scale diagnostics for community-initiated research in underserved populations during a pandemic. The success of RADx-UP shows that community-based research to increase access and uptake of testing among underserved populations can be accomplished during a pandemic with tools, resources, and multidisciplinary expertise provided by a centralized testing-specific coordinating center. We developed adaptive tools to support individual testing strategies and frameworks for these diverse studies and ensured continuous monitoring of testing strategies and use of study data. In a rapidly evolving setting of tremendous uncertainty, the TC provided essential and real-time technical expertise to support safe, effective, and adaptive testing. The lessons learned go beyond this pandemic and can serve as a framework for rapid deployment of testing in response to future crises, especially when populations are affected inequitably.

The multi-tissue landscape of somatic mtDNA mutations indicates tissue-specific accumulation and removal in aging.

Accumulation of somatic mutations in the mitochondrial genome (mtDNA) has long been proposed as a possible mechanism of mitochondrial and tissue dysfunction that occurs during aging. A thorough characterization of age-associated mtDNA somatic mutations has been hampered by the limited ability to detect low-frequency mutations. Here, we used Duplex Sequencing on eight tissues of an aged mouse cohort to detect >89,000 independent somatic mtDNA mutations and show significant tissue-specific increases during aging across all tissues examined which did not correlate with mitochondrial content and tissue function. G→A/C→T substitutions, indicative of replication errors and/or cytidine deamination, were the predominant mutation type across all tissues and increased with age, whereas G→T/C→A substitutions, indicative of oxidative damage, were the second most common mutation type, but did not increase with age regardless of tissue. We also show that clonal expansions of mtDNA mutations with age is tissue- and mutation type-dependent. Unexpectedly, mutations associated with oxidative damage rarely formed clones in any tissue and were significantly reduced in the hearts and kidneys of aged mice treated at late age with elamipretide or nicotinamide mononucleotide. Thus, the lack of accumulation of oxidative damage-linked mutations with age suggests a life-long dynamic clearance of either the oxidative lesions or mtDNA genomes harboring oxidative damage.

Hypoxia induces transgenerational epigenetic inheritance of small RNAs.

Animals sense and adapt to decreased oxygen availability, but whether and how hypoxia exposure in ancestors can elicit phenotypic consequences in normoxia-reared descendants are unclear. We show that hypoxia educes an intergenerational reduction in lipids and a transgenerational reduction in fertility in the nematode Caenorhabditis elegans. The transmission of these epigenetic phenotypes is dependent on repressive histone-modifying enzymes and the argonaute HRDE-1. Feeding naive C. elegans small RNAs extracted from hypoxia-treated worms is sufficient to induce a fertility defect. Furthermore, the endogenous small interfering RNA F44E5.4/5 is upregulated intergenerationally in response to hypoxia, and soaking naive normoxia-reared C. elegans with F44E5.4/5 double-stranded RNA (dsRNA) is sufficient to induce an intergenerational fertility defect. Finally, we demonstrate that labeled F44E5.4/5 dsRNA is itself transmitted from parents to children. Our results suggest that small RNAs respond to the environment and are sufficient to transmit non-genetic information from parents to their naive children.

An atypical RNA quadruplex marks RNAs as vectors for gene silencing.

The addition of poly(UG) ('pUG') repeats to 3' termini of mRNAs drives gene silencing and transgenerational epigenetic inheritance in the metazoan Caenorhabditis elegans. pUG tails promote silencing by recruiting an RNA-dependent RNA polymerase (RdRP) that synthesizes small interfering RNAs. Here we show that active pUG tails require a minimum of 11.5 repeats and adopt a quadruplex (G4) structure we term the pUG fold. The pUG fold differs from known G4s in that it has a left-handed backbone similar to Z-RNA, no consecutive guanosines in its sequence, and three G quartets and one U quartet stacked non-sequentially. The compact pUG fold binds six potassium ions and brings the RNA ends into close proximity. The biological importance of the pUG fold is emphasized by our observations that porphyrin molecules bind to the pUG fold and inhibit both gene silencing and binding of RdRP. Moreover, specific 7-deaza substitutions that disrupt the pUG fold neither bind RdRP nor induce RNA silencing. These data define the pUG fold as a previously unrecognized RNA structural motif that drives gene silencing. The pUG fold can also form internally within larger RNA molecules. Approximately 20,000 pUG-fold sequences are found in noncoding regions of human RNAs, suggesting that the fold probably has biological roles beyond gene silencing.

Resistance of mitochondrial DNA to cadmium and Aflatoxin B1 damage-induced germline mutation accumulation in C. elegans.

Mitochondrial DNA (mtDNA) is prone to mutation in aging and over evolutionary time, yet the processes that regulate the accumulation of de novo mtDNA mutations and modulate mtDNA heteroplasmy are not fully elucidated. Mitochondria lack certain DNA repair processes, which could contribute to polymerase error-induced mutations and increase susceptibility to chemical-induced mtDNA mutagenesis. We conducted error-corrected, ultra-sensitive Duplex Sequencing to investigate the effects of two known nuclear genome mutagens, cadmium and Aflatoxin B1, on germline mtDNA mutagenesis in Caenorhabditis elegans. Detection of thousands of mtDNA mutations revealed pervasive heteroplasmy in C. elegans and that mtDNA mutagenesis is dominated by C:G → A:T mutations generally attributed to oxidative damage. However, there was no effect of either exposure on mtDNA mutation frequency, spectrum, or trinucleotide context signature despite a significant increase in nuclear mutation rate after aflatoxin B1 exposure. Mitophagy-deficient mutants pink-1 and dct-1 accumulated significantly higher levels of mtDNA damage compared to wild-type C. elegans after exposures. However, there were only small differences in mtDNA mutation frequency, spectrum, or trinucleotide context signature compared to wild-type after 3050 generations, across all treatments. These findings suggest mitochondria harbor additional previously uncharacterized mechanisms that regulate mtDNA mutational processes across generations.

Myometrial oxidative stress drives MED12 mutations in leiomyoma.

BACKGROUND: More than 70% of leiomyomas (LM) harbor MED12 mutations, primarily in exon 2 at c.130-131(GG). The cause of MED12 mutations in myometrial cells remains largely unknown. We hypothesized that increased ROS promotes MED12 mutations in myometrial cells through the oxidation of guanine nucleotides followed by misrepair. METHODS: Genomic oxidative burden (8-OHdG) was evaluated in vitro and in vivo by immunohistochemistry. MED12 mutations were examined by Sanger sequencing and deep sequencing. Transcriptome examined by RNA-seq was performed in myometrium with and without LM, in primary myometrial cells treated with ROS. 8-OHdG mediated misrepair was analyzed by CRISPR/Cas9. RESULTS: Uteri with high LM burden had a significantly higher rate of MED12 mutations than uteri with low LM burden. Compelling data suggest that the uterus normally produces reactive oxidative species (ROS) in response to stress, and ROS levels in LM are elevated due to metabolic defects. We demonstrated that genomic oxidized guanine (8-OHdG) was found at a significantly higher level in the myometrium of uteri that had multiple LM compared to myometrium without LM. Transcriptome and pathway analyses detected ROS stress in myometrium with LM. Targeted replacement of guanine with 8-OHdG at MED12 c.130 by CRISPR/Cas9 significantly increased the misrepair of G>T. Exposure of primary myometrial cells to oxidative stress in vitro increased misrepair/mutations as detected by duplex sequencing. CONCLUSIONS: Together, our data identified a clear connection between increased myometrial oxidative stress and a high rate of MED12 mutations that may underlie the risk of LM development and severity in women of reproductive age.

RNA polymerase inaccuracy underlies SARS-CoV-2 variants and vaccine heterogeneity.

Both the SARS-CoV-2 virus and its mRNA vaccines depend on RNA polymerases (RNAP)1,2; however, these enzymes are inherently error-prone and can introduce variants into the RNA3. To understand SARS-CoV-2 evolution and vaccine efficacy, it is critical to identify the extent and distribution of errors introduced by the RNAPs involved in each process. Current methods lack the sensitivity and specificity to measure de novo RNA variants in low input samples like viral isolates3. Here, we determine the frequency and nature of RNA errors in both SARS-CoV-2 and its vaccine using a targeted Accurate RNA Consensus sequencing method (tARC-seq). We found that the viral RNA-dependent RNAP (RdRp) makes ~1 error every 10,000 nucleotides - higher than previous estimates4. We also observed that RNA variants are not randomly distributed across the genome but are associated with certain genomic features and genes, such as S (Spike). tARC-seq captured a number of large insertions, deletions and complex mutations that can be modeled through non-programmed RdRp template switching. This template switching feature of RdRp explains many key genetic changes observed during the evolution of different lineages worldwide, including Omicron. Further sequencing of the Pfizer-BioNTech COVID-19 vaccine revealed an RNA variant frequency of ~1 in 5,000, meaning most of the vaccine transcripts produced in vitro by T7 phage RNAP harbor a variant. These results demonstrate the extraordinary genetic diversity of viral populations and the heterogeneous nature of an mRNA vaccine fueled by RNAP inaccuracy. Along with functional studies and pandemic data, tARC-seq variant spectra can inform models to predict how SARS-CoV-2 may evolve. Finally, our results may help improve future vaccine development and study design as mRNA therapies continue to gain traction.

The Complicated Nature of Somatic mtDNA Mutations in Aging.

Mitochondria are the main source of energy used to maintain cellular homeostasis. This aspect of mitochondrial biology underlies their putative role in age-associated tissue dysfunction. Proper functioning of the electron transport chain (ETC), which is partially encoded by the extra-nuclear mitochondrial genome (mtDNA), is key to maintaining this energy production. The acquisition of de novo somatic mutations that interrupt the function of the ETC have long been associated with aging and common diseases of the elderly. Yet, despite over 30 years of study, the exact role(s) mtDNA mutations play in driving aging and its associated pathologies remains under considerable debate. Furthermore, even fundamental aspects of age-related mtDNA mutagenesis, such as when mutations arise during aging, where and how often they occur across tissues, and the specific mechanisms that give rise to them, remain poorly understood. In this review, we address the current understanding of the somatic mtDNA mutations, with an emphasis of when, where, and how these mutations arise during aging. Additionally, we highlight current limitations in our knowledge and critically evaluate the controversies stemming from these limitations. Lastly, we highlight new and emerging technologies that offer potential ways forward in increasing our understanding of somatic mtDNA mutagenesis in the aging process.

Secondary structure prediction for RNA sequences including N6-methyladenosine.

There is increasing interest in the roles of covalently modified nucleotides in RNA. There has been, however, an inability to account for modifications in secondary structure prediction because of a lack of software and thermodynamic parameters. We report the solution for these issues for N6-methyladenosine (m6A), allowing secondary structure prediction for an alphabet of A, C, G, U, and m6A. The RNAstructure software now works with user-defined nucleotide alphabets of any size. We also report a set of nearest neighbor parameters for helices and loops containing m6A, using experiments. Interestingly, N6-methylation decreases folding stability for adenosines in the middle of a helix, has little effect on folding stability for adenosines at the ends of helices, and increases folding stability for unpaired adenosines stacked on a helix. We demonstrate predictions for an N6-methylation-activated protein recognition site from MALAT1 and human transcriptome-wide effects of N6-methylation on the probability of adenosine being buried in a helix.

Nuclear magnetic resonance reveals a two hairpin equilibrium near the 3'-splice site of influenza A segment 7 mRNA that can be shifted by oligonucleotides.

Influenza A kills hundreds of thousands of people globally every year and has the potential to generate more severe pandemics. Influenza A's RNA genome and transcriptome provide many potential therapeutic targets. Here, nuclear magnetic resonance (NMR) experiments suggest that one such target could be a hairpin loop of 8 nucleotides in a pseudoknot that sequesters a 3' splice site in canonical pairs until a conformational change releases it into a dynamic 2 × 2-nt internal loop. NMR experiments reveal that the hairpin loop is dynamic and able to bind oligonucleotides as short as pentamers. A 3D NMR structure of the complex contains 4 and likely 5 bp between pentamer and loop. Moreover, a hairpin sequence was discovered that mimics the equilibrium of the influenza hairpin between its structure in the pseudoknot and upon release of the splice site. Oligonucleotide binding shifts the equilibrium completely to the hairpin secondary structure required for pseudoknot folding. The results suggest this hairpin can be used to screen for compounds that stabilize the pseudoknot and potentially reduce splicing.

Nuclear Magnetic Resonance Spectra and AMBER OL3 and ROC-RNA Simulations of UCUCGU Reveal Force Field Strengths and Weaknesses for Single-Stranded RNA.

Single-stranded regions of RNA are important for folding of sequences into 3D structures and for design of therapeutics targeting RNA. Prediction of ensembles of 3D structures for single-stranded regions often involves classical mechanical approximations of interactions defined by quantum mechanical calculations on small model systems. Nuclear magnetic resonance (NMR) spectra and molecular dynamics (MD) simulations of short single strands provide tests for how well the approximations model many of the interactions. Here, the NMR spectra for UCUCGU at 2, 15, and 30 °C are compared to simulations with the AMBER force fields, OL3 and ROC-RNA. This is the first such comparison to an oligoribonucleotide containing an internal guanosine nucleotide (G). G is particularly interesting because of its many H-bonding groups, large dipole moment, and proclivity for both syn and anti conformations. Results reveal formation of a G amino to phosphate non-bridging oxygen H-bond. The results also demonstrate dramatic differences in details of the predicted structures. The variations emphasize the dependence of predictions on individual parameters and their balance with the rest of the force field. The NMR data can serve as a benchmark for future force fields.

Probing the mechanisms of two exonuclease domain mutators of DNA polymerase ϵ.

We report the properties of two mutations in the exonuclease domain of the Saccharomyces cerevisiae DNA polymerase ϵ. One, pol2-Y473F, increases the mutation rate by about 20-fold, similar to the catalytically dead pol2-D290A/E290A mutant. The other, pol2-N378K, is a stronger mutator. Both retain the ability to excise a nucleotide from double-stranded DNA, but with impaired activity. pol2-Y473F degrades DNA poorly, while pol2-N378K degrades single-stranded DNA at an elevated rate relative to double-stranded DNA. These data suggest that pol2-Y473F reduces the capacity of the enzyme to perform catalysis in the exonuclease active site, while pol2-N378K impairs partitioning to the exonuclease active site. Relative to wild-type Pol ϵ, both variants decrease the dNTP concentration required to elicit a switch between proofreading and polymerization by more than an order of magnitude. While neither mutation appears to alter the sequence specificity of polymerization, the N378K mutation stimulates polymerase activity, increasing the probability of incorporation and extension of a mismatch. Considered together, these data indicate that impairing the primer strand transfer pathway required for proofreading increases the probability of common mutations by Pol ϵ, elucidating the association of homologous mutations in human DNA polymerase ϵ with cancer.

A replication-linked mutational gradient drives somatic mutation accumulation and influences germline polymorphisms and genome composition in mitochondrial DNA.

Mutations in mitochondrial DNA (mtDNA) cause maternally inherited diseases, while somatic mutations are linked to common diseases of aging. Although mtDNA mutations impact health, the processes that give rise to them are under considerable debate. To investigate the mechanism by which de novo mutations arise, we analyzed the distribution of naturally occurring somatic mutations across the mouse and human mtDNA obtained by Duplex Sequencing. We observe distinct mutational gradients in G→A and T→C transitions delimited by the light-strand origin and the mitochondrial Control Region (mCR). The gradient increases unequally across the mtDNA with age and is lost in the absence of DNA polymerase γ proofreading activity. In addition, high-resolution analysis of the mCR shows that important regulatory elements exhibit considerable variability in mutation frequency, consistent with them being mutational 'hot-spots' or 'cold-spots'. Collectively, these patterns support genome replication via a deamination prone asymmetric strand-displacement mechanism as the fundamental driver of mutagenesis in mammalian DNA. Moreover, the distribution of mtDNA single nucleotide polymorphisms in humans and the distribution of bases in the mtDNA across vertebrate species mirror this gradient, indicating that replication-linked mutations are likely the primary source of inherited polymorphisms that, over evolutionary timescales, influences genome composition during speciation.

piRNAs coordinate poly(UG) tailing to prevent aberrant and perpetual gene silencing.

Noncoding RNAs have emerged as mediators of transgenerational epigenetic inheritance (TEI) in a number of organisms. A robust example of such RNA-directed TEI is the inheritance of gene-silencing states following RNA interference (RNAi) in the metazoan C. elegans. During RNAi inheritance, gene silencing is transmitted by a self-perpetuating cascade of siRNA-directed poly(UG) tailing of mRNA fragments (pUGylation), followed by siRNA synthesis from poly(UG)-tailed mRNA templates (termed pUG RNA/siRNA cycling). Despite the self-perpetuating nature of pUG RNA/siRNA cycling, RNAi inheritance is finite, suggesting that systems likely exist to prevent indefinite RNAi-triggered gene silencing. Here we show that, in the absence of Piwi-interacting RNAs (piRNAs), an animal-specific class of small noncoding RNA, RNAi-based gene silencing can become essentially permanent, lasting at near 100% penetrance for more than 5 years and hundreds of generations. This perpetual gene silencing is mediated by continuous pUG RNA/siRNA cycling. Further, we find that piRNAs coordinate endogenous RNAi pathways to prevent germline-expressed genes, which are not normally subjected to TEI, from entering a state of continual and irreversible epigenetic silencing also mediated by continuous maintenance of pUG RNA/siRNA cycling. Together, our results show that one function of C. elegans piRNAs is to insulate germline-expressed genes from aberrant and runaway inactivation by the pUG RNA/siRNA epigenetic inheritance system.

PolyG-DS: An ultrasensitive polyguanine tract-profiling method to detect clonal expansions and trace cell lineage.

Polyguanine tracts (PolyGs) are short guanine homopolymer repeats that are prone to accumulating mutations when cells divide. This feature makes them especially suitable for cell lineage tracing, which has been exploited to detect and characterize precancerous and cancerous somatic evolution. PolyG genotyping, however, is challenging because of the inherent biochemical difficulties in amplifying and sequencing repetitive regions. To overcome this limitation, we developed PolyG-DS, a next-generation sequencing (NGS) method that combines the error-correction capabilities of duplex sequencing (DS) with enrichment of PolyG loci using CRISPR-Cas9-targeted genomic fragmentation. PolyG-DS markedly reduces technical artifacts by comparing the sequences derived from the complementary strands of each original DNA molecule. We demonstrate that PolyG-DS genotyping is accurate, reproducible, and highly sensitive, enabling the detection of low-frequency alleles (<0.01) in spike-in samples using a panel of only 19 PolyG markers. PolyG-DS replicated prior results based on PolyG fragment length analysis by capillary electrophoresis, and exhibited higher sensitivity for identifying clonal expansions in the nondysplastic colon of patients with ulcerative colitis. We illustrate the utility of this method for resolving the phylogenetic relationship among precancerous lesions in ulcerative colitis and for tracing the metastatic dissemination of ovarian cancer. PolyG-DS enables the study of tumor evolution without prior knowledge of tumor driver mutations and provides a tool to perform cost-effective and easily scalable ultra-accurate NGS-based PolyG genotyping for multiple applications in biology, genetics, and cancer research.

Translational precision medicine: an industry perspective.

In the era of precision medicine, digital technologies and artificial intelligence, drug discovery and development face unprecedented opportunities for product and business model innovation, fundamentally changing the traditional approach of how drugs are discovered, developed and marketed. Critical to this transformation is the adoption of new technologies in the drug development process, catalyzing the transition from serendipity-driven to data-driven medicine. This paradigm shift comes with a need for both translation and precision, leading to a modern Translational Precision Medicine approach to drug discovery and development. Key components of Translational Precision Medicine are multi-omics profiling, digital biomarkers, model-based data integration, artificial intelligence, biomarker-guided trial designs and patient-centric companion diagnostics. In this review, we summarize and critically discuss the potential and challenges of Translational Precision Medicine from a cross-industry perspective.