Engineering complex synthetic transcriptional programs with CRISPR RNA scaffolds.

Eukaryotic cells execute complex transcriptional programs in which specific loci throughout the genome are regulated in distinct ways by targeted regulatory assemblies. We have applied this principle to generate synthetic CRISPR-based transcriptional programs in yeast and human cells. By extending guide RNAs to include effector protein recruitment sites, we construct modular scaffold RNAs that encode both target locus and regulatory action. Sets of scaffold RNAs can be used to generate synthetic multigene transcriptional programs in which some genes are activated and others are repressed. We apply this approach to flexibly redirect flux through a complex branched metabolic pathway in yeast. Moreover, these programs can be executed by inducing expression of the dCas9 protein, which acts as a single master regulatory control point. CRISPR-associated RNA scaffolds provide a powerful way to construct synthetic gene expression programs for a wide range of applications, including rewiring cell fates or engineering metabolic pathways.

Elastomeric electrolytes for high-energy solid-state lithium batteries.

The use of lithium metal anodes in solid-state batteries has emerged as one of the most promising technologies for replacing conventional lithium-ion batteries1,2. Solid-state electrolytes are a key enabling technology for the safe operation of lithium metal batteries as they suppress the uncontrolled growth of lithium dendrites. However, the mechanical properties and electrochemical performance of current solid-state electrolytes do not meet the requirements for practical applications of lithium metal batteries. Here we report a class of elastomeric solid-state electrolytes with a three-dimensional interconnected plastic crystal phase. The elastomeric electrolytes show a combination of mechanical robustness, high ionic conductivity, low interfacial resistance and high lithium-ion transference number. The in situ-formed elastomer electrolyte on copper foils accommodates volume changes for prolonged lithium plating and stripping processes with a Coulombic efficiency of 100.0 per cent. Moreover, the elastomer electrolytes enable stable operation of the full cells under constrained conditions of a limited lithium source, a thin electrolyte and a high-loading LiNi0.83Mn0.06Co0.11O2 cathode at a high voltage of 4.5 volts at ambient temperature, delivering a high specific energy exceeding 410 watt-hours per kilogram of electrode plus electrolyte. The elastomeric electrolyte system presents a powerful strategy for enabling stable operation of high-energy, solid-state lithium batteries.

Sequential application of anticancer drugs enhances cell death by rewiring apoptotic signaling networks.

Crosstalk and complexity within signaling pathways and their perturbation by oncogenes limit component-by-component approaches to understanding human disease. Network analysis of how normal and oncogenic signaling can be rewired by drugs may provide opportunities to target tumors with high specificity and efficacy. Using targeted inhibition of oncogenic signaling pathways, combined with DNA-damaging chemotherapy, we report that time-staggered EGFR inhibition, but not simultaneous coadministration, dramatically sensitizes a subset of triple-negative breast cancer cells to genotoxic drugs. Systems-level analysis-using high-density time-dependent measurements of signaling networks, gene expression profiles, and cell phenotypic responses in combination with mathematical modeling-revealed an approach for altering the intrinsic state of the cell through dynamic rewiring of oncogenic signaling pathways. This process converts these cells to a less tumorigenic state that is more susceptible to DNA damage-induced cell death by reactivation of an extrinsic apoptotic pathway whose function is suppressed in the oncogene-addicted state.

STING controls nociception via type I interferon signalling in sensory neurons.

The innate immune regulator STING is a critical sensor of self- and pathogen-derived DNA. DNA sensing by STING leads to the induction of type-I interferons (IFN-I) and other cytokines, which promote immune-cell-mediated eradication of pathogens and neoplastic cells1,2. STING is also a robust driver of antitumour immunity, which has led to the development of STING activators and small-molecule agonists as adjuvants for cancer immunotherapy3. Pain, transmitted by peripheral nociceptive sensory neurons (nociceptors), also aids in host defence by alerting organisms to the presence of potentially damaging stimuli, including pathogens and cancer cells4,5. Here we demonstrate that STING is a critical regulator of nociception through IFN-I signalling in peripheral nociceptors. We show that mice lacking STING or IFN-I signalling exhibit hypersensitivity to nociceptive stimuli and heightened nociceptor excitability. Conversely, intrathecal activation of STING produces robust antinociception in mice and non-human primates. STING-mediated antinociception is governed by IFN-Is, which rapidly suppress excitability of mouse, monkey and human nociceptors. Our findings establish the STING-IFN-I signalling axis as a critical regulator of physiological nociception and a promising new target for treating chronic pain.

