A hybridization-based approach for quantitative and low-bias single-stranded DNA ligation

Rapid discovery of functional RNA domains

Many strategies have been implemented to enrich an RNA population for a selectable function, but demarcation of the optimal functional motifs or minimal structures within longer libraries remains a lengthy and tedious process. To overcome this problem, we have developed a technique that isolates minimal active segments from complex heterogeneous pools of RNAs. This method allows for truncations to occur at both 5' and 3' ends of functional domains and introduces independent primer-binding sequences, thereby removing sequence and structure bias introduced by constant-sequence regions. We show examples of minimization for genomic and synthetic aptamers and demonstrate that the method can directly reveal an active RNA assembled from multiple strands, facilitating the development of heterodimeric structures used in cellular sensors. This approach provides a pipeline to experimentally define the boundaries of active domains and accelerate the discovery of functional RNAs.

A review on the impact of single-stranded library preparation on plasma cell-free diversity for cancer detection

In clinical oncology, cell-free DNA (cfDNA) has shown immense potential in its ability to noninvasively detect cancer at various stages and monitor the progression of therapy. Despite the rapid improvements in cfDNA liquid biopsy approaches, achieving the required sensitivity to detect rare tumor-derived cfDNA still remains a challenge. For next-generation sequencing, the perceived presentation of cfDNA is strongly linked to the extraction and library preparation protocols. Conventional double-stranded DNA library preparation (dsDNA-LP) focuses on assessing ~167bp double-stranded mononucleosomal (mncfDNA) and its other oligonucleosomal cell-free DNA counterparts in plasma. However, dsDNA-LP methods fail to include short, single-stranded, or nicked DNA in the final library preparation, biasing the representation of the actual cfDNA populations in plasma. The emergence of single-stranded library preparation (ssDNA-LP) strategies over the past decade has now allowed these other populations of cfDNA to be studied from plasma. With the use of ssDNA-LP, single-stranded, nicked, and ultrashort cfDNA can be comprehensively assessed for its molecular characteristics and clinical potential. In this review, we overview the current literature on applications of ssDNA-LP on plasma cfDNA from a potential cancer liquid biopsy perspective. To this end, we discuss the molecular principles of single-stranded DNA adapter ligation, how library preparation contributes to the understanding of native cfDNA characteristics, and the potential for ssDNA-LP to improve the sensitivity of circulating tumor DNA detection. Additionally, we review the current literature on the newly reported species of plasma ultrashort single-stranded cell-free DNA plasma, which appear biologically distinct from mncfDNA. We conclude with a discussion of future perspectives of ssDNA-LP for liquid biopsy endeavors.

Observation of coordinated RNA folding events by systematic cotranscriptional RNA structure probing

RNA begins to fold as it is transcribed by an RNA polymerase. Consequently, RNA folding is constrained by the direction and rate of transcription. Understanding how RNA folds into secondary and tertiary structures therefore requires methods for determining the structure of cotranscriptional folding intermediates. Cotranscriptional RNA chemical probing methods accomplish this by systematically probing the structure of nascent RNA that is displayed from an RNA polymerase. Here, we describe a concise, high-resolution cotranscriptional RNA chemical probing procedure called variable length Transcription Elongation Complex RNA structure probing (TECprobe-VL). We demonstrate the accuracy and resolution of TECprobe-VL by replicating and extending previous analyses of ZTP and fluoride riboswitch folding and mapping the folding pathway of a ppGpp-sensing riboswitch. In each system, we show that TECprobe-VL identifies coordinated cotranscriptional folding events that mediate transcription antitermination. Our findings establish TECprobe-VL as an accessible method for mapping cotranscriptional RNA folding pathways.

In vivo structure probing of RNA in Archaea: novel insights into the ribosome structure of Methanosarcina acetivorans

Structure probing combined with next-generation sequencing (NGS) has provided novel insights into RNA structure-function relationships. To date, such studies have focused largely on bacteria and eukaryotes, with little attention given to the third domain of life, archaea. Furthermore, functional RNAs have not been extensively studied in archaea, leaving open questions about RNA structure and function within this domain of life. With archaeal species being diverse and having many similarities to both bacteria and eukaryotes, the archaea domain has the potential to be an evolutionary bridge. In this study, we introduce a method for probing RNA structure in vivo in the archaea domain of life. We investigated the structure of ribosomal RNA (rRNA) from Methanosarcina acetivorans, a well-studied anaerobic archaeal species, grown with either methanol or acetate. After probing the RNA in vivo with dimethyl sulfate (DMS), Structure-seq2 libraries were generated, sequenced, and analyzed. We mapped the reactivity of DMS onto the secondary structure of the ribosome, which we determined independently with comparative analysis, and confirmed the accuracy of DMS probing in M. acetivorans Accessibility of the rRNA to DMS in the two carbon sources was found to be quite similar, although some differences were found. Overall, this study establishes the Structure-seq2 pipeline in the archaea domain of life and informs about ribosomal structure within M. acetivorans.

