Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • shiny_cuffdiff: On-the-fly interrogation and plotting of gene expression results

    Hi all. I'd like to introduce the community to shiny_cuffdiff

    http://nzgl.github.io/shiny_cuffdiff/

    Put simply, this code uses the outputs of the Tuxedo workflow (Tophat2 > Cufflinks > Cuffdiff) which is a cuffdiff SQL database file called cuffdiff.db. Then, using the R package cummeRbund, various plots and outputs are created on-the fly.

    shiny_cuffdiff consists of two Shiny apps:

    gene_expression_panel_cuffdiff:

    Takes a cuffdiff database (rebuilding if required) and displays information on a single named gene (as defined by the gene short name or the XLOC number). After the gene information is defined, gene expression plots can be restricted by sample.

    gene_comparison_panel_cuffdiff:

    Takes a cuffdiff database (rebuilding if required) and displays information on a set of named genes (as defined by the gene short name or the XLOC number). Displays a heatmap and a barplot. After the gene information is defined, gene expression plots can be restricted by sample.

    This project is at an early phase of development, but is under active development.

    Participants (either users, testers, or code developers) are encouraged!


    We are also developing Shiny code to interrogate the outputs of other common gene expression / differential expression workflows (initially starting with DESeq2), which will be published as a separate repository.

    Dan Jones
    School of Biological Sciences
    The University of Auckland

    Please feel free to contact me through SeqAnswers or Github.
    Last edited by danwiththeplan; 01-24-2016, 07:42 PM.

  • #2
    Hi

    I was just testing the program out. This seems to be a great way to help non-Bioinformaticians to visualize their differential expression data.

    There is an extra ")" at the end of the code to run it through shiny in R.

    Have you considered something like this: QuickRNASeq http://sourceforge.net/projects/quickrnaseq/ or Degust: http://vicbioinformatics.com/degust/...?code=example#
    Last edited by Zapages; 01-25-2016, 10:36 AM.

    Comment


    • #3
      Hi, thanks for your reply. I'm not able to replicate this error. If you mean this error:
      Code:
      error in evaluating the argument 'object' in selecting a method for function 'expressionPlot': Error in sqliteSendQuery(con, statement, bind.data) : 
        error in statement: near ")": syntax error
      ..then this occurs when you aren't giving the program a valid cuffdiff directory (not the cuffdiff.db file, the cuffdiff directory), or the gene name or XLOC number does not exist. Note that the plotting isn't instant, if you can see cummeRbund getting the gene information:

      Code:
      Getting gene information:
      	FPKM
      	Differential Expression Data
      	Annotation Data
      	Replicate FPKMs
      	Counts
      .then it's working and the plots will pop up soon.

      This is not so much an error in the code, more an example of error handling that could be more clear and informative. As is noted, this is a program in early development and I do note that error handling is the next thing I'm working on. Additionally, the plots are not that nice, so I will be using ggplot2 to make nicer plots in the very near future.

      I'm aware of Degust and QuickRNAseq. Degust does something very similar but does not handle cuffdiff outputs at all. Valid comparison, but until they add capability to handle cuffdiff, I think this work has a place. Additionally Degust requires that you upload your data (EDIT: you can download the code and run it on your computer, so I'm wrong here)
      Not sure the comparison to QuickRNAseq is valid, since it's an entire workflow that starts right from the FASTQ file (as I understand it).
      Last edited by danwiththeplan; 01-25-2016, 01:01 PM.

      Comment


      • #4
        Really interesting project. Looking forward to contributing to it's development. Thanks for picking up this mantle.

        -Loyal

        Comment


        • #5
          Yes, I agree this work definitely is useful and has a place. As this could be specifically designed for cuffdiff and other R based tool sets.

          Currently, using gplot2 to create heatmaps takes a bit of time with R for DeSeq and EdgeR data.

          I was able to fix the error, by going into the load_packages.R files and fixing up the URL addresses. My R was not accepting the 's' in https://...

          I look forward to future versions of your tool.

          Best Regards.

          Comment


          • #6
            Originally posted by Zapages View Post
            Yes, I agree this work definitely is useful and has a place. As this could be specifically designed for cuffdiff and other R based tool sets.

            Currently, using gplot2 to create heatmaps takes a bit of time with R for DeSeq and EdgeR data.

            I was able to fix the error, by going into the load_packages.R files and fixing up the URL addresses. My R was not accepting the 's' in https://...

            I look forward to future versions of your tool.

            Best Regards.
            Great, thanks, this is now fixed. Prettier plots (using ggplot2) are definitely a plan for the near future. Hoping to get publication-level plots that can be generated on the fly. To avoid spamming Seqanswers, it might be better to keep bug reports as personal messages or in github.

            Questions and clarifications on how to use the software are fine

            Comment


            • #7
              Hi

              I just tested the program out again with all of my RNA-Seq datasets. It works really great. Thank you for creating this application.

              Maybe to help users create publication like graphs/charts, maybe there could be an option to save the images as *.svg files. This could allow the user to edit svg file (with svg editors) and have publication like results.

              Also if I were to use this for publication. How should we cite the application?

              Many thanks in advance.

              -Zain

              Comment


              • #8
                Originally posted by Zapages View Post
                Hi

                I just tested the program out again with all of my RNA-Seq datasets. It works really great. Thank you for creating this application.

                Maybe to help users create publication like graphs/charts, maybe there could be an option to save the images as *.svg files. This could allow the user to edit svg file (with svg editors) and have publication like results.

                Also if I were to use this for publication. How should we cite the application?

                Many thanks in advance.

                -Zain
                Thanks for the feedback! To answer your questions:

                1) For future releases, we will include citation information, but for now, please cite:

                Jones,D., Stuckey, A., and Fan, V. (2016). shiny_cuffdiff (version 0.1): On-the-fly interrogation and plotting of gene expression results. http://nzgl.github.io/shiny_cuffdiff/

                Remember to also cite the underlying R code (cummeRbund) and the Tuxedo workflow.

                2) Excellent suggestion regarding .svg files. I will add this to the list of features we are implementing. The ability to save graphics in various formats is definitely on the to-do list.

                Following the project on Github is a good idea, as we are actively working on new features and better error handling. For example, we have already implemented better methods of selecting the database file and are working on making the graphics more visually pleasing and customisable. We expect to release version 0.2 early next week, so keep an eye out!

                Again, thanks for the feedback, please Star us on Github!
                Last edited by danwiththeplan; 02-03-2016, 05:39 PM.

                Comment


                • #9
                  nice a bit scary to have it in a browser...
                  what the "include replicates" option means?

                  Comment

                  Latest Articles

                  Collapse

                  • seqadmin
                    Strategies for Sequencing Challenging Samples
                    by seqadmin


                    Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
                    03-22-2024, 06:39 AM
                  • seqadmin
                    Techniques and Challenges in Conservation Genomics
                    by seqadmin



                    The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

                    Avian Conservation
                    Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
                    03-08-2024, 10:41 AM

                  ad_right_rmr

                  Collapse

                  News

                  Collapse

                  Topics Statistics Last Post
                  Started by seqadmin, Yesterday, 06:37 PM
                  0 responses
                  10 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, Yesterday, 06:07 PM
                  0 responses
                  9 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 03-22-2024, 10:03 AM
                  0 responses
                  51 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 03-21-2024, 07:32 AM
                  0 responses
                  67 views
                  0 likes
                  Last Post seqadmin  
                  Working...
                  X