Single Base Pair Research Articles

Rescue of the endogenous FVIII expression in hemophilia A mice using CRISPR Cas9 mRNA LNPs

Gene editing provides a promising alternative approach that may achieve long-term FVIII expression for hemophilia A (HemA) treatment. In this study, we investigated in vivo correction of a mutant factor VIII (FVIII) gene in HemA mice. We first developed MC3-based LNPs for efficient mRNA delivery into liver sinusoidal endothelial cells (LSECs), the major site of FVIII biosynthesis. To target a five base pair deletion in FVIII exon 1 in a specific HemA mouse strain, we injected LNPs encapsulating Cas9 mRNA and specifically designed sgRNAs intravenously for in vivo gene editing of the mutant FVIII. Indel variants generated at the mutant site contained mostly a single base pair deletion , resulting in frame-shift correction of FVIII gene. Sustained endogenous FVIII activity up to 6% was achieved over 26 weeks in treated HemA mice. Sequencing data indicated an average gene editing rate of 15.3% in LSECs. Our study suggests that optimized MC3 LNP formulations, combined with CRISPR/Cas9 technology, can effectively correct the mutant FVIII gene in LSECs and restore FVIII activity for therapeutic treatment of HemA.

The Thermodynamic and Kinetic Properties of the dA-rU DNA-RNA Hybrid Base Pair Investigated via Molecular Dynamics Simulations.

DNA-RNA hybrid duplexes play essential roles during the reverse transcription of RNA viruses and DNA replication. The opening and conformation changes of individual base pairs are critical to their biological functions. However, the microscopic mechanisms governing base pair closing and opening at the atomic level remain poorly understood. In this study, we investigated the thermodynamic and kinetic parameters of the dA-rU base pair in a DNA-RNA hybrid duplex using 4 μs all-atom molecular dynamics (MD) simulations at different temperatures. Our results showed that the thermodynamic parameters of the dA-rU base pair aligned with the predictions of the nearest-neighbor model and were close to those of the AU base pair in RNA. The temperature dependence of the average lifetimes of both the open and the closed states, as well as the transition path times, were obtained. The free-energy barrier for a single base pair opening and closing arises from an increase in enthalpy due to the disruption of the base-stacking interactions and hydrogen bonding, along with an entropy loss attributed to the accompanying restrictions, such as torsional angle constraints and solvent viscosity.

CRISPR-Mediated Gene Editing Tool Development for Patients with Usher Syndrome Type II

Usher syndrome, a rare genetic disease affecting hearing and vision, is an autosomal recessive condition affecting 4 to 17 per 100,000 people and accounting for about 50 percent of hereditary deaf-blindness cases. Usher syndrome type II is caused by mutations of the USH2A gene, a 15.7 kb gene encoding the protein usherin, which localizes to photoreceptor cilium and cochlear hair cells. This subtype of Usher syndrome includes hearing loss from birth and progressive loss of vision, prompting retinitis pigmentosa (RP). Prime editing is a viable option for addressing USH2A mutations and restoring usherin expression and function. To begin, a HEK293T cell line needs to be established for subsequent testing of prime editors. Since prime editing to make the cell line was unsuccessful, homology-directed repair was done and successfully created a cell population with the mutation of interest (single base-pair deletion). Now that this cell population has been created, prime editors will be tested and optimized to address the mutation in a HEK293T cell line and potentially patient-derived induced pluripotent stem cells that are differentiated into photoreceptor precursor cells.

Aminoacyl-tRNA synthetases

Aminoacyl-tRNA synthetases hold the key to the genetic code and assign nucleic acid-based codons to amino acids, the building blocks of proteins. In their ability to recognize identity elements on transfer RNAs (tRNAs), some as simple as a single base pair, they ensure that the same proteins are formed each time information embedded in DNA is transcribed into messenger RNA (mRNA) and translated into proteins (Figure 1A). Thus, aminoacyl-tRNA synthetase active sites are conserved; however, since their evolutionary origin, their functions have been co-opted, expanded on and played novel roles during evolution. Below, we provide an overview of the many functions of aminoacyl-tRNA synthetases - from their role in translation, one of the most fundamental processes of all life, to newly discovered, diverse functions.

Prediction and functional interpretation of inter-chromosomal genome architecture from DNA sequence with TwinC.

