Comments (1)
Hey Matteo,
Thanks for using QTLseqr in your work! Happy to hear you were able to get some results.
Yes in theory you can use a reduced representation sequencing system with QTLseqr but you might need to be a little creative.
So first of all my question is do you have GBS data for the bulk as a bulk (ie one pooled DNA sample per bulk and then sequenced on GBS) or do you have GBS data for every single individual in the bulk (ie each of say n=15 individuals was GBSed separately for each of the bulks).
If you have 1 sample per bulk, you should be able to just call SNPs as if you were working with regular WGS and though i've personally never done this (and not sure about the biases and read depth issues with GBS and this approach) it should be able to run like the regular pipeline.
If you have the second case where you are essentially doing an in-silico pooling of individuals that were all GBSed separately then this can be done but there is not really precedence for how to do this statistically. I'm thinking about the best way to develop this perhaps in the future.
But that being said, I have recently done this by manually averaging the SNP-indeces for all the individuals in the Bulk and then setting up the dataframes to work with the QTLseqr scripts to compare the bulks. This worked for my case (single dominant gene) and we were able to map the gene. I did have to use the QTLseqr scripts source code a bit differently and not directly using the commands developed for the package to wrangle the data to work. Maybe this can help you
preprint here: https://www.biorxiv.org/content/10.1101/2022.04.13.487913v1.abstract
scripts here: https://github.com/bmansfeld/CMD2_project/blob/main/CMD2_mapping_and_phenotype_scripts.R
This might help you get were you want to go until I have time to develop this more in depth for the future.
Hope this helps let me know if you have questions hopefully I can help you!
Ben
from qtlseqr.
Related Issues (20)
- Help HOT 3
- G' Distribution Fitting HOT 7
- How to deal with the multi-sample ? HOT 2
- getQTLTable
- getQTLTable without running runGprimeAnalysis HOT 1
- Installation failed: error in running command when use devtools::install_github("bmansfeld/QTLseqr") to Install HOT 3
- #REEOR# importFromGATK# Error in file(file, "rt") : cannot open the connection HOT 1
- The script running normal and no error produced, why no plot generated? just only produce a table of CSV HOT 1
- Tag releases? HOT 1
- BUG? HOT 2
- Error: Problem with `mutate()` column `tricubeDeltaSNP` HOT 1
- argument 'x' contains missing values HOT 1
- windowSize change, popStruc change HOT 3
- Help me please... HOT 8
- df <- importFromTable() troubles HOT 3
- No QTL signal when running runQTLseqAnalysis() but a clear signal when running runGprimeAnalysis() HOT 3
- cant install this R package HOT 2
- runQTLseqAnalysis error HOT 4
- install error!! HOT 1
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from qtlseqr.