# PeakPrime: a peak-guided primer design pipeline for target enrichment in 3′-end RNA-seq

**Authors:** Franco Poma-Soto, Brecht Soulliaert, Hanne Van Droogenbroeck, Pieter Mestdagh, Jo Vandesompele

PMC · DOI: 10.1093/bioadv/vbag080 · Bioinformatics Advances · 2026-03-19

## TL;DR

PeakPrime is a pipeline that improves 3′-end RNA-seq by designing primers to enrich for specific transcripts, increasing read allocation and detection of low-abundance genes.

## Contribution

PeakPrime introduces a reproducible, peak-guided primer design pipeline for targeted 3′-end RNA-seq using MACS2, Primer3, and Bowtie2.

## Key findings

- Targeted libraries increased reads assigned to target genes from <1% to ~25%–42%.
- Gene-level expression correlations with random-primed RNA-seq were high (Pearson r≈0.88).
- Low-abundance transcript detection was significantly improved with the new method.

## Abstract

Targeted enrichment can offset the bias and depth requirements of random-primed second-strand synthesis in 3'-end RNA-seq by reallocating reads to transcripts of interest. We present PeakPrime, a reproducible Nextflow pipeline that (i) identifies high-coverage 3 ′ RNA-seq regions via MACS2, (ii) selects exonic windows around significant peaks, (iii) designs strand-appropriate cDNA primers via Primer3, and (iv) screens specificity with Bowtie2, producing publication-ready plots and quality control (QC) summaries.

We demonstrate coverage-informed primer design using a custom 3 ′ protocol that employs a barcoded oligo(dT)-VN primer, optimize parameters on IMR-32 neuroblastoma and Universal Human Reference RNA data, and validate primers by quantitative PCR (qPCR) and targeted 3 ′ RNA-seq. Targeted libraries preserve gene-level expression correlations with random-primed 3 ′ RNA-seq (Pearson r≈0.88 across conditions), while boosting the fraction of reads assigned to target genes from <1% to ∼25%–42% and increasing detection of the lowest-abundance transcripts. While empirically tested in bulk RNA-seq, PeakPrime enables extensions to single-cell and spatial workflows where oligo(dT) capture is followed by in situ reverse transcription and second-strand synthesis.

PeakPrime source code, example configs, and documentation are available at https://github.com/OncoRNALab/PeakPrime.

## Linked entities

- **Diseases:** neuroblastoma (MONDO:0005072)
- **Species:** Homo sapiens (taxon 9606)

## Full-text entities

- **Species:** Homo sapiens (human, species) [taxon 9606]

## Full text

_Full body text omitted from this summary view._ Fetch the complete paper as Markdown: https://tomesphere.com/paper/PMC13034549/full.md

## Figures

5 figures with captions in the complete paper: https://tomesphere.com/paper/PMC13034549/full.md

## References

27 references — full list in the complete paper: https://tomesphere.com/paper/PMC13034549/full.md

---
Source: https://tomesphere.com/paper/PMC13034549