Single Nucleotide Change Research Articles

BIOM-61. CSF CTDNA ANALYSIS OF PATIENTS WITH DIFFUSE GLIOMAS REVEALS DIAGNOSTIC MOLECULAR ALTERATIONS

Abstract There is a critical need for better strategies to diagnose and monitor patients with diffuse gliomas. Liquid biopsies are less invasive than tissue biopsies and are amenable to serial collection at multiple timepoints. However, analysis of circulating tumor DNA (ctDNA) in cerebrospinal fluid (CSF) in patients with diffuse gliomas remains underutilized in clinical practice due to technical challenges and a scarcity of studies demonstrating clinical utility. We evaluated ctDNA in 49 CSF samples from patients with various types of diffuse gliomas with a next generation sequencing (NGS) panel interrogating 613 cancer-related genes. Cell-free DNA concentration obtained from CSF (0.4-5mL) ranged from 0.009 to 60ng/µL. Input DNA concentration ranged from 0.12-30ng. Mutations were detected in 29/49 (59.2%) samples including 22/49 (44.9%) with single nucleotide changes, 14/49 (28.6%) with copy number changes, and 1 fusion (2.0%). The majority of samples without mutations detected (n=20) had suboptimal DNA input (<30ng): 17/20 (85%). We identified diagnostic alterations in CSF including mutations in IDH1 (3/5 IDH-mutant gliomas), TERTp and EGFR amplification (7/34 glioblastomas), H3-3A (4/6 pediatric high-grade gliomas). Some of these mutations, in combination with other clinical parameters, allowed us to identify specific CNS tumor types based on the 2021 WHO classification of CNS tumors, resulting in the reclassification of 2 patients. Matched tumor tissue sequencing results were available for 23/49 (46.9%) patients and the results showed good concordance with CSF ctDNA analysis. In summary, our results suggest that the NGS analysis of CSF-ctDNA can provide clinically relevant information in patients with diffuse gliomas.

Cytosine base editors with increased PAM and deaminase motif flexibility for gene editing in zebrafish

Cytosine base editing is a powerful tool for making precise single nucleotide changes in cells and model organisms like zebrafish, which are valuable for studying human diseases. However, current base editors struggle to edit cytosines in certain DNA contexts, particularly those with GC and CC pairs, limiting their use in modelling disease-related mutations. Here we show the development of zevoCDA1, an optimized cytosine base editor for zebrafish that improves editing efficiency across various DNA contexts and reduces restrictions imposed by the protospacer adjacent motif. We also create zevoCDA1-198, a more precise editor with a narrower editing window of five nucleotides, minimizing off-target effects. Using these advanced tools, we successfully generate zebrafish models of diseases that were previously challenging to create due to sequence limitations. This work enhances the ability to introduce human pathogenic mutations in zebrafish, broadening the scope for genomic research with improved precision and efficiency.

Environment by environment interactions (ExE) differ across genetic backgrounds (ExExG).

While the terms "gene-by-gene interaction" (GxG) and "gene-by-environment interaction" (GxE) are widely recognized in the fields of quantitative and evolutionary genetics, "environment-byenvironment interaction" (ExE) is a term used less often. In this study, we find that environmentby-environment interactions are a meaningful driver of phenotypes, and moreover, that they differ across different genotypes (suggestive of ExExG). To support this conclusion, we analyzed a large dataset of roughly 1,000 mutant yeast strains with varying degrees of resistance to different antifungal drugs. Our findings reveal that the effectiveness of a drug combination, relative to single drugs, often differs across drug resistant mutants. Remarkably, even mutants that differ by only a single nucleotide change can have dramatically different drug × drug (ExE) interactions. We also introduce a new framework that more accurately predicts the direction and magnitude of ExE interactions for some mutants. Understanding how ExE interactions change across genotypes (ExExG) is crucial not only for modeling the evolution of pathogenic microbes, but also for enhancing our knowledge of the underlying cell biology and the sources of phenotypic variance within populations. While the significance of ExExG interactions has been overlooked in evolutionary and population genetics, these fields and others stand to benefit from understanding how these interactions shape the complex behavior of living systems.

