Hello,
We are relatively new to RNAseq analysis and are seeing unusual results in our cuffdiff output, where a number of our genes are reporting “-Nan” for one treatment FPKM in the final output, while the FPKM values are seemingly fine at all other steps of the analysis. We are concerned that this may be an error, as these Nan values are not associated with the same animal across a number of different analyses (animal A is problem in comparison 1, but is fine in comparison 2 where animal 2 is now a problem) nor are they consistent when the analysis is repeated, despite the constant presence (approximately 2/3 of genes reporting Nan are consistent between repeated analyses, while 1/3 are unique). I’ve provided some additional details of our analysis below.
This study involves 65 free-ranging animals from multiple locations, seasons, sex, etc. Sequencing was performed by Illumnia HiSeq, approximately 25M 100 base single-end reads per animal.
Analyses have been performed using the Tuxedo Suite of tools in a local instance of Galaxy.
Following grooming and trimming, TopHat v 2.0.9 (Galaxy tool v 0.7) was used to map the reads from each individual to the species-specific Ensembl genome.
Cufflinks v 2.2.1 (Galaxy tool v 2.2.1.0) was then run using reference annotation, bias correction, multi read correct and cufflinks effective length correction. The FPKM values for each animal at this point in the analysis seem fine, no “-Nan” values are present in the output.
Cufflinks gtfs from all individuals were combined using Cuffmerge v 2.2.1 (Galaxy tool v 2.2.1.0) using reference annotation and sequence data.
Multiple separate cuffdiff v 2.2.1 (Galaxy tool v 2.2.1.0) analyses were run (location, season, sex, etc.) using the cuffmerge gtf file and tophat generated bam files for each animal. The cuffdiff parameters used included, geometric library normalization, pooled dispersion estimation, 0.05 FDR, Min alignment count 10, multi-read correct, bias correction and cufflinks effective length correction.
In the differential expression testing and FPKM tracking output files (both gene and transcript level), one treatment will have “-Nan” for a number (5-10%) of features, resulting in a no test status for DE. From the read group tracking file we can see that this is due to a single animal from the treatment group reporting “-Nan” FPKM values; however when looking at the cufflinks output, the FPKM value there seems fine and is well within the range of FPKMs reported for other animals.
All 65 animals are used in each cuffdiff analysis, but are assigned to different groups depending on the variable tested. As such the n for each group is different, but is always >10. The cuffmerge gtf and tophat bam files used in each cuffdiff comparison are the same. However, different animals report –Nan in the various cufdiff analyses and if we repeat an analysis, the Nan values are still present and in the same animal, but the genes reporting these are different.
I apologize for the longwinded explanation. If anyone has seen similar results and can provide any insight to help us determine if these are random errors or legitimate biological results we appreciate your response.
Thanks!
Seqanswers Leaderboard Ad
Collapse
Latest Articles
Collapse
-
by seqadmin
The COVID-19 pandemic highlighted the need for proactive pathogen surveillance systems. As ongoing threats like avian influenza and newly emerging infections continue to pose risks, researchers are working to improve how quickly and accurately pathogens can be identified and tracked. In a recent SEQanswers webinar, two experts discussed how next-generation sequencing (NGS) and machine learning are shaping efforts to monitor viral variation and trace the origins of infectious...-
Channel: Articles
03-24-2025, 11:48 AM -
-
by seqadmin
This year’s Advances in Genome Biology and Technology (AGBT) General Meeting commemorated the 25th anniversary of the event at its original venue on Marco Island, Florida. While this year’s event didn’t include high-profile musical performances, the industry announcements and cutting-edge research still drew the attention of leading scientists.
The Headliner
The biggest announcement was Roche stepping back into the sequencing platform market. In the years since...-
Channel: Articles
03-03-2025, 01:39 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 03-20-2025, 05:03 AM
|
0 responses
49 views
0 reactions
|
Last Post
by seqadmin
03-20-2025, 05:03 AM
|
||
Started by seqadmin, 03-19-2025, 07:27 AM
|
0 responses
57 views
0 reactions
|
Last Post
by seqadmin
03-19-2025, 07:27 AM
|
||
Started by seqadmin, 03-18-2025, 12:50 PM
|
0 responses
49 views
0 reactions
|
Last Post
by seqadmin
03-18-2025, 12:50 PM
|
||
Started by seqadmin, 03-03-2025, 01:15 PM
|
0 responses
200 views
0 reactions
|
Last Post
by seqadmin
03-03-2025, 01:15 PM
|