Seqanswers Leaderboard Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • Aholton
    Member
    • Jul 2012
    • 23

    New to Bioinformatics

    Hello, I am new to Bioinformatics and it's software. I was wondering if I could ask a sample question and see what discussions came up so I could learn a bit more about it.

    Let's say I have 2 RNA sequences from humans in FASTQ format and wish to do a comparison of the two through gene expression to figure out what the differences are between the two sequences.

    I have looked at programs like the Tuxedo Suite, Galaxy, and BLAST. I would like to use our own 12 core processor to do the calculations, some website interface could be used as well.

    Thank you for all your help!

    PS: Is there anyway to get as specific as what membrane proteins are expressed in one but not the other?
  • Krish_143
    Member
    • Jan 2012
    • 45

    #2
    first you have to find which membrane proteins i.e which is mem prot or not and specific to type .. later you have to find gene expression which one was expressed or not. then by combing these two you can get result whether type of membrane proteins expressed are not.. it is one way i know.
    Krishna

    Comment

    • Krish_143
      Member
      • Jan 2012
      • 45

      #3
      you have to change subject ( what you want ) and post again then you can get more replies.. i guess.
      Krishna

      Comment

      • Aholton
        Member
        • Jul 2012
        • 23

        #4
        Thanks for the fast response, think I'll edit the title and try again. Mainly just wanted to know which program would be able to compare the two sequences and tell me what was different about them.

        Comment

        • Krish_143
          Member
          • Jan 2012
          • 45

          #5
          do you have Fastq format (raw reads) or fasta format(sequence).. ?
          Krishna

          Comment

          • Aholton
            Member
            • Jul 2012
            • 23

            #6
            They are in Fastq format, I ran them through Clean_reads to cut out some of the lesser quality parts. First time doing all this so hopefully I'm on the right track.

            Currently trying to get the Tuxedo Suite (Tophat, bowtie, etc.) on our server and see if I can't get some results from there.

            Kinda dove head first into the project. hahaha

            Comment

            • Krish_143
              Member
              • Jan 2012
              • 45

              #7
              ok thats fine.. your correct.. after that finally for expression you have to use cummeRbund, EdgeR or any other..

              for finding whether what type of membrane protein. you can paste fasta sequence http://www.imtech.res.in/raghava/pslpred/

              mix this two results.
              Krishna

              Comment

              • Aholton
                Member
                • Jul 2012
                • 23

                #8
                Ok so if I was to write a guide for this, what is the Tophat/Bowtie part actually doing to the sequence? Is it lining them up to where they are easier to read?

                Comment

                • Krish_143
                  Member
                  • Jan 2012
                  • 45

                  #9
                  simply Tophat is a tool used for mapping the short reads against the genome or transcriptiome. bowtie and bowtie2 are internally used in tophat for old and new versions respectively.

                  Yeah it was easy you can search in google scholar for articles you can understand after reading.
                  how to use why to use specifically..
                  Krishna

                  Comment

                  • Aholton
                    Member
                    • Jul 2012
                    • 23

                    #10
                    Thanks! I tried to Google it to understand more, but I was quickly lost in all the jargon that was being used.

                    Comment

                    • billstevens
                      Senior Member
                      • Mar 2012
                      • 120

                      #11
                      Hey, I know what its like.

                      I did what you are doing a few months ago. Let me try and outline it for you:

                      So Bowtie will take your reads and map them to your genome. Eukaryotic proteins use splicing, so you need a program that will figure out where your junctions, etc. are. That's where Tophat comes in. It takes all the reads that didn't map contiguously to the genome and figures out where it was spliced. At the end of this, you'll have a file called acceptedhits.bam. You are going to want to run Samtools on this, "samtools flagstat" which will give you your stats. It will tell you how many reads you have, and how many reads mapped. Obviously, the higher the number, the better.

                      You are then going to run Cufflinks on your hits.

                      Haha, on second thought, you can get almost all of this by reading the nature protocol paper by trapnell et al.



                      Enjoy

                      Comment

                      • Aholton
                        Member
                        • Jul 2012
                        • 23

                        #12
                        Thank you for the quick explanation, I'll be scouring that paper during lunch today. Hopefully it can put some things into perspective for me.

                        Comment

                        Latest Articles

                        Collapse

                        • seqadmin
                          New Genomics Tools and Methods Shared at AGBT 2025
                          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...
                          03-03-2025, 01:39 PM
                        • seqadmin
                          Investigating the Gut Microbiome Through Diet and Spatial Biology
                          by seqadmin




                          The human gut contains trillions of microorganisms that impact digestion, immune functions, and overall health1. Despite major breakthroughs, we’re only beginning to understand the full extent of the microbiome’s influence on health and disease. Advances in next-generation sequencing and spatial biology have opened new windows into this complex environment, yet many questions remain. This article highlights two recent studies exploring how diet influences microbial...
                          02-24-2025, 06:31 AM

                        ad_right_rmr

                        Collapse

                        News

                        Collapse

                        Topics Statistics Last Post
                        Started by seqadmin, 03-20-2025, 05:03 AM
                        0 responses
                        17 views
                        0 reactions
                        Last Post seqadmin  
                        Started by seqadmin, 03-19-2025, 07:27 AM
                        0 responses
                        18 views
                        0 reactions
                        Last Post seqadmin  
                        Started by seqadmin, 03-18-2025, 12:50 PM
                        0 responses
                        19 views
                        0 reactions
                        Last Post seqadmin  
                        Started by seqadmin, 03-03-2025, 01:15 PM
                        0 responses
                        186 views
                        0 reactions
                        Last Post seqadmin  
                        Working...