Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Differential gene expression analysis without reference

    I have to find the differential gene expression between two genotypes of a plant species sequenced using 454 pyrosequencer.

    There is no reference genome and the closest species -glycine max aligns poorly with the reads.

    How does one go about DE analysis in this case?

    Should i combine the reads of both genotypes, assemble them, use that as a reference genome.Then map the reads of each genotype to this reference and continue the analysis?

    Thank you.

  • #2
    Hi cerebralrust,

    I was curious if you've tried M. truncatula? This a close relative of G. max.
    Have you had any luck with genomes off of phytozome?
    Last edited by phoss; 05-03-2012, 05:11 PM.

    Comment


    • #3
      You might try one of the "no genome" assemblers like Trinity or Abyss to build a "gene" library from your data. I think those put together some consensus set of sequences assembled from your reads. The you could build a bowtie reference from those FASTA sequences and align your reads to it with bowtie. Finally you can count reads aligned to each one and compare samples using something like DESeq.

      You'll need some major computer power to run Trinity, from what I hear. That process of assembling sequences from reads is much more reasource consuming than the bowtie alignment stage.
      /* Shawn Driscoll, Gene Expression Laboratory, Pfaff
      Salk Institute for Biological Studies, La Jolla, CA, USA */

      Comment


      • #4
        De novo assembly, as sdriscoll mentioned, is the typical solution when no decent reference genome exists.

        Given that you're dealing with gene expression, I assume you have transcriptome reads, in which case you could look into Trans-ABySS, which is the transcriptome-specific version of ABySS. It is a single software pipeline that aims to assemble reads into transcripts and quantify transcript abundance, all without a reference genome. In theory you would end up with two sets of transcripts and expression levels, after which standard DE analysis could be conducted. Although finding corresponding transcripts between the two sets could be tricky...

        Software link is here and paper is here
        Last edited by jujubix; 05-03-2012, 01:26 PM.

        Comment


        • #5
          indeed. do you think that one would have to engage in a massive pairwise BLAST session between assemblies in order to match them up?

          Maybe, for that reason, it would be easiest to pool all reads into a massive FASTQ and run them through ABySS at once to get a master list of transcripts and then perform quantification through other means.
          /* Shawn Driscoll, Gene Expression Laboratory, Pfaff
          Salk Institute for Biological Studies, La Jolla, CA, USA */

          Comment


          • #6
            Yeah, at this point building a common reference via assembly is looking mighty tempting. This of course, assuming cerebralrust has the major computer power to run everything

            Comment


            • #7
              yeah. i'd be a little nervous to try it myself. but that's why i have more than one computer.
              /* Shawn Driscoll, Gene Expression Laboratory, Pfaff
              Salk Institute for Biological Studies, La Jolla, CA, USA */

              Comment


              • #8
                Hello members. Thank you for your valuable advice.

                I've run assemblies on my data using Trinity,Newbler,MIRA,velvet on my HP laptop which has 4GB RAM and i3 processsor. About 800k reads with both genotypes pooled together.
                No, i have not tried M.truncatula, phoss.I will,thanks!
                I've pooled all the reads and assembled using MIRA + CAP3.Trinity, although a really good assembler is quite bad for plant genomes.(poor annotation, poor N50 etc)
                Yes it is transcriptome. Now i suppose i will map the reads back to this 'reference', quantify and continue with the analyses.

                Thanks for the Abyss suggestion and paper, jujubix & sdriscoll. I will try it out.

                Comment

                Latest Articles

                Collapse

                • 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
                • seqadmin
                  The Impact of AI in Genomic Medicine
                  by seqadmin



                  Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
                  02-26-2024, 02:07 PM

                ad_right_rmr

                Collapse

                News

                Collapse

                Topics Statistics Last Post
                Started by seqadmin, 03-14-2024, 06:13 AM
                0 responses
                34 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-08-2024, 08:03 AM
                0 responses
                72 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-07-2024, 08:13 AM
                0 responses
                81 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-06-2024, 09:51 AM
                0 responses
                68 views
                0 likes
                Last Post seqadmin  
                Working...
                X