ShorkieΒΆ
Predicting dynamic expression patterns in budding yeast with a fungal DNA language model.
Shorkie is a sequence-to-expression model for Saccharomyces cerevisiae: a masked DNA language model pretrained on hundreds of fungal genomes, then fine-tuned on 5,215 epigenomic and transcriptomic tracks β including a large set of transcriptional-regulator induction RNA-seq experiments generated for this study β to predict RNA-seq coverage and variant effects.
What you can do with ShorkieΒΆ
Run the 8-fold ensemble over any 16 kb window and get predicted coverage across all 5,215 tracks.
Compute logSED for a SNP β the metric used in the paperβs eQTL and MPRA benchmarks β from the command line or a notebook.
Masked-token prediction and self-attention embeddings from the fungal LM that Shorkie is built on.
Transfer the LM trunk to your own RNA-seq / ChIP data, including a small demo you can actually run.
A figure-by-figure tour of the paper, each linked to an executed notebook that regenerates it and re-checks the published numbers.
Every dataset curated for this study β models, four pretraining corpora, reference genome, and the eQTL/MPRA benchmarks β with sizes and commands.
Quick startΒΆ
git clone --recurse-submodules https://github.com/calico/shorkie-paper.git
cd shorkie-paper
conda env create -f environment.yml && conda activate yeast_ml
pip install -e external/baskerville-yeast -e external/westminster -e .
cp config/paths.example.yaml config/paths.yaml
bash data/download.sh --minimal # 8 Shorkie folds (~0.46 GB)
bash data/download.sh --genome -u <your-gcp-project> # R64 FASTA + GTF
python minimal_example/run_shorkie_variant.py \
--model_dir ./my_shorkie
See Installation for details and Quick start for a walkthrough of that commandβs output.
The three released modelsΒΆ
Model |
What it is |
Use it for |
|---|---|---|
Shorkie_LM |
Masked DNA language model, pretrained on 165 Saccharomycetales genomes |
Masked-base prediction, sequence embeddings, and as the starting point for fine-tuning |
Shorkie |
8-fold supervised ensemble, fine-tuned from Shorkie_LM on 5,215 tracks |
Coverage prediction and variant-effect scoring β this is the main model |
Shorkie_Random_Init |
Same architecture and data, trained from scratch (no LM pretraining, lr 5e-4) |
The ablation that isolates what LM pretraining contributes |
All three are public and catalogued with checksums β see Released models.
Getting started
Using Shorkie
Background
- Analysis gallery
- Figure 1 β The fungal corpus and the model
- Figure 2 β What the language model learned
- Figure 3 β Predicting RNA-seq coverage
- Figure 4 β Promoter and splicing motifs
- Figure 5 β Dynamic induction time courses
- Figure 6 β MPRA variant effects
- Figure 7 β cis-eQTL variant effects
- Running them yourself
- How Shorkie was trained
- Datasets
- Reproducing the paper figures
Reference