COVID-19 mRNA vaccines as hypothetical epigenetic players: Results from an in silico analysis, considerations and perspectives

Objectives: To investigate in silico the occurrence of epigenetic crosstalk by nucleotide sequence complementarity between the BNT162b2 mRNA vaccine and whole human genome, including coding and noncoding (nc)RNA genes. To correlate these results with those obtained with the original spike (S) gene of Severe Acute Respiratory Syndrome CoronaVirus-2 (SARS-CoV-2). Methods: The publicly available FASTA sequence of the BNT162b2 mRNA vaccine and the SARS-CoV-2 isolate Wuhan-Hu-1 S gene (NC_045512.2) were used separately as key input to the Ensembl.org library to evaluate base pair match to human GRCh38 genome. Human coding and noncoding genes harboring hits were assessed for functional activity and health effects using bioinformatics tools and GWAS databases. Results: The BLAT analysis against the human GRCh38 genome revealed a total of 37 hits for BNT162b2 mRNA and no hits for the SARS-CoV-2 S gene. More specifically, BNT162b2 mRNA matched 19 human genes whose protein products are variously involved in enzyme reactions, nucleotide or cation binding, signaling, and carrier functions. In BLASTN analysis of ncRNA genes, BNT162b2 mRNA and SARS-CoV-2 S gene matched 17 and 24 different human genomic regions, respectively. Overall, characterization of the matched noncoding sequences revealed stronger interference with epigenetic pathways for BNT162b2 mRNA compared with the original S gene. Conclusion: This pivotal in silico analysis shows that SARS-CoV-2 S gene and the BNT162b2 mRNA vaccine exhibit Watson-Crick nucleotide complementarity with human coding or noncoding genes. Although they do not share the same complementarity pattern, both may disrupt epigenetic mechanisms in target cells, potentially leading to long-term complications.