Single-molecule mapping of replisome progression.

Fundamental aspects of DNA replication, such as the anatomy of replication stall sites, how replisomes are influenced by gene transcription, and whether the progression of sister replisomes is coordinated, are poorly understood. Available techniques do not allow the precise mapping of the positions of individual replisomes on chromatin. We have developed a method called Replicon-seq that entails the excision of full-length replicons by controlled nuclease cleavage at replication forks. Replicons are sequenced using Nanopore, which provides a single-molecule readout of long DNA. Using Replicon-seq, we found that sister replisomes function autonomously and yet progress through chromatin with remarkable consistency. Replication forks that encounter obstacles pause for a short duration but rapidly resume synthesis. The helicase Rrm3 plays a critical role both in mitigating the effect of protein barriers and with facilitating efficient termination. Replicon-seq provides a high-resolution means of defining how individual replisomes move across the genome.

The SARS-CoV-2 subgenome landscape and its novel regulatory features.

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is currently a global pandemic. CoVs are known to generate negative subgenomes (subgenomic RNAs [sgRNAs]) through transcription-regulating sequence (TRS)-dependent template switching, but the global dynamic landscapes of coronaviral subgenomes and regulatory rules remain unclear. Here, using next-generation sequencing (NGS) short-read and Nanopore long-read poly(A) RNA sequencing in two cell types at multiple time points after infection with SARS-CoV-2, we identified hundreds of template switches and constructed the dynamic landscapes of SARS-CoV-2 subgenomes. Interestingly, template switching could occur in a bidirectional manner, with diverse SARS-CoV-2 subgenomes generated from successive template-switching events. The majority of template switches result from RNA-RNA interactions, including seed and compensatory modes, with terminal pairing status as a key determinant. Two TRS-independent template switch modes are also responsible for subgenome biogenesis. Our findings reveal the subgenome landscape of SARS-CoV-2 and its regulatory features, providing a molecular basis for understanding subgenome biogenesis and developing novel anti-viral strategies.

Targeting the m6A RNA modification pathway blocks SARS-CoV-2 and HCoV-OC43 replication.

N6-methyladenosine (m6A) is an abundant internal RNA modification, influencing transcript fate and function in uninfected and virus-infected cells. Installation of m6A by the nuclear RNA methyltransferase METTL3 occurs cotranscriptionally; however, the genomes of some cytoplasmic RNA viruses are also m6A-modified. How the cellular m6A modification machinery impacts coronavirus replication, which occurs exclusively in the cytoplasm, is unknown. Here we show that replication of SARS-CoV-2, the agent responsible for the COVID-19 pandemic, and a seasonal human ß-coronavirus HCoV-OC43, can be suppressed by depletion of METTL3 or cytoplasmic m6A reader proteins YTHDF1 and YTHDF3 and by a highly specific small molecule METTL3 inhibitor. Reduction of infectious titer correlates with decreased synthesis of viral RNAs and the essential nucleocapsid (N) protein. Sites of m6A modification on genomic and subgenomic RNAs of both viruses were mapped by methylated RNA immunoprecipitation sequencing (meRIP-seq). Levels of host factors involved in m6A installation, removal, and recognition were unchanged by HCoV-OC43 infection; however, nuclear localization of METTL3 and cytoplasmic m6A readers YTHDF1 and YTHDF2 increased. This establishes that coronavirus RNAs are m6A-modified and host m6A pathway components control ß-coronavirus replication. Moreover, it illustrates the therapeutic potential of targeting the m6A pathway to restrict coronavirus reproduction.

Splicing Kinetics and Coordination Revealed by Direct Nascent RNA Sequencing through Nanopores.

Understanding how splicing events are coordinated across numerous introns in metazoan RNA transcripts requires quantitative analyses of transient RNA processing events in living cells. We developed nanopore analysis of co-transcriptional processing (nano-COP), in which nascent RNAs are directly sequenced through nanopores, exposing the dynamics and patterns of RNA splicing without biases introduced by amplification. Long nano-COP reads reveal that, in human and Drosophila cells, splicing occurs after RNA polymerase II transcribes several kilobases of pre-mRNA, suggesting that metazoan splicing transpires distally from the transcription machinery. Inhibition of the branch-site recognition complex SF3B rapidly diminished global co-transcriptional splicing. We found that splicing order does not strictly follow the order of transcription and is associated with cis-acting elements, alternative splicing, and RNA-binding factors. Further, neighboring introns in human cells tend to be spliced concurrently, implying that splicing of these introns occurs cooperatively. Thus, nano-COP unveils the organizational complexity of RNA processing.