Functional genomic screens with death rate analyses reveal mechanisms of drug action.

A common approach for understanding how drugs induce their therapeutic effects is to identify the genetic determinants of drug sensitivity. Because 'chemo-genetic profiles' are performed in a pooled format, inference of gene function is subject to several confounding influences related to variation in growth rates between clones. In this study, we developed Method for Evaluating Death Using a Simulation-assisted Approach (MEDUSA), which uses time-resolved measurements, along with model-driven constraints, to reveal the combination of growth and death rates that generated the observed drug response. MEDUSA is uniquely effective at identifying death regulatory genes. We apply MEDUSA to characterize DNA damage-induced lethality in the presence and absence of p53. Loss of p53 switches the mechanism of DNA damage-induced death from apoptosis to a non-apoptotic death that requires high respiration. These findings demonstrate the utility of MEDUSA both for determining the genetic dependencies of lethality and for revealing opportunities to potentiate chemo-efficacy in a cancer-specific manner.

Author Correction: Elpistostege and the origin of the vertebrate hand.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Elpistostege and the origin of the vertebrate hand.

The evolution of fishes to tetrapods (four-limbed vertebrates) was one of the most important transformations in vertebrate evolution. Hypotheses of tetrapod origins rely heavily on the anatomy of a few tetrapod-like fish fossils from the Middle and Late Devonian period (393-359 million years ago)1. These taxa-known as elpistostegalians-include Panderichthys2, Elpistostege3,4 and Tiktaalik1,5, none of which has yet revealed the complete skeletal anatomy of the pectoral fin. Here we report a 1.57-metre-long articulated specimen of Elpistostege watsoni from the Upper Devonian period of Canada, which represents-to our knowledge-the most complete elpistostegalian yet found. High-energy computed tomography reveals that the skeleton of the pectoral fin has four proximodistal rows of radials (two of which include branched carpals) as well as two distal rows that are organized as digits and putative digits. Despite this skeletal pattern (which represents the most tetrapod-like arrangement of bones found in a pectoral fin to date), the fin retains lepidotrichia (fin rays) distal to the radials. We suggest that the vertebrate hand arose primarily from a skeletal pattern buried within the fairly typical aquatic pectoral fin of elpistostegalians. Elpistostege is potentially the sister taxon of all other tetrapods, and its appendages further blur the line between fish and land vertebrates.

Cells of the adult human heart.

Cardiovascular disease is the leading cause of death worldwide. Advanced insights into disease mechanisms and therapeutic strategies require a deeper understanding of the molecular processes involved in the healthy heart. Knowledge of the full repertoire of cardiac cells and their gene expression profiles is a fundamental first step in this endeavour. Here, using state-of-the-art analyses of large-scale single-cell and single-nucleus transcriptomes, we characterize six anatomical adult heart regions. Our results highlight the cellular heterogeneity of cardiomyocytes, pericytes and fibroblasts, and reveal distinct atrial and ventricular subsets of cells with diverse developmental origins and specialized properties. We define the complexity of the cardiac vasculature and its changes along the arterio-venous axis. In the immune compartment, we identify cardiac-resident macrophages with inflammatory and protective transcriptional signatures. Furthermore, analyses of cell-to-cell interactions highlight different networks of macrophages, fibroblasts and cardiomyocytes between atria and ventricles that are distinct from those of skeletal muscle. Our human cardiac cell atlas improves our understanding of the human heart and provides a valuable reference for future studies.

A periplasmic phospholipase that maintains outer membrane lipid asymmetry in Pseudomonas aeruginosa.

The outer membrane of Gram-negative bacteria is unique in both structure and function. The surface-exposed outer leaflet is composed of lipopolysaccharide, while the inner leaflet is composed of glycerophospholipids. This lipid asymmetry creates mechanical strength, lowers membrane permeability, and is necessary for virulence in many pathogens. Glycerophospholipids that mislocalize to the outer leaflet are removed by the Mla pathway, which consists of the outer membrane channel MlaA, the periplasmic lipid carrier MlaC, and the inner membrane transporter MlaBDEF. The opportunistic pathogen Pseudomonas aeruginosa has two proteins of the MlaA family: PA2800 and PA3239. Here, we show that PA2800 is part of a canonical Mla pathway, while PA3239 functions with the putative lipase PA3238. While loss of either pathway individually has little to no effect on outer membrane integrity, loss of both pathways weakens the outer membrane permeability barrier and increases production of the secondary metabolite pyocyanin. We propose that mislocalized glycerophospholipids are removed from the outer leaflet by PA3239 (renamed MlaZ), transferred to PA3238 (renamed MlaY), and degraded. This pathway streamlines recycling of glycerophospholipid degradation products by removing glycerophospholipids from the outer leaflet prior to degradation.

