Accurate expression quantification from nanopore direct RNA sequencing with NanoCount.

Accurately quantifying gene and isoform expression changes is essential to understanding cell functions, differentiation and disease. Sequencing full-length native RNAs using long-read direct RNA sequencing (DRS) has the potential to overcome many limitations of short and long-read sequencing methods that require RNA fragmentation, cDNA synthesis or PCR. However, there are a lack of tools specifically designed for DRS and its ability to identify differential expression in complex organisms is poorly characterised. We developed NanoCount for fast, accurate transcript isoform quantification in DRS and demonstrate it outperforms similar methods. Using synthetic controls and human SH-SY5Y cell differentiation into neuron-like cells, we show that DRS accurately quantifies RNA expression and identifies differential expression of genes and isoforms. Differential expression of 231 genes, 333 isoforms, plus 27 isoform switches were detected between undifferentiated and differentiated SH-SY5Y cells and samples clustered by differentiation state at the gene and isoform level. Genes upregulated in neuron-like cells were associated with neurogenesis. NanoCount quantification of thousands of novel isoforms discovered with DRS likewise enabled identification of their differential expression. Our results demonstrate enhanced DRS isoform quantification with NanoCount and establish the ability of DRS to identify biologically relevant differential expression of genes and isoforms.

Nanopore ReCappable sequencing maps SARS-CoV-2 5' capping sites and provides new insights into the structure of sgRNAs.

The SARS-CoV-2 virus has a complex transcriptome characterised by multiple, nested subgenomic RNAsused to express structural and accessory proteins. Long-read sequencing technologies such as nanopore direct RNA sequencing can recover full-length transcripts, greatly simplifying the assembly of structurally complex RNAs. However, these techniques do not detect the 5' cap, thus preventing reliable identification and quantification of full-length, coding transcript models. Here we used Nanopore ReCappable Sequencing (NRCeq), a new technique that can identify capped full-length RNAs, to assemble a complete annotation of SARS-CoV-2 sgRNAs and annotate the location of capping sites across the viral genome. We obtained robust estimates of sgRNA expression across cell lines and viral isolates and identified novel canonical and non-canonical sgRNAs, including one that uses a previously un-annotated leader-to-body junction site. The data generated in this work constitute a useful resource for the scientific community and provide important insights into the mechanisms that regulate the transcription of SARS-CoV-2 sgRNAs.

Chemical and photochemical error rates in light-directed synthesis of complex DNA libraries.

Nucleic acid microarrays are the only tools that can supply very large oligonucleotide libraries, cornerstones of the nascent fields of de novo gene assembly and DNA data storage. Although the chemical synthesis of oligonucleotides is highly developed and robust, it is not error free, requiring the design of methods that can correct or compensate for errors, or select for high-fidelity oligomers. However, outside the realm of array manufacturers, little is known about the sources of errors and their extent. In this study, we look at the error rate of DNA libraries synthesized by photolithography and dissect the proportion of deletion, insertion and substitution errors. We find that the deletion rate is governed by the photolysis yield. We identify the most important substitution error and correlate it to phosphoramidite coupling. Besides synthetic failures originating from the coupling cycle, we uncover the role of imperfections and limitations related to optics, highlight the importance of absorbing UV light to avoid internal reflections and chart the dependence of error rate on both position on the array and position within individual oligonucleotides. Being able to precisely quantify all types of errors will allow for optimal choice of fabrication parameters and array design.

pycoMeth: a toolbox for differential methylation testing from Nanopore methylation calls.

We present pycoMeth, a toolbox to store, manage and analyze DNA methylation calls from long-read sequencing data obtained using the Oxford Nanopore Technologies sequencing platform. Building on a novel, rapid-access, read-level and reference-anchored methylation storage format MetH5, we propose efficient algorithms for haplotype aware, multi-sample consensus segmentation and differential methylation testing. We show that MetH5 is more efficient than existing solutions for storing Oxford Nanopore Technologies methylation calls, and carry out benchmarking for pycoMeth segmentation and differential methylation testing, demonstrating increased performance and sensitivity compared to existing solutions designed for short-read methylation data.

Natural genetic variation quantitatively regulates heart rate and dimension.