A critical spotlight on the paradigms of FFPE-DNA sequencing

In the late 19th century, formalin fixation with paraffin-embedding (FFPE) of tissues was developed as a fixation and conservation method and is still used to this day in routine clinical and pathological practice. The implementation of state-of-the-art nucleic acid sequencing technologies has sparked much interest for using historical FFPE samples stored in biobanks as they hold promise in extracting new information from these valuable samples. However, formalin fixation chemically modifies DNA, which potentially leads to incorrect sequences or misinterpretations in downstream processing and data analysis. Many publications have concentrated on one type of DNA damage, but few have addressed the complete spectrum of FFPE-DNA damage. Here, we review mitigation strategies in (I) pre-analytical sample quality control, (II) DNA repair treatments, (III) analytical sample preparation and (IV) bioinformatic analysis of FFPE-DNA. We then provide recommendations that are tested and illustrated with DNA from 13-year-old liver specimens, one FFPE preserved and one fresh frozen, applying target-enriched sequencing. Thus, we show how DNA damage can be compensated, even when using low quantities (50 ng) of fragmented FFPE-DNA (DNA integrity number 2.0) that cannot be amplified well (Q129 bp/Q41 bp = 5%). Finally, we provide a checklist called 'ERROR-FFPE-DNA' that summarises recommendations for the minimal information in publications required for assessing fitness-for-purpose and inter-study comparison when using FFPE samples.

Genome-wide measurement of DNA replication fork directionality and quantification of DNA replication initiation and termination with Okazaki fragment sequencing

Studying the dynamics of genome replication in mammalian cells has been historically challenging. To reveal the location of replication initiation and termination in the human genome, we developed Okazaki fragment sequencing (OK-seq), a quantitative approach based on the isolation and strand-specific sequencing of Okazaki fragments, the lagging strand replication intermediates. OK-seq quantitates the proportion of leftward- and rightward-oriented forks at every genomic locus and reveals the location and efficiency of replication initiation and termination events. Here we provide the detailed experimental procedures for performing OK-seq in unperturbed cultured human cells and budding yeast and the bioinformatics pipelines for data processing and computation of replication fork directionality. Furthermore, we present the analytical approach based on a hidden Markov model, which allows automated detection of ascending, descending and flat replication fork directionality segments revealing the zones of replication initiation, termination and unidirectional fork movement across the entire genome. These tools are essential for the accurate interpretation of human and yeast replication programs. The experiments and the data processing can be accomplished within six days. Besides revealing the genome replication program in fine detail, OK-seq has been instrumental in numerous studies unravelling mechanisms of genome stability, epigenome maintenance and genome evolution.

Pseudouridylation of chloroplast ribosomal RNA contributes to low temperature acclimation in rice

Ribosomal RNAs (rRNAs) undergo many modifications during transcription and maturation; homeostasis of rRNA modifications is essential for chloroplast biogenesis in plants. The chloroplast acts as a hub to sense environmental signals, such as cold temperature. However, how RNA modifications contribute to low temperature responses remains unknown. Here we reveal that pseudouridine (Ψ) modification of rice chloroplast rRNAs mediated by the pseudouridine synthase (OsPUS1) contributes to cold tolerance at seedling stage. Loss-function of OsPUS1 leads to abnormal chloroplast development and albino seedling phenotype at low temperature. We find that OsPUS1 is accumulated upon cold and binds to chloroplast precursor rRNAs (pre-rRNAs) to catalyse the pseudouridylation on rRNA. These modifications on chloroplast rRNAs could be required for their processing, as the reduction of mature chloroplast rRNAs and accumulation of pre-rRNAs are observed in ospus1-1 at low temperature. Therefore, the ribosome activity and translation in chloroplasts is disturbed in ospus1-1. Furthermore, transcriptome and translatome analysis reveals that OsPUS1 balances growth and stress-responsive state, preventing excess reactive oxygen species accumulation. Taken together, our findings unveil a crucial function of Ψ in chloroplast ribosome biogenesis and cold tolerance in rice, with potential applications in crop improvement.

Enhanced transcriptome-wide RNA G-quadruplex sequencing for low RNA input samples with rG4-seq 2.0

Background

RNA G-quadruplexes (rG4s) are non-canonical structural motifs that have diverse functional and regulatory roles, for instance in transcription termination, alternative splicing, mRNA localization and stabilization, and translational process. We recently developed the RNA G-quadruplex structure sequencing (rG4-seq) technique and described rG4s in both eukaryotic and prokaryotic transcriptomes. However, rG4-seq suffers from a complicated gel purification step and limited PCR product yield, thus requiring a high amount of RNA input, which limits its applicability in more physiologically or clinically relevant studies often characterized by the limited availability of biological material and low RNA abundance. Here, we redesign and enhance the workflow of rG4-seq to address this issue.

Results