Blockade of trans PD-L1 interaction with CD80 augments antitumor immunity.

PD-L1 has two receptors: PD-1 and CD80. Previous reports assumed that PD-L1 and CD80 interacted in trans, but recent reports showed that only cis PD-L1/CD80 interactions existed, and prevention of cis PD-L1/CD80 interactions on antigen-presenting cells (APCs) reduced antitumor immunity via augmenting PD-L1/PD-1 and CD80/CTLA4 interactions between T and APCs. Here, using tumor-bearing mice capable of cis and trans or trans only PD-L1/CD80 interactions, we show that trans PD-L1/CD80 interactions do exist between tumor and T cells, and the effects of trans PD-L1/CD80 interactions require tumor cell expression of MHC-I and T cell expression of CD28. The blockade of PD-L1/CD80 interactions in mice with both cis and trans interactions or with only trans interactions augments antitumor immunity by expanding IFN-γ-producing CD8+ T cells and IFN-γ-dependent NOS2-expressing tumor-associated macrophages. Our studies indicate that although cis and trans PD-L1/CD80 interactions may have opposite effects on antitumor immunity, the net effect of blocking PD-L1/CD80 interactions in vivo augments CD8+ T cell-mediated antitumor immunity.

Dual thermal ecotypes coexist within a nearly genetically identical population of the unicellular marine cyanobacterium Synechococcus.

The extent and ecological significance of intraspecific functional diversity within marine microbial populations is still poorly understood, and it remains unclear if such strain-level microdiversity will affect fitness and persistence in a rapidly changing ocean environment. In this study, we cultured 11 sympatric strains of the ubiquitous marine picocyanobacterium Synechococcus isolated from a Narragansett Bay (RI) phytoplankton community thermal selection experiment. Thermal performance curves revealed selection at cool and warm temperatures had subdivided the initial population into thermotypes with pronounced differences in maximum growth temperatures. Curiously, the genomes of all 11 isolates were almost identical (average nucleotide identities of >99.99%, with >99% of the genome aligning) and no differences in gene content or single nucleotide variants were associated with either cool or warm temperature phenotypes. Despite a very high level of genomic similarity, sequenced epigenomes for two strains showed differences in methylation on genes associated with photosynthesis. These corresponded to measured differences in photophysiology, suggesting a potential pathway for future mechanistic research into thermal microdiversity. Our study demonstrates that present-day marine microbial populations can harbor cryptic but environmentally relevant thermotypes which may increase their resilience to future rising temperatures.

Genetic subgroups inform on pathobiology in adult and pediatric Burkitt lymphoma.

Burkitt lymphoma (BL) accounts for most pediatric non-Hodgkin lymphomas, being less common but significantly more lethal when diagnosed in adults. Much of the knowledge of the genetics of BL thus far has originated from the study of pediatric BL (pBL), leaving its relationship to adult BL (aBL) and other adult lymphomas not fully explored. We sought to more thoroughly identify the somatic changes that underlie lymphomagenesis in aBL and any molecular features that associate with clinical disparities within and between pBL and aBL. Through comprehensive whole-genome sequencing of 230 BL and 295 diffuse large B-cell lymphoma (DLBCL) tumors, we identified additional significantly mutated genes, including more genetic features that associate with tumor Epstein-Barr virus status, and unraveled new distinct subgroupings within BL and DLBCL with 3 predominantly comprising BLs: DGG-BL (DDX3X, GNA13, and GNAI2), IC-BL (ID3 and CCND3), and Q53-BL (quiet TP53). Each BL subgroup is characterized by combinations of common driver and noncoding mutations caused by aberrant somatic hypermutation. The largest subgroups of BL cases, IC-BL and DGG-BL, are further characterized by distinct biological and gene expression differences. IC-BL and DGG-BL and their prototypical genetic features (ID3 and TP53) had significant associations with patient outcomes that were different among aBL and pBL cohorts. These findings highlight shared pathogenesis between aBL and pBL, and establish genetic subtypes within BL that serve to delineate tumors with distinct molecular features, providing a new framework for epidemiologic, diagnostic, and therapeutic strategies.

