Full-length nanopore sequencing of circular RNA landscape in peripheral blood cells following sequential BNT162b2 mRNA vaccination

Circular RNAs (circRNA) lack 5' or 3' ends; their unique covalently closed structures prevent RNA degradation by exonucleases. These characteristics provide circRNAs with high pharmaceutical stability and biostability relative to current standard-of-care linear mRNAs. CircRNA levels are reportedly associated with certain human diseases, making them novel disease biomarkers and a noncanonical class of therapeutic targets. In this study, the endogenous circRNAs underlying the response to BNT162b2 mRNA vaccination were evaluated. To this end, peripheral blood samples were subjected to full-length sequencing of circRNAs via nanopore sequencing and transcriptome sequencing. Fifteen samples, comprising pre-, first, and second vaccination cohorts, were obtained from five healthcare workers with no history of SARS-CoV-2 infection or previous vaccination. A total of 4706 circRNAs were detected; following full-length sequencing, 4217 novel circRNAs were identified as being specifically expressed during vaccination. These circRNAs were enriched in the binding motifs of stress granule assemblies and SARS-CoV-2 RNA binding proteins, namely poly(A) binding protein cytoplasmic 1 (PABPC1), pumilio RNA binding family member 1 (PUM1), and Y box binding protein 1 (YBX1). Moreover, 489 circRNAs were identified as previously reported miRNA sponges. The differentially expressed circRNAs putatively originated from plasma B cells compared to circRNAs reported in human blood single-cell RNA sequencing datasets. The pre- and post-vaccination differences observed in the circRNA expression landscape in response to the SARS-CoV-2 BNT162b2 mRNA vaccine.