Nanopore strand-specific mismatch enables de novo detection of bacterial DNA modifications.

DNA modifications in bacteria present diverse types and distributions, playing crucial functional roles. Current methods for detecting bacterial DNA modifications via nanopore sequencing typically involve comparing raw current signals to a methylation-free control. In this study, we found that bacterial DNA modification induces errors in nanopore reads. And these errors are found only in one strand but not the other, showing a strand-specific bias. Leveraging this discovery, we developed Hammerhead, a pioneering pipeline designed for de novo methylation discovery that circumvents the necessity of raw signal inference and a methylation-free control. The majority (14 out of 16) of the identified motifs can be validated by raw signal comparison methods or by identifying corresponding methyltransferases in bacteria. Additionally, we included a novel polishing strategy employing duplex reads to correct modification-induced errors in bacterial genome assemblies, achieving a reduction of over 85% in such errors. In summary, Hammerhead enables users to effectively locate bacterial DNA methylation sites from nanopore FASTQ/FASTA reads, thus holds promise as a routine pipeline for a wide range of nanopore sequencing applications, such as genome assembly, metagenomic binning, decontaminating eukaryotic genome assembly, and functional analysis for DNA modifications.

Genomics in the long-read sequencing era.

Long-read sequencing (LRS) technologies have provided extremely powerful tools to explore genomes. While in the early years these methods suffered technical limitations, they have recently made significant progress in terms of read length, throughput, and accuracy and bioinformatics tools have strongly improved. Here, we aim to review the current status of LRS technologies, the development of novel methods, and the impact on genomics research. We will explore the most impactful recent findings made possible by these technologies focusing on high-resolution sequencing of genomes and transcriptomes and the direct detection of DNA and RNA modifications. We will also discuss how LRS methods promise a more comprehensive understanding of human genetic variation, transcriptomics, and epigenetics for the coming years.

Beyond the exome: What's next in diagnostic testing for Mendelian conditions.

Despite advances in clinical genetic testing, including the introduction of exome sequencing (ES), more than 50% of individuals with a suspected Mendelian condition lack a precise molecular diagnosis. Clinical evaluation is increasingly undertaken by specialists outside of clinical genetics, often occurring in a tiered fashion and typically ending after ES. The current diagnostic rate reflects multiple factors, including technical limitations, incomplete understanding of variant pathogenicity, missing genotype-phenotype associations, complex gene-environment interactions, and reporting differences between clinical labs. Maintaining a clear understanding of the rapidly evolving landscape of diagnostic tests beyond ES, and their limitations, presents a challenge for non-genetics professionals. Newer tests, such as short-read genome or RNA sequencing, can be challenging to order, and emerging technologies, such as optical genome mapping and long-read DNA sequencing, are not available clinically. Furthermore, there is no clear guidance on the next best steps after inconclusive evaluation. Here, we review why a clinical genetic evaluation may be negative, discuss questions to be asked in this setting, and provide a framework for further investigation, including the advantages and disadvantages of new approaches that are nascent in the clinical sphere. We present a guide for the next best steps after inconclusive molecular testing based upon phenotype and prior evaluation, including when to consider referral to research consortia focused on elucidating the underlying cause of rare unsolved genetic disorders.

The functional impact of 1,570 individual amino acid substitutions in human OTC.

Deleterious mutations in the X-linked gene encoding ornithine transcarbamylase (OTC) cause the most common urea cycle disorder, OTC deficiency. This rare but highly actionable disease can present with severe neonatal onset in males or with later onset in either sex. Individuals with neonatal onset appear normal at birth but rapidly develop hyperammonemia, which can progress to cerebral edema, coma, and death, outcomes ameliorated by rapid diagnosis and treatment. Here, we develop a high-throughput functional assay for human OTC and individually measure the impact of 1,570 variants, 84% of all SNV-accessible missense mutations. Comparison to existing clinical significance calls, demonstrated that our assay distinguishes known benign from pathogenic variants and variants with neonatal onset from late-onset disease presentation. This functional stratification allowed us to identify score ranges corresponding to clinically relevant levels of impairment of OTC activity. Examining the results of our assay in the context of protein structure further allowed us to identify a 13 amino acid domain, the SMG loop, whose function appears to be required in human cells but not in yeast. Finally, inclusion of our data as PS3 evidence under the current ACMG guidelines, in a pilot reclassification of 34 variants with complete loss of activity, would change the classification of 22 from variants of unknown significance to clinically actionable likely pathogenic variants. These results illustrate how large-scale functional assays are especially powerful when applied to rare genetic diseases.

