Kate (Kathryn) Lawrence
@itskatelawrence.bsky.social
59 followers 85 following 29 posts
Stanford Genetics PhD Candidate MIT Physics 2020 I also ride bikes
Posts Media Videos Starter Packs
itskatelawrence.bsky.social
Yes you're right - I think I need to look into our data a bit more and see how this is playing out. Even for the larger clusters, the PC is just a linear combination of the expression values so it would just be a QTL on a weighted sum/difference over more variables
itskatelawrence.bsky.social
so if something is a PC2 QTL, then it would be pulling expression away from that main axis of correlated variation, and then to a larger difference in expression levels for any given individual.
itskatelawrence.bsky.social
Thank you, I think you are right - I've been conceptualizing (at least for a correlated 2 gene cluster) PC1 capturing the shared main axis of variation, and PC2 capturing the spread away from that main axis (so for positively correlated variables I think you're right that would be the difference)
itskatelawrence.bsky.social
Just brainstorming, but do you think multi-ancestry GWAS/QTLs or comparing between ancestries could help to at least give us different LD confounders and tease apart those relationships? I'm not as familiar with the literature in that space though
itskatelawrence.bsky.social
I think establishing a gold-standard set of single-gene effect variants would be difficult - there are just so many strange mechanisms in biology! For instance, NMD of mutated mRNAs can cause up-regulation of paralogs. So even a PTVs can have multi-gene effects www.nature.com/articles/s41...
Genetic compensation triggered by mutant mRNA degradation - Nature
Transcriptional adaptation, a genetic compensation process by which organisms respond to mutations by upregulating related genes, is triggered by mRNA decay and involves a sequence-dependent mechanism...
www.nature.com
itskatelawrence.bsky.social
How to get at this in an unbiased way is further complicated by the fact that its fairly common for paralogs or functionally related genes to be near each other, like those IFNa/b and OAS1/OAS2/OASL loci
itskatelawrence.bsky.social
I don't, I think it would be really interesting to think of a way to systematically access how often it is one or the other. I like @kauralasoo.bsky.social 's point about comparing putative proxitropic regulatory QTLs with missense or splicing QTLs which are more likely to be single-gene effects.
itskatelawrence.bsky.social
It was a delight to work with you on this Tami! 🫶 hope for more collaboration in the future :)
itskatelawrence.bsky.social
That's an interesting point and definitely worth thinking about. I think a case with multiple causal variants in LD with different marginal effects on each of the genes near them would require orthogonal approaches (perhaps MPRAs?) to disentangle.
Multiple causal variants underlie genetic associations in humans
Some human quantitative traits are determined by multiple linked genetic variants.
www.science.org
itskatelawrence.bsky.social
I hadn't seen your paper, thank you for sharing! I was wondering, in your analysis, did you see if cis-eQTLs did worse in the prioritization task when the pQTL gene they attempted to pick shared regulation with nearby genes or overlapped nearby genes?
itskatelawrence.bsky.social
Ribeiro et. al. has a nice analysis showing eSNPs shared between multiple eGenes are more indeed pleiotropic. We also see this in our data www.nature.com/articles/s41...
itskatelawrence.bsky.social
I think that both of these occur. Sometimes maybe a pcQTL effects multiple genes and we need the power boost from looking at all the genes to detect it, but only one gene actually impacts the trait, but I believe sometimes it is a multi-gene or even second order gene co-expression effect
itskatelawrence.bsky.social
Yes I think you are exactly right! They tend to be weaker effects. We also have some preliminary results not in this preprint indicating that as you increase sample size, you end up discovering many of the loci as eQTLs that were only pcQTLs at the smaller sample size
itskatelawrence.bsky.social
Very exciting you'll be coming to Stanford, I would love to chat when you're here! Do you think extending COWAS to larger groups of proteins with PCs would be feasible?
itskatelawrence.bsky.social
Thank you! No worries, I only just made the account so it would have been impossible to find!
itskatelawrence.bsky.social
If you're interested in a bit more detail about what we did, check out our tweetorial here! bsky.app/profile/itsk...
itskatelawrence.bsky.social
If you're interested in a bit more detail about what we did, check out our tweetorial here! bsky.app/profile/itsk...
itskatelawrence.bsky.social
Excited to share my first PhD paper in the @sbmontgom.bsky.social lab with @tamigj.bsky.social (www.biorxiv.org/content/10.1...)! Standard QTL methods treat each gene independently. But what if a single variant regulates multiple nearby genes at once - what we call “allelic proxitropy”? 🧵 ⬇️
Standard methods are equivalent to a flashlight, looking at each gene independently. We combine signals from multiple genes, turning a floodlight onto the genome.
itskatelawrence.bsky.social
Thank you for the shoutout! I just posted a tweetorial (bsky.app/profile/itsk...) if anyone is interested in a bit more about what we did!
itskatelawrence.bsky.social
Excited to share my first PhD paper in the @sbmontgom.bsky.social lab with @tamigj.bsky.social (www.biorxiv.org/content/10.1...)! Standard QTL methods treat each gene independently. But what if a single variant regulates multiple nearby genes at once - what we call “allelic proxitropy”? 🧵 ⬇️
Standard methods are equivalent to a flashlight, looking at each gene independently. We combine signals from multiple genes, turning a floodlight onto the genome.
itskatelawrence.bsky.social
🧵5/5 Make allelic proxitropy work for you! Even this simple approach substantially increased colocalizations. We’re excited to see future refinements as we move beyond the “one variant, one gene” paradigm to get a more complete picture of the complexities of local gene regulation 🔬🧬✨
itskatelawrence.bsky.social
🧵4/5 But are new pcQTLs relevant? We examined associations with GWAS for complex traits and diseases and found over 500 new colocalizations – a 34% increase over single-gene-based analyses! In fact, variants with distributed effects are more likely to be GWAS hits than those affecting just one gene