Threat management priorities for conserving Antarctic biodiversity.

Antarctic terrestrial biodiversity faces multiple threats, from invasive species to climate change. Yet no large-scale assessments of threat management strategies exist. Applying a structured participatory approach, we demonstrate that existing conservation efforts are insufficient in a changing world, estimating that 65% (at best 37%, at worst 97%) of native terrestrial taxa and land-associated seabirds are likely to decline by 2100 under current trajectories. Emperor penguins are identified as the most vulnerable taxon, followed by other seabirds and dry soil nematodes. We find that implementing 10 key threat management strategies in parallel, at an estimated present-day equivalent annual cost of US$23 million, could benefit up to 84% of Antarctic taxa. Climate change is identified as the most pervasive threat to Antarctic biodiversity and influencing global policy to effectively limit climate change is the most beneficial conservation strategy. However, minimising impacts of human activities and improved planning and management of new infrastructure projects are cost-effective and will help to minimise regional threats. Simultaneous global and regional efforts are critical to secure Antarctic biodiversity for future generations.

Regulation and function of transposable elements in cancer genomes.

Over half of human genomic DNA is composed of repetitive sequences generated throughout evolution by prolific mobile genetic parasites called transposable elements (TEs). Long disregarded as "junk" or "selfish" DNA, TEs are increasingly recognized as formative elements in genome evolution, wired intimately into the structure and function of the human genome. Advances in sequencing technologies and computational methods have ushered in an era of unprecedented insight into how TE activity impacts human biology in health and disease. Here we discuss the current views on how TEs have shaped the regulatory landscape of the human genome, how TE activity is implicated in human cancers, and how recent findings motivate novel strategies to leverage TE activity for improved cancer therapy. Given the crucial role of methodological advances in TE biology, we pair our conceptual discussions with an in-depth review of the inherent technical challenges in studying repeats, specifically related to structural variation, expression analyses, and chromatin regulation. Lastly, we provide a catalog of existing and emerging assays and bioinformatic software that altogether are enabling the most sophisticated and comprehensive investigations yet into the regulation and function of interspersed repeats in cancer genomes.

Evidence of a genetic background predisposing to complex regional pain syndrome type 1.