TERRA ONTseq: a long-read-based sequencing pipeline to study the human telomeric transcriptome.

The long noncoding RNA TERRA is transcribed from telomeres in virtually all eukaryotes with linear chromosomes. In humans, TERRA transcription is driven in part by promoters comprising CpG dinucleotide-rich repeats of 29 bp repeats, believed to be present in half of the subtelomeres. Thus far, TERRA expression has been analyzed mainly using molecular biology-based approaches that only generate partial and somehow biased results. Here, we present a novel experimental pipeline to study human TERRA based on long-read sequencing (TERRA ONTseq). By applying TERRA ONTseq to different cell lines, we show that the vast majority of human telomeres produce TERRA and that the cellular levels of TERRA transcripts vary according to their chromosomes of origin. Using TERRA ONTseq, we also identified regions containing TERRA transcription start sites (TSSs) in more than half of human subtelomeres. TERRA TSS regions are generally found immediately downstream from 29 bp repeat-related sequences, which appear to be more widespread than previously estimated. Finally, we isolated a novel TERRA promoter from the highly expressed subtelomere of the long arm of Chromosome 7. With the development of TERRA ONTseq, we provide a refined picture of human TERRA biogenesis and expression and we equip the scientific community with an invaluable tool for future studies.

ReadCurrent: a VDCNN-based tool for fast and accurate nanopore selective sequencing.

Nanopore selective sequencing allows the targeted sequencing of DNA of interest using computational approaches rather than experimental methods such as targeted multiplex polymerase chain reaction or hybridization capture. Compared to sequence-alignment strategies, deep learning (DL) models for classifying target and nontarget DNA provide large speed advantages. However, the relatively low accuracy of these DL-based tools hinders their application in nanopore selective sequencing. Here, we present a DL-based tool named ReadCurrent for nanopore selective sequencing, which takes electric currents as inputs. ReadCurrent employs a modified very deep convolutional neural network (VDCNN) architecture, enabling significantly lower computational costs for training and quicker inference compared to conventional VDCNN. We evaluated the performance of ReadCurrent across 10 nanopore sequencing datasets spanning human, yeasts, bacteria, and viruses. We observed that ReadCurrent achieved a mean accuracy of 98.57% for classification, outperforming four other DL-based selective sequencing methods. In experimental validation that selectively sequenced microbial DNA from human DNA, ReadCurrent achieved an enrichment ratio of 2.85, which was higher than the 2.7 ratio achieved by MinKNOW using the sequence-alignment strategy. In summary, ReadCurrent can rapidly classify target and nontarget DNA with high accuracy, providing an alternative in the toolbox for nanopore selective sequencing. ReadCurrent is available at https://github.com/Ming-Ni-Group/ReadCurrent.

NASTRA: accurate analysis of short tandem repeat markers by nanopore sequencing with repeat-structure-aware algorithm.

Short-tandem repeats (STRs) are the type of genetic markers extensively utilized in biomedical and forensic applications. Due to sequencing noise in nanopore sequencing, accurate analysis methods are lacking. We developed NASTRA, an innovative tool for Nanopore Autosomal Short Tandem Repeat Analysis, which overcomes traditional database-based methods' limitations and provides a precise germline analysis of STR genetic markers without the need for allele sequence reference. Demonstrating high accuracy in cell line authentication testing and paternity testing, NASTRA significantly surpasses existing methods in both speed and accuracy. This advancement makes it a promising solution for rapid cell line authentication and kinship testing, highlighting the potential of nanopore sequencing for in-field applications.

Replicon-seq: seeing is believing.

Claussin et al. introduce Replicon-seq, a new genome-wide DNA sequencing technology that monitors the progression of individual replisomes at high resolution in vivo.

Beyond sequencing: machine learning algorithms extract biology hidden in Nanopore signal data.

Nanopore sequencing provides signal data corresponding to the nucleotide motifs sequenced. Through machine learning-based methods, these signals are translated into long-read sequences that overcome the read size limit of short-read sequencing. However, analyzing the raw nanopore signal data provides many more opportunities beyond just sequencing genomes and transcriptomes: algorithms that use machine learning approaches to extract biological information from these signals allow the detection of DNA and RNA modifications, the estimation of poly(A) tail length, and the prediction of RNA secondary structures. In this review, we discuss how developments in machine learning methodologies contributed to more accurate basecalling and lower error rates, and how these methods enable new biological discoveries. We argue that direct nanopore sequencing of DNA and RNA provides a new dimensionality for genomics experiments and highlight challenges and future directions for computational approaches to extract the additional information provided by nanopore signal data.