The polygenic contribution to heart development and function along the health-disease continuum remains unresolved. To gain insight into the genetic basis of quantitative cardiac phenotypes, we utilize highly inbred Japanese rice fish models, Oryzias latipes, and Oryzias sakaizumii. Employing automated quantification of embryonic heart rates as core metric, we profiled phenotype variability across five inbred strains. We observed maximal phenotypic contrast between individuals of the HO5 and the HdrR strain. HO5 showed elevated heart rates associated with embryonic ventricular hypoplasia and impaired adult cardiac function. This contrast served as the basis for genome-wide mapping. In a segregation population of 1192 HO5 x HdrR F2 embryos, we mapped 59 loci (173 genes) associated with heart rate. Experimental validation of the top 12 candidate genes in loss-of-function models revealed their causal and distinct impact on heart rate, development, ventricle size, and arrhythmia. Our study uncovers new diagnostic and therapeutic targets for developmental and electrophysiological cardiac diseases and provides a novel scalable approach to investigate the intricate genetic architecture of the vertebrate heart.

Long-read sequencing of diagnosis and post-therapy medulloblastoma reveals complex rearrangement patterns and epigenetic signatures.

Cancer genomes harbor a broad spectrum of structural variants (SVs) driving tumorigenesis, a relevant subset of which escape discovery using short-read sequencing. We employed Oxford Nanopore Technologies (ONT) long-read sequencing in a paired diagnostic and post-therapy medulloblastoma to unravel the haplotype-resolved somatic genetic and epigenetic landscape. We assembled complex rearrangements, including a 1.55-Mbp chromothripsis event, and we uncover a complex SV pattern termed templated insertion (TI) thread, characterized by short (mostly <1 kb) insertions showing prevalent self-concatenation into highly amplified structures of up to 50 kbp in size. TI threads occur in 3% of cancers, with a prevalence up to 74% in liposarcoma, and frequent colocalization with chromothripsis. We also perform long-read-based methylome profiling and discover allele-specific methylation (ASM) effects, complex rearrangements exhibiting differential methylation, and differential promoter methylation in cancer-driver genes. Our study shows the advantage of long-read sequencing in the discovery and characterization of complex somatic rearrangements.

The Medaka Inbred Kiyosu-Karlsruhe (MIKK) panel.

BACKGROUND: Unraveling the relationship between genetic variation and phenotypic traits remains a fundamental challenge in biology. Mapping variants underlying complex traits while controlling for confounding environmental factors is often problematic. To address this, we establish a vertebrate genetic resource specifically to allow for robust genotype-to-phenotype investigations. The teleost medaka (Oryzias latipes) is an established genetic model system with a long history of genetic research and a high tolerance to inbreeding from the wild. RESULTS: Here we present the Medaka Inbred Kiyosu-Karlsruhe (MIKK) panel: the first near-isogenic panel of 80 inbred lines in a vertebrate model derived from a wild founder population. Inbred lines provide fixed genomes that are a prerequisite for the replication of studies, studies which vary both the genetics and environment in a controlled manner, and functional testing. The MIKK panel will therefore enable phenotype-to-genotype association studies of complex genetic traits while allowing for careful control of interacting factors, with numerous applications in genetic research, human health, drug development, and fundamental biology. CONCLUSIONS: Here we present a detailed characterization of the genetic variation across the MIKK panel, which provides a rich and unique genetic resource to the community by enabling large-scale experiments for mapping complex traits.

Genomic variations and epigenomic landscape of the Medaka Inbred Kiyosu-Karlsruhe (MIKK) panel.

BACKGROUND: The teleost medaka (Oryzias latipes) is a well-established vertebrate model system, with a long history of genetic research, and multiple high-quality reference genomes available for several inbred strains. Medaka has a high tolerance to inbreeding from the wild, thus allowing one to establish inbred lines from wild founder individuals. RESULTS: We exploit this feature to create an inbred panel resource: the Medaka Inbred Kiyosu-Karlsruhe (MIKK) panel. This panel of 80 near-isogenic inbred lines contains a large amount of genetic variation inherited from the original wild population. We use Oxford Nanopore Technologies (ONT) long read data to further investigate the genomic and epigenomic landscapes of a subset of the MIKK panel. Nanopore sequencing allows us to identify a large variety of high-quality structural variants, and we present results and methods using a pan-genome graph representation of 12 individual medaka lines. This graph-based reference MIKK panel genome reveals novel differences between the MIKK panel lines and standard linear reference genomes. We find additional MIKK panel-specific genomic content that would be missing from linear reference alignment approaches. We are also able to identify and quantify the presence of repeat elements in each of the lines. Finally, we investigate line-specific CpG methylation and performed differential DNA methylation analysis across these 12 lines. CONCLUSIONS: We present a detailed analysis of the MIKK panel genomes using long and short read sequence technologies, creating a MIKK panel-specific pan genome reference dataset allowing for investigation of novel variation types that would be elusive using standard approaches.

Nanopore RNA Sequencing Analysis.

