Unraveling C-to-U RNA editing events from direct RNA sequencing.

In mammals, RNA editing events involve the conversion of adenosine (A) in inosine (I) by ADAR enzymes or the hydrolytic deamination of cytosine (C) in uracil (U) by the APOBEC family of enzymes, mostly APOBEC1. RNA editing has a plethora of biological functions, and its deregulation has been associated with various human disorders. While the large-scale detection of A-to-I is quite straightforward using the Illumina RNAseq technology, the identification of C-to-U events is a non-trivial task. This difficulty arises from the rarity of such events in eukaryotic genomes and the challenge of distinguishing them from background noise. Direct RNA sequencing by Oxford Nanopore Technology (ONT) permits the direct detection of Us on sequenced RNA reads. Surprisingly, using ONT reads from wild-type (WT) and APOBEC1-knock-out (KO) murine cell lines as well as in vitro synthesized RNA without any modification, we identified a systematic error affecting the accuracy of the Cs call, thereby leading to incorrect identifications of C-to-U events. To overcome this issue in direct RNA reads, here we introduce a novel machine learning strategy based on the isolation Forest (iForest) algorithm in which C-to-U editing events are considered as sequencing anomalies. Using in vitro synthesized and human ONT reads, our model optimizes the signal-to-noise ratio improving the detection of C-to-U editing sites with high accuracy, over 90% in all samples tested. Our results suggest that iForest, known for its rapid implementation and minimal memory requirements, is a promising tool to denoise ONT reads and reliably identify RNA modifications.

COVID-19 annual update: a narrative review.

Three and a half years after the pandemic outbreak, now that WHO has formally declared that the emergency is over, COVID-19 is still a significant global issue. Here, we focus on recent developments in genetic and genomic research on COVID-19, and we give an outlook on state-of-the-art therapeutical approaches, as the pandemic is gradually transitioning to an endemic situation. The sequencing and characterization of rare alleles in different populations has made it possible to identify numerous genes that affect either susceptibility to COVID-19 or the severity of the disease. These findings provide a beginning to new avenues and pan-ethnic therapeutic approaches, as well as to potential genetic screening protocols. The causative virus, SARS-CoV-2, is still in the spotlight, but novel threatening virus could appear anywhere at any time. Therefore, continued vigilance and further research is warranted. We also note emphatically that to prevent future pandemics and other world-wide health crises, it is imperative to capitalize on what we have learnt from COVID-19: specifically, regarding its origins, the world's response, and insufficient preparedness. This requires unprecedented international collaboration and timely data sharing for the coordination of effective response and the rapid implementation of containment measures.

DNA Deamination Is Required for Human APOBEC3A-Driven Hepatocellular Carcinoma In Vivo.

Although the APOBEC3 family of single-stranded DNA cytosine deaminases is well-known for its antiviral factors, these enzymes are rapidly gaining attention as prominent sources of mutation in cancer. APOBEC3's signature single-base substitutions, C-to-T and C-to-G in TCA and TCT motifs, are evident in over 70% of human malignancies and dominate the mutational landscape of numerous individual tumors. Recent murine studies have established cause-and-effect relationships, with both human APOBEC3A and APOBEC3B proving capable of promoting tumor formation in vivo. Here, we investigate the molecular mechanism of APOBEC3A-driven tumor development using the murine Fah liver complementation and regeneration system. First, we show that APOBEC3A alone is capable of driving tumor development (without Tp53 knockdown as utilized in prior studies). Second, we show that the catalytic glutamic acid residue of APOBEC3A (E72) is required for tumor formation. Third, we show that an APOBEC3A separation-of-function mutant with compromised DNA deamination activity and wildtype RNA-editing activity is defective in promoting tumor formation. Collectively, these results demonstrate that APOBEC3A is a "master driver" that fuels tumor formation through a DNA deamination-dependent mechanism.

Identification and characterization of novel ETV4 splice variants in prostate cancer.

