I was following the tutorial “DNA Methylation data analysis” and the MethylDackel tool that it says to use for Methylation bias and metric extraction can not be found. I tried to download the tutorial’s workflow and it says it is missing tools from the toolshed. Is there anyway to still use that tool? I have all my files in BAM format from using the bwameth tool but without the next tool I am stuck
Hi @ntillquist
The tool may not be installed at the public Galaxy server you are working at currently. Or, you own Galaxy if that is where you are working.
Try working at one of the servers (or the Docker instance, if available) listed for the tutorial under the pull-down menu “Available on these Galaxies”.
GTN Tutorial DNA Methylation data analysis
GTN Topic: Epigenetics Galaxy Training!
Screenshot with that menu expanded. UseGalaxy.eu is one good choice as GTN tutorials are routinely tested there. This tutorial happens to also be associated with a pre-configured Docker image – navigate back up to the topic level to find the details if that interests you.
Hope that helps!
Hi @jennaj . I am trying to run this tutorial on use galaxy.eu.
In the “Hands-on: Methylation extraction with MethylDackel section”
I cannot find the option to " * Set the parameter Extract fractional methylation (only) at each position. This is mutually exclusive with --counts, --logit, and --methylKit
to Yes
." , as mentioned in the tutorial.
Please suggest
Try clicking to open the tool directly. That points to the EU server, and opens the version of the tool supported in the tutorial. Be sure to look under advanced options too for options.
Let us know if that is not enough and can look closer. These tutorials were just updated so are current.
Hi,
@jennaj . I opened the tool from usegalaxy.eu server and executed a job.
I looked in advanced options as well but could not find the specified function to set the parameter.
The tool can be opened directly from the tutorial. That points to the specific version of the tool that matches the tutorial.
Or, you can navigate using the version’s menu Changing the tool version
I’m explaining since you might need to do this across tutorials and tools, so is good to learn how Give that a try and if still stuck, I’ll explain the steps with screenshots.
@jennaj . Indeed an older version of the tool , had the missing parameter as advanced settings that helped me .
Another question based on the course tutorial : for compute matrix it uses imported file CpGislands.bed. As per my understanding, this is for human genome . I am using mm10 methylation datasets. Please could you help me where can I find the CpGislands. bed to carry out the plotProfile and differential methylation.
UCSC https://genome.ucsc.edu/ has a track for this data associated with mm10. See their Table Browser. You can go there first and check the box for output to Galaxy, or open it from within Galaxy with the tool UCSC main.
Hi @jennaj
Thanks . I will do that . I reckon in the tutorial at no point it is mentioned to import CpGislands.bed as input and can be confusing . Perhaps adding it as data to upload prior might be useful.
On a separate note I checked the workflow and would just like to clarify few parameters.
In the tutorial :
MethylDackel tool with the following parameters:
- Choose at the first option
Load reference genome from
Local cache
and forUsing reference genome
the valuehg38
. - Select for the option
sorted_alignments.bam
the computed bam file of step 4 of thebwameth
alignment. - Use for
What do you want to do?
the valueDetermine the position-dependent methylation bias in the dataset, producing diagnostic SVG images
. - Set the parameters
By default, if only one read in a pair aligns (a singleton) then it's ignored.
andBy default, paired-end alignments with the properly-paired bit unset in the FLAG field are ignored. Note that the definition of concordant and discordant is based on your aligner settings.
toYes
.
Should we change the default settings of the above mentioned parameters and set it to “yes”.
Hello @akanksha_bafna
If your reads are paired-end, then choose the pair-end end option with this tool, too. The concordant alignments flag in the BAM will be important to avoid over counting. SAMtools stats can count those up, and Filter BAM can remove the excess reads if you want to reduce the BAM size for some reason (example: to reduce the load on downstream tools).