Three-dimensional nuclear DNA architecture comprises well-studied intra-chromosomal (cis) folding and less characterized inter-chromosomal (trans) interfaces. Current predictive models of 3D genome folding can effectively infer pairwise cis-chromatin interactions from the primary DNA sequence but generally ignore trans contacts. There is an unmet need for robust models of trans-genome organization that provide insights into their underlying principles and functional relevance. We present TwinC, an interpretable convolutional neural network model that reliably predicts trans contacts measurable through genome-wide chromatin conformation capture (Hi-C). TwinC uses a paired sequence design from replicate Hi-C experiments to learn single base pair relevance in trans interactions across two stretches of DNA. The method achieves high predictive accuracy (AUROC=0.80) on a cross-chromosomal test set from Hi-C experiments in heart tissue. Mechanistically, the neural network learns the importance of compartments, chromatin accessibility, clustered transcription factor binding and G-quadruplexes in forming trans contacts. In summary, TwinC models and interprets trans genome architecture, shedding light on this poorly understood aspect of gene regulation.

Genetic diagnostic outcomes from a 10-year research programme in autism in Aotearoa New Zealand

ABSTRACT Autism is a relatively common neurodevelopmental difference with considerable phenotypic heterogeneity impacting cognitive, sensory, and social processing, and often co-occurs with other conditions. Therefore, there is not a one-size-fits-all clinical support pathway for autistic individuals following diagnosis. DNA sequencing technology has enabled the discovery of genes causative of, or associated with, autism. Unsurprisingly, genetic heterogeneity goes hand-in-hand with the phenotypic heterogeneity for this condition; with causative genetic variation ranging from single base pair changes to complex chromosomal rearrangements in more than 100 different genes. This study captures a snapshot (201 individuals) of the autistic population (both clinically referred and self-referred) in Aotearoa New Zealand and documents a decade’s research effort to refine diagnosis using a flexible and customised genome-wide sequencing approach. The diagnostic yield in this phenotypically disparate cohort was 12.9%, with an additional 15.9% of individuals harbouring ‘likely causal’ variants, providing the groundwork to tailor clinical, social, and educational care. Importantly, this study reveals the diagnostic utility of customised genetic screening for autism across a phenotypically diverse autistic population.

Identification of streptomycin-resistant Erwinia amylovora in Iowa.

Erwinia amylovora is a bacterial pathogen that causes fire blight, an important disease in apples and pears. Applying the antibiotic streptomycin during the phenological bloom stage is considered the most effective management tactic for fire blight. Though streptomycin-resistant (SmR) E. amylovora populations have emerged in major U.S. apple-producing regions, antibiotic resistance data for medium to small-sized apple-producing regions like the Midwest is still lacking. This short communication collected symptomatic fire blight samples from Iowa apple orchards during 2022 and 2023, where recent fire blight outbreaks persisted despite streptomycin use. Among E. amylovora isolates from seven counties in central and eastern Iowa, around 90% of them were SmR. All SmR isolates exhibited a single base pair mutation in codon 43 of the rpsL gene, conferring resistance to streptomycin levels exceeding 1,000 μg/mL. Through clustered regularly interspaced short palindromic repeat (CRISPR) analysis, we characterized two E. amylovora genotypes unique to our region. Whole genome sequencing of one representative SmR isolate, IA01, confirmed its CRISPR genotype and subsequent phylogenetic analysis suggested that IA01 is genetically similar to Michigan isolates and distinct from those in eastern and western regions of North America. Furthermore, the disease-causing ability of IA01 was comparable to that of the highly virulent Ea110 strain, a streptomycin sensitive strain isolated from Michigan, in immature pears. Overall, this study underscores the urgent need for regional strategies to address antibiotic resistance and provide insights into its genetic basis and geographic distribution which are crucial for sustainable orchard management.

The effect of the loop on the thermodynamic and kinetic of single base pair in pseudoknot.

RNA pseudoknots are RNA molecules with specialized three-dimensional structures that play important roles in various biological processes. To understand the functions and mechanisms of pseudoknots, it is essential to elucidate their structures and folding pathways. The most fundamental step in RNA folding is the opening and closing of a base pair. The effect of flexible loops on the base pair in pseudoknots remains unclear. In this work, we use molecular dynamics simulations and Markov state model to study the configurations, thermodynamic and kinetic of single base pair in pseudoknots. We find that the presence of the loop leads to a trap state. In addition, the rate-limiting step for the formation of base pair is the disruption of the trap state, rather than the open state to the closed state, which is quite different from the previous studies on non-pseudoknot RNA. For the thermodynamic parameters in pseudoknots, we find that the entropy difference upon opening the base pair between this simulation and the nearest-neighbor model results from the different entropy of different lengths of loop in solution. The thermodynamic parameters of the stack in pseudoknot are close to the nearest-neighbor parameters. The bases on the loop have different distribution patterns in different states, and the slow transition states of the loop are determined by the orientation of the bases.