ETV4, one of ETS proteins overexpressed in prostate cancer, promotes migration, invasion, and proliferation in prostate cells. This study identifies a series of previously unknown ETV4 alternatively spliced transcripts in human prostate cell lines. Their expression has been validated using several unbiased techniques, including Nanopore sequencing. Most of these transcripts originate from an in-frame exon skipping and, thus, are expected to be translated into ETV4 protein isoforms. Functional analysis of the most abundant among these isoforms shows that they still bear an activity, namely a reduced ability to promote proliferation and a residual ability to regulate the transcription of ETV4 target genes. Alternatively spliced genes are common in cancer cells: an analysis of the TCGA dataset confirms the abundance of these novel ETV4 transcripts in prostate tumors, in contrast to peritumoral tissues. Since none of their translated isoforms have acquired a higher oncogenic potential, such abundance is likely to reflect the tumor deranged splicing machinery. However, it is also possible that their interaction with the canonical variants may contribute to the biology and the clinics of prostate cancer. Further investigations are needed to elucidate the biological role of these ETV4 transcripts and of their putative isoforms.

Detecting cell-of-origin and cancer-specific methylation features of cell-free DNA from Nanopore sequencing.

The Oxford Nanopore (ONT) platform provides portable and rapid genome sequencing, and its ability to natively profile DNA methylation without complex sample processing is attractive for point-of-care real-time sequencing. We recently demonstrated ONT shallow whole-genome sequencing to detect copy number alterations (CNAs) from the circulating tumor DNA (ctDNA) of cancer patients. Here, we show that cell type and cancer-specific methylation changes can also be detected, as well as cancer-associated fragmentation signatures. This feasibility study suggests that ONT shallow WGS could be a powerful tool for liquid biopsy.

MODOMICS: a database of RNA modification pathways. 2021 update.

The MODOMICS database has been, since 2006, a manually curated and centralized resource, storing and distributing comprehensive information about modified ribonucleosides. Originally, it only contained data on the chemical structures of modified ribonucleosides, their biosynthetic pathways, the location of modified residues in RNA sequences, and RNA-modifying enzymes. Over the years, prompted by the accumulation of new knowledge and new types of data, it has been updated with new information and functionalities. In this new release, we have created a catalog of RNA modifications linked to human diseases, e.g., due to mutations in genes encoding modification enzymes. MODOMICS has been linked extensively to RCSB Protein Data Bank, and sequences of experimentally determined RNA structures with modified residues have been added. This expansion was accompanied by including nucleotide 5'-monophosphate residues. We redesigned the web interface and upgraded the database backend. In addition, a search engine for chemically similar modified residues has been included that can be queried by SMILES codes or by drawing chemical molecules. Finally, previously available datasets of modified residues, biosynthetic pathways, and RNA-modifying enzymes have been updated. Overall, we provide users with a new, enhanced, and restyled tool for research on RNA modification. MODOMICS is available at https://iimcb.genesilico.pl/modomics/.

Fam72a enforces error-prone DNA repair during antibody diversification.

Efficient humoral responses rely on DNA damage, mutagenesis and error-prone DNA repair. Diversification of B cell receptors through somatic hypermutation and class-switch recombination are initiated by cytidine deamination in DNA mediated by activation-induced cytidine deaminase (AID)1 and by the subsequent excision of the resulting uracils by uracil DNA glycosylase (UNG) and by mismatch repair proteins1-3. Although uracils arising in DNA are accurately repaired1-4, how these pathways are co-opted to generate mutations and double-strand DNA breaks in the context of somatic hypermutation and class-switch recombination is unknown1-3. Here we performed a genome-wide CRISPR-Cas9 knockout screen for genes involved in class-switch recombination and identified FAM72A, a protein that interacts with the nuclear isoform of UNG (UNG2)5 and is overexpressed in several cancers5. We show that the FAM72A-UNG2 interaction controls the levels of UNG2 and that class-switch recombination is defective in Fam72a-/- B cells due to the upregulation of UNG2. Moreover, we show that somatic hypermutation is reduced in Fam72a-/- B cells and that its pattern is skewed upon upregulation of UNG2. Our results are consistent with a model in which FAM72A interacts with UNG2 to control its physiological level by triggering its degradation, regulating the level of uracil excision and thus the balance between error-prone and error-free DNA repair. Our findings have potential implications for tumorigenesis, as reduced levels of UNG2 mediated by overexpression of Fam72a would shift the balance towards mutagenic DNA repair, rendering cells more prone to acquire mutations.

