Ribosome Profiling data analysis pipeline: preprocessing, differential expression analysis, GO analysis and codon occupancy analysis based on Dr. Puneet Sharma workflow.
script 1-10 were executed on the cluster script 11-15 were executed locally on R R session information is available as R_session_info.txt
Ribo-seQC package used has been modified from the original ( https://github.com/Tim-Yu/RiboseQC ) Please check script 11 to see how to modify the code or use the modified package uploaded as RiboseQC
(3 sources were used to obtain sno-, sn- and r-RNA; tRNA; and rRNA):
FASTA file containing rRNA, snRNA and snoRNA sequences can be generated from Ensembl biomart:
- Database: Ensembl Genes (Version)
- Dataset: Human genes (Version)
- Click “Filters”
- Expand the “Gene” tab
- Select Transcript type: rRNA, snRNA, snoRNA
- Click “Attributes”
- Select “Sequences” (Radio button) and expand the “Sequences” tab
- Select “Unspliced (gene)”
- Expand the “Header information” tab
- Select “Gene stable ID”
- Click “Count”
- Click “Results”
- Download the FASTA file
tRNA FASTA sequences can be obtained from Genomic tRNA database:
- Under “Links to Most Viewed Genomes”, select “Homo sapiens (GRCh38/hg38)”
- On the left hand side, select “FASTA Seqs”
- Download “High confidence tRNA sequences: hg38-tRNAs.fa”
NCBI Nucleotide module to procure 45S, 28S, 18S, 5.8S, 5S rRNA sequences.
- Search: biomol_rRNA[prop] AND “Homo sapiens”[Organism].
- On the top right, click “Send to”
- Select “File”
- Select Format: FASTA
- Click “Create File”
Concatenate the 3 files into a single FASTA file to form undesired RNA file named GRCh38_p13_r_t_sno_sn_RNA_ENSEMBL_NCBI_GtRNAdb.fa.
The reference genome annotation can be procured from Ensembl by clicking on the DNA(FASTA) of Homo sapiens. Upon clicking the link you will be redirected to the Ensembl FTP server where you will encounter a bunch of files. The file you are interested in will be named “Homo_sapiens.GRCh38.dna.primary_assembly.fa.gz” (download and rename to GRCh38.dna.primary_assembly.fa.gz)
Also download the relevant GTF file from Ensembl under “Gene sets” named Homo_sapiens.GRCh38.108.gtf.gz (at the time of writing).
A FASTA file containing sequences of principle splice isoforms that are extended by 18 nt in 5’- and 3’- ends can be generated from Ensembl biomart by clicking-on/selecting the following options in order:
- Database: Ensembl Genes (Version)
- Dataset: Human genes (Version)
- Click “Filters”
- Expand the “Gene” tab
- Select Transcript type: protein_coding
- Check APPRIS annotation
- Click “Attributes”
- Select “Sequences” (Radio button) and expand the “Sequences” tab
- Select “Coding sequence”
- Check Upstream flank and provide the value “18”
- Check the Downstream flank and provide the value “18”
- Expand the “Header information” tab
- Select “Transcript stable ID” and “Gene name”
- Click “Count”
- Click “Results”
- Download the FASTA file (called GRCh38_APPRIS_CDS_18.fa in this module)