Identification of DNA Termini in Sequencing Data through Combined Analysis of End Capture and Local Strand Bias
Creators: name: Wang, William
name: Fire, Andrew
Publisher: name: Stanford Digital Repository
Year: 2023
Subjects: subject: Genetics
subject: FOS: Biological sciences
subject: Molecular genetics
subject: Genomics
subject: Cancer
subject: DNA
subject: DNA termini
subject: Extrachromosomal DNA
subject: Linear extrachromosomal DNA
subject: Next-generation sequencing (NGS)
subject: High-throughput sequencing (HTS)
subject: Illumina
subject: HIV
subject: Retrovirus
subject: Bacteriophage
subject: Phage
subject: Bacteriophage lambda
subject: Restriction enzyme
subject: HIV replication
subject: HIV reverse transcription
subject: Double-stranded DNA
subject: Single-stranded DNA
subject: RNA
subject: Purine-rich sequences
subject: Central polypurine tract (cPPT)
subject: Alternative polypurine tracts (altPPT)
subject: Adenovirus
subject: Mitochondrial DNA
subject: Mitochondria
subject: Caenorhabditis elegans
subject: Chlamydomonas reinhardtii
subject: Innate immunity
subject: Biology

description: Detecting DNA termini, such as ends of linear extrachromosomal DNA, plays an essential role in understanding the structure and functions of DNA molecules. Here we describe an approach combining direct and indirect computational methods to detect DNA termini from next-generation short-read sequencing. While a direct inference of ends can come from mapping the specific capture points of DNA fragments, this approach is insufficient for analytical pipelines where the DNA termini are not captured. Thus, we add an indirect detection of ends based on strand bias, the difference in sequence representation between the plus and minus strands of DNA in a dataset. Termini are reflected by a strong strand bias, with inward-facing reads greatly enriched over outward-facing reads in the immediate proximity of any end. Applying this analysis to negative control regions (where DNA is continuous and with no known termini), we observe no strong end capture peaks or strand bias. Applying to positive control regions where known DNA termini are present yields strong strand bias signals even in cases where blocked termini prevent end capture (for a protein-blocked adenovirus), or where ends are not explicitly captured (tagmentation of restriction digested lambda DNA). Analysis of a more complex situation (HIV replication) produces a picture that includes both the known termini of the reverse-transcribed genome (the PBS [primer binding site] on the negative strand and the PPT [3’ polypurine tract] on the positive strand) as well as a signal corresponding to a previously described additional initiation site for second strand synthesis (cPPT [central polypurine tract]). These results confirm the ability to detect DNA structural discontinuities in a pooled sample where high throughput shotgun sequence data is available. In addition to the known initiation sequence in the HIV genome, we detect a signal of positive strand DNA termini at several positions on the plus strand sequence. These sites share several characteristics with the previously characterized second strand initiation sites (the cPPT and 3’ PPT sites): (i) observed spike in directly captured cDNA ends, (ii) an indirect terminus signal evident in localized strand bias, (iii) a strong preference for forward-facing termini, (iv) an upstream purine-rich motif, and (v) a decrease in terminus signal at late time points after infection. These characteristics are consistent in duplicate samples in two different genotypes (wild type and integrase-lacking HIV).  The observation of distinct internal termini associated with multiple purine rich regions suggests the possibility that multiple internal initiations of second strand synthesis might contribute to HIV replication through acceleration of second strand synthesis and/or strand displacement at the HIV 3’ end.