Nanopore sequencing from liquid biopsy: analysis of copy number variations from cell-free DNA of lung cancer patients.

In the "precision oncology" era the characterization of tumor genetic features is a pivotal step in cancer patients' management. Liquid biopsy approaches, such as analysis of cell-free DNA from plasma, represent a powerful and noninvasive strategy to obtain information about the genomic status of the tumor. Sequencing-based analyses of cell-free DNA, currently performed with second generation sequencers, are extremely powerful but poorly scalable and not always accessible also due to instrumentation costs. Third generation sequencing platforms, such as Nanopore sequencers, aim at overcoming these obstacles but, unfortunately, are not designed for cell-free DNA analysis.Here we present a customized workflow to exploit low-coverage Nanopore sequencing for the detection of copy number variations from plasma of cancer patients. Whole genome molecular karyotypes of 6 lung cancer patients and 4 healthy subjects were successfully produced with as few as 2 million reads, and common lung-related copy number alterations were readily detected.This is the first successful use of Nanopore sequencing for copy number profiling from plasma DNA. In this context, Nanopore represents a reliable alternative to Illumina sequencing, with the advantages of minute instrumentation costs and extremely short analysis time.The availability of protocols for Nanopore-based cell-free DNA analysis will make this analysis finally accessible, exploiting the full potential of liquid biopsy both for research and clinical purposes.

A mark of disease: how mRNA modifications shape genetic and acquired pathologies.

RNA modifications have recently emerged as a widespread and complex facet of gene expression regulation. Counting more than 170 distinct chemical modifications with far-reaching implications for RNA fate, they are collectively referred to as the epitranscriptome. These modifications can occur in all RNA species, including messenger RNAs (mRNAs) and noncoding RNAs (ncRNAs). In mRNAs the deposition, removal, and recognition of chemical marks by writers, erasers and readers influence their structure, localization, stability, and translation. In turn, this modulates key molecular and cellular processes such as RNA metabolism, cell cycle, apoptosis, and others. Unsurprisingly, given their relevance for cellular and organismal functions, alterations of epitranscriptomic marks have been observed in a broad range of human diseases, including cancer, neurological and metabolic disorders. Here, we will review the major types of mRNA modifications and editing processes in conjunction with the enzymes involved in their metabolism and describe their impact on human diseases. We present the current knowledge in an updated catalog. We will also discuss the emerging evidence on the crosstalk of epitranscriptomic marks and what this interplay could imply for the dynamics of mRNA modifications. Understanding how this complex regulatory layer can affect the course of human pathologies will ultimately lead to its exploitation toward novel epitranscriptomic therapeutic strategies.

Live-Cell Quantification of APOBEC1-Mediated RNA Editing: A Comparison of RNA Editing Assays.

APOBEC1 is a member of the AID/APOBECs, a group of deaminases responsible for the editing of C>U in both DNA and RNA. APOBEC1 is physiologically involved in C>U RNA editing: while hundreds of targets have been discovered in mice, in humans the only well-characterized target of APOBEC1 is the apolipoprotein B (ApoB) transcript. APOBEC1 edits a CAA codon into a stop codon, which causes the translation of a truncated form of ApoB. A number of assays have been developed to investigate this process. Early assays, poisoned primer extension and Sanger sequencing, have focused on accuracy and sensitivity but rely on extraction of the RNA from tissues and cells. More recently, the need to visualize the RNA editing process directly in live cells have led to the development of fluorescence-based tools. These assays detect RNA editing through reporters whose editing causes a change in cellular localization or a change in fluorescent properties. Here we review the available assays to quantify RNA editing, and we present the protocol for cytofluorimetric analysis using a double-fluorescent reporter.

Evidence for host-dependent RNA editing in the transcriptome of SARS-CoV-2.