The recent advent of Nanopore sequencing allows for the sequencing of full-length RNA or cDNA molecules. This new type of data introduces new challenges from the computational point of view, and requires new software as well as dedicated analysis pipelines. In this chapter, we guide the reader through the typical analysis steps required to process the raw data produced by the instrument into a table of counts suitable for downstream analyses. We first describe the procedure to convert raw direct RNA-Seq and cDNA-Seq data into sequences in fastq format. We then outline how to perform quality control and filtering steps and how to map the filtered long reads to a reference transcriptome or genome.

AAV8 locoregional delivery induces long-term expression of an immunogenic transgene in macaques despite persisting local inflammation.

Adeno-associated virus (AAV) vectors are considered efficient vectors for gene transfer, as illustrated by recent successful clinical trials targeting retinal or neurodegenerative disorders. However, limitations as host immune responses to AAV capsid or transduction of limited regions must still be overcome. Here, we focused on locoregional (LR) intravenous perfusion vector delivery that allows transduction of large muscular areas and is considered to be less immunogenic than intramuscular (IM) injection. To confirm this hypothesis, we injected 6 cynomolgus monkeys with an AAV serotype 8 (AAV8) vector encoding for the highly immunogenic GFP driven by either a muscle-specific promoter (n = 3) or a cytomegalovirus (CMV) promoter (n = 3). We report that LR delivery allows long-term GFP expression in the perfused limb (up to 1 year) despite the initiation of a peripheral transgene-specific immune response. The analysis of the immune status of the perfused limb shows that LR delivery induces persisting inflammation. However, this inflammation is not sufficient to result in transgene clearance and is balanced by resident regulatory T cells. Overall, our results suggest that LR delivery promotes persisting transgene expression by induction of Treg cells in situ and might be a safe alternative to IM route to target large muscle territories for the expression of secreted therapeutic factors.

RNA modifications detection by comparative Nanopore direct RNA sequencing.

RNA molecules undergo a vast array of chemical post-transcriptional modifications (PTMs) that can affect their structure and interaction properties. In recent years, a growing number of PTMs have been successfully mapped to the transcriptome using experimental approaches relying on high-throughput sequencing. Oxford Nanopore direct-RNA sequencing has been shown to be sensitive to RNA modifications. We developed and validated Nanocompore, a robust analytical framework that identifies modifications from these data. Our strategy compares an RNA sample of interest against a non-modified control sample, not requiring a training set and allowing the use of replicates. We show that Nanocompore can detect different RNA modifications with position accuracy in vitro, and we apply it to profile m6A in vivo in yeast and human RNAs, as well as in targeted non-coding RNAs. We confirm our results with orthogonal methods and provide novel insights on the co-occurrence of multiple modified residues on individual RNA molecules.

Single-Stranded DNA Virus Sequencing (SSV-Seq) for Characterization of Residual DNA and AAV Vector Genomes.

With the success of clinical trials using recombinant adeno-associated viral vectors (rAAV), regulatory agencies ask for a more comprehensive characterization of process- and product- related impurities found in rAAV stocks in order to assess the potential risks for patients. During production, rAAV capsids are known to internalize illegitimate DNA fragments in addition to their recombinant genome. These contaminants can come from plasmid or helper virus DNA as well as from the producer host cell. Here, we describe a method based on high-throughput sequencing to identify and quantify residual DNA in rAAV vector lots. Contrary to qPCR, SSV-Seq (Single-Stranded DNA Virus Sequencing) offers a nonselective approach to determine the percentage of each DNA contaminant and analyze rAAV vector genome identity.

Genetics instability of wtAAV2 genome and AAV promoter activities in the Baculovirus/Sf9 cells system.

The human Adeno-Associated Virus serotype 2 (wtAAV2) is a common non-pathological virus and its recombinant form (rAAV) is widely used as gene therapy vector. Although rAAVs are routinely produced in the Baculovirus/Sf9 cell system, wtAAV2 has never been studied in this context. We tried to produce wtAAV2 in the baculovirus/Sf9 cell system hypothesizing that the wtAAV2 may be considered as a normal recombinant AAV transgene. Through our attempts to produce wtAAV2 in Baculovirus/Sf9, we found that wtAAV2 p5 promoter, which controls the expression of large Rep proteins in mammalian cells, was active in this system. p5 promoter activity in the baculovirus/Sf9 cell system leads to the expression of Rep78 that finally excises wtAAV2 genome from the baculovirus genome during the earliest phases of baculovirus stock production. Via p5 promoter expression kinetics and strand specific RNA-Seq analysis of wtAAV2, rAAV and Rep2/Cap2 cassettes in the baculovirus context we could demonstrate that wtAAV2 native promoters, p5, p19 and p40 are all active in the context of the baculovirus system and lead to the expression of different proteins and peptides. In addition, this study has proven that the baculovirus brings at least some of the helper functions needed in the AAV replication/life cycle.

