Many disease-associated variants are thought to be regulatory but are not present in existing catalogues of expression quantitative trait loci (eQTL). We hypothesise that these variants may regulate expression in specific biological contexts, such as stimulated immune cells. Here, we used human iPSC-derived macrophages to map eQTLs across 24 cellular conditions. We found that 76% of eQTLs detected in at least one stimulated condition were also found in naive cells. The percentage of response eQTLs (reQTLs) varied widely across conditions (3.7% - 28.4%), with reQTLs specific to a single condition being rare (1.11%). Despite their relative rarity, reQTLs were overrepresented among disease-colocalizing eQTLs. We nominated an additional 21.7% of disease effector genes at GWAS loci via colocalization of reQTLs, with 38.6% of these not found in the Genotype-Tissue Expression (GTEx) catalogue. Our study highlights the diversity of genetic effects on expression and demonstrates how condition-specific regulatory variation can enhance our understanding of common disease risk alleles.

The majority of immune-mediated disease (IMD) risk loci are located in non-coding regions of the genome, making it difficult to decipher their functional effects. To assess the extent to which alternative splicing contributes to IMD risk, we mapped genetic variants associated with alternative splicing (splicing quantitative trait loci or sQTL) in macrophages exposed to 24 cellular conditions. We found that genes involved in innate immune response pathways undergo extensive differential splicing in response to stimulation and detected significant sQTL effects for 5,734 genes across all conditions. We colocalised sQTL signals for over 700 genes with IMD-associated risk loci from 21 IMDs with high confidence (PP4 ≥ 0.75). Approximately half of the colocalisations implicate lowly-used splice junctions (mean usage ratio < 0.1). Finally, we demonstrate how an inflammatory bowel disease (IBD) risk allele increases the usage of a lowly-used isoform of PTPN2, a negative regulator of inflammation. Together, our findings highlight the role alternative splicing plays in IMD risk, and suggest that lowly-used splicing events significantly contribute to complex disease risk.

Common genetic variants across individuals modulate the cellular response to pathogens and are implicated in diverse immune pathologies, yet how they dynamically alter the response upon infection is not well understood. Here, we triggered antiviral responses in human fibroblasts from 68 healthy donors, and profiled tens of thousands of cells using single-cell RNA-sequencing. We developed GASPACHO (GAuSsian Processes for Association mapping leveraging Cell HeterOgeneity), a statistical approach designed to identify nonlinear dynamic genetic effects across transcriptional trajectories of cells. This approach identified 1,275 expression quantitative trait loci (local false discovery rate 10%) that manifested during the responses, many of which were colocalized with susceptibility loci identified by genome-wide association studies of infectious and autoimmune diseases, including the OAS1 splicing quantitative trait locus in a COVID-19 susceptibility locus. In summary, our analytical approach provides a unique framework for delineation of the genetic variants that shape a wide spectrum of transcriptional responses at single-cell resolution.

Mutational signatures are imprints of pathophysiological processes arising through tumorigenesis. We generated isogenic CRISPR–Cas9 knockouts (∆) of 43 genes in human induced pluripotent stem cells, cultured them in the absence of added DNA damage and performed whole-genome sequencing of 173 subclones. ∆OGG1, ∆UNG, ∆EXO1, ∆RNF168, ∆MLH1, ∆MSH2, ∆MSH6, ∆PMS1 and ∆PMS2 produced marked mutational signatures indicative of them being critical mitigators of endogenous DNA modifications. Detailed analyses revealed mutational mechanistic insights, including how 8-oxo-2′-deoxyguanosine elimination is sequence context specific while uracil clearance is sequence context independent. Mismatch repair (MMR) deficiency signatures are engendered by oxidative damage (C > A transversions) and differential misincorporation by replicative polymerases (T > C and C > T transitions), and we propose a reverse template slippage model for T > A transversions. ∆MLH1, ∆MSH6 and ∆MSH2 signatures were similar to each other but distinct from ∆PMS2. Finally, we developed a classifier, MMRDetect, where application to 7,695 whole-genome-sequenced cancers showed enhanced detection of MMR-deficient tumors, with implications for responsiveness to immunotherapies. Nik-Zainal and colleagues leverage CRISPR–Cas9 and whole-genome sequencing to examine mutational patterns following knockout of 42 human DNA repair genes. They further develop and validate a clinically relevant tool to detect mismatch repair-deficient tumors.