Case studies

Projects

Seven bioinformatics projects, each written up as a case study: the question, the approach, the results and their limits. Code for all of them is on GitHub.

Forest plot of GPNMB Mendelian randomisation odds ratios for Parkinson’s disease across SomaScan and Olink instrument sets

MSc dissertation · Statistical genetics

GPNMB and Parkinson’s disease: Mendelian randomisation and colocalisation

Two-sample MR, LD-aware sensitivity analysis and Bayesian colocalisation to test whether genetically predicted GPNMB protein levels are linked to Parkinson’s disease risk. The evidence supports prioritisation, not causality.

  • Mendelian randomisation
  • Colocalisation
  • SuSiE
  • pQTL / eQTL
  • R
  • OpenGWAS
Read case study →
Regional plot of the GPNMB locus from the pipeline: UKB-PPP Olink pQTL association above, Parkinson’s disease GWAS below

Human genetics · Target discovery · Nextflow

Proteome-wide cis-MR and colocalisation scan for Parkinson’s disease

Scaled the dissertation from one protein to a Nextflow pipeline built to test each of the ~2,900 UKB-PPP plasma proteins with cis-pQTL MR and colocalisation, then apply FDR control, evidence tiers and an HTML report.

  • Nextflow
  • Mendelian randomisation
  • Colocalisation
  • UKB-PPP
  • Python
  • pytest
Read case study →
UMAP of 83,484 substantia nigra nuclei coloured by author cell-type annotation

Single-nucleus transcriptomics · Parkinson’s disease

GPNMB in Parkinson’s disease: snRNA-seq re-analysis

Re-analysed 83,484 substantia nigra nuclei from 29 donors to localise GPNMB and test PD–control differences with donor-level pseudobulk models. GPNMB was highest in microglia and about two-fold higher in PD microglia (P = 0.009), a nominal effect that did not survive transcriptome-wide FDR.

  • snRNA-seq
  • Scanpy
  • Pseudobulk
  • PyDESeq2
  • Microglia
  • Propeller
Read case study →
Forest plot of GPNMB hazard ratios: 2.00 unadjusted falling to 1.15 after adjustment for grade and IDH status

Biostatistics · Survival analysis

GPNMB and glioma survival: separating a prognostic signal from confounding

Pre-specified Cox analysis of 624 TCGA diffuse gliomas. GPNMB’s hazard ratio fell from 2.00 per SD unadjusted to 1.15 (95% CI 0.96–1.38) after adjusting for grade and IDH/1p19q status, with an optimism-corrected C-index gain of 0.0009.

  • Cox regression
  • Kaplan–Meier
  • Bootstrap
  • lifelines
  • TCGA
  • pytest
Read case study →
Mean MDS-UPDRS and MoCA trajectories over 15 years in the PPMI Parkinson's disease cohort, observed and model-based

Biostatistics · Longitudinal and survival analysis

Parkinson’s disease progression in PPMI: trajectories and clinical milestones

Mixed-effects and Cox models of 1,789 PPMI participants. Medication-adjusted motor progression slowed from 2.0 to 1.2 MDS-UPDRS III points a year, MoCA held for five years before declining, and age was the strongest predictor of milestones.

  • Mixed-effects models
  • Cox regression
  • Kaplan–Meier
  • statsmodels
  • lifelines
  • PPMI
Read case study →
Galaxy workflow graph for M. tuberculosis WGS variant analysis

Microbial genomics · WGS

M. tuberculosis WGS variant analysis (Galaxy)

End-to-end Galaxy workflow from SRA reads to annotated resistance variants, run on 6 clinical isolates. It called rpoB and katG mutations, confirmed in IGV, and classified 4 isolates as genotypically MDR-TB.

  • WGS
  • Galaxy
  • BWA-MEM2
  • bcftools
  • SnpEff
  • IGV
Read case study →
Volcano plot of dexamethasone versus control in airway smooth muscle cells with known response genes labelled

Transcriptomics · Pipelines

RNA-seq pipeline (Nextflow + Docker)

Samplesheet-driven Nextflow DSL2 pipeline from raw FASTQ to differential expression in one command, with every tool in a versioned container. Validated on real human data: 8 of 8 known dexamethasone-response genes recovered.

  • Nextflow
  • Docker
  • RNA-seq
  • STAR
  • Salmon
  • DESeq2
Read case study →
Saturation mutagenesis of the TERT promoter with the C228T and C250T cancer mutations highlighted

Regulatory genomics · Machine learning

RegulonML: regulatory variant effects in unseen promoters and enhancers

Predicted the effects of 41,724 single-base variants in 21 regulatory elements from DNA sequence alone, scoring binding-site gain and loss for 1,019 JASPAR motifs and testing only on loci the model had never seen.

  • MPRA
  • Regulatory genomics
  • JASPAR
  • scikit-learn
  • Gradient boosting
  • Grouped CV
Read case study →

Other work

Machine learning, LLMs and data engineering

Mostly course capstones, kept here because they show range. Each links straight to its repository.

Get in touch

I’m looking for bioinformatics roles in statistical genetics, transcriptomics and NGS analysis, especially where clinical or cell and gene therapy experience helps. Based in London, open to hybrid and remote.