RNA-Seq data preprocessing

NO N. A. Omelyanchuk
DW D. S. Wiebe
DN D. D. Novikova
VL V. G. Levitsky
NK N. Klimova
VG V. Gorelova
CW C. Weinholdt
GV G. V. Vasiliev
EZ E. V. Zemlyanskaya
NK N. A. Kolchanov
AK A. V. Kochetov
IG I. Grosse
VM V. V. Mironova
request Request a Protocol
ask Ask a question
Favorite

Raw reads were mapped onto the Arabidopsis genome (TAIR10) using SHRiMP (v2.2.3)27. Uniquely mapped reads were counted using the htseq-count program with the annotation file from TAIR1028. The DEGs were detected using the Bioconductor package DESeq29. As we investigated specificity in response to auxin, we filtered out the transcripts whose fold changes tended to zero or infinity. The genes were defined as DEGs if the Benjamini–Hochberg adjusted p-value was smaller than 0.05. Further functional annotation analysis was performed with 789 genes significantly upregulated and 659 genes downregulated by auxin.

Do you have any questions about this protocol?

Post your question to gather feedback from the community. We will also invite the authors of this article to respond.

post Post a Question
0 Q&A