Non-canonical start codons confer context-dependent advantages in carbohydrate utilization for commensal E. coli in the murine gut.

Resource competition is a driver of gut microbiota composition. Bacteria can outcompete metabolically similar rivals through the limitation of shared growth-fuelling nutrients. The mechanisms underlying this remain unclear for bacteria with identical sets of metabolic genes. Here we analysed the lactose utilization operon in the murine commensal Escherichia coli 8178. Using in vitro and in vivo approaches, we showed that translation of the lactose utilization repressor gene lacI from its native non-canonical GTG start codon increases the basal expression of the lactose utilization cluster, enhancing adaptation to lactose consumption. Consequently, a strain carrying the wild type lacI GTG start codon outperformed the lacI ATG start codon mutant in the mouse intestine. This advantage was attenuated upon limiting host lactose intake through diet shift or altering the mutant frequency, emphasizing the context-dependent effect of a single nucleotide change on the bacterial fitness of a common member of the gut microbiota. Coupled with a genomic analysis highlighting the selection of non-ATG start codons in sugar utilization regulator genes across the Enterobacteriaceae family, our data exposed an unsuspected function of non-canonical start codons in metabolic competition.

Hyperactive Nickase Activity Improves Adenine Base Editing.

Base editing technologies enable programmable single-nucleotide changes in target DNA without double-stranded DNA breaks. Adenine base editors (ABEs) allow precise conversion of adenine (A) to guanine (G). However, limited availability of optimized deaminases as well as their variable efficiencies across different target sequences can limit the ability of ABEs to achieve effective adenine editing. Here, we explored the use of a TurboCas9 nickase in an ABE to improve its genome editing activity. The resulting TurboABE exhibits amplified editing efficiency on a variety of adenine target sites without increasing off-target editing in DNA and RNA. An interesting feature of TurboABE is its ability to significantly improve the editing frequency at bases with normally inefficient editing rates in the editing window of each target DNA. Development of improved ABEs provides new possibilities for precise genetic modification of genes in living cells.

Enhanced catabolism of glycine betaine and derivatives provides improved osmotic stress protection in Methylorubrum extorquens PA1.

Integration of metabolites into the overall metabolic network of a cell requires careful coordination dependent upon the ultimate usage of the metabolite. Different stoichiometric needs, and thus pathway fluxes, must exist for compounds destined for diverse uses, such as carbon sources, nitrogen sources, or stress-protective agents. Herein, we expand upon our previous work that highlighted the nature of glycine betaine (GB) metabolism in Methylobacteria to examine the utilization of GB-derivative compounds dimethylglycine (DMG) and sarcosine into Methylorubrum extorquens in different metabolic capacities, including as sole nitrogen and/or carbon sources. We isolated gain-of-function mutations that allowed M. extorquens PA1 to utilize dimethylglycine as a carbon source and dimethylglycine and sarcosine as nitrogen source. Characterization of mutants demonstrated selection for variants of the AraC-like regulator Mext_3735 that confer constitutive expression of the GB metabolic gene cluster, allowing direct utilization of the downstream GB derivatives. Finally, among the distinct isolates examined, we found that catabolism of the osmoprotectant used for selection (GB or dimethylglycine) enhanced osmotic stress resistance provided in the presence of that particular osmolyte. Thus, access to the carbon and nitrogen and osmoprotective effects of GB and DMG are made readily accessible through adaptive mutations. In M. extorquens PA1, the limitations to exploiting this group of compounds appear to exist predominantly at the levels of gene regulation and functional activity, rather than being constrained by transport or toxicity.IMPORTANCEOsmotic stress is a common challenge for bacteria colonizing the phyllosphere, where glycine betaine (GB) can be found as a prevalent osmoprotectant. Though Methylorubrum extorquens PA1 cannot use GB or its demethylation products, dimethylglycine (DMG) and sarcosine, as a sole carbon source, utilization is highly selectable via single nucleotide changes for both GB and DMG growth. The innate inability to use these compounds is due to limited flux through steps in the pathway and regulatory constraints. Herein, the characterization of the transcriptional regulator, Mext_3735 (GbdR), expands our understanding of the various roles in which GB derivatives can be used in M. extorquens PA1. Interestingly, increased catabolism of GB and derivatives does not interfere with, but rather improves, the ability of cells to thrive under increased salt stress conditions, suggesting that metabolic flux improves stress tolerance rather than providing a distinct tension between uses.

