Improved Nanopore full-length cDNA sequencing by PCR-suppression.

Full-length transcript sequencing remains a main goal of RNA sequencing. However, even the application of long-read sequencing technologies such as Oxford Nanopore Technologies still fail to yield full-length transcript sequencing for a significant portion of sequenced reads. Since these technologies can sequence reads that are far longer than the longest known processed transcripts, the lack of efficiency to obtain full-length transcripts from good quality RNAs stems from library preparation inefficiency rather than the presence of degraded RNA molecules. It has previously been shown that addition of inverted terminal repeats in cDNA during reverse transcription followed by single-primer PCR creates a PCR suppression effect that prevents amplification of short molecules thus enriching the library for longer transcripts. We adapted this method for Nanopore cDNA library preparation and show that not only is PCR efficiency increased but gene body coverage is dramatically improved. The results show that implementation of this simple strategy will result in better quality full-length RNA sequencing data and make full-length transcript sequencing possible for most of sequenced reads.

Evolution of naturally arising SARS-CoV-2 defective interfering particles.

Defective interfering (DI) particles arise during virus propagation, are conditional on parental virus for replication and packaging, and interfere with viral expansion. There is much interest in developing DIs as anti-viral agents. Here we characterize DI particles that arose following serial passaging of SARS-CoV-2 at high multiplicity of infection. The prominent DIs identified have lost ~84% of the SARS-CoV-2 genome and are capable of attenuating parental viral titers. Synthetic variants of the DI genomes also interfere with infection and can be used as conditional, gene delivery vehicles. In addition, the DI genomes encode an Nsp1-10 fusion protein capable of attenuating viral replication. These results identify naturally selected defective viral genomes that emerged and stably propagated in the presence of parental virus.

Nanopore long-read RNA-seq and absolute quantification delineate transcription dynamics in early embryo development of an insect pest.

The olive fruit fly, Bactrocera oleae, is the most important pest for the olive fruit but lacks adequate transcriptomic characterization that could aid in molecular control approaches. We apply nanopore long-read RNA-seq with internal RNA standards allowing absolute transcript quantification to analyze transcription dynamics during early embryo development for the first time in this organism. Sequencing on the MinION platform generated over 31 million reads. Over 50% of the expressed genes had at least one read covering its entire length validating our full-length approach. We generated a de novo transcriptome assembly and identified 1768 new genes and a total of 79,810 isoforms; a fourfold increase in transcriptome diversity compared to the current NCBI predicted transcriptome. Absolute transcript quantification per embryo allowed an insight into the dramatic re-organization of maternal transcripts. We further identified Zelda as a possible regulator of early zygotic genome activation in B. oleae and provide further insights into the maternal-to-zygotic transition. These data show the utility of long-read RNA in improving characterization of non-model organisms that lack a fully annotated genome, provide potential targets for sterile insect technic approaches, and provide the first insight into the transcriptome landscape of the developing olive fruit fly embryo.

Publisher Correction: Regulation of cellular sterol homeostasis by the oxygen responsive noncoding RNA lincNORS.

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

Regulation of cellular sterol homeostasis by the oxygen responsive noncoding RNA lincNORS.

We hereby provide the initial portrait of lincNORS, a spliced lincRNA generated by the MIR193BHG locus, entirely distinct from the previously described miR-193b-365a tandem. While inducible by low O2 in a variety of cells and associated with hypoxia in vivo, our studies show that lincNORS is subject to multiple regulatory inputs, including estrogen signals. Biochemically, this lincRNA fine-tunes cellular sterol/steroid biosynthesis by repressing the expression of multiple pathway components. Mechanistically, the function of lincNORS requires the presence of RALY, an RNA-binding protein recently found to be implicated in cholesterol homeostasis. We also noticed the proximity between this locus and naturally occurring genetic variations highly significant for sterol/steroid-related phenotypes, in particular the age of sexual maturation. An integrative analysis of these variants provided a more formal link between these phenotypes and lincNORS, further strengthening the case for its biological relevance.