RNA-Seq Analysis of an Antisense Sequence Optimized for Exon Skipping in Duchenne Patients Reveals No Off-Target Effect.

Non-coding uridine-rich small nuclear RNAs (UsnRNAs) have emerged in recent years as effective tools for exon skipping for the treatment of Duchenne muscular dystrophy (DMD), a degenerative muscular genetic disorder. We recently showed the high capacity of a recombinant adeno-associated virus (rAAV)-U7snRNA vector to restore the reading frame of the DMD mRNA in the muscles of DMD dogs. We are now moving toward a phase I/II clinical trial with an rAAV-U7snRNA-E53, carrying an antisense sequence designed to hybridize exon 53 of the human DMD messenger. As observed for genome-editing tools, antisense sequences present a risk of off-target effects, reflecting partial hybridization onto unintended transcripts. To characterize the clinical antisense sequence, we studied its expression and explored the occurrence of its off-target effects in human in vitro models of skeletal muscle and liver. We presented a comprehensive methodology combining RNA sequencing and in silico filtering to analyze off-targets. We showed that U7snRNA-E53 induced the effective exon skipping of the DMD transcript without inducing the notable deregulation of transcripts in human cells, neither at gene expression nor at the mRNA splicing level. Altogether, these results suggest that the use of the rAAV-U7snRNA-E53 vector for exon skipping could be safe in eligible DMD patients.

Accurate Identification and Quantification of DNA Species by Next-Generation Sequencing in Adeno-Associated Viral Vectors Produced in Insect Cells.

Recombinant adeno-associated viral (rAAV) vectors have proven excellent tools for the treatment of many genetic diseases and other complex diseases. However, the illegitimate encapsidation of DNA contaminants within viral particles constitutes a major safety concern for rAAV-based therapies. Moreover, the development of rAAV vectors for early-phase clinical trials has revealed the limited accuracy of the analytical tools used to characterize these new and complex drugs. Although most published data concerning residual DNA in rAAV preparations have been generated by quantitative PCR, we have developed a novel single-strand virus sequencing (SSV-Seq) method for quantification of DNA contaminants in AAV vectors produced in mammalian cells by next-generation sequencing (NGS). Here, we describe the adaptation of SSV-Seq for the accurate identification and quantification of DNA species in rAAV stocks produced in insect cells. We found that baculoviral DNA was the most abundant contaminant, representing less than 2.1% of NGS reads regardless of serotype (2, 8, or rh10). Sf9 producer cell DNA was detected at low frequency (≤0.03%) in rAAV lots. Advanced computational analyses revealed that (1) baculoviral sequences close to the inverted terminal repeats preferentially underwent illegitimate encapsidation, and (2) single-nucleotide variants were absent from the rAAV genome. The high-throughput sequencing protocol described here enables effective DNA quality control of rAAV vectors produced in insect cells, and is adapted to conform with regulatory agency safety requirements.

Identification of miRNAs Involved in Reprogramming Acinar Cells into Insulin Producing Cells.

Reprogramming acinar cells into insulin producing cells using adenoviral (Ad)-mediated delivery of Pdx1, Ngn3 and MafA (PNM) is an innovative approach for the treatment of diabetes. Here, we aimed to investigate the molecular mechanisms involved in this process and in particular, the role of microRNAs. To this end, we performed a comparative study of acinar-to-ß cell reprogramming efficiency in the rat acinar cell line AR42J and its subclone B13 after transduction with Ad-PNM. B13 cells were more efficiently reprogrammed than AR42J cells, which was demonstrated by a strong activation of ß cell markers (Ins1, Ins2, IAPP, NeuroD1 and Pax4). miRNome panels were used to analyze differentially expressed miRNAs in acinar cells under four experimental conditions (i) non-transduced AR42J cells, (ii) non-transduced B13 cells, (iii) B13 cells transduced with Ad-GFP vectors and (iv) B13 cells transduced with Ad-PNM vectors. A total of 59 miRNAs were found to be differentially expressed between non-transduced AR42J and B13 cells. Specifically, the miR-200 family was completely repressed in B13 cells, suggesting that these cells exist in a less differentiated state than AR42J cells and as a consequence they present a greater plasticity. Adenoviral transduction per se induced dedifferentiation of acinar cells and 11 miRNAs were putatively involved in this process, whereas 8 miRNAs were found to be associated with PNM expression. Of note, Ad-PNM reprogrammed B13 cells presented the same levels of miR-137-3p, miR-135a-5p, miR-204-5p and miR-210-3p of those detected in islets, highlighting their role in the process. In conclusion, this study led to the identification of miRNAs that might be of compelling importance to improve acinar-to-ß cell conversion for the future treatment of diabetes.