BACKGROUND: Complex regional pain syndrome type 1 (CRPS-1) is a rare, disabling and sometimes chronic disorder usually arising after a trauma. This exploratory study examined whether patients with chronic CRPS-1 have a different genetic profile compared with those who do not have the condition. METHODS: Exome sequencing was performed to seek altered non-synonymous SNP allele frequencies in a discovery cohort of well-characterised patients with chronic CRPS-1 (n=34) compared with population databases. Identified SNP alleles were confirmed by Sanger sequencing and sought in a replication cohort (n=50). Gene expression of peripheral blood macrophages was assessed. RESULTS: In the discovery cohort, the rare allele frequencies of four non-synonymous SNPs were statistically increased. The replication cohort confirmed this finding. In a chronic pain cohort, these alleles were not overexpressed. In total, 25 out of 84 (29.8%) patients with CRPS-1 expressed a rare allele. The SNPs were rs41289586 in ANO10, rs28360457 in P2RX7, rs1126930 in PRKAG1 and rs80308281 in SLC12A9. Males were more likely than females to have a rare SNP allele, 8 out of 14 (57.1%) vs 17 out of 70 (24.3%) (Fisher's p=0.023). ANO10, P2RX7, PRKAG1 and SLC12A9 were all expressed in macrophages from healthy human controls. CONCLUSION: A single SNP in each of the genes ANO10, P2RX7, PRKAG1 and SLC12A9 was associated with developing chronic CRPS-1, with more males than females expressing these rare alleles. Our work suggests the possibility that a permissive genetic background is an important factor in the development of CRPS-1.

Persistent transcriptional changes in cardiac adaptive immune cells following myocardial infarction: New evidence from the re-analysis of publicly available single cell and nuclei RNA-sequencing data sets.

INTRODUCTION: Chronic immunopathology contributes to the development of heart failure after a myocardial infarction. Both T and B cells of the adaptive immune system are present in the myocardium and have been suggested to be involved in post-MI immunopathology. METHODS: We analyzed the B and T cell populations isolated from previously published single cell RNA-sequencing data sets (PMID: 32130914, PMID: 35948637, PMID: 32971526 and PMID: 35926050), of the mouse and human heart, using differential expression analysis, functional enrichment analysis, gene regulatory inferences, and integration with autoimmune and cardiovascular GWAS. RESULTS: Already at baseline, mature effector B and T cells are present in the human and mouse heart, having increased activity in transcription factors maintaining tolerance (e.g. DEAF1, JDP2, SPI-B). Following MI, T cells upregulate pro-inflammatory transcript levels (e.g. Cd11, Gzmk, Prf1), while B cells upregulate activation markers (e.g. Il6, Il1rn, Ccl6) and collagen (e.g. Col5a2, Col4a1, Col1a2). Importantly, pro-inflammatory and fibrotic transcription factors (e.g. NFKB1, CREM, REL) remain active in T cells, while B cells maintain elevated activity in transcription factors related to immunoglobulin production (e.g. ERG, REL) in both mouse and human post-MI hearts. Notably, genes differentially expressed in post-MI T and B cells are associated with cardiovascular and autoimmune disease. CONCLUSION: These findings highlight the varied and time-dependent dynamic roles of post-MI T and B cells. They appear ready-to-go and are activated immediately after MI, thus participate in the acute wound healing response. However, they subsequently remain in a state of pro-inflammatory activation contributing to persistent immunopathology.

Preclinical characterization of a non-peptidomimetic HIV protease inhibitor with improved metabolic stability.

Protease inhibitors (PIs) remain an important component of antiretroviral therapy for the treatment of HIV-1 infection due to their high genetic barrier to resistance development. Nevertheless, the two most commonly prescribed HIV PIs, atazanavir and darunavir, still require co-administration with a pharmacokinetic boosting agent to maintain sufficient drug plasma levels which can lead to undesirable drug-drug interactions. Herein, we describe GS-9770, a novel investigational non-peptidomimetic HIV PI with unboosted once-daily oral dosing potential due to improvements in its metabolic stability and its pharmacokinetic properties in preclinical animal species. This compound demonstrates potent inhibitory activity and high on-target selectivity for recombinant HIV-1 protease versus other aspartic proteases tested. In cell culture, GS-9770 inhibits Gag polyprotein cleavage and shows nanomolar anti-HIV-1 potency in primary human cells permissive to HIV-1 infection and against a broad range of HIV subtypes. GS-9770 demonstrates an improved resistance profile against a panel of patient-derived HIV-1 isolates with resistance to atazanavir and darunavir. In resistance selection experiments, GS-9770 prevented the emergence of breakthrough HIV-1 variants at all fixed drug concentrations tested and required multiple protease substitutions to enable outgrowth of virus exposed to escalating concentrations of GS-9770. This compound also remained fully active against viruses resistant to drugs from other antiviral classes and showed no in vitro antagonism when combined pairwise with drugs from other antiretroviral classes. Collectively, these preclinical data identify GS-9770 as a potent, non-peptidomimetic once-daily oral HIV PI with potential to overcome the persistent requirement for pharmacological boosting with this class of antiretroviral agents.

A phase II trial of UGT1A1 genotype-guided FOLFIRI plus bevacizumab as first-line therapy for advanced, unresectable colorectal cancer.

BACKGROUND: FOLFIRI is a standard regimen for metastatic colorectal cancer (mCRC). We hypothesized that a pharmacogenomic-directed strategy where more efficient irinotecan metabolizers (UGT1A1 *1/*1 homozygotes and *1/*28 heterozygotes) receive higher-than-standard irinotecan doses would improve progression-free survival (PFS) compared to non-genotype selected historical controls with acceptable toxicity. METHODS: In this phase II multicenter study irinotecan dosing in first-line FOLFIRI and bevacizumab for mCRC was based on UGT1A1 genotype with *1/*1, *1/*28, and *28/*28 patients receiving 310 mg/m2, 260 mg/m2, and 180 mg/m2, respectively. Primary endpoint was PFS. Secondary endpoints were investigator and patient-reported adverse events, and estimation of overall survival (OS). RESULTS: One-hundred patients were enrolled with 91 evaluable for PFS and 83 evaluable for best response. Median PFS was 12.5 months (90% CI 10.9, 15.4), shorter than the anticipated alternative hypothesis of 14 months. PFS by genotype was 12.5 months (90% CI 10.9, 17.4) for *1/*1, 14.6 months (90% CI 11.8, 17.5) for *1/*28, and 6 months (90% CI 2.3, 7.7) for *28/28, respectively. OS was 24.5 months (90% CI 19.1, 30.7) and by genotype was 26.5 (90% CI 19.1, 32.9), 25.9 (90% CI 17.6, 37.7), and 13.4 (90% CI 2.3, 20.5) months for *1/*1, *1/*28, and *28/*28, respectively. G3/4 toxicity was similar between all subgroups, including diarrhea and neutropenia. CONCLUSIONS: A pharmacogenomic-directed irinotecan strategy improved PFS in the *1/*1 and *1/*28 genotypes with higher rates of neutropenia and similar rates of diarrhea compared to expected with standard FOLFIRI dosing. However, improvements in response rate and PFS were modest. This strategy should not change standard practice for mCRC patients in the first-line setting.

The role of alternative autophagy in cell viability and response to paclitaxel treatment in v-Ha-ras-transformed NIH 3T3 cells.

In confluent v-Ha-ras-transformed NIH 3T3 fibroblasts (Ras-NIH 3T3), LC3 downregulation may precede a decrease in canonical autophagy, thus contributing to cell survival. Herein, we aimed to investigate the role of alternative autophagy in the viability of long-term cultures of Ras-NIH 3T3 cells and their parental NIH 3T3 cells. As cell confluence increased with the culture period, the level of alternative autophagy, as assessed through Lamp2-Rab9 co-localization, gradually decreased in both cell lines. However, Ras-NIH 3T3 cells maintained higher levels of alternative autophagy than the parental cells did. Rab9 knockdown minimally affected NIH 3T3 cells while drastically reducing the viability of Ras-NIH 3T3 cells, which suggested that alternative autophagy plays a critical role in Ras-NIH 3T3 cells. In contrast, reactive oxygen species (ROS) production in Ras-NIH 3T3 cells was higher than that in NIH 3T3 cells during long-term culture. Moreover, NIH 3T3 cells exhibited a continual decrease in mitochondrial mass, whereas Ras-NIH 3T3 cells maintained high mitochondrial mass. Immunofluorescence analysis of mitochondrial membrane marker proteins and mitochondrial membrane potential (MMP) also demonstrated a temporal pattern of changes similar to those of mitochondrial mass. This finding could be attributed to the relatively higher level of alternative autophagy in Ras-NIH 3T3 cells facilitating the removal of damaged mitochondria. Paclitaxel treatment in Ras-NIH 3T3 cells induced an increase in canonical autophagy rates along with suppression of alternative autophagy. Ras-NIH 3T3 cells showed high sensitivity to paclitaxel at the early stage of culture, but as cell confluence increased, resistance to paclitaxel increased, showing a similar level of cell viability to the vehicle control group. The study findings suggest that alternative autophagy is more important than canonical autophagy for maintaining cell survival in response to an unfavorable environment, such as during high cell confluence and exposure to anticancer agents.

SARS-CoV-2 genomic diversity and the implications for qRT-PCR diagnostics and transmission.

The COVID-19 pandemic has sparked an urgent need to uncover the underlying biology of this devastating disease. Though RNA viruses mutate more rapidly than DNA viruses, there are a relatively small number of single nucleotide polymorphisms (SNPs) that differentiate the main SARS-CoV-2 lineages that have spread throughout the world. In this study, we investigated 129 RNA-seq data sets and 6928 consensus genomes to contrast the intra-host and inter-host diversity of SARS-CoV-2. Our analyses yielded three major observations. First, the mutational profile of SARS-CoV-2 highlights intra-host single nucleotide variant (iSNV) and SNP similarity, albeit with differences in C > U changes. Second, iSNV and SNP patterns in SARS-CoV-2 are more similar to MERS-CoV than SARS-CoV-1. Third, a significant fraction of insertions and deletions contribute to the genetic diversity of SARS-CoV-2. Altogether, our findings provide insight into SARS-CoV-2 genomic diversity, inform the design of detection tests, and highlight the potential of iSNVs for tracking the transmission of SARS-CoV-2.