Methodologies for Transcript Profiling Using Long-Read Technologies.

RNA sequencing using next-generation sequencing technologies (NGS) is currently the standard approach for gene expression profiling, particularly for large-scale high-throughput studies. NGS technologies comprise high throughput, cost efficient short-read RNA-Seq, while emerging single molecule, long-read RNA-Seq technologies have enabled new approaches to study the transcriptome and its function. The emerging single molecule, long-read technologies are currently commercially available by Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT), while new methodologies based on short-read sequencing approaches are also being developed in order to provide long range single molecule level information-for example, the ones represented by the 10x Genomics linked read methodology. The shift toward long-read sequencing technologies for transcriptome characterization is based on current increases in throughput and decreases in cost, making these attractive for de novo transcriptome assembly, isoform expression quantification, and in-depth RNA species analysis. These types of analyses were challenging with standard short sequencing approaches, due to the complex nature of the transcriptome, which consists of variable lengths of transcripts and multiple alternatively spliced isoforms for most genes, as well as the high sequence similarity of highly abundant species of RNA, such as rRNAs. Here we aim to focus on single molecule level sequencing technologies and single-cell technologies that, combined with perturbation tools, allow the analysis of complete RNA species, whether short or long, at high resolution. In parallel, these tools have opened new ways in understanding gene functions at the tissue, network, and pathway levels, as well as their detailed functional characterization. Analysis of the epi-transcriptome, including RNA methylation and modification and the effects of such modifications on biological systems is now enabled through direct RNA sequencing instead of classical indirect approaches. However, many difficulties and challenges remain, such as methodologies to generate full-length RNA or cDNA libraries from all different species of RNAs, not only poly-A containing transcripts, and the identification of allele-specific transcripts due to current error rates of single molecule technologies, while the bioinformatics analysis on long-read data for accurate identification of 5' and 3' UTRs is still in development.

Transcript Profiling Using Long-Read Sequencing Technologies.

RNA sequencing using next-generation sequencing (NGS, RNA-Seq) technologies is currently the standard approach for gene expression profiling, particularly for large-scale high-throughput studies. NGS technologies comprise short-read RNA-Seq (dominated by Illumina) and long-read RNA-Seq technologies provided by Pacific Bioscience (PacBio) and Oxford Nanopore Technologies (ONT). Although short-read sequencing technologies are the most widely used, long-read technologies are increasingly becoming the standard approach for de novo transcriptome assembly and isoform expression quantification due to the complex nature of the transcriptome which consists of variable lengths of transcripts and multiple alternatively spliced isoforms for most genes. In this chapter, we describe experimental procedures for library preparation, sequencing, and associated data analysis approaches for PacBio and ONT with a major focus on full length cDNA synthesis, de novo transcriptome assembly, and isoform quantification.

Current and Future Methods for mRNA Analysis: A Drive Toward Single Molecule Sequencing.

The transcriptome encompasses a range of species including messenger RNA, and other noncoding RNA such as rRNA, tRNA, and short and long noncoding RNAs. Due to the huge role played by mRNA in development and disease, several methods have been developed to sequence and characterize mRNA, with RNA sequencing (RNA-Seq) emerging as the current method of choice particularly for large high-throughput studies. Short-read RNA-Seq which involves sequencing of short cDNA fragments and computationally assembling them to reconstruct the transcriptome, or aligning them to a reference is the most widely used approach. However, due to inherent limitations of this approach in de novo transcriptome assembly and isoform quantification, long-read RNA-Seq approaches, which also happen to be single molecule sequencing approaches, are increasingly becoming the standard for de novo transcriptome assembly and isoform quantification. In this chapter, we review the technical aspects of the current methods of RNA-Seq, both short and long-read approaches, and data analysis methods available. We discuss recent advances in single-cell RNA-Seq and direct RNA-Seq approaches, which perhaps will dominate the future of RNA-Seq.

