Wow, you really illustrate the workflow so systematically. As a beginner, I only knew a little about RNA-seq analysis, so this helps a lot! Besides, you mention and explain lots of practical details, unlike some other video providers. I gain so much from your videos! Many thanks to you!!!
@kitdordkhar49642 жыл бұрын
KZfaq is a library. You are one of the books that I frequently read. Thanks! :-)
@Bioinformagician2 жыл бұрын
Flattered and grateful :)
@UT924 күн бұрын
You are brilliant! Thank you for this crystal clear explanation!!!
@aspenkayceedowning65883 ай бұрын
Thank you so much for this. I am in my honors and I am the only one on my team that can do the bioinformatics and my classes haven't caught up yet. Thank you so much for this tutorial. I am liking, subscribing and commenting to show support so you can keep doing this
@shrishteekandoi1202 жыл бұрын
You’re the best !
@chienlee4695 Жыл бұрын
I really like your videos. This tutorial help me go over the bumpy learning curve! I am so excited that I can get the read count by myself. Thanks for putting your time to make these videos.
@josefinagarciadiaz478 Жыл бұрын
I just want to say thank you very much for doing this, im just starting out in my PhD studies and you have been really helpful!
@wolveshowlhigh5 ай бұрын
Thank you!!! You've been a lifesaver for me as a newbie to the lab. Would love to see the ATAC-seq version of this !
@Ice84letters3 ай бұрын
Thank you so much for this videos i am studying bioinformatics but we made little practice so i have many empty knowledge in some parts of this pipelines...your videos are helping me so much!! 💜
@aamirmalik77402 жыл бұрын
Awesome, thank you for doing this video [much awaited] -- so nice of you Khushbu -- Please keep making such great stuff -- Many thanks
@Luvinlife4112 жыл бұрын
Thank you, thank you for these wonderful videos! Perfect for beginners and they've helped immensely for completing my thesis!
@sivanim383010 ай бұрын
Great explanation for a beginner to get started with RNA seq! Cheers to you!
@edepoisdoenem86343 ай бұрын
Nice vídeo. I use R in Windows-based systems but your video was useful to understand the steps and the logic of this analysis.
@codewithme_19882 жыл бұрын
You r doing better than other. Keep it up
@akshayaupadhyay8063 Жыл бұрын
Hi, great videos! You mentioned that you would be detailing on how to install tools(HISAT etc) for the pre-processing steps. Could you elaborate on that, probably if there is a script for running on Mac terminal/Shell then it would be great! Thank you!
@kingsleyibeabuchi3461 Жыл бұрын
Thank you very much for this wonderful video. Keep it up. I was also able to perform most of these tasks in Galaxy
@hozifaelgadal6233 ай бұрын
thank you , this was such a help , I really loved your videos
@rahultanpure13382 жыл бұрын
Great informational video. Thanks a ton.🙏 If possible please make a video on de novo Trinity pipeline.
@Bioinformagician2 жыл бұрын
Sure, will definitely plan on creating a video covering de novo transcriptomic assembly.
@nk129875 ай бұрын
Thank you so much for this video it's very helpful. Can you please make a video of further steps like count matrix and finding deg's .
@apedike Жыл бұрын
This is so clear--love your content
@kingsleyibeabuchi3461 Жыл бұрын
Great job. Your vidoes have been very helpful. I am sincerely grateful
@bastiwala08084 ай бұрын
Hello there. It is a very informative video kindly make a tutorial video of miRNA processing from fastq to readcount. Thanks
@debashruti49845 ай бұрын
Hi! Wonderful video..makes things so easy to understand. I had one question though. After alignment when do we use Cufflinks and when to use featureCounts? What differences are there in the two tools? First one being cufflinks-cuffmerge-cuffdiff and the latter being featureCounts-Deseq/EdgeR
@kiplimosimon142911 ай бұрын
Thanks for the nice and informative session.
@felipevilicich980 Жыл бұрын
Brilliant! Thank you very much for the tutorial!
@nayanbhowmik501011 ай бұрын
Thanks for your beautiful videos! Can you make a similar video for PE reads and include the step of generating genome index files using HISAT2? That would be so helpful.
@louisepitcher61802 жыл бұрын
Amazing easy to follow video! I was Just wondering if you could clarify how the lines of code would change for paired ends reads, specifically when using the TruSeq Stranded mRNA kit for bulk RNA-seq library prep. Thanks!!
@Bioinformagician2 жыл бұрын
I am glad you liked my video! ☺️ Starting at the trimming step, Trimmomatic will be required to run on "PE" mode, followed by HISAT (change -U parameter for paired end reads). There is a strandedness parameter in both - HISAT (--rna-strandness parameter) and featurecounts (-S parameter). You can test strandedness of your protocol by running infer_experiment.py script from RSeQC package (rseqc.sourceforge.net/#infer-experiment-py).
@christianahoyindamolajemiy5531 Жыл бұрын
Thank you so much. Very understandable and easy to follow.
@Gbemi78 Жыл бұрын
You are such a blessing.
@juniorkika2545 ай бұрын
Very helpful video.Thank you!
@yaarobaltali6076 Жыл бұрын
Thank you so much for this helpful videos please keep making such a great content and I have one question, when are you going to upload a video covering FastQC analysis
@Bioinformagician Жыл бұрын
I shall make a video on FastQC soon, once I finish the videos I have currently lined up. Thanks!
@user-ns7ng2hz9w10 ай бұрын
Great and clear… Thanks. I got two questions though: 1. what’s the difference between the count matrix and TPM/FPKB/RPKB? Noting that TPM/FPKB/RPKB can’t be use for DEG, while count matrix can. 2. Why Does bowtie can’t be used for alignment to genome?
@alvaroruiztabas56272 жыл бұрын
Amazing helpful video. Congratulations!!!
@rashmitan6867 Жыл бұрын
You are doing such a great job ! Kudos! Would be willing to show the analysis of 10x spatial transcriptomics in the future ?
@histephenson007 Жыл бұрын
Brilliantly done!
@escastorage74272 жыл бұрын
cannot be happier
@ashishbhaladhare Жыл бұрын
thank you, its really informative for me.
@parnakanodia22623 ай бұрын
Hi thank you for your fabulous lecture, could you please post a lecture (on urgent basis if possible) on how to remove duplicates from fastq files? Using oxford academic bioinformatics' minirmd?
@tushardhyani39312 жыл бұрын
Thank you for this video !!
@truehighs78453 ай бұрын
Love it, thanks!!! I have a question, indeed not only RNA will yield a sequence that can be difficult to reconnect the exome itself because expression of a given mRNA, the problem that persists is correlating what expression for what mRNA, otherwise whatever counting operation you make, it will not distinguish between various RNA outputs, it will only be a composite signal of all the RNA present in your sample - even if you drill down to a single cell assay, you still have a composite signal hard to get any significance from?
@adeelabbasi7375 Жыл бұрын
please make a video on De Novo RNA seq with paired end reads - its a hot topic
@ritikabassi1998 Жыл бұрын
It was a wonderful video.
@acramulhaquekabir58522 жыл бұрын
Great video. I was wondering when are you planning to show the next steps. From count matrix to DEA
@Bioinformagician2 жыл бұрын
I have covered those steps in detail in two of my videos - 1. kzfaq.info/get/bejne/ZshiZ9Ch37PDZYE.html 2. kzfaq.info/get/bejne/heB-rbJp1t3boGQ.html
@subarnathakur6973 Жыл бұрын
Please make a video of setting up a pipeline using Trinity for denovo sequence assembly
@grsbiosciences Жыл бұрын
Hi madam, how to convert raw data.gz file to fpkm file. You have used in your other video fpkn file and proceseed later using deseq2 library. There you have taken directly tab file. Buti have fastp.no.tab file. How i should convert yhis to fpkm tab file
@arpanparichha6784 Жыл бұрын
This is really helpful
@ishakyusuf46162 жыл бұрын
Thanks for your video. It is amazing.
@ishakyusuf46162 жыл бұрын
I wander if the purpose of --rna-strandness R and -S 2 are to accelerate the hisat2 and featurecounts works?
@Bioinformagician2 жыл бұрын
--rna-strandedness provides information on orientation of sequencing reads. This video explains this concept really well: kzfaq.info/get/bejne/ldR7iraJuamdkWQ.html
@ishakyusuf4616 Жыл бұрын
@@Bioinformagician I see. thanks🤞
@sdmasroor11 ай бұрын
Thank you so much!
@ignaciomazzitelli41823 ай бұрын
Brilliant!
@amitrupani98982 жыл бұрын
Thanks for this great tutorial! Question - If you have high duplicates rate (say 70-80% in case of FFPE) pre alignment, would you remove them? There are some wet lab QC metrics like DV200, RIN, etc that also determines the effectiveness of downstream analysis might be worth considering/including for your tutorial on QC. Keep it up and Cheers!
@Bioinformagician2 жыл бұрын
A general consensus about dealing with duplicate reads is to NOT remove them. The level of duplication could be technical (FFPE for example) or could be due to higher coverage or expression levels of certain genes. My concern with duplicate reads are 1. how do we distinguish reads based on source of duplication and 2. In the process of removing those reads, we could risk losing out on true expression levels. Thank you for bringing in the wet lab aspects that determine RNA quality to my notice. Definitely very important, will surely talk about it in my RNA-Seq QC video. Thanks again!
@amitrupani98982 жыл бұрын
@@Bioinformagician Thanks for sharing your thoughts. Looking forward to upcoming content. Cheers!
@sutrishakundu21475 ай бұрын
I don't have the genome sequencing data...I need to do de novo assembly of my transcriptome data...can you please elaborate how can I do this? What additional steps I need to follow before hisat2?
@chrisdoan32102 жыл бұрын
Thank you for your videos! When I run this code: hisat2 -q --rna-strandness R -x HISAT2/grch38/genome -U data/demo_trimmed.fastq. I got error message: hisat2: command not found. Would you tell me how to fix this? Does the way we install hitsat2 matter? What is location of hisat2 should be?
@Bioinformagician2 жыл бұрын
Error is telling you it cannot find hisat2 executable. You need to provide the complete path to hisat2 (path to where hisat2 is downloaded).
@chrisdoan32102 жыл бұрын
Hi@@Bioinformagician, I tried to follow your suggestion: /Users/Chris/Bash_script/HISAT2/hisat2 -q --rna-strandness R -x HISAT2/grch38/genome -U data/demo_trimmed.fastq | samtools sort -o HISAT2/demo_trimmed.bam. I got this message: dyld[1437]: Library not loaded: @rpath/libcrypto.1.0.0.dylib. I looked for the answer on the Internet but still couldn't fix that. I stucked at this code. Would you tell me how to fix it? Thank you so much!
@rushupatel2 жыл бұрын
Good job!!!
@HopefulHopelessly Жыл бұрын
Hi, wonderful video so far (I'm only a little over halfway through) I was wondering why you chose to trim it by 10 base pairs? Was that just an arbitrary decision or is it recommended? Thank you!!
@deepanshugarg4257 Жыл бұрын
How to use featureCounts command for paired end reads? and hisat2 as well?
@kajalpanchal8239 Жыл бұрын
Thank You !
@freezingtolerance7493 Жыл бұрын
Thank you for your video. it was very helpful to learn about rna-seq workflow... just quick questions: after hisat2 (alignment) is it essential to do variant calling stage before read counts? or I can just do readcounts using featurecount or htseq-count without process of SNP.
@renqiuguoli1258 Жыл бұрын
Good show, thx a lot
@sonalidey14547 ай бұрын
Maam why raw sequence number and after adapter removal sequence number varies ... as all raw read have adapter
@Rita-t1eАй бұрын
Is there any way, i can get an hour tutorial? I am Dr. Rita Bhui working as data scientist, needs some guidance on one of my project. I am stuck at calculating FPKM values after a proper alignment. I would be happy to pay you for your time.
@rafaelamachadotugores47692 жыл бұрын
thank you, girl
@CostanzoPadovano2 жыл бұрын
Thank you for this tutorial! I hope that you will create a patreon to support your channel :)
@Bioinformagician2 жыл бұрын
Working on creating one. Until then, you could support me by buying me a coffee - www.buymeacoffee.com/bioinformagic
@anishanna51252 жыл бұрын
Thank you for the very useful tutorial. Could you please advise how to exclude the multi-mapped reads that were mapped by HISAT2? I knocked out a gene (and confirmed the KO by qPCR and western blot before sequencing) but the KO gene was not differentially expressed in my RNAseq. The only explanation that I could think of it that the KO gene sequence is very close to multiple other genes, thus multi-mapped sequences from other similar genes could masked the gene deletion. Thank you!
@Bioinformagician2 жыл бұрын
Did you quantify your reads using featureCounts or any quantification tool that excludes multi-overlap reads? Also, can you check how many reads map uniquely to the gene in your BAM files (in both WT and KO) using sam flags?
@chrisdoan32102 жыл бұрын
Hi Bioinformatician. Would you please tell me why we have many pipelines for RNA-Seq? So these steps below are not considered pipeline (workflow): 1. Check quality with FastQC 2. Trim reads with Trimmomatic 3. Align reads to the reference with STAR 4. Calculate gene hit counts with FeatureCounts 5. Compare hit counts between groups with DESeq2
@tarahsify Жыл бұрын
You're awsomee❤❤
@ShivangiSinghMBB2 жыл бұрын
Can please post a tutorial for vcf file analysis with R
@Bioinformagician2 жыл бұрын
Definitely have plans on covering this topic. Thanks for the suggestion :)
@RaviKumar-cb7tw2 ай бұрын
Can you please tell me how to install Hisat2 in my computer? I am using Macbook AIr M2 chip and Sonoma 14.5 OS
@user-uq3qh2cy9v Жыл бұрын
Hi, this video is very concise. Thanks for sharing. I'm still struggling install HISAT2 onto M1 mac. How did you do it?
@quinattasneemrafique5368 ай бұрын
The link to the raw reads is not attached in the description. Kindly provide please
@kathirvelm7552 жыл бұрын
Great mam. .I need help for building cnv pipeline.. will be a great help if u tell anyway for contacting you .
@Bioinformagician2 жыл бұрын
You can find my email ID under every video.
@sMr_Borgov10 ай бұрын
Hi, quick question, the count table display gene ID's and their counts, how do I get more gene information? (such as gene name, symbol, and annotation added ). Thanks
@mohdfahad62245 ай бұрын
if i have replicates, at which step i will have to merge the two rna seq data?
@Tekofilic Жыл бұрын
Hi can someone guide me about the error: Failed to open the annotation file ../celegan/caenorhabditis.gtf, or it’s format is incorrect or it contains no ‘exon’ features
@sofiapuvogel3663 Жыл бұрын
Thank you very much for this video! I have a question for paired end data. After doing samtools sort, the reverse read and the forward read are not consequetive reads in the sorted bam anymore. Is that okey?
@tolga1292 Жыл бұрын
When will you release the video about post alignment qc? :)
@chrisdoan32102 жыл бұрын
Would you please guide me on how to set up the tools in this tutorial for a Mac M1 user because I got a lot of difficulties and errors for about a month? Thank you so much!
@Bioinformagician2 жыл бұрын
I set this pipeline up in Mac M1 pro, you should be able to reproduce this without major difficulties. Can you share what errors/difficulties you faced?
@chrisdoan32102 жыл бұрын
@@Bioinformagician would you please tell me how to update bash on an M1 mac? I tried to do it and asked the Internet but still couldn't figure out how. I got this error: Error: Cannot install in Homebrew on ARM processor in Intel default prefix (/usr/local)! when I try: brew install bash. The default bash version on Mac is old. Thank you so much!
@chrisdoan32102 жыл бұрын
@@Bioinformagician Would you please tell me how to set up bash on Mac M1? The default bash is obsolete and I got error when trying to upgrade bash.
@Bioinformagician2 жыл бұрын
@@chrisdoan3210 You are not required to update bash. It should work fine on an M1 Mac. What error do you see when you try to run this pipeline in your bash?
@chrisdoan32102 жыл бұрын
@@Bioinformagician I don't have an error now but I think it is better to use bash in the future than zsh.
@maytelopez-cascales6113 Жыл бұрын
Hi, great tutorial, i have a question. I have aligned with star, bulk RNAseq from voxels clones, and i have 52% uniquely mapped reads, 30% of multimapped, and the rest of not aligned. Can you explain different aspects of the output of star, if they have good quality or not? with FASTQC i cant find any adapter.
@user-mn4gu7kv2d Жыл бұрын
Hi, thanks for the video! But not all the commands are working. I am using the terminal in macOS and changed to the bash shell. What's the problem and how to fix it?
@user-vk8bd1re8c3 ай бұрын
❤❤❤
@juanete692 ай бұрын
What happens with reads aligned multiple times?
@gcbicca Жыл бұрын
I run another dataset of rnaseq and the results of fastqc for this dataset is the same but with Sequence Duplication levels with a flag, this is significant in experiment of Rnaseq that analyze control vs treatment?
@ynpan910 Жыл бұрын
Hi I have watched couple of your videos and they are so informative! For this pipeline, can I ask if there are R packages that can do the same thing? bc I am not familiar with linux or bash scripts but I will be learning them. Thanks for all the efforts you put in the videos.
@shreyaslabhsetwar90559 ай бұрын
Hey, I wish to quantify human micro RNAs from double stranded fastq files. Which reference shall I use?
@serrasonmez7912 Жыл бұрын
hello, thank you for your videos, and can you please make video about the github, so uploading and presenting the projects on github?
@linkandash802510 ай бұрын
@Bioinformagician can you do an Orthofinder tutorial ?
@JaskaranSingh-om5mv Жыл бұрын
how to download the tools required for the pipeline on a M2 MacBook?
@anamikapandey4769 Жыл бұрын
first of all i would like to thank you for your youtube videos which are very interesting and understandable.please keep posting more videos like this. i have one issue while i am preparing my counts table for chip seq analysis. as when i my giving the last command of features counts then its gives error all the time and the error is "failed to open the annotation file , format is incorrect or it contains no exon features" . i am just brainstorming with it from past couple of months. i will be heartily obliged to you if you could help me regarding this .
@shivanikushwaha38398 ай бұрын
I get the same issue. Did you find a solution?
@nourarifi2642 Жыл бұрын
thank you for your video. I have a question : Can we use exactly the same pipeline to analyse single cell RNA sequencing data or it's not the same ?? thank u in advance
@mostafakarami76878 ай бұрын
I have the same question! I'm curious what would be different from bulk sequencing.
@francescosilvestro2092 Жыл бұрын
Thanks a lot for this tutorial, it helps me a lot to start with. What's your math background? Have you studied advanced statistic or something related to it?
@recordinghoney Жыл бұрын
Hi! Great video. I was wondering if there is a similar quality check for CEL files. Do you happen to know?
@Tekofilic Жыл бұрын
This is an awesome tutorial. I really appreciate it. While analysing my own RNA-seq data, I am facing these two issues, kindly someone enlighten me. (1) I hv download the trimmomatic jar file from the official website, then I don’t know why it pops the error: invalid or corrupt jar file. (2) I skipped the trimmomatic step and I jumped to step 3 for aligning it with the genome, although I have a folder containing the grch38, it throws an error: (ERR): "HISAT2/genome/grch38" does not exist Exiting now .. santools sort: failed to read header from "."
@quinattasneemrafique5368 ай бұрын
Can you please make a detailed video on Metagenomics
@Bioinformagician8 ай бұрын
Thanks for the suggestion. I am currently working on some metagenomics content. Please stay tuned!
@quinattasneemrafique5368 ай бұрын
@@Bioinformagicianthat's great news, waiting eagerly 😍
@kruthiiirao Жыл бұрын
❤
@vargas4762 Жыл бұрын
Thanks for this tutorial! Could u help me on how can i create an index based on other genome like Aedes aegypti. i want to make RNA-seq analysis of infected mosquitoes, wich file should i use as reference genome ?
@Bioinformagician Жыл бұрын
If a reference genome is not available, you could assemble the genome first using the sequence assembly from NCBI. Once assembled, you could use it as a reference genome.
@anamikapandey4769 Жыл бұрын
hello ma'am , i just wanna say that when i am running your provided pipelines with code then at the time of creating features counts , it is showing error "failed to open the annotation file , format is incorrect or it contains no exon features" . please correct me , as i have tried many ways but didn't get the answer.
@seatballmub2 ай бұрын
I figured it out. I extracted the 'Homo_sapiens.GRCh38.106.gtf.gz' file to my working directory under the filename 'homo38.106.gtf' for simplicity. Then, I entered the following command exactly into the terminal: featureCounts -S 2 -a homo38.106.gtf -o quants/demo_featurecounts.txt HISAT2/demo_trimmed.bam I don't know why this worked/was even necessary for me. You'll notice this omits the ../ preceding the filename in her example code that I've copied below. featureCounts -S 2 -a ../hg38/Homo_sapiens.GRCh38.106.gtf -o quants/demo_featurecounts.txt HISAT2/demo_trimmed.bam
@juliangrandvallet5359 Жыл бұрын
So... after this, how to analyze it on R?
@zulfakhan2655 Жыл бұрын
can we convert fasta instead of fastq to count data?
@anxu4239 Жыл бұрын
Could you please make a video of how to trans counts data into FPKM or RPKM or TPM?
@Bioinformagician Жыл бұрын
You could use tools like RSEM for read quantification which gives you FPKM and TPM values along with expected counts.