Structural and biophysical analysis of cytochrome P450 2C9*14 and *27 variants in complex with losartan

The human cytochrome P450 (CYP) 1, 2 and 3 families of enzymes are responsible for the biotransformation of a majority of the currently available pharmaceutical drugs. The highly polymorphic CYP2C9 predominantly metabolizes many drugs including anticoagulant S-warfarin, anti-hypertensive losartan, anti-diabetic tolbutamide, analgesic ibuprofen, etc. There are >80 single nucleotide changes identified in CYP2C9, many of which significantly alter the clearance of important drugs. Here we report the structural and biophysical analysis of two polymorphic variants, CYP2C9*14 (Arg125His) and CYP2C9*27 (Arg150Leu) complexed with losartan. The X-ray crystal structures of the CYP2C9*14 and *27 illustrate the binding of two losartan molecules, one in the active site near heme and another on the periphery. Both losartan molecules are bound in an identical conformation to that observed in the previously solved CYP2C9 wild-type complex, however, the number of losartan differs from the wild-type structure, which showed binding of three molecules. Additionally, isothermal titration calorimetry experiments reveal a lower binding affinity of losartan with *14 and *27 variants when compared to the wild-type. Overall, the results provide new insights into the effects of these genetic polymorphisms and suggests a possible mechanism contributing to reduced metabolic activity in patients carrying these alleles.

The 3'UTR Polymorphisms in the NLRP3 Gene Associated with the Risk of COPD and Their Putative Effects on the microRNA Mechanism.

Aims: Evaluating the association between a single nucleotide polymorphism in the 3' untranslated region (3'UTR) of the miRNA binding site of the NLRP3 gene and the occurrence and development of chronic obstructive pulmonary disease (COPD) and providing information to aid in the early detection and treatment of COPD. Materials and Methods: The regulatory single nuclear polymorphisms (SNPs) located in NLRP3 3'UTR were searched by using the dbSNP database and miRNA binding site prediction database. Meanwhile, samples from COPD patients and healthy controls in the same period were used for verification. The clinical baseline information of all subjects was collected, and the transcription level and protein expression level of NLRP3 and the expression level of inflammatory factors downstream of NLRP3 were detected. The effects of SNPs' single nucleotide changes on the transcription and expression of inflammatory factors were analyzed. Results: The study included 418 participants (249 in the COPD group and 169 in the control group). NLRP3 SNPs with miRNA binding sites include rs10754558 (G > C), rs1664774076 (ATAT > del), and rs1664775106 (C > G). Furthermore, two genotypes, GCG and GCA, were discovered to have a linkage mutation at 3'UTR 459-461. COPD susceptibility is tightly associated with the expression of the rs1664774076 del/del genotype, and the risk of COPD increased by 2.770 times (p = 0.003). Type 459-461 GCA was substantially related to the likelihood of developing COPD at various stages (p < 0.05). Except for rs10754558, all homozygous mutants increased NLRP3 mRNA and protein levels. NLRP3 had the greatest area under the receiver operating characteristic (ROC) curve for predicting the development and diagnosis of COPD when compared with its downstream inflammatory variables (AUC = 0.9291). Conclusions: The NLRP3 rs1664774076 del/del genotype is a COPD susceptibility gene, and the GCA genotype at 459-461 can be used as an early predictor of COPD exacerbation. The NLRP3 3'UTR polymorphism may alter the loss of miRNA binding sites, leading to an increase in NLRP3 expression. In the development of COPD, NLRP3 has a better diagnostic value than traditional inflammatory factors. The Clinical Trials Registration number Z: protocol KY01-2020-11-06.

A novel variant in the FLNB gene associated with spondylocarpotarsal synostosis syndrome.