We developed rG4-seq 2.0 by introducing a new ssDNA adapter containing deoxyuridine during library preparation to enhance library quality with no gel purification step, less PCR amplification cycles and higher yield of PCR products. We demonstrate that rG4-seq 2.0 produces high-quality cDNA libraries that support reliable and reproducible rG4 identification at varying RNA inputs, including RNA mounts as low as 10 ng. rG4-seq 2.0 also improved the rG4-seq calling outcome and nucleotide bias in rG4 detection persistent in rG4-seq 1.0. We further provide in vitro mapping of rG4 in the HEK293T cell line, and recommendations for assessing RNA input and sequencing depth for individual rG4 studies based on transcript abundance.

Conclusions

rG4-seq 2.0 can improve the identification and study of rG4s in low abundance transcripts, and our findings can provide insights to optimize cDNA library preparation in other related methods.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01448-3.

Improvements to Hybridization-Ligation Enzyme-Linked Immunosorbent Assay Methods to Overcome Bioanalytical Challenges Posed by Novel Oligonucleotide Therapeutics

As oligonucleotides (ONs) and similar nucleic acid therapeutic modalities enter development pipelines, there is continual need to develop bioanalytical methodologies addressing unique challenges they pose. Novel ONs back bone chemistries, especially those enabling stereochemical control, and base modifications are being exploited to improve pharmacological properties, potency, and increase half-lives. These changes have strained established methods, oftentimes precluding development of assays sensitive and specific enough to meet the needs of preclinical programs. For stereopure ONs representing a single molecular species, nontrivial presence of chain-shortened metabolites in biological samples necessitate assays with high specificity. To meet these needs, this report presents a toolbox of novel techniques, easy to implement for existing hybridization-ligation enzyme-linked immunosorbent assay formats, which address this challenge and yield significant sensitivity and specificity enhancements. Ligation efficiency was improved up to 61-fold through addition of polyethylene glycol, betaine, or dimethylsulfoxide, mitigating major differences among sequence-matched ONs of varying stereopurity, enabling sensitivities below 0.100 ng/mL for quantitation. These improvements enabled further refinement of capture probe designs engendering sufficient specificity to discriminate N-1 chain-shortened metabolites at both the 5′ and 3′ end of the ONs. These generalizable methods advance the performance of mainstay bioanalytical assays, facilitating research and development of innovative ONs therapeutics.

Identifying transcript 5' capped ends in Plasmodium falciparum

Background

The genome of the human malaria parasite Plasmodium falciparum is poorly annotated, in particular, the 5' capped ends of its mRNA transcripts. New approaches are needed to fully catalog P. falciparum transcripts for understanding gene function and regulation in this organism.

Methods

We developed a transcriptomic method based on next-generation sequencing of complementary DNA (cDNA) enriched for full-length fragments using eIF4E, a 5' cap-binding protein, and an unenriched control. DNA sequencing adapter was added after enrichment of full-length cDNA using two different ligation protocols. From the mapped sequence reads, enrichment scores were calculated for all transcribed nucleotides and used to calculate P-values of 5' capped nucleotide enrichment. Sensitivity and accuracy were increased by combining P-values from replicate experiments. Data were obtained for P. falciparum ring, trophozoite and schizont stages of intra-erythrocytic development.

Results

5' capped nucleotide signals were mapped to 17,961 non-overlapping P. falciparum genomic intervals. Analysis of the dominant 5' capped nucleotide in these genomic intervals revealed the presence of two groups with distinctive epigenetic features and sequence patterns. A total of 4,512 transcripts were annotated as 5' capped based on the correspondence of 5' end with 5' capped nucleotide annotated from full-length cDNA data.

Discussion

The presence of two groups of 5' capped nucleotides suggests that alternative mechanisms may exist for producing 5' capped transcript ends in P. falciparum. The 5' capped transcripts that are antisense, outside of, or partially overlapping coding regions may be important regulators of gene function in P. falciparum.

A Fast and Efficient Single-stranded Genomic Library Preparation Method Optimized for Ancient DNA

We present a protocol to prepare extracted DNA for sequencing on the Illumina sequencing platform that has been optimized for ancient and degraded DNA. Our approach, the Santa Cruz Reaction or SCR, uses directional splinted ligation of Illumina’s P5 and P7 adapters to convert natively single-stranded DNA and heat denatured double-stranded DNA into sequencing libraries in a single enzymatic reaction. To demonstrate its efficacy in converting degraded DNA molecules, we prepare 5 ancient DNA extracts into sequencing libraries using the SCR and 2 of the most commonly used approaches for preparing degraded DNA for sequencing: BEST, which targets and converts double-stranded DNA, and ssDNA2.0, which targets and converts single-stranded DNA. We then compare the efficiency with which each approach recovers unique molecules, or library complexity, given a standard amount of DNA input. We find that the SCR consistently outperforms the BEST protocol in recovering unique molecules and, despite its relative simplicity to perform and low cost per library, has similar performance to ssDNA2.0 across a wide range of DNA inputs. The SCR is a cost- and time-efficient approach that minimizes the loss of unique molecules and makes accessible a taxonomically, geographically, and a temporally broader sample of preserved remains for genomic analysis.