The COVID-19 outbreak has become a global health risk, and understanding the response of the host to the SARS-CoV-2 virus will help to combat the disease. RNA editing by host deaminases is an innate restriction process to counter virus infection, but it is not yet known whether this process operates against coronaviruses. Here, we analyze RNA sequences from bronchoalveolar lavage fluids obtained from coronavirus-infected patients. We identify nucleotide changes that may be signatures of RNA editing: adenosine-to-inosine changes from ADAR deaminases and cytosine-to-uracil changes from APOBEC deaminases. Mutational analysis of genomes from different strains of Coronaviridae from human hosts reveals mutational patterns consistent with those observed in the transcriptomic data. However, the reduced ADAR signature in these data raises the possibility that ADARs might be more effective than APOBECs in restricting viral propagation. Our results thus suggest that both APOBECs and ADARs are involved in coronavirus genome editing, a process that may shape the fate of both virus and patient.

AAV-mediated FOXG1 gene editing in human Rett primary cells.

Variations in the Forkhead Box G1 (FOXG1) gene cause FOXG1 syndrome spectrum, including the congenital variant of Rett syndrome, characterized by early onset of regression, Rett-like and jerky movements, and cortical visual impairment. Due to the largely unknown pathophysiological mechanisms downstream the impairment of this transcriptional regulator, a specific treatment is not yet available. Since both haploinsufficiency and hyper-expression of FOXG1 cause diseases in humans, we reasoned that adding a gene under nonnative regulatory sequences would be a risky strategy as opposed to a genome editing approach where the mutated gene is reversed into wild-type. Here, we demonstrate that an adeno-associated viruses (AAVs)-coupled CRISPR/Cas9 system is able to target and correct FOXG1 variants in patient-derived fibroblasts, induced Pluripotent Stem Cells (iPSCs) and iPSC-derived neurons. Variant-specific single-guide RNAs (sgRNAs) and donor DNAs have been selected and cloned together with a mCherry/EGFP reporter system. Specific sgRNA recognition sequences were inserted upstream and downstream Cas9 CDS to allow self-cleavage and inactivation. We demonstrated that AAV serotypes vary in transduction efficiency depending on the target cell type, the best being AAV9 in fibroblasts and iPSC-derived neurons, and AAV2 in iPSCs. Next-generation sequencing (NGS) of mCherry+/EGFP+ transfected cells demonstrated that the mutated alleles were repaired with high efficiency (20-35% reversion) and precision both in terms of allelic discrimination and off-target activity. The genome editing strategy tested in this study has proven to precisely repair FOXG1 and delivery through an AAV9-based system represents a step forward toward the development of a therapy for Rett syndrome.

High rate of HDR in gene editing of p.(Thr158Met) MECP2 mutational hotspot.

Rett syndrome is a progressive neurodevelopmental disorder which affects almost exclusively girls, caused by variants in MECP2 gene. Effective therapies for this devastating disorder are not yet available and the need for tight regulation of MECP2 expression for brain to properly function makes gene replacement therapy risky. For this reason, gene editing with CRISPR/Cas9 technology appears as a preferable option for the development of new therapies. To study the disease, we developed and characterized a human neuronal model obtained by genetic reprogramming of patient-derived primary fibroblasts into induced Pluripotent Stem Cells. This cellular model represents an important source for our studies, aiming to correct MECP2 variants in neurons which represent the primarily affected cell type. We engineered a gene editing toolkit composed by a two-plasmid system to correct a hotspot missense variant in MECP2, c.473 C > T (p.(Thr158Met)). The first construct expresses the variant-specific sgRNA and the Donor DNA along with a fluorescent reporter system. The second construct brings Cas9 and targets for auto-cleaving, to avoid long-term Cas9 expression. NGS analysis on sorted cells from four independent patients demonstrated an exceptionally high editing efficiency, with up to 80% of HDR and less than 1% of indels in all patients, outlining the relevant potentiality of the approach for Rett syndrome therapy.

New frontiers to cure Alport syndrome: COL4A3 and COL4A5 gene editing in podocyte-lineage cells.