Genetic disorders involved in skeleton system arise due to the disturbance in skeletal development, growth and homeostasis. Filamin B is an actin binding protein which is large dimeric protein which cross link actin cytoskeleton filaments into dynamic structure. A single nucleotide changes in the FLNB gene causes spondylocarpotarsal synostosis syndrome, a rare bone disorder due to which the fusion of carpels and tarsals synostosis occurred along with fused vertebrae. In the current study we investigated a family residing in north-western areas of Pakistan. The whole exome sequencing of proband was performed followed by Sanger sequencing of all family members of the subject to validate the variant segregation within the family. Bioinformatics tools were utilized to assess the pathogenicity of the variant. Whole Exome Sequencing revealed a novel variant (NM_001457: c.209C>T and p.Pro70Leu) in the FLNB gene which was homozygous missense mutation in the FLNB gene. The variant was further validated and visualized by Sanger sequencing and protein structure studies respectively as mentioned before. The findings have highlighted the importance of the molecular diagnosis in SCT (spondylocarpotarsal synostosis syndrome) for genetic risk counselling in consanguineous families.

Abstract PO5-25-09: rs2242652 POLYMORPHISM OF THE hTERT GENE IN WOMEN WITH BREAST CANCER

Abstract Introduction: Breast cancer is the most common neoplasm in women around the world. There are several risk factors for this cancer: sex, age, body mass index (BMI), hormone therapy, family history, genetic mutations, among several others. The vast majority of breast tumors are characterized as carcinomas and may originate from ducts or lobules. Immunohistochemistry is related to the expressions of the estrogen receptor (ER) and progesterone receptor (PR), as well as the expression of the human epidermal growth factor receptor 2 (HER2), cell proliferation or the absence of any expression from receptors on the tumor membrane. The telomeres protect the chromosomes natural endings from suffering DNA loss. Telomerase is a ribonuclear protein complex (RNA and protein) that neutralizes telomeres' shortening. Recurrent point mutations in the TERT telomerase promoter have been identified in several cancers, and the rs2242652 polymorphism is a single nucleotide change, and has been described as associated with breast cancer risk. Objective: The objective of this study is to analyze the TERT rs2242652 polymorphism in breast cancer women and to test the correlation of such data with the prognosis and various clinical variables. Methods: The study was conducted as an observational clinical study of women with breast cancer. A blood sample was collected from the patient and PCR was performed to extract the DNA and to analyze the telomerase polymorphisms in this population. To ensure a more descriptive analysis, the data was used in the construction of tables and graphs, that were then studied in a qualitative manner. The associations were analyzed through direct correlation tests and were followed by comparisons of categorization groups. Results and discussion: Fifty-four women presented polymorphism of the telomerase hTERT gene rs2242652. They were divided into two groups according to the A or G allele that confer the GA/AA and GG gene variants. It was found that the GA/AA group had a trend towards more aggressive tumors, presenting 33% triple negative and 22% of tumors with HER2 overexpression, different from the GG group: 19% triple negative and 10% HER2 overexpression. More aggressive T4 tumors, were more expressive in the GA/AA gene variant group representing 20% of the patients, whereas in the GG group they represented only 8%. Conclusion: The data found in this study on the telomerase TERT gene polymorphism rs2242652 with breast cancer refers to a tendency found in tumors that have a more severe prognosis, to associate with the gene variant linked to the A allele of this gene, however a larger sample number is required for statistical confirmation. Citation Format: Miriam Alvares, Elisa Mascarenhas, Joao Matos. rs2242652 POLYMORPHISM OF THE hTERT GENE IN WOMEN WITH BREAST CANCER [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 PO5-25-09.

Validation of Single-Nucleotide Mosaic Variants Through Droplet Digital PCR.