7-Deaza-2'-deoxyisoguanosine, a Noncanonical Nucleoside for Nucleic Acid Code Expansion and New DNA Constructs: Nucleobase Functionalization of Inverse Watson-Crick and Purine-Purine Base Pairs.

7-Deaza-2'-deoxyisoguanosine forms stable inverse Watson-Crick base pairs with 5-methyl-2'-deoxyisocytidine and purine-purine base pairs with 2'-deoxyguanosine or 5-aza-7-deaza-2'-deoxyguanosine. Both base pairs expand the genetic coding system. The manuscript reports on the functionalization of these base pairs with halogen atoms and clickable side chains introduced at 7-position of the 7-deazapurine base. Oligonucleotides containing the functionalized base pairs were prepared by solid-phase synthesis. To this end, a series of phosphoramidites were synthesized and clickable side chains with short and long linkers were incorporated in oligonucleotides. Fluorescent pyrene conjugates were obtained by postmodification. Functionalization of DNA with a single inverse Watson-Crick base pair by halogens or clickable residues has only a minor impact on duplex stability. Pyrene click adducts increase (long linker) or decrease (short linker) the double helix stability. Stable hybrid duplexes were constructed containing three consecutive purine-purine pairs of 7-functionalized 7-deaza-2'-deoxyisoguanine with guanine or 5-aza-7-deazaguanine in the center and Watson-Crick pairs at both ends. The incorporation of a hybrid base pair tract of 7-deaza-2'-deoxyisoguanosine/5-aza-7-deaza-2'-deoxyguanosine pairs stabilizes the double helix strongly. Fluorescence intensity of pyrene short linker adducts increased when the 7-deazapurine base was positioned opposite to 5-methylisocytosine (inverse base pair) compared to purine-purine base pairs with guanine or 5-aza-7-deazaguanine in opposite positions. For long liker adducts, the situation is more complex. Circular dichroism (CD) spectra of purine DNA differ to those of Watson-Crick double helices and are indicative for the new DNA constructs. The impact of 7-deaza-2'-deoxyisoguanine base pair functionalization is studied for the first time and all experimental details are reported to prepare DNA functionalized at the 7-deazaisoguanine site. The influence of single and multiple incorporations on DNA structure and stability is shown. Clickable residues introduced at the 7-position of the 7-deazaisoguanine base provide handles for Huisgen-Sharpless-Meldal click cycloadditions without harming the stability of purine-pyrimidine and purine-purine base pairs. Other chemistries might be used for bioconjugation. Our investigation paves the way for the functionalization of a new DNA related recognition system expanding the common Watson-Crick regime.

Genotyping SNPs and Indels: A method to improve the scope and sensitivity of High-Resolution melt (HRM) analysis based applications

High-resolution melt (HRM) analysis is a closed-tube technique for detecting single nucleotide polymorphisms (SNPs). However, it has limited use in high-resolution melting devices, even those with high thermal accuracy (HTA). In addition to the cost of switching to these specialized devices, the presence of nearest neighbour neutral changes (class III, IV SNPs and small indels) made HRM-based assays a challenging task due to reduced sensitivity. This study aimed to design a common modified competitive amplification of differently melting amplicons (CADMA)-based assay to address these challenges by generating allele-specific qPCR products that are detectable on most qPCR platforms.For this study, SNPs were selected from all four classes of SNPs (class I: C/T or G/A mutation; class II: C/A or G/T mutation; class III: G/C mutation; class IV: A/T mutation). A single base pair and 19 bp indels were also chosen to simulate how CADMA primers could be designed for indels of varying lengths. The melting temperatures (Tm) were determined using IDT oligoAnalyzer. qPCR and melt data acquisition were performed on the CFX96 qPCR platform, and the melt curve data were analyzed using Precision Melt software (Bio-Rad, USA). The clusters for different genotypes were successfully identified with the aid of the control samples, and Tm predictions were carried out using the uMelt batch and Tm online tools for comparison.Using HRM-qPCR assays based on the modified CADMA method, genotyping of various SNPs was successfully carried out. For some SNPs, similarly shaped melt curves were observed for homozygotes and heterozygotes, making shape-based genotype prediction difficult. The Tm values calculated via the Blake and Delcourts (1998) method were the closest to the experimental Tm values after adjusting for the salt concentration.Since HRM assays usually depend on the ΔTm caused by mutations, they are prone to a high error rate due to nearest neighbour neutral changes. The technique developed in this study significantly reduces the failure rates in HRM-based genotyping and could be applied to any SNP or indel in any platform. It is crucial to have a deep understanding of the melt instrument, its accuracy and the nature of the target (SNP class or indel length and GC content of the flanking region). Furthermore, the availability of controls is essential for a high success rate.