Novel target capture DNA library preparation method using CircLigase-mediated hook ligation

Target capture DNA library preparation methods primarily use the biotin-streptavidin interaction to enrich target DNA, requiring complex operations and a time-consuming workflow with a series of wash steps at regulated temperatures. Here, a new method is presented termed 'hook ligation' for target DNA library preparation, using CircLigase to capture target DNA. The concept was validated by library preparation, sequencing, and acquisition of target DNA fragment information through bioinformatics analysis. An efficient reference point for single-strand DNA ligation to a hook sequence using CircLigase is also provided and it is shown that the efficiency can be influenced by the length and the position of the complementary area on the target DNA.

Anomalous Reverse Transcription through Chemical Modifications in Polyadenosine Stretches

Thermostable reverse transcriptases are workhorse enzymes underlying nearly all modern techniques for RNA structure mapping and for the transcriptome-wide discovery of RNA chemical modifications. Despite their wide use, these enzymes' behaviors at chemical modified nucleotides remain poorly understood. Wellington-Oguri et al. recently reported an apparent loss of chemical modification within putatively unstructured polyadenosine stretches modified by dimethyl sulfate or 2' hydroxyl acylation, as probed by reverse transcription. Here, reanalysis of these and other publicly available data, capillary electrophoresis experiments on chemically modified RNAs, and nuclear magnetic resonance spectroscopy on (A)₁₂ and variants show that this effect is unlikely to arise from an unusual structure of polyadenosine. Instead, tests of different reverse transcriptases on chemically modified RNAs and molecules synthesized with single 1-methyladenosines implicate a previously uncharacterized reverse transcriptase behavior: near-quantitative bypass through chemical modifications within polyadenosine stretches. All tested natural and engineered reverse transcriptases (MMLV; SuperScript II, III, and IV; TGIRT-III; and MarathonRT) exhibit this anomalous bypass behavior. Accurate DMS-guided structure modeling of the polyadenylated HIV-1 3' untranslated region requires taking into account this anomaly. Our results suggest that poly(rA-dT) hybrid duplexes can trigger an unexpectedly effective reverse transcriptase bypass and that chemical modifications in mRNA poly(A) tails may be generally undercounted.

Probing RNA structure in vivo

RNA structure underpins many essential functions in biology. New chemical reagents and techniques for probing RNA structure in living cells have emerged in recent years. High-throughput, genome-wide techniques such as Structure-seq2 and DMS-MaPseq exploit nucleobase modification by dimethylsulfate (DMS) to obtain complete structuromes, and are applicable to multiple domains of life and conditions. New reagents such as 1-ethyl-3-(3-dimethylaminopropyl)carbodiimide (EDC), glyoxal, and nicotinoyl azide (NAz) greatly expand the capabilities of nucleobase probing in cells. Additionally, ribose-targeting reagents in selective 2'-hydroxyl acylation and primer extension (SHAPE) detect RNA flexibility in vivo. These techniques, coupled with crosslinking nucleobases in psoralen analysis of RNA interactions and structures (PARIS), provide new and diverse ways to elucidate RNA secondary and tertiary structure in vivo and genome-wide.

Personalized circulating tumor DNA analysis to detect residual disease after neoadjuvant therapy in breast cancer

Longitudinal analysis of circulating tumor DNA (ctDNA) has shown promise for monitoring treatment response. However, most current methods lack adequate sensitivity for residual disease detection during or after completion of treatment in patients with nonmetastatic cancer. To address this gap and to improve sensitivity for minute quantities of residual tumor DNA in plasma, we have developed targeted digital sequencing (TARDIS) for multiplexed analysis of patient-specific cancer mutations. In reference samples, by simultaneously analyzing 8 to 16 known mutations, TARDIS achieved 91 and 53% sensitivity at mutant allele fractions (AFs) of 3 in 104 and 3 in 105, respectively, with 96% specificity, using input DNA equivalent to a single tube of blood. We successfully analyzed up to 115 mutations per patient in 80 plasma samples from 33 women with stage I to III breast cancer. Before treatment, TARDIS detected ctDNA in all patients with 0.11% median AF. After completion of neoadjuvant therapy, ctDNA concentrations were lower in patients who achieved pathological complete response (pathCR) compared to patients with residual disease (median AFs, 0.003 and 0.017%, respectively, P = 0.0057, AUC = 0.83). In addition, patients with pathCR showed a larger decrease in ctDNA concentrations during neoadjuvant therapy. These results demonstrate high accuracy for assessment of molecular response and residual disease during neoadjuvant therapy using ctDNA analysis. TARDIS has achieved up to 100-fold improvement beyond the current limit of ctDNA detection using clinically relevant blood volumes, demonstrating that personalized ctDNA tracking could enable individualized clinical management of patients with cancer treated with curative intent.