Alport syndrome (AS) is an inherited genetic disorder characterized by range of alterations from glomerular basement membrane abnormalities up to end-stage renal disease. Pathogenic variants in the collagen α3, α4, and α5 encoding genes are causative both of the autosomal dominant and of the X-linked forms of AS. Podocytes are the only renal cells that are able to produce the COL(IV)a3-a4a5 heterotrimer. We have previously demonstrated how it is possible to isolate podocyte-lineage cells from urine of patients, providing an easily accessible cellular model closer to the podocytes' physiological conditions. Taking advantage of disease-relevant cell lines, we employed a two-plasmid approach in order to achieve a beneficial and stable variant-specific correction using CRISPR/Cas9 genome editing. One plasmid carries a Donor DNA and a reporter system mCherry/GFP to track the activity of Cas9 in cells. The other plasmid carries a self-cleaving SpCas9 and the variant-specific sgRNA. We have analyzed two stable podocyte-lineage cell lines, harboring a variant in the X-linked COL4A5 (p.(Gly624Asp)) and in the autosomal COL4A3 gene (p.(Gly856Glu)). We have achieved reversion of variants greater than 40% with undesired insertions/deletions lower than 15%. Overall, we have demonstrated a new gene therapy approach directly on patients' cells, key players of Alport pathogenesis, and we have reverted COL4 causative variants towards the wild type state. These results, in combination with preclinical models, could open new frontiers in the management and the treatment of the disorder.

A fluorescent reporter for quantification and enrichment of DNA editing by APOBEC-Cas9 or cleavage by Cas9 in living cells.

Base editing is an exciting new genome engineering technology. C-to-T mutations in genomic DNA have been achieved using ribonucleoprotein complexes comprised of rat APOBEC1 single-stranded DNA deaminase, Cas9 nickase (Cas9n), uracil DNA glycosylase inhibitor (UGI), and guide (g)RNA. Here, we report the first real-time reporter system for quantification of APOBEC-mediated base editing activity in living mammalian cells. The reporter expresses eGFP constitutively as a marker for transfection or transduction, and editing restores functionality of an upstream mCherry cassette through the simultaneous processing of two gRNA binding regions that each contain an APOBEC-preferred 5'TCA target site. Using this system as both an episomal and a chromosomal editing reporter, we show that human APOBEC3A-Cas9n-UGI and APOBEC3B-Cas9n-UGI base editing complexes are more efficient than the original rat APOBEC1-Cas9n-UGI construct. We also demonstrate coincident enrichment of editing events at a heterologous chromosomal locus in reporter-edited, mCherry-positive cells. The mCherry reporter also quantifies the double-stranded DNA cleavage activity of Cas9, and may therefore be adaptable for use with many different CRISPR systems. The combination of a rapid, fluorescence-based editing reporter system and more efficient, structurally defined DNA editing enzymes broadens the versatility of the rapidly expanding toolbox of genome editing and engineering technologies.

An efficient method to enrich for knock-out and knock-in cellular clones using the CRISPR/Cas9 system.

Clustered Regularly Interspaced Short Palindromic Repeats-associated protein 9 nuclease (CRISPR/Cas9) and Transcription Activator-Like Effector Nucleases (TALENs) are versatile tools for genome editing. Here we report a method to increase the frequency of Cas9-targeted cellular clones. Our method is based on a chimeric construct with a Blasticidin S Resistance gene (bsr) placed out-of-frame by a surrogate target sequence. End joining of the CRISPR/Cas9-induced double-strand break on the surrogate target can place the bsr in frame, thus providing temporary resistance to Blasticidin S: this is used to enrich for cells where Cas9 is active. By this approach, in a real experimental setting, we disrupted the Aicda gene in ~70% of clones from CH12F3 lymphoma cells (>40% biallelically). With the same approach we knocked in a single nucleotide to reconstruct the frame of Aicda in these null cells, restoring the function in ~37% of the clones (less than 10% by the standard approach). Targeting of single nucleotide changes in other genes yielded analogous results. These results support our enrichment method as an efficient tool in genome editing.