Intrinsic DNA topology as a prioritization metric in genomic fine-mapping studies

In genomic fine-mapping studies, some approaches leverage annotation data to prioritize likely functional polymorphisms. However, existing annotation sources often present challenges as many: lack data for novel variants, offer no context for noncoding regions, and/or are confounded with linkage disequilibrium. We propose a novel annotation source - sequence-dependent DNA topology - as a prioritization metric for fine-mapping. DNA topology and function are well-intertwined, and as an intrinsic DNA property, it is readily applicable to any genomic region. Here, we constructed and applied, Minor Groove Width (MGW), as a prioritization metric. Using an established MGW-prediction method, we generated an MGW census for 199,038,197 SNPs across the human genome. Summarizing a SNPs change in MGW ({Delta}MGW) as a Euclidean distance, {Delta}MGW exhibited a strongly right-skewed distribution, highlighting the infrequency of SNPs that generate dissimilar shape profiles. We hypothesized that phenotypically-associated SNPs can be prioritized by {Delta}MGW. We applied Bayesian and frequentist MGW-prioritization approaches to three non-coding regions associated with System Lupus Erythematosus in multiple ancestries. In two regions, including {Delta}MGW resolved the association to a single, trans-ancestral, SNP, corroborated by external functional data. Together, this study presents the first usage of sequence-dependent DNA topology as a prioritization metric in genomic association studies. Graphical AbstractWe hypothesize that SNPs imposing dissimilar minor groove width profiles ({Delta}MGW) are more likely to alter function. {Delta}MGW was interrogated genome-wide and then used as a weighting metric for fine-mapping associations. O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=153 SRC="FIGDIR/small/837245v2_ufig1.gif" ALT="Figure 1"> View larger version (42K): org.highwire.dtl.DTLVardef@d350adorg.highwire.dtl.DTLVardef@155a2cborg.highwire.dtl.DTLVardef@132e30aorg.highwire.dtl.DTLVardef@1d7b065_HPS_FORMAT_FIGEXP M_FIG C_FIG