High-throughput determination of RNA structures

RNA performs and regulates a diverse range of cellular processes, with new functional roles being uncovered at a rapid pace. Interest is growing in how these functions are linked to RNA structures that form in the complex cellular environment. A growing suite of technologies that use advances in RNA structural probes, high-throughput sequencing and new computational approaches to interrogate RNA structure at unprecedented throughput are beginning to provide insights into RNA structures at new spatial, temporal and cellular scales.

Terminator-free template-independent enzymatic DNA synthesis for digital information storage

DNA is an emerging medium for digital data and its adoption can be accelerated by synthesis processes specialized for storage applications. Here, we describe a de novo enzymatic synthesis strategy designed for data storage which harnesses the template-independent polymerase terminal deoxynucleotidyl transferase (TdT) in kinetically controlled conditions. Information is stored in transitions between non-identical nucleotides of DNA strands. To produce strands representing user-defined content, nucleotide substrates are added iteratively, yielding short homopolymeric extensions whose lengths are controlled by apyrase-mediated substrate degradation. With this scheme, we synthesize DNA strands carrying 144 bits, including addressing, and demonstrate retrieval with streaming nanopore sequencing. We further devise a digital codec to reduce requirements for synthesis accuracy and sequencing coverage, and experimentally show robust data retrieval from imperfectly synthesized strands. This work provides distributive enzymatic synthesis and information-theoretic approaches to advance digital information storage in DNA.

Adoption of DNA as a data storage medium could be accelerated with specialized synthesis processes and codecs. The authors describe TdT-mediated DNA synthesis in which data is stored in transitions between non-identical nucleotides and the use of synchronization markers to provide error tolerance.

Systematic evaluation and optimization of the experimental steps in RNA G-quadruplex structure sequencing

cDNA library preparation is important for many high-throughput sequencing applications, such as RNA G-quadruplex structure sequencing (rG4-seq). A systematic evaluation of the procedures of the experimental pipeline, however, is lacking. Herein, we perform a comprehensive assessment of the 5 key experimental steps involved in the cDNA library preparation of rG4-seq, and identify better reaction conditions and/or enzymes to carry out each of these key steps. Notably, we apply the improved methods to fragmented cellular RNA, and show reduced RNA input requirement, lower transcript abundance variations between biological replicates, as well as lower transcript coverage bias when compared to prior arts. In addition, the time to perform these steps is substantially reduced to hours. Our method and results can be directly applied in protocols that require cDNA library preparation, and provide insights to the further development of simple and efficient cDNA library preparation for different biological applications.

Sensitive detection of structural features and rearrangements in long, structured RNA molecules

Technical innovations in structural probing have drastically advanced the field of RNA structure analysis. These advances have led to parallel approaches developed in separate labs for analyzing RNA structure and dynamics. With the wealth of methodologies available, it can be difficult to determine which is best suited for a given application. Here, using a long, highly structured viral RNA as an example (the positive strand genome of Hepatitis C Virus), we present a semi-comprehensive analysis and describe the major approaches for analyzing the architecture of RNA that is modified with structure-sensitive probes. Additionally, we present an updated method for generating in vitro transcribed and folded RNA that maintains native secondary structures in long RNA molecules. We anticipate that the methods described here will streamline the use of current approaches and help investigators who are unfamiliar with structure probing, obviating the need for time-consuming and expensive optimization.

Insights from resolving protein-DNA interactions at near base-pair resolution

One of the central goals in molecular biology is to understand how cell-type-specific expression patterns arise through selective recruitment of RNA polymerase II (Pol II) to a subset of gene promoters. Pol II needs to be recruited to a precise genomic position at the proper time to produce messenger RNA from a DNA template. Ostensibly, transcription is a relatively simple cellular process; yet, experimentally measuring and then understanding the combinatorial possibilities of transcriptional regulators remain a daunting task. Since its introduction in 1985, chromatin immunoprecipitation (ChIP) has remained a key tool for investigating protein-DNA contacts in vivo. Over 30 years of intensive research using ChIP have provided numerous insights into mechanisms of gene regulation. As functional genomic technologies improve, they present new opportunities to address key biological questions. ChIP-exo is a refined version of ChIP-seq that significantly reduces background signal, while providing near base-pair mapping resolution for protein-DNA interactions. This review discusses the evolution of the ChIP assay over the years; the methodological differences between ChIP-seq, ChIP-exo and ChIP-nexus; and highlight new insights into epigenetic and transcriptional mechanisms that were uniquely enabled with the near base-pair resolution of ChIP-exo.

Small RNA-seq: The RNA 5'-end adapter ligation problem and how to circumvent it