Benchmarking of the Oxford Nanopore MinION sequencing for quantitative and qualitative assessment of cDNA populations.

To assess the performance of the Oxford Nanopore Technologies MinION sequencing platform, cDNAs from the External RNA Controls Consortium (ERCC) RNA Spike-In mix were sequenced. This mix mimics mammalian mRNA species and consists of 92 polyadenylated transcripts with known concentration. cDNA libraries were generated using a template switching protocol to facilitate the direct comparison between different sequencing platforms. The MinION performance was assessed for its ability to sequence the cDNAs directly with good accuracy in terms of abundance and full length. The abundance of the ERCC cDNA molecules sequenced by MinION agreed with their expected concentration. No length or GC content bias was observed. The majority of cDNAs were sequenced as full length. Additionally, a complex cDNA population derived from a human HEK-293 cell line was sequenced on an Illumina HiSeq 2500, PacBio RS II and ONT MinION platforms. We observed that there was a good agreement in the measured cDNA abundance between PacBio RS II and ONT MinION (rpearson = 0.82, isoforms with length more than 700bp) and between Illumina HiSeq 2500 and ONT MinION (rpearson = 0.75). This indicates that the ONT MinION can sequence quantitatively both long and short full length cDNA molecules.

IRF5:RelA interaction targets inflammatory genes in macrophages.

Interferon Regulatory Factor 5 (IRF5) plays a major role in setting up an inflammatory macrophage phenotype, but the molecular basis of its transcriptional activity is not fully understood. In this study, we conduct a comprehensive genome-wide analysis of IRF5 recruitment in macrophages stimulated with bacterial lipopolysaccharide and discover that IRF5 binds to regulatory elements of highly transcribed genes. Analysis of protein:DNA microarrays demonstrates that IRF5 recognizes the canonical IRF-binding (interferon-stimulated response element [ISRE]) motif in vitro. However, IRF5 binding in vivo appears to rely on its interactions with other proteins. IRF5 binds to a noncanonical composite PU.1:ISRE motif, and its recruitment is aided by RelA. Global gene expression analysis in macrophages deficient in IRF5 and RelA highlights the direct role of the RelA:IRF5 cistrome in regulation of a subset of key inflammatory genes. We map the RelA:IRF5 interaction domain and suggest that interfering with it would offer selective targeting of macrophage inflammatory activities.

Extensive regulation of the non-coding transcriptome by hypoxia: role of HIF in releasing paused RNApol2.

Hypoxia is central to both ischaemic and neoplastic diseases. However, the non-coding transcriptional response to hypoxia is largely uncharacterized. We undertook integrated genomic analyses of both non-coding and coding transcripts using massively parallel sequencing and interfaced this data with pan-genomic analyses of hypoxia-inducible factor (HIF) and RNApol2 binding in hypoxic cells. These analyses revealed that all classes of RNA are profoundly regulated by hypoxia and implicated HIF as a major direct regulator of both the non-coding and coding transcriptome, acting predominantly through release of pre-bound promoter-paused RNApol2. These findings indicate that the transcriptional response to hypoxia is substantially more extensive than previously considered.

TGF-ß/Smad2/3 signaling directly regulates several miRNAs in mouse ES cells and early embryos.