Nanopore sequencing of internal 2'-PO4 modifications installed by RNA repair.

Ligation by plant and fungal RNA ligases yields an internal 2'-phosphate group on each RNA ligation product. In budding yeast, this covalent mark occurs at the splice junction of two targets of ligation: intron-containing tRNAs and the messenger RNA HAC1 The repertoire of RNA molecules repaired by RNA ligation has not been explored due to a lack of unbiased approaches for identifying RNA ligation products. Here, we define several unique signals produced by 2'-phosphorylated RNAs during nanopore sequencing. A 2'-phosphate at the splice junction of HAC1 mRNA inhibits 5' → 3' degradation, enabling detection of decay intermediates in yeast RNA repair mutants by nanopore sequencing. During direct RNA sequencing, intact 2'-phosphorylated RNAs on HAC1 and tRNAs produce diagnostic changes in nanopore current properties and base calling features, including stalls produced as the modified RNA translocates through the nanopore motor protein. These approaches enable directed studies to identify novel RNA repair events in other contexts.

Nanopore-based RNA sequencing deciphers the formation, processing, and modification steps of rRNA intermediates in archaea.

Ribosomal RNA (rRNA) maturation in archaea is a complex multistep process that requires well-defined endo- and exoribonuclease activities to generate fully mature linear rRNAs. However, technical challenges prevented detailed mapping of rRNA processing steps and a systematic analysis of rRNA maturation pathways across the tree of life. In this study, we used long-read (PCR)-cDNA and direct RNA nanopore-based sequencing to study rRNA maturation in three archaeal model organisms, namely the Euryarchaea Haloferax volcanii and Pyrococcus furiosus and the Crenarchaeon Sulfolobus acidocaldarius Compared to standard short-read protocols, nanopore sequencing facilitates simultaneous readout of 5'- and 3'-positions, which is required for the classification of rRNA processing intermediates. More specifically, we (i) accurately detect and describe rRNA maturation stages by analysis of terminal read positions of cDNA reads and thereupon (ii) explore the stage-dependent installation of the KsgA-mediated dimethylations in H. volcanii using base-calling and signal characteristics of direct RNA reads. Due to the single-molecule sequencing capacity of nanopore sequencing, we could detect hitherto unknown intermediates with high confidence, revealing details about the maturation of archaea-specific circular rRNA intermediates. Taken together, our study delineates common principles and unique features of rRNA processing in euryarchaeal and crenarchaeal representatives, thereby significantly expanding our understanding of rRNA maturation pathways in archaea.

Nanoblot: an R-package for visualization of RNA isoforms from long-read RNA-sequencing data.

RT-PCR and northern blots have long been used to study RNA isoforms usage for single genes. Recent advancements in long-read sequencing have yielded unprecedented information about the usage and abundance of these RNA isoforms. However, visualization of long-read sequencing data remains challenging due to the high information density. To alleviate these issues, we have developed NanoBlot, an open-source R-package that generates northern blot and RT-PCR-like images from long-read sequencing data. NanoBlot requires aligned, positionally sorted and indexed BAM files. Plotting is based around ggplot2 and is easily customizable. Advantages of NanoBlot include a robust system for designing probes to visualize isoforms including excluding reads based on the presence or absence of a specified region, an elegant solution to representing isoforms with continuous variations in length, and the ability to overlay multiple genes in the same plot using different colors. We present examples of nanoblots compared to actual northern blot data. In addition to traditional gel-like images, the NanoBlot package can also output other visualizations such as violin plots and 3'-RACE-like plots focused on 3'-end isoforms visualization. The use of the NanoBlot package should provide a simple answer to some of the challenges of visualizing long-read RNA-sequencing data.

NanoDeep: a deep learning framework for nanopore adaptive sampling on microbial sequencing.