The preparation of small RNA cDNA sequencing libraries depends on the unbiased ligation of adapters to the RNA ends. Small RNA with 5' recessed ends are poor substrates for enzymatic adapter ligation, but this 5' adapter ligation problem can go undetected if the library preparation steps are not monitored. Here we illustrate the severity of the 5' RNA end ligation problem using several pre-miRNA-like hairpins that allow us to expand the definition of the problem to include 5' ends close to a hairpin stem, whether recessed or in a short extension. The ribosome profiling method can avoid a difficult 5' adapter ligation, but the enzyme typically used to circularize the cDNA has been reported to be biased, calling into question the benefit of this workaround. Using the TS2126 RNA ligase 1 (a.k.a. CircLigase) as the circularizing enzyme, we devised a bias test for the circularization of first strand cDNA. All possible dinucleotides were circle-ligated with similar efficiency. To re-linearize the first strand cDNA in the ribosome profiling approach, we introduce an improved method wherein a single ribonucleotide is placed between the sequencing primer binding sites in the reverse transcriptase primer, which later serves as the point of re-linearization by RNase A. We incorporate this step into the ribosomal profiling method and describe a complete improved library preparation method, Coligo-seq, for the sequencing of small RNA with secondary structure close to the 5' end. This method accepts a variety of 5' modified RNA, including 5' monophosphorylated RNA, as demonstrated by the construction of a HeLa cell microRNA cDNA library.

Determining 3'-Termini and Sequences of Nascent Single-Stranded Viral DNA Molecules during HIV-1 Reverse Transcription in Infected Cells

Monitoring of nucleic acid intermediates during virus replication provides insights into the effects and mechanisms of action of antiviral compounds and host cell proteins on viral DNA synthesis. Here we address the lack of a cell-based, high-coverage, and high-resolution assay that is capable of defining retroviral reverse transcription intermediates within the physiological context of virus infection. The described method captures the 3'-termini of nascent complementary DNA (cDNA) molecules within HIV-1 infected cells at single nucleotide resolution. The protocol involves harvesting of whole cell DNA, targeted enrichment of viral DNA via hybrid capture, adaptor ligation, size fractionation by gel purification, PCR amplification, deep sequencing, and data analysis. A key step is the efficient and unbiased ligation of adaptor molecules to open 3'-DNA termini. Application of the described method determines the abundance of reverse transcripts of each particular length in a given sample. It also provides information about the (internal) sequence variation in reverse transcripts and thereby any potential mutations. In general, the assay is suitable for any questions relating to DNA 3'-extension, provided that the template sequence is known.

A Method for Single-Stranded Ancient DNA Library Preparation

Genomic library preparation from highly degraded DNA is more efficient when library molecules are prepared separately from the complementary strands of DNA fragments. We describe a protocol in which libraries are constructed from single DNA strands in a three-step procedure: single-stranded ligation of the first adapter with T4 DNA ligase in the presence of a splinter oligonucleotide, copying of the DNA strand with a proofreading polymerase, and blunt-end ligation of the second double-stranded adapter with T4 DNA ligase.

Simplified ChIP-exo assays

ChIP-seq and ChIP-exo identify where proteins bind along any genome in vivo. Although ChIP-seq is widely adopted in academic research, it has inherently high noise. In contrast, ChIP-exo has relatively low noise and achieves near-base pair resolution. Consequently, and unlike other genomic assays, ChIP-exo provides structural information on genome-wide binding proteins. Construction of ChIP-exo libraries is technically difficult. Here we describe greatly simplified ChIP-exo methods, each with use-specific advantages. This is achieved through assay optimization and use of Tn5 tagmentation and/or single-stranded DNA ligation. Greater library yields, lower processing time, and lower costs are achieved. In comparing assays, we reveal substantial limitations in other ChIP-based assays. Importantly, the new ChIP-exo assays allow high-resolution detection of some protein-DNA interactions in organs and in as few as 27,000 cells. It is suitable for high-throughput parallelization. The simplicity of ChIP-exo now makes it a highly appropriate substitute for ChIP-seq, and for broader adoption.

While ChIP-exo is low noise and highly informative regarding genome-wide binding proteins, libraries are difficult to construct. Here the authors present a simplified ChIP-exo method for high-resolution detection of interactions.

Detecting RNA G-Quadruplexes (rG4s) in the Transcriptome

RNA G-quadruplex (rG4) secondary structures are proposed to play key roles in fundamental biological processes that include the modulation of transcriptional, co-transcriptional, and posttranscriptional events. Recent methodological developments that include predictive algorithms and structure-based sequencing have enabled the detection and mapping of rG4 structures on a transcriptome-wide scale at high sensitivity and resolution. The data generated by these studies provide valuable insights into the potentially diverse roles of rG4s in biology and open up a number of mechanistic hypotheses. Herein we highlight these methodologies and discuss the associated findings in relation to rG4-related biological mechanisms.

Structure-seq2: sensitive and accurate genome-wide profiling of RNA structure in vivo