The Transforming Growth Factor-ß (TGF-ß) signaling pathway is one of the major pathways essential for normal embryonic development and tissue homeostasis, with anti-tumor but also pro-metastatic properties in cancer. This pathway directly regulates several target genes that mediate its downstream functions, however very few microRNAs (miRNAs) have been identified as targets. miRNAs are modulators of gene expression with essential roles in development and a clear association with diseases including cancer. Little is known about the transcriptional regulation of the primary transcripts (pri-miRNA, pri-miR) from which several mature miRNAs are often derived. Here we present the identification of miRNAs regulated by TGF-ß signaling in mouse embryonic stem (ES) cells and early embryos. We used an inducible ES cell system to maintain high levels of the TGF-ß activated/phosphorylated Smad2/3 effectors, which are the transcription factors of the pathway, and a specific inhibitor that blocks their activation. By performing short RNA deep-sequencing after 12 hours Smad2/3 activation and after 16 hours inhibition, we generated a database of responsive miRNAs. Promoter/enhancer analysis of a subset of these miRNAs revealed that the transcription of pri-miR-181c/d and the pri-miR-341∼3072 cluster were found to depend on activated Smad2/3. Several of these miRNAs are expressed in early mouse embryos, when the pathway is known to play an essential role. Treatment of embryos with TGF-ß inhibitor caused a reduction of their levels confirming that they are targets of this pathway in vivo. Furthermore, we showed that pri-miR-341∼3072 transcription also depends on FoxH1, a known Smad2/3 transcription partner during early development. Together, our data show that miRNAs are regulated directly by the TGF-ß/Smad2/3 pathway in ES cells and early embryos. As somatic abnormalities in functions known to be regulated by the TGF-ß/Smad2/3 pathway underlie tumor suppression and metastasis, this research also provides a resource for miRNAs involved in cancer.

Extensive characterization of NF-κB binding uncovers non-canonical motifs and advances the interpretation of genetic functional traits.

BACKGROUND: Genetic studies have provided ample evidence of the influence of non-coding DNA polymorphisms on trait variance, particularly those occurring within transcription factor binding sites. Protein binding microarrays and other platforms that can map these sites with great precision have enhanced our understanding of how a single nucleotide polymorphism can alter binding potential within an in vitro setting, allowing for greater predictive capability of its effect on a transcription factor binding site. RESULTS: We have used protein binding microarrays and electrophoretic mobility shift assay-sequencing (EMSA-Seq), a deep sequencing based method we developed to analyze nine distinct human NF-κB dimers. This family of transcription factors is one of the most extensively studied, but our understanding of its DNA binding preferences has been limited to the originally described consensus motif, GGRRNNYYCC. We highlight differences between NF-κB family members and also put under the spotlight non-canonical motifs that have so far received little attention. We utilize our data to interpret the binding of transcription factors between individuals across 1,405 genomic regions laden with single nucleotide polymorphisms. We also associated binding correlations made using our data with risk alleles of disease and demonstrate its utility as a tool for functional studies of single nucleotide polymorphisms in regulatory regions. CONCLUSIONS: NF-κB dimers bind specifically to non-canonical motifs and these can be found within genomic regions in which a canonical motif is not evident. Binding affinity data generated with these different motifs can be used in conjunction with data from chromatin immunoprecipitation-sequencing (ChIP-Seq) to enable allele-specific analyses of expression and transcription factor-DNA interactions on a genome-wide scale.

High-resolution genome-wide mapping of HIF-binding sites by ChIP-seq.

Hypoxia-inducible factor (HIF) regulates the major transcriptional cascade central to the response of all mammalian cells to alterations in oxygen tension. Expression arrays indicate that many hundreds of genes are regulated by this pathway, controlling diverse processes that in turn orchestrate both oxygen delivery and utilization. However, the extent to which HIF exerts direct versus indirect control over gene expression together with the factors dictating the range of HIF-regulated genes remains unclear. Using chromatin immunoprecipitation linked to high throughput sequencing, we identify HIF-binding sites across the genome, independently of gene architecture. Using gene set enrichment analysis, we demonstrate robust associations with the regulation of gene expression by HIF, indicating that these sites operate over long genomic intervals. Analysis of HIF-binding motifs demonstrates sequence preferences outside of the core RCGTG-binding motif but does not reveal any additional absolute sequence requirements. Across the entire genome, only a small proportion of these potential binding sites are bound by HIF, although occupancy of potential sites was enhanced approximately 20-fold at normoxic DNAse1 hypersensitivity sites (irrespective of distance from promoters), suggesting that epigenetic regulation of chromatin may have an important role in defining the response to hypoxia.