Abdul Muntakim Rafi

Research

Towards mapping the cis-regulatory landscape.

We have gained substantial ability to predict how much a DNA sequence is expressed in different cell types. But we do not fully understand the cis-regulatory landscape, or how it varies as trans factors change. I believe it is possible to map the entire cis-regulatory landscape, in every cell type, by generating the right data and using predictive modelling. I work on developing technologies that generate that data for the purpose of training models, and on methods to train, evaluate, interpret, and reliably apply those models.

The questions I try to address in my work.

How do we generate the data?

Most data used to train models of gene regulation were generated to characterise biology, not to train models. I work on developing technologies to synthesize sequence libraries and to design experiments whose purpose is to train models.

How do we train the best models?

Model performance depends on the architecture, the training strategy, and how the training data are processed and different data points are weighted.

What have the models learned?

Standard evaluations can overstate what a model has learned, and models can learn spurious correlations and dataset structure that are predictive without being causal. I work on detecting these and accounting for them, and on evaluations robust enough to probe the different aspects of the underlying problem. Separating learned regulatory logic from memorization is a step towards causal modelling.

Can we trust a single prediction?

Model performance is usually reported over a whole held-out test set, capturing different aspects of the underlying problem. That does not tell us how much to trust an individual prediction, such as a single variant effect. I work on estimating the reliability of individual predictions.

Can we trust how we read them?

Most biological conclusions from these models are drawn through interpretation methods such as attribution and perturbation, which are themselves not extensively tested. I work on assessing how faithful different interpretation methods are.

Projects

Projects, published and in progress.

Chromosome-scale sequence from outside the host

In progress

Short oligos sit in one fixed context, so they cannot report on promoter–gene distance, chromatin or real transcripts. Sequence carried on yeast artificial chromosomes behaves like an extra chromosome and has never been under selection in the organism reading it. We annotated a single YAC in Luthra et al., then increased the data tenfold to train models on it, which became Yorzoi, built with Timon Schneider and Tom Ellis at Imperial College London. We are now scaling the data tenfold again for Yakformer.

High-information-content libraries

In progress

Designing sequence synthesis approaches that create high-information-content sequence libraries.

SynBIOM

In progress

A biochemical simulator of cis-regulation. Users specify the transcription factors at play, with their binding sites or PWMs, concentrations, activity levels and accessibility potentiation, and SynBIOM simulates expression from sequence. We use it to benchmark several aspects of sequence-to-expression modeling.

pairFrag — genome-wide homology mapping

In preparation

Making homology-aware evaluation something any group can do without repeating the computation.

hashFrag — homology, leakage and memorization

Preprint

Neither chromosomal nor random train/test splits account for homology within a species, so standard evaluations of genome-trained models are inflated. I measured how far, showed that the dependence on training-set similarity is not monotonic, and released hashFrag, which detects homology and partitions data at roughly a hundredth of the compute of exhaustive alignment. Its recommendation is to stratify a test set rather than build a fully orthogonal one, because an orthogonal split hides the bias instead of exposing it.

nextFrag — benchmarking active learning

Preprint

If every sequence has to be paid for, each one should be chosen to be informative. We benchmarked six selection strategies across architectures, datasets and configurations, simulated on pools that had already been measured, so the benchmark itself needed no new experiment. All beat random sampling, uncertainty-based methods did best while being cheapest to compute, and most of the gain from many small acquisition rounds survives with fewer, larger ones — which is what makes lab-in-the-loop practical. Selected sequences look distinctive, but selecting directly on those properties never matched active learning: informativeness is a property of the model’s ignorance, not of the sequence. Building on this, we are extending the work to large-scale experimental data, to report how active learning is best done in genomics.

gRely — reliability of individual predictions

Preprint

Aggregate benchmark scores cannot tell a user whether a single variant-effect prediction is correct. I built a meta-model that estimates the probability that it is, from features of the variant, gene, tissue and model. Its top-scoring fifth reaches 97% sign concordance against 54% in the bottom fifth, and it stays discriminative among the low-magnitude variants that effect-size filtering discards, which is where most GWAS signal is expected to act. It transfers zero-shot to other architectures, so reliability looks like a property of the locus rather than of the model. Begun during an internship at Genentech.

Random Promoter DREAM Challenge

Published

Random sequence removes the homology and sample-size ceilings at once. I designed and ran an open challenge on 6.7 million random promoters measured in yeast, with more than 110 teams and 28 final models; nineteen beat the previous state of the art. The winner had the fewest parameters and three of the top five used no transformer, so I built a framework that recombines the entrants’ modules across architectures and trainers to find out which choice carried the gain. It was the trainer. Models tuned on random yeast sequence then transferred to other species and assays.

Interested in collaborating on any of these? Reach out .

Funding

Grants and compute.

Four Year Doctoral Fellowship (4YF)

University of British Columbia

96,000 CAD over four years. UBC’s flagship doctoral fellowship, awarded with tuition to top PhD students.

2021 – 2025

Continual improvement of gene regulatory models

The Digital Research Alliance of Canada · PI: Carl de Boer

Priority access to GPUs. Co-wrote the proposal.

2025 – 2026

Evaluation, optimization and continual improvement of sequence-based cis-regulatory models

Advanced Research Computing, UBC · PI: Carl de Boer

20,000 CAD in Microsoft Azure credit. Co-applicant; wrote the proposal.

Aug 2023 – Jun 2024

Lossless preprocessing of the sequence and expression space to improve sequence-based gene regulatory models

School of Biomedical Engineering, UBC · PI: Carl de Boer

6,000 CAD. Co-applicant; wrote the proposal and hired a Co-op student through the grant.

May – Aug 2023

Random Promoter DREAM Challenge 2022

TPU Research Cloud, Google · PI: Carl de Boer (UBC), Pablo Meyer (IBM Research), Jake Albrecht (Sage Bionetworks)

50 TPU quotas. Sole graduate student on the organising committee.

May – Jul 2022

Identifying selection on human gene expression with gene regulatory models

The Digital Research Alliance of Canada · PI: Carl de Boer

Priority access to GPUs. Co-wrote the proposal.

2022 – 2024

Efficient edge inference benchmarking for AI-driven applications

Mitacs Accelerate · PI: Jonathan Wu

15,000 CAD. Sole co-applicant; wrote the proposal.

Nov 2020 – Mar 2021

Spatio-temporal human activity recognition on manufacturing floors

Mitacs Accelerate · PI: Jonathan Wu

22,500 CAD. Co-applicant; assisted with proposal writing.

Oct 2019 – Apr 2020