The detection, validation, and subsequent interpretation of potentially mosaic single-nucleotide variants (SNV) within next-generation sequencing data remains a challenge in both research and clinical laboratory settings. The ability to identify mosaic variants in high genome coverage sequencing data at levels of ≤1% underscores the necessity for developing guidelines and best practices to verify these variants orthogonally. Droplet digital PCR (ddPCR) has proven to be a powerful and precise method that allows for the determination of low-level variant fractions within a given sample. Herein we describe two precise ddPCR methods using either a fluorescent TaqMan hydrolysis probe approach or an EvaGreen fluorescent dye protocol. The TaqMan approach relies on two different fluorescent probes (FAM and HEX/VIC), each designed to amplify selectively only in the presence of a single nucleotide change denoting the variant or reference position. The fractional abundance is then calculated to determine the relative quantities of both alleles in the final sample. The EvaGreen protocol relies on two independent reactions with oligonucleotide primers designed with the single nucleotide change denoting the variant at the penultimate position of the primer. The relative amplification efficiency of both primer sets (reference and variant) can be compared to determine the mosaic level of a given variant. As the cost of high-coverage sequencing continues to decrease, the identification of potentially mosaic variants will also increase. The approaches outlined will allow clinicians and researchers a more precise determination of the true mosaic level of a given variant allowing them to better assess not only its potential pathogenicity but also its possible recurrence risk when offering genetic counseling to families. © 2024 Wiley Periodicals LLC. Basic Protocol: Droplet digital PCR (ddPCR) with TaqMan hydrolysis probes Alternate Protocol: EvaGreen oligonucleotide-specific ddPCR.

Convergent genetic adaptation of Escherichia coli in minimal media leads to pleiotropic divergence.

Adaptation in an environment can either be beneficial, neutral or disadvantageous in another. To test the genetic basis of pleiotropic behaviour, we evolved six lines of E. coli independently in environments where glucose and galactose were the sole carbon sources, for 300 generations. All six lines in each environment exhibit convergent adaptation in the environment in which they were evolved. However, pleiotropic behaviour was observed in several environmental contexts, including other carbon environments. Genome sequencing reveals that mutations in global regulators rpoB and rpoC cause this pleiotropy. We report three new alleles of the rpoB gene, and one new allele of the rpoC gene. The novel rpoB alleles confer resistance to Rifampicin, and alter motility. Our results show how single nucleotide changes in the process of adaptation in minimal media can lead to wide-scale pleiotropy, resulting in changes in traits that are not under direct selection.

Optimized allele-specific silencing of the dominant-negative COL6A1 G293R substitution causing collagen VI-related dystrophy

Collagen VI-related dystrophies (COL6-RDs) are a group of severe, congenital-onset muscular dystrophies for which there is no effective causative treatment. Dominant-negative mutations are common in COL6A1, COL6A2, and COL6A3 genes, encoding the collagen α1, α2, and α3 (VI) chains. They act by incorporating into the hierarchical assembly of the three α (VI) chains and consequently produce a dysfunctional collagen VI extracellular matrix, while haploinsufficiency for any of the COL6 genes is not associated with disease. Hence, allele-specific transcript inactivation is a valid therapeutic strategy, although selectively targeting a pathogenic single nucleotide variant is challenging. Here, we develop a small interfering RNA (siRNA) that robustly, and in an allele-specific manner, silences a common glycine substitution (G293R) caused by a single nucleotide change in COL6A1 gene. By intentionally introducing an additional mismatch into the siRNA design, we achieved enhanced specificity toward the mutant allele. Treatment of patient-derived fibroblasts effectively reduced the levels of mutant transcripts while maintaining unaltered wild-type transcript levels, rescuing the secretion and assembly of collagen VI matrix by reducing the dominant-negative effect of mutant chains. Our findings establish a promising treatment approach for patients with the recurrent dominantly negative acting G293R glycine substitution.

Abstract A028: Subclonal lineage recording in glioblastoma using CRISPR base editing