Nanopore sequencers can enrich or deplete the targeted DNA molecules in a library by reversing the voltage across individual nanopores. However, it requires substantial computational resources to achieve rapid operations in parallel at read-time sequencing. We present a deep learning framework, NanoDeep, to overcome these limitations by incorporating convolutional neural network and squeeze and excitation. We first showed that the raw squiggle derived from native DNA sequences determines the origin of microbial and human genomes. Then, we demonstrated that NanoDeep successfully classified bacterial reads from the pooled library with human sequence and showed enrichment for bacterial sequence compared with routine nanopore sequencing setting. Further, we showed that NanoDeep improves the sequencing efficiency and preserves the fidelity of bacterial genomes in the mock sample. In addition, NanoDeep performs well in the enrichment of metagenome sequences of gut samples, showing its potential applications in the enrichment of unknown microbiota. Our toolkit is available at https://github.com/lysovosyl/NanoDeep.

Uncovering the dynamics and consequences of RNA isoform changes during neuronal differentiation.

Static gene expression programs have been extensively characterized in stem cells and mature human cells. However, the dynamics of RNA isoform changes upon cell-state-transitions during cell differentiation, the determinants and functional consequences have largely remained unclear. Here, we established an improved model for human neurogenesis in vitro that is amenable for systems-wide analyses of gene expression. Our multi-omics analysis reveals that the pronounced alterations in cell morphology correlate strongly with widespread changes in RNA isoform expression. Our approach identifies thousands of new RNA isoforms that are expressed at distinct differentiation stages. RNA isoforms mainly arise from exon skipping and the alternative usage of transcription start and polyadenylation sites during human neurogenesis. The transcript isoform changes can remodel the identity and functions of protein isoforms. Finally, our study identifies a set of RNA binding proteins as a potential determinant of differentiation stage-specific global isoform changes. This work supports the view of regulated isoform changes that underlie state-transitions during neurogenesis.

Analyzing RNA posttranscriptional modifications to decipher the epitranscriptomic code.

The discovery of RNA silencing has revealed that non-protein-coding sequences (ncRNAs) can cover essential roles in regulatory networks and their malfunction may result in severe consequences on human health. These findings have prompted a general reassessment of the significance of RNA as a key player in cellular processes. This reassessment, however, will not be complete without a greater understanding of the distribution and function of the over 170 variants of the canonical ribonucleotides, which contribute to the breathtaking structural diversity of natural RNA. This review surveys the analytical approaches employed for the identification, characterization, and detection of RNA posttranscriptional modifications (rPTMs). The merits of analyzing individual units after exhaustive hydrolysis of the initial biopolymer are outlined together with those of identifying their position in the sequence of parent strands. Approaches based on next generation sequencing and mass spectrometry technologies are covered in depth to provide a comprehensive view of their respective merits. Deciphering the epitranscriptomic code will require not only mapping the location of rPTMs in the various classes of RNAs, but also assessing the variations of expression levels under different experimental conditions. The fact that no individual platform is currently capable of meeting all such demands implies that it will be essential to capitalize on complementary approaches to obtain the desired information. For this reason, the review strived to cover the broadest possible range of techniques to provide readers with the fundamental elements necessary to make informed choices and design the most effective possible strategy to accomplish the task at hand.

Accurate prenatal diagnosis of facioscapulohumeral muscular dystrophy 1 using nanopore sequencing.

BACKGROUND: Facioscapulohumeral muscular dystrophy 1 (FSHD1) is an autosomal dominant muscular disorder mainly caused by the contraction and hypomethylation of the D4Z4 repeat array in chromosome 4q35. Prenatal diagnosis of FSHD1 is challenging due to the highly repetitive and long genomic structure. In this study, a pregnant woman diagnosed with FSHD1 using optical genome mapping sought assistance for a healthy offspring. METHODS: At the 17th week of gestation, she underwent amniocentesis, and genomic DNA (gDNA) was extracted from amniocytes. Whole-genome sequencing of the gDNA was performed using the nanopore MinION platform. RESULTS: Despite a sequencing depth of only 7.3×, bioinformatic analyses revealed that the fetus inherited four D4Z4 repeat units with the permissive 4qA from the mother and the eight D4Z4 repeat units with the non-permissive 4qB from the father. To validate the results, SNP-based linkage analyses were conducted with gDNA from the proband, the proband's father and proband's amniocytes. Results indicated that the fetus inherited the maternal pathogenic haplotype based on 144 informative SNPs. Linkage analysis was consistent with the nanopore sequencing. CONCLUSION: Nanopore sequencing proves to be an accurate and direct method for genetic testing of monogenic diseases at the single-nucleotide level. This study represents the first application of nanopore sequencing in the prenatal diagnosis of FSHD1, providing a significant advantage for patients with de novo mutations.