RNA serves many functions in biology such as splicing, temperature sensing, and innate immunity. These functions are often determined by the structure of RNA. There is thus a pressing need to understand RNA structure and how it changes during diverse biological processes both in vivo and genome-wide. Here, we present Structure-seq2, which provides nucleotide-resolution RNA structural information in vivo and genome-wide. This optimized version of our original Structure-seq method increases sensitivity by at least 4-fold and improves data quality by minimizing formation of a deleterious by-product, reducing ligation bias, and improving read coverage. We also present a variation of Structure-seq2 in which a biotinylated nucleotide is incorporated during reverse transcription, which greatly facilitates the protocol by eliminating two PAGE purification steps. We benchmark Structure-seq2 on both mRNA and rRNA structure in rice (Oryza sativa). We demonstrate that Structure-seq2 can lead to new biological insights. Our Structure-seq2 datasets uncover hidden breaks in chloroplast rRNA and identify a previously unreported N¹-methyladenosine (m¹A) in a nuclear-encoded Oryza sativa rRNA. Overall, Structure-seq2 is a rapid, sensitive, and unbiased method to probe RNA in vivo and genome-wide that facilitates new insights into RNA biology.

Single-stranded DNA library preparation from highly degraded DNA using T4 DNA ligase

DNA library preparation for high-throughput sequencing of genomic DNA usually involves ligation of adapters to double-stranded DNA fragments. However, for highly degraded DNA, especially ancient DNA, library preparation has been found to be more efficient if each of the two DNA strands are converted into library molecules separately. We present a new method for single-stranded library preparation, ssDNA2.0, which is based on single-stranded DNA ligation with T4 DNA ligase utilizing a splinter oligonucleotide with a stretch of random bases hybridized to a 3΄ biotinylated donor oligonucleotide. A thorough evaluation of this ligation scheme shows that single-stranded DNA can be ligated to adapter oligonucleotides in higher concentration than with CircLigase (an RNA ligase that was previously chosen for end-to-end ligation in single-stranded library preparation) and that biases in ligation can be minimized when choosing splinters with 7 or 8 random nucleotides. We show that ssDNA2.0 tolerates higher quantities of input DNA than CircLigase-based library preparation, is less costly and better compatible with automation. We also provide an in-depth comparison of library preparation methods on degraded DNA from various sources. Most strikingly, we find that single-stranded library preparation increases library yields from tissues stored in formalin for many years by several orders of magnitude.

Interpreting Reverse Transcriptase Termination and Mutation Events for Greater Insight into the Chemical Probing of RNA

Chemical probing has the power to provide insight into RNA conformation in vivo and in vitro, but interpreting the results depends on methods to detect the chemically modified nucleotides. Traditionally, the presence of modified bases was inferred from their ability to halt reverse transcriptase during primer extension and the locations of termination sites observed by electrophoresis or sequencing. More recently, modification-induced mutations have been used as a readout for chemical probing data. Given the variable propensity for mismatch incorporation and read-through with different reverse transcriptases, we examined how termination and mutation events compare to each other in the same chemical probing experiments. We found that mutations and terminations induced by dimethyl sulfate probing are both specific for methylated bases, but these two measures have surprisingly little correlation and represent largely nonoverlapping indicators of chemical modification data. We also show that specific biases for modified bases depend partly on local sequence context and that different reverse transcriptases show different biases toward reading a modification as a stop or a mutation. These results support approaches that incorporate analysis of both termination and mutation events into RNA probing experiments.

Facile single-stranded DNA sequencing of human plasma DNA via thermostable group II intron reverse transcriptase template switching

High-throughput single-stranded DNA sequencing (ssDNA-seq) of cell-free DNA from plasma and other bodily fluids is a powerful method for non-invasive prenatal testing, and diagnosis of cancers and other diseases. Here, we developed a facile ssDNA-seq method, which exploits a novel template-switching activity of thermostable group II intron reverse transcriptases (TGIRTs) for DNA-seq library construction. This activity enables TGIRT enzymes to initiate DNA synthesis directly at the 3′ end of a DNA strand while simultaneously attaching a DNA-seq adapter without end repair, tailing, or ligation. Initial experiments using this method to sequence E. coli genomic DNA showed that the TGIRT enzyme has surprisingly robust DNA polymerase activity. Further experiments showed that TGIRT-seq of plasma DNA from a healthy individual enables analysis of nucleosome positioning, transcription factor-binding sites, DNA methylation sites, and tissues-of-origin comparably to established methods, but with a simpler workflow that captures precise DNA ends.

RNA-seq of human reference RNA samples using a thermostable group II intron reverse transcriptase

