Comments (1)
Think we fixed this in Slack, feel free to re-open! https://nfcore.slack.com/archives/CFXHZD9TR/p1694596296008679
Edmund
6 months ago
That's a good question, it definitely doesn't have out of the box PINTS support, but we can tweak the params. Looks like it should be similar to GRO-seq or PRO-seq?
I'll get back to you on it when I have some more time, could you make a GitHub issue for it? https://github.com/nf-core/nascent/issues/new?template=feature_request.yml
Edmund
6 months ago
If you could find a repo with an example of processing it, that would be a huge help or a publication on how the processing is done(Might be in the one you sent, just link in the issue for the section of it please!)
東島佳毅
6 months ago
Thank you for your response. Below is the processing method that my previous colleague used. This work is published this year (PMID: 37024579). I want to analyze EU-seq data by myself with nf-core!
東島佳毅
6 months ago
EU RNA-seq differential gene expression analysis
Previously published (GSE146328)21 EU RNA-seq data of ESCs treated with DMSO or A-485 (10 µM, 30, 60 and 120 min) were reprocessed. Trimming of adapters and low-quality sequences (Phred quality score < 20) was performed using Cutadapt v.4.2 (https://doi.org/10.14806/ej.17.1.200). Read sequences were aligned to mm10 (mouse) using the Burrows–Wheeler Aligner (BWA) with default parameters (BWA v.0.7.10)59. Multi-mapped reads and reads with more than three mismatches were removed using SAMtools (v.1.4)60. Reads mapped to the ribosomal and transfer RNA region obtained from the UCSC genome browser were removed using BEDTools (v.2.23)61. For 30, 60 and 120 min time point data, based on the polymerase II elongation rates, maximum 30, 90 and 240 kb gene body regions from the TSS were used for differential gene expression analysis. The number of reads mapped to the defined regions was counted using HTSeq (v.0.11.1)62. log2 fold change and P values were calculated at individual time points using DEseq2 (v.1.32.0)63 with the default scaling method (median of relative abundance). Expressed genes were determined as described in the reference genome annotation section; low-expressed genes were filtered further if the mean EU RNA-seq read count of control and A-485 treatment condition was less than 20. 東島佳毅
6 months ago
It looks like GRO-seq or PRO-seq but I don’t think same. I will also make a GitHUB issue for this.
東島佳毅
6 months ago
I have made the comments on a GitHUB as well. Thank you very much!
Edmund
6 months ago
Awesome thank you so much!
東島佳毅
6 months ago
Hi again, I am tentatively setting the assay_type as “GROseq” and analyzing my EU-seq data. Then, I have the following error just after starting the pipeline. Do you have any idea how to solve this issue? Sorry to ask such a basic question. Thank you for your help in advance! file:///Users/higashijimayoshiki/Library/CloudStorage/Dropbox/Bioinformatics/230920_EU-seq_analysis/execution_report_2023-09-20_18-20-09.html
from nascent.
Related Issues (20)
- Bedtools sort output isn't deterministic
- BEDTOOLS_INTERSECT_FILTER is not invoked before accessing the output attribute in release version HOT 3
- Error executing process > 'NFCORE_NASCENT:NASCENT:ALIGN_BWA:BWA_MEM (1)'
- Clean up Multiqc channels
- test_full profile is broken
- "ERROR ~ fromIndex = -1" during samtools sort step HOT 6
- Add in nf-validation for samplesheet HOT 2
- Flipped bigwig strands in PROseq samples HOT 2
- Issues of creating BWA index with custom fasta and GTF
- Test results say chr22 when the ref is chr21?
- Discrepancy between bigwigs generated by bedtools and deeptools
- Saving BWA MEM log
- Consider Scatter Gather for processing speed
- Issue with concatenation of technical replicates
- Support STAR
- Example of integration with ATAC-seq pipeline output
- motif production for region
- Example of 3D genome integration
- Improve reporting when no transcripts are found
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 nascent.