Analyzing single-cell bisulfite sequencing data with MethSCAn

Single-cell bisulfite sequencing (scBS) is a technique that enables the assessment of DNA methylation at single-base pair and single-cell resolution. The analysis of large datasets obtained from scBS requires preprocessing to reduce the data size, improve the signal-to-noise ratio and provide interpretability. Typically, this is achieved by dividing the genome into large tiles and averaging the methylation signals within each tile. Here we demonstrate that this coarse-graining approach can lead to signal dilution. We propose improved strategies to identify more informative regions for methylation quantification and a more accurate quantitation method than simple averaging. Our approach enables better discrimination of cell types and other features of interest and reduces the need for large numbers of cells. We also present an approach to detect differentially methylated regions between groups of cells and demonstrate its ability to identify biologically meaningful regions that are associated with genes involved in the core functions of specific cell types. Finally, we present the software tool MethSCAn for scBS data analysis (https://anders-biostat.github.io/MethSCAn).

Influence of empirically derived filtering parameters, amplicon sequence variant, and operational taxonomic unit pipelines on assessing rumen microbial diversity

Microbes play an important role in human and animal health, as well as animal productivity. The host microbial interactions within ruminants play a critical role in animal health and productivity and provide up to 70% of the animal's energy needs in the form of fermentation products. As such, many studies have investigated microbial community composition to understand the microbial community changes and factors that affect microbial colonization and persistence. Advances in next-generation sequencing technologies and the low cost of sequencing have led many studies to use 16S rDNA-based analysis tools for interrogation of microbiomes at a much finer scale than traditional culturing. However, methods that rely on single base pair differences for bacterial taxa clustering may inflate or underestimate diversity, leading to inaccurate identification of bacterial diversity. Therefore, in this study, we sequenced mock communities of known membership and abundance to establish filtration parameters to reduce the inflation of microbial diversity due to PCR and sequencing errors. Additionally, we evaluated the effect of the resulting filtering parameters proposed using established bioinformatic pipelines on a study consisting of Holstein and Jersey cattle to identify breed and treatment effects on the bacterial community composition and the impact of filtering on global microbial community structure analysis and results. Filtration resulted in a sharp reduction in bacterial taxa identified, yet retain most sequencing data (retaining >79% of sequencing reads) when analyzed using 3 different microbial analysis pipelines (DADA2, Mothur, USEARCH). After filtration, conclusions from α-diversity and β-diversity tests showed very similar results across all analysis methods. The mock community-based filtering parameters proposed in this study help provide a more realistic estimation of bacterial diversity. Additionally, filtration reduced the variation between microbiome analysis methods and helped to identify microbial community differences that could have been missed due to the large animal-to-animal variation observed in the unfiltered data. As such, we believe the new filtering parameters described in this study will help to obtain diversity estimates that are closer to realistic values, improve the ability to detecting microbial community differences, and help to better understand microbial community changes in 16S rDNA-based studies.

A single base pair duplication in the SLC33A1 gene is associated with fetal losses and neonatal lethality in Manech Tête Rousse dairy sheep.

We recently discovered that the Manech Tête Rousse (MTR) deficient homozygous haplotype 2 (MTRDHH2) probably carries a recessive lethal mutation in sheep. In this study, we fine-mapped this region through whole-genome sequencing of five MTRDHH2 heterozygous carriers and 95 non-carriers from various ovine breeds. We identified a single base pair duplication within the SLC33A1 gene, leading to a frameshift mutation and a premature stop codon (p.Arg246Alafs*3). SLC33A1 encodes a transmembrane transporter of acetyl-coenzyme A that is crucial for cellular metabolism. To investigate the lethality of this mutation in homozygous MTR sheep, we performed at-risk matings using artificial insemination (AI) between heterozygous SLC33A1 variant carriers (SLC33A1_dupG). Pregnancy was confirmed 15 days post-AI using a blood test measuring interferon Tau-stimulated MX1 gene expression. Ultrasonography between 45 and 60 days post-AI revealed a 12% reduction in AI success compared with safe matings, indicating embryonic/fetal loss. This was supported by the MX1 differential expression test suggesting fetal losses between 15 and 60 days of gestation. We also observed a 34.7% pre-weaning mortality rate in 49 lambs born from at-risk matings. Homozygous SLC33A1_dupG lambs accounted for 47% of this mortality, with deaths occurring mostly within the first 5 days without visible clinical signs. Therefore, appropriate management of SLC33A1_dupG with an allele frequency of 0.04 in the MTR selection scheme would help increase overall fertility and lamb survival.

Abstract IA020: Functional studies of genetic variation using precision genome editing

Abstract Tumor genomes often harbor a complex spectrum of single nucleotide mutations, small indels, and large chromosome rearrangements that can perturb coding and non-coding regions of the genome in ways that remain poorly understood. The mutational processes responsible for the genesis of these events are also not static; instead, they continue to operate throughout disease evolution and further diversify the genetic and phenotypic landscape of cancer cells. Mounting evidence suggests that tumor genotype can be an important determinant of disease progression and therapy response, such as by modulating the sensitivity or resistance of cancer cells to mutant-specific drugs. These types of observations have motivated efforts to treat cancers with genome-informed therapies, highlighting the need to understand how different genetic variants precisely affect gene function and overall tumor phenotypes. Variant-function maps that provide a mechanistic understanding of the biology driven by cancer-associated mutations are needed to design these types of treatment strategies. Precision genome editing technologies like base and prime editing are uniquely suited to tackle this problem. Nevertheless, deploying these methods for systematic variant-function studies and disease modeling in vivo has not been straightforward due to lack of robust and scalable platforms capable of assessing editing efficiency and precision, particularly at endogenous loci. With this goal in mind, we previously developed and applied high-throughput base editing ‘sensor’ approaches that link endogenous genome editing outcomes with synthetic DNA-based readouts and cellular fitness measurements. Using these approaches, we showed that several uncharacterized mutant p53 alleles drive cancer cell proliferation and in vivo tumor development. Building upon this work, we recently developed new prime editing guide RNA design tools and sensor-based approaches that similarly couple quantitative editing outcomes to cellular fitness, allowing us to significantly expand the breadth of cancer-associated mutations that can be interrogated using these precision genome editing technologies. In this talk, I will describe ongoing work using base and prime editing sensor libraries to probe the biological impact of thousands of functionally-distinct genetic variants in diverse types of protein-coding genes and families to learn how these influence various cancer phenotypes. I will also discuss how the implementation of these modular technologies will allow researchers to functionally disentangle the impact of endogenous and exogenous mutational processes on functional selection and tumor evolution with close to single base pair resolution. This generalizable precision genome editing framework will facilitate the functional interrogation of genetic variants across diverse biological contexts, providing much needed insight into cancer variant-function relationships that could be leveraged to develop more precise cancer treatment paradigms, including synthetic lethal therapies that exploit tumor genotype. Citation Format: Francisco J. Sánchez-Rivera, Samuel I. Gould, Alexandra N. Wuest, Kexin Dong, Grace A. Johnson, Alvin Hsu, Varun K. Narendra, Ondine Atwa, Stuart S. Levine, David R. Liu. Functional studies of genetic variation using precision genome editing [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: Expanding and Translating Cancer Synthetic Vulnerabilities; 2024 Jun 10-13; Montreal, Quebec, Canada. Philadelphia (PA): AACR; Mol Cancer Ther 2024;23(6 Suppl):Abstract nr IA020.

Prediction of Deleterious Single Amino Acid Polymorphisms with a Consensus Holdout Sampler.

Single Amino Acid Polymorphisms (SAPs) or nonsynonymous Single Nucleotide Variants (nsSNVs) are the most common genetic variations. They result from missense mutations where a single base pair substitution changes the genetic code in such a way that the triplet of bases (codon) at a given position is coding a different amino acid. Since genetic mutations sometimes cause genetic diseases, it is important to comprehend and foresee which variations are harmful and which ones are neutral (not causing changes in the phenotype). This can be posed as a classification problem. Computational methods using machine intelligence are gradually replacing repetitive and exceedingly overpriced mutagenic tests. By and large, uneven quality, deficiencies, and irregularities of nsSNVs datasets debase the convenience of artificial intelligence-based methods. Subsequently, strong and more exact approaches are needed to address these problems. In the present work paper, we show a consensus classifier built on the holdout sampler, which appears strong and precise and outflanks all other popular methods. We produced 100 holdouts to test the structures and diverse classification variables of diverse classifiers during the training phase. The finest performing holdouts were chosen to develop a consensus classifier and tested using a k-fold (1 ≤ k ≤5) cross-validation method. We also examined which protein properties have the biggest impact on the precise prediction of the effects of nsSNVs. Our Consensus Holdout Sampler outflanks other popular algorithms, and gives excellent results, highly accurate with low standard deviation. The advantage of our method emerges from using a tree of holdouts, where diverse LM/AI-based programs are sampled in diverse ways.

Isotope-Edited Variable Temperature Infrared Spectroscopy for Measuring Transition Temperatures of Single A-T Watson-Crick Base Pairs in DNA Duplexes.

Experimental methods to determine transition temperatures for individual base pair melting events in DNA duplexes are lacking despite intense interest in these thermodynamic parameters. Here, we determine the dimensions of the thymine (T) C2═O stretching vibration when it is within the DNA duplex via isotopic substitutions at other atomic positions in the structure. First, we determined that this stretching state was localized enough to specific atoms in the molecule to make submolecular scale measurements of local structure and stability in high molecular weight complexes. Next, we develop a new isotope-edited variable temperature infrared method to measure melting transitions at various locations in a DNA structure. As an initial test of this "sub-molecular scale thermometer", we applied our T13C2 difference infrared signal to measure location-dependent melting temperatures (TmL) in a DNA duplex via variable temperature attenuated total reflectance Fourier transform infrared (VT-ATR-FTIR) spectroscopy. We report that the TmL of a single Watson-Crick A-T base pair near the end of an A-T rich sequence (poly T) is ∼34.9 ± 0.7°C. This is slightly lower than the TmL of a single base pair near the middle position of the poly T sequence (TmL ∼35.6±0.2°C). In addition, we also report that the TmL of a single Watson-Crick A-T base pair near the end of a 50% G-C sequence (12-mer) is ∼52.5 ± 0.3°C, which is slightly lower than the global melting Tm of the 12-mer sequence (TmL ∼54.0±0.9°C). Our results provide direct physical evidence for end fraying in DNA sequences with our novel spectroscopic methods.

Localized Surface Plasmon Resonance Optical Biosensor for Simple Detection of Deoxyribonucleic Acid Mismatches

Optical biosensors are optical technologies that evaluate changes in the refractive index as they monitor non‐covalent molecular interactions in real time. These make use of unsophisticated, label‐free analytical approaches, which do not require dyes to produce a visible signal. In this study, the efficiency of localized surface plasmon resonance (LSPR) biosensor in detecting a single nucleotide mismatch in deoxyribonucleic acid is examined. The detection is based on the hybridization of a target DNA at 100 ng μL−1 with a complementary biotinylated probe as well as a partially complementary biotinylated with one nucleotide mismatch probe on a gold‐coated surface. Both probes are used at a concentration of 0.1 μm. The LSPR exhibited sensitivity by differentiating sample M+ from sample C+ through varying transmission intensities of 0.28 and 0.26 μA, respectively. Based on these findings, this approach demonstrates a great potential due to its ability to distinguish samples that differ with a single base pair, and its efficiency will be explored in the development of a point‐of‐care device as a simpler and cost‐effective approach for detection of various biologically and medically significant mutations such as antimicrobial resistance mutations. More work is underway to determine the robustness of the LSPR biosensor using the biotin–neutravidin approach.

Abstract GS03-09: Characterization and proposed therapeutic exploitation of fusion RNAs in metastatic breast cancers

Abstract BACKGROUND: Large-scale genomic studies such as The Cancer Genome Atlas (TCGA) and Pan-Cancer Analysis of Whole Genomes (PCAWG) show that the breast cancer (BrCa) genome is dominated by structural variation (SV) rather than single base pair mutations, producing a fertile environment for gene fusions. In this study, we implement a rigorous, expression-based approach to create a comprehensive landscape of fusion RNAs in metastatic breast cancer (MBC). We find fusion RNAs—many of which are novel involving known oncogenes—are surprisingly common in the advanced setting and credential their use as base-editing therapeutic targets. METHODS: Two retrospective cohorts of MBC RNA-sequencing data were analyzed— Dana-Farber Cancer Institute CCPM (n = 252 cases, 276 specimens), MichiganCSER (n = 171 cases, 190 specimens)—with a Fusion MetaCaller that integrates 5 unsupervised fusion-finding algorithms. Fusion RNAs identified in at least 2 callers (High-Confidence) and absent in RNA-seq from normal tissue (Cancer-Specific) were classified as HCCS-Fusions. Further removal of common artifactual fusions was performed using public databases. Expression of each HCCS-Fusion was quantified using a supervised method (FusionInspector)—expressed HCCS-Fusions were defined as having a Fusion Fragment Per Million (FFPM) value > 0.1 and at least 10% of read counts mapping to the fusion breakpoint versus the flanking 5’ or 3’ partners’ exons. Normalized gene-level expression abundances were calculated to correlate transcriptomic features (gene expression, PAM50) with fusion RNAs. Recurrent, potentially pathogenic fusion RNAs were annotated using OncoKB and outlier expressed fusions (Q3 FFPM + [1.5 X IQR]) were interrogated. RESULTS: The frequency of HCCS-Fusions differed between subtypes with basal BrCa harboring the most per tumor followed by Her2, LumB and LumA—with a median HCCS-Fusion count of 13, 12, 7, 4 respectively. 64.5% of cases harbored a HCCS-Fusion in an OncoKB cancer-related gene. The most recurrent fusions involving a cancer-related gene were 5’ ESR1 fusions (14 cases)—all with in-frame breakpoints near exon6/7, disrupting ESR1’s ligand binding domain. 13 of 14 ESR1 fusions were called in Luminal B (LumB) metastases (Fisher’s exact enrichment p < 0.005 vs other subtypes)—defining an ESR1 fusion frequency of 6.5% in LumB disease. Beyond ESR1, we identify recurrent, low-frequency (2-4 cases) in-frame kinase fusions involving FGFR2, ADK, TLK2, PRKCA, BRAF, CHKA, CSNK1D, NEK11, TNIK—some potentially targetable with FDA-approved small molecule inhibitors—as well as recurrent, predicted loss-of-function fusion RNAs in NF1, MSI2, USP32, PTEN, and CDH1. Lastly, 33.6% of cases harbored at least one outlier expressed fusion RNA; including highly expressed in-frame fusions involving known BrCa mediators such as ERBB2, BRCA1, ARID1B, RPS6KB1/2, PIK3R3, AXIN1, TGFB1/2, FOXP1, PAK1, and CREBBP. CONCLUSIONS: Taken together, these results demonstrate that fusion RNAs in MBC—some recurrent, many highly expressed and unique to individual tumors—are common. We create the most comprehensive catalog of ESR1 fusions in MBC, better define their frequency, discover their enrichment in LumB-like tumors, and will discuss clinicopathologic and transcriptomic features associated with ESR1 fusion positive disease. We identify druggable fusions that would likely be missed by current testing standards, find recurrent loss-of-function fusion RNAs, and show that over one-third of metastatic cases harbor at least one outlier expressed fusion—many of which involve BrCa-related genes. In summary, we propose that fusion RNAs are a driving and perhaps overlooked mechanism of tumor evolution in therapy-resistant disease and postulate fusion transcripts present a compelling therapeutic opportunity in MBC. Preliminary data targeting fusion RNA breakpoints using a novel RNA base-editing approach will be discussed. Citation Format: Nolan Priedigkeit, Alinés Lebrón-Torres, Janny Liao, Jean-Baptiste Alberge, Stefania Morganti, Jakob Weiss, Jorge Gomez Tejeda Zanudo, Albert Grinshpun, Melissa Hughes, Karla Helvie, Kerry Sendrick, Kyleen Nguyen, Sarah Strauss, Janet Files, Maxwell Lloyd, Nikhil Wagle, Chip Stewart, Eric Winer, Bruce Johnson, Yvonne Li, Rinath Jeselsohn, Sara Tolaney, Daniel Abravanel, Nancy Lin, Heather Parsons, Gad Getz, Steffi Oesterreich, Adrian Lee, Todd Golub. Characterization and proposed therapeutic exploitation of fusion RNAs in metastatic breast cancers [abstract]. In: Proceedings of the 2023 San Antonio Breast Cancer Symposium; 2023 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2024;84(9 Suppl):Abstract nr GS03-09.

Pathophysiology of Carnitine Palmitoyltransferase 2 (CPT2) Deficiency in the Dahl Salt-sensitive R at

Carnitine palmitoyltransferase 2 (CPT2) is an essential enzyme of the carnitine shuttle, which facilitates fatty acid -oxidation in mitochondria. This study aims to characterize the various effects of CPT2 deficiency on the Dahl salt-sensitive (SS) rat metabolome and its concomitant influence on renal homeostasis. CRISPR/Cas9 was used to insert a single base pair into the SS Cpt2 gene, causing a frameshift mutation found to be lethal in homozygotes (SSCpt2−/−). Western blotting and qPCR analysis of renal cortex tissue confirmed a reduced expression of CPT2 in heterozygous (SSCpt2+/-) rats compared to wild-type (WT) littermates; therefore, male and female heterozygotes were used for all experiments. Lipid metabolomics of renal cortex samples collected from SSCpt2+/- animals at 8 weeks old on a baseline 0.4% NaCl diet revealed a significant upregulation in inflammatory lipids such as 11-HETE (p<0.01) and 12-HHT (p<0.05) when compared to WT littermates. Urine analysis from the SSCpt2+/- group demonstrated a severe reduction in acetoacetic acid production (p<0.001), indicating ketogenesis is likely impacted by the loss of CPT2 in -oxidation. Additionally, TCA cycle metabolites fumaric acid (p<0.05) and oxaloacetic acid (p<0.01) are significantly downregulated in urine from SSCpt2+/- animals. Radio telemeters were implanted into the rats at 8 weeks of age to continuously record mean arterial pressure (MAP) and monitor the development of SS hypertension over the course of several dietary treatment protocols. Three weeks of either a 4% NaCl high-salt (HS) diet or a HS + ketogenic (HSK) diet with equivalent caloric content (4:1 caloric ratio, fat: carbohydrate) were administered to the animals. No difference in MAP elevation occurred between groups after the HS challenge alone, but SSCpt2+/- animals on the HSK diet had a significantly attenuated development of hypertension (p=0.03). Further metabolic analysis after the HSK diet illustrated that both groups had similar increases in long-chain fatty acids, liver enzymes, and cholesterol levels with a decrease in short-chain fatty acids. These experiments show that the partial loss of CPT2 in the Dahl SS rat leads to dysregulation across multiple metabolic pathways. Moreover, a ketogenic diet appears to protect heterozygous animals from the development of salt-sensitive hypertension. Funding: NIH R35 HL135749 and VA I01 BX004024. This is the full abstract presented at the American Physiology Summit 2024 meeting and is only available in HTML format. There are no additional versions or additional content available for this abstract. Physiology was not involved in the peer review process.

First report of Perkinsus marinus occurrence associated with wild Pacific oysters Crassostrea gigas from the west coast of Korea

This study reports the occurrence of Perkinsus marinus associated with wild Pacific oyster (Crassostrea gigas) specimens collected along the west coast of Korea. Confirmation of P. marinus presence was achieved by conventional PCR using World Organization of Animal Health (WOAH)-recommended primers that specifically targeted regions of the rDNA locus (ITS1, 5.8S, and ITS2). Sequencing of 10 samples revealed two distinct sequences differing by a single base pair, indicating potential haplotype variability. One sequence closely resembled the P. marinus strain found in Maryland, USA, whereas the other exhibited divergence, indicative of species diversity in the Korean strain, as was evident from the haplotype network analysis. Further validation involved the Ray’s Fluid Thioglycollate Medium (RFTM) assay, which initially yielded inconclusive results, possibly due to low infection intensity. Subsequently, RFTM and 2 M NaOH assays conducted on the isolates in the present study, cultured P. marinus cells in standard DMEM/F12 medium, and a positive P. marinus strain (ATCC 50509), revealed characteristic hypnospores of P. marinus upon Lugol’s iodine staining. These comprehensive investigations underscore the conclusive confirmation of P. marinus in Korean waters and mark a significant milestone in our understanding of the distribution and characteristics of this parasite in previously unreported regions.