Next-generation RNA sequencing (RNA-seq) has revolutionized our ability to analyze transcriptomes. Current RNA-seq methods are highly reproducible, but each has biases resulting from different modes of RNA sample preparation, reverse transcription, and adapter addition, leading to variability between methods. Moreover, the transcriptome cannot be profiled comprehensively because highly structured RNAs, such as tRNAs and snoRNAs, are refractory to conventional RNA-seq methods. Recently, we developed a new method for strand-specific RNA-seq using thermostable group II intron reverse transcriptases (TGIRTs). TGIRT enzymes have higher processivity and fidelity than conventional retroviral reverse transcriptases plus a novel template-switching activity that enables RNA-seq adapter addition during cDNA synthesis without using RNA ligase. Here, we obtained TGIRT-seq data sets for well-characterized human RNA reference samples and compared them to previous data sets obtained for these RNAs by the Illumina TruSeq v2 and v3 methods. We find that TGIRT-seq recapitulates the relative abundance of human transcripts and RNA spike-ins in ribo-depleted, fragmented RNA samples comparably to non-strand-specific TruSeq v2 and better than strand-specific TruSeq v3. Moreover, TGIRT-seq is more strand specific than TruSeq v3 and eliminates sampling biases from random hexamer priming, which are inherent to TruSeq. The TGIRT-seq data sets also show more uniform 5' to 3' gene coverage and identify more splice junctions, particularly near the 5' ends of mRNAs, than do the TruSeq data sets. Finally, TGIRT-seq enables the simultaneous profiling of mRNAs and lncRNAs in the same RNA-seq experiment as structured small ncRNAs, including tRNAs, which are essentially absent with TruSeq.

Protein-DNA binding in high-resolution

Recent advances in experimental and computational methodologies are enabling ultra-high resolution genome-wide profiles of protein-DNA binding events. For example, the ChIP-exo protocol precisely characterizes protein-DNA cross-linking patterns by combining chromatin immunoprecipitation (ChIP) with 5' → 3' exonuclease digestion. Similarly, deeply sequenced chromatin accessibility assays (e.g. DNase-seq and ATAC-seq) enable the detection of protected footprints at protein-DNA binding sites. With these techniques and others, we have the potential to characterize the individual nucleotides that interact with transcription factors, nucleosomes, RNA polymerases and other regulatory proteins in a particular cellular context. In this review, we explain the experimental assays and computational analysis methods that enable high-resolution profiling of protein-DNA binding events. We discuss the challenges and opportunities associated with such approaches.

Targeted Detection of G-Quadruplexes in Cellular RNAs

The G-quadruplex (G4) is a non-canonical nucleic acid structure which regulates important cellular processes. RNA G4s have recently been shown to exist in human cells and be biologically significant. Described herein is a new approach to detect and map RNA G4s in cellular transcripts. This method exploits the specific control of RNA G4-cation and RNA G4-ligand interactions during reverse transcription, by using a selective reverse transcriptase to monitor RNA G4-mediated reverse transcriptase stalling (RTS) events. Importantly, a ligation-amplification strategy is coupled with RTS, and enables detection and mapping of G4s in important, low-abundance cellular RNAs. Strong evidence is provided for G4 formation in full-length cellular human telomerase RNA, offering important insights into its cellular function.

An optimized kit-free method for making strand-specific deep sequencing libraries from RNA fragments

Deep sequencing of strand-specific cDNA libraries is now a ubiquitous tool for identifying and quantifying RNAs in diverse sample types. The accuracy of conclusions drawn from these analyses depends on precise and quantitative conversion of the RNA sample into a DNA library suitable for sequencing. Here, we describe an optimized method of preparing strand-specific RNA deep sequencing libraries from small RNAs and variably sized RNA fragments obtained from ribonucleoprotein particle footprinting experiments or fragmentation of long RNAs. Our approach works across a wide range of input amounts (400 pg to 200 ng), is easy to follow and produces a library in 2–3 days at relatively low reagent cost, all while giving the user complete control over every step. Because all enzymatic reactions were optimized and driven to apparent completion, sequence diversity and species abundance in the input sample are well preserved.

SHAPE-Seq 2.0: systematic optimization and extension of high-throughput chemical probing of RNA secondary structure with next generation sequencing

RNA structure is a primary determinant of its function, and methods that merge chemical probing with next generation sequencing have created breakthroughs in the throughput and scale of RNA structure characterization. However, little work has been done to examine the effects of library preparation and sequencing on the measured chemical probe reactivities that encode RNA structural information. Here, we present the first analysis and optimization of these effects for selective 2′-hydroxyl acylation analyzed by primer extension sequencing (SHAPE-Seq). We first optimize SHAPE-Seq, and show that it provides highly reproducible reactivity data over a wide range of RNA structural contexts with no apparent biases. As part of this optimization, we present SHAPE-Seq v2.0, a ‘universal’ method that can obtain reactivity information for every nucleotide of an RNA without having to use or introduce a specific reverse transcriptase priming site within the RNA. We show that SHAPE-Seq v2.0 is highly reproducible, with reactivity data that can be used as constraints in RNA folding algorithms to predict structures on par with those generated using data from other SHAPE methods. We anticipate SHAPE-Seq v2.0 to be broadly applicable to understanding the RNA sequence–structure relationship at the heart of some of life's most fundamental processes.