Abstract Novel lineage tracing approaches are needed to undercover lineage relationships between diverse cells in complex tumors such as glioblastoma (GBM). Intratumoral heterogeneity in GBM drives treatment resistance and recurrence. This intratumoral heterogeneity is reflected by multiple cellular states seen within tumors and across patients. Unfortunately, little is known about their lineage relationships and shared origins. Here we’ve developed a novel lineage tracing technology to track the clonal and subclonal evolution in patient-derived GBM models. Our adenine base editor lineage tracing (ABELT) system uses a Cas9 adenine base editor (ABE) to induce heritable single-nucleotide changes in the 3’ untranslated region of endogenous transcripts. These heritable marks distinguish individual clonal populations and can be used to trace the relationship between progeny within each clone. We demonstrate our system's in vitro lineage tracing capacity and show that the ABELT system faithfully records lineage, which is recoverable with unmodified single-cell RNA sequencing approaches. We then track the relationships between evolving cellular states across multiple divisions in an evolving GBM model. Our ABELT technology enables the rapid engineering and lineage tracing of cancer models incompatible with previous approaches. This novel approach will enable our group and others to explore GBM intratumoral heterogeneity with subclonal resolution to determine transcriptional plasticity that ultimately leads to treatment resistance. Citation Format: Abigail C. Marshall, Aaron McKenna. Subclonal lineage recording in glioblastoma using CRISPR base editing [abstract]. In: Proceedings of the AACR Special Conference on Brain Cancer; 2023 Oct 19-22; Minneapolis, Minnesota. Philadelphia (PA): AACR; Cancer Res 2024;84(5 Suppl_1):Abstract nr A028.

Full-length sequence of the novel allele, HLA-B*58:01:40, identified in a Chinese individual.

HLA-B*58:01:40 differs from HLA-B*58:01:01 by a single nucleotide change in exon 3, 507 C- > T (codon 145.3 CGC- > CGT).

Genotype-by-environment interactions govern fitness changes associated with adaptive mutations in two-component response systems.

Introduction: Two-component response systems (TCRS) are the main mechanism by which prokaryotes acclimate to changing environments. These systems are composed of a membrane bound histidine kinase (HK) that senses external signals and a response regulator (RR) that activates transcription of response genes. Despite their known role in acclimation, little is known about the role TCRS play in environmental adaptation. Several experimental evolution studies have shown the acquisition of mutations in TCRS during adaptation, therefore here we set out to characterize the adaptive mechanism resulting from these mutations and evaluate whether single nucleotide changes in one gene could induce variable genotype-by-environment (GxE) interactions. Methods: To do this, we assessed fitness changes and differential gene expression for four adaptive mutations in cusS, the gene that encodes the HK CusS, acquired by Escherichia coli during silver adaptation. Results: Fitness assays showed that as the environment changed, each mutant displayed a unique fitness profile with greatest fitness in the original selection environment. RNAseq then indicated that, in ± silver nitrate, each mutant induces a primary response that upregulates cusS, its RR cusR, and constitutively expresses the target response genes cusCFBA. This then induces a secondary response via differential expression of genes regulated by the CusR through TCRS crosstalk. Finally, each mutant undergoes fitness tuning through unique tertiary responses that result in gene expression patterns specific for the genotype, the environment and optimized for the original selection conditions. Discussion: This three-step response shows that different mutations in a single gene leads to individualized phenotypes governed by unique GxE interactions that not only contribute to transcriptional divergence but also to phenotypic plasticity.

Green, yellow, or cyan? Introduction of color change mutations into a green thermostable fluorescent protein and characterization during an introduction to biochemistry lab course.

Green fluorescent protein has long been a favorite protein for demonstrating protein purification in the biochemistry lab course. The protein's vivid green color helps demonstrate to students the concept(s) behind affinity or ion exchange chromatography. We designed a series of introduction to biochemistry labs utilizing a thermostable green protein (TGP-E) engineered to have unusually high thermostability. This protein allows students to proceed through purification and characterization without the need to keep protein samples on ice. The 5-week lab series begins with an introduction to molecular biology techniques during weeks 1 and 2, where site-directed mutagenesis is used introduce, a single nucleotide change that shifts the fluorescent spectra of TGP-E to either cyan (CTP-E) or yellow (YTP-E). Students identify successful mutagenesis reaction by the color of a small expression sample after induction with IPTG. Next, students purify either the TGP-E (control-typically one group volunteers), YTP-E, or CTP-E protein as a 1-week lab. During the following week's lab, students run SDS-PAGE to verify protein purity, bicinchoninic acid assay to quantify protein yield, and absorbance and fluorescence spectra to characterize their protein's fluorescent character. The final lab in the series investigates the thermostability of YTP-E and CTP-E compared with TGP-E using a fluorescence plate reader. This 5-week series of experiments provide students with experience in several key biochemistry techniques and allows the students to compare properties of mutations. At the end of the course, the students will write a research report and give a short presentation over their results.

