diff --git a/README.md b/README.md
index 94825d4..5cdb0a6 100644
--- a/README.md
+++ b/README.md
@@ -15,3 +15,7 @@ Example scripts for different tasks:
 - Preparing data: `data/mlm/Snakefile`
 - Training: `analysis/mlm/train_512_convnet_only_athaliana.sh`
 - Variant effect prediction: `analysis/mlm/run_vep.py`
+
+## Citation
+Gonzalo Benegas, Sanjit Singh Batra and Yun S. Song "DNA language models are powerful zero-shot predictors of non-coding variant effects" bioRxiv (2022)  
+DOI: [10.1101/2022.08.22.504706](https://doi.org/10.1101/2022.08.22.504706)