Advanced Characterization of DNA Molecules in rAAV Vector Preparations by Single-stranded Virus Next-generation Sequencing.

Recent successful clinical trials with recombinant adeno-associated viral vectors (rAAVs) have led to a renewed interest in gene therapy. However, despite extensive developments to improve vector-manufacturing processes, undesirable DNA contaminants in rAAV preparations remain a major safety concern. Indeed, the presence of DNA fragments containing antibiotic resistance genes, wild-type AAV, and packaging cell genomes has been found in previous studies using quantitative polymerase chain reaction (qPCR) analyses. However, because qPCR only provides a partial view of the DNA molecules in rAAV preparations, we developed a method based on next-generation sequencing (NGS) to extensively characterize single-stranded DNA virus preparations (SSV-Seq). In order to validate SSV-Seq, we analyzed three rAAV vector preparations produced by transient transfection of mammalian cells. Our data were consistent with qPCR results and showed a quasi-random distribution of contaminants originating from the packaging cells genome. Finally, we found single-nucleotide variants (SNVs) along the vector genome but no evidence of large deletions. Altogether, SSV-Seq could provide a characterization of DNA contaminants and a map of the rAAV genome with unprecedented resolution and exhaustiveness. We expect SSV-Seq to pave the way for a new generation of quality controls, guiding process development toward rAAV preparations of higher potency and with improved safety profiles.

Short-lived recombinant adeno-associated virus transgene expression in dystrophic muscle is associated with oxidative damage to transgene mRNA.

Preclinical gene therapy strategies using recombinant adeno-associated virus (AAV) vectors in animal models of Duchenne muscular dystrophy have shown dramatic phenotype improvements, but long-lasting efficacy remains questionable. It is believed that in dystrophic muscles, transgene persistence is hampered, notably by the progressive loss of therapeutic vector genomes resulting from muscle fibers degeneration. Intracellular metabolic perturbations resulting from dystrophin deficiency could also be additional factors impacting on rAAV genomes and transgene mRNA molecular fate. In this study, we showed that rAAV genome loss is not the only cause of reduced transgene mRNA level and we assessed the contribution of transcriptional and post-transcriptional factors. We ruled out the implication of transgene silencing by epigenetic mechanisms and demonstrated that rAAV inhibition occurred mostly at the post-transcriptional level. Since Duchenne muscular dystrophy (DMD) physiopathology involves an elevated oxidative stress, we hypothesized that in dystrophic muscles, transgene mRNA could be damaged by oxidative stress. In the mouse and dog dystrophic models, we found that rAAV-derived mRNA oxidation was increased. Interestingly, when a high expression level of a therapeutic transgene is achieved, oxidation is less pronounced. These findings provide new insights into rAAV transductions in dystrophic muscles, which ultimately may help in the design of more effective clinical trials.

Adeno-associated viral vector-mediated transgene expression is independent of DNA methylation in primate liver and skeletal muscle.

Recombinant adeno-associated viral (rAAV) vectors can support long-term transgene expression in quiescent tissues. Intramuscular (i.m.) administration of a single-stranded AAV vector (ssAAV) in the nonhuman primate (NHP) results in a peak protein level at 2-3 months, followed by a decrease over several months before reaching a steady-state. To investigate transgene expression and vector genome persistence, we previously demonstrated that rAAV vector genomes associate with histones and form a chromatin structure in NHP skeletal muscle more than one year after injection. In the mammalian nucleus, chromatin remodeling via epigenetic modifications plays key role in transcriptional regulation. Among those, CpG hyper-methylation of promoters is a known hallmark of gene silencing. To assess the involvement of DNA methylation on the transgene expression, we injected NHP via the i.m. or the intravenous (i.v.) route with a recombinant ssAAV2/1 vector. The expression cassette contains the transgene under the transcriptional control of the constitutive Rous Sarcoma Virus promoter (RSVp). Total DNA isolated from NHP muscle and liver biopsies from 1 to 37 months post-injection was treated with sodium bisulfite and subsequently analyzed by pyrosequencing. No significant CpG methylation of the RSVp was found in rAAV virions or in vector DNA isolated from NHP transduced tissues. Direct de novo DNA methylation appears not to be involved in repressing transgene expression in NHP after gene transfer mediated by ssAAV vectors. The study presented here examines host/vector interactions and the impact on transgene expression in a clinically relevant model.