Characterization of the novel HLA-B*35:08:17 allele.

HLA-B*35:08:17 differs from B*35:08:01:01 by a single nucleotide change in codon 105 (CCC > CCG).

Efficient detection of single nucleotide variants in targeted genomic loci.

Single nucleotide variants (SNVs), including single nucleotide polymorphisms, are often associated with morphological and/or physiological abnormalities in various organisms. Targeted genomic DNA can be amplified and directly sequenced to detect these mutations, but this method is relatively time consuming and expensive. We recently established the heteroduplex mobility assay to detect genetic mutations as an easy, low-cost method in genome editing, but detecting such small genetic differences remains difficult. Here, we developed a new, simple method to detect single nucleotide changes in the zebrafish genome by polymerase chain reaction (PCR) and electrophoresis. We first designed a specific single stranded DNA with four tandem guanine nucleotides inserted beside the mutation site, called guanine-inserted primer (GIP). When reannealing, hybridized complexes of GIP and PCR amplicons with or without 1-bp-mutated alleles form different bulge structures, presumably leading to different mobilities on a polyacrylamide gel. This GIP-interacting mobility assay is easy to use; therefore, it could contribute to the detection of SNVs in any organism.

Adaptation of SARS-CoV-2 to ACE2H353K mice reveals new spike residues that drive mouse infection.

The Coronavirus Disease 2019 (COVID-19) pandemic continues to cause extraordinary loss of life and economic damage. Animal models of severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) infection are needed to better understand disease pathogenesis and evaluate preventive measures and therapies. While mice are widely used to model human disease, mouse angiotensin converting enzyme 2 (ACE2) does not bind the ancestral SARS-CoV-2 spike protein to mediate viral entry. To overcome this limitation, we "humanized" mouse Ace2 using CRISPR gene editing to introduce a single amino acid substitution, H353K, predicted to facilitate S protein binding. While H353K knockin Ace2 (mACE2H353K) mice supported SARS-CoV-2 infection and replication, they exhibited minimal disease manifestations. Following 30 serial passages of ancestral SARS-CoV-2 in mACE2H353K mice, we generated and cloned a more virulent virus. A single isolate (SARS2MA-H353K) was prepared for detailed studies. In 7-11-month-old mACE2H353K mice, a 104 PFU inocula resulted in diffuse alveolar disease manifested as edema, hyaline membrane formation, and interstitial cellular infiltration/thickening. Unexpectedly, the mouse-adapted virus also infected standard BALB/c and C57BL/6 mice and caused severe disease. The mouse-adapted virus acquired five new missense mutations including two in spike (K417E, Q493K), one each in nsp4, nsp9, and M and a single nucleotide change in the 5' untranslated region. The Q493K spike mutation arose early in serial passage and is predicted to provide affinity-enhancing molecular interactions with mACE2 and further increase the stability and affinity to the receptor. This new model and mouse-adapted virus will be useful to evaluate COVID-19 disease and prophylactic and therapeutic interventions.IMPORTANCEWe developed a new mouse model with a humanized angiotensin converting enzyme 2 (ACE2) locus that preserves native regulatory elements. A single point mutation in mouse ACE2 (H353K) was sufficient to confer in vivo infection with ancestral severe acute respiratory syndrome-coronavirus-2 virus. Through in vivo serial passage, a virulent mouse-adapted strain was obtained. In aged mACE2H353K mice, the mouse-adapted strain caused diffuse alveolar disease. The mouse-adapted virus also infected standard BALB/c and C57BL/6 mice, causing severe disease. The mouse-adapted virus acquired five new missense mutations including two in spike (K417E, Q493K), one each in nsp4, nsp9, and M and a single nucleotide change in the 5' untranslated region. The Q493K spike mutation arose early in serial passage and is predicted to provide affinity-enhancing molecular interactions with mACE2 and further increase the stability and affinity to the receptor. This new model and mouse-adapted virus will be useful to evaluate COVID-19 disease and prophylactic and therapeutic interventions.