Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • The DGE result can't confirm the micorarray data

    Hi,

    It is my first time to use the house-make reagents to prepare DGE-tag profiling library for sequencing. The sequencing result seem to not perfectly comfirm the microarray result and they are from the same samples (seeClick image for larger version

Name:	sample1.JPG
Views:	1
Size:	62.7 KB
ID:	307415). The eighty percent reads can map to the arabidopsis cDNA sequences, but only about fifty percent reads can map to tag table from the nearest restriction site to 3' teminal of cDNA. The tag table made by myself. Is it normal? Could you give me some advice on sample preparation?

    Maybe the insufficiency of enzyme digestion result in the low proportion of reads mapping to tag table. But I don't think it's the main reason of the low consistency with microarray data.

    Thanks.

  • #2
    Could you please precise how the "tag count" values were obtained? It seems that you are working with SAGE tags, is that correct?
    Is there a published study where SAGE-based DGE is compared with microarray data (and with a good correlation)? It could help to identify what could differ from your experiment.

    Comment


    • #3
      Originally posted by steven View Post
      Could you please precise how the "tag count" values were obtained? It seems that you are working with SAGE tags, is that correct?
      Is there a published study where SAGE-based DGE is compared with microarray data (and with a good correlation)? It could help to identify what could differ from your experiment.
      Thanks.
      We followed the protocol of "Preparing Samples for Digital Gene Expression-Tag Profiling with NlaIII" from illumina to generates a unique 16 bp for each transcript, anchored with the recognition site by the restriction enzyme NlaIII. The 16 bp sequencing data combined with the known NlaIII restriction site ("CATG") data generates the unique 20 bp tag used for annotation. The quantitative expression level of the unique transcripts is demonstrated by the number of times the sequence is detected. The tags come from the nearest restriction site to 3' teminal of arabidopsis cDNA.

      There is an official introduction in the attachmentrnaDGETagProfiling.pdf.

      Comment


      • #4
        Thanks for the info!
        Indeed, a partial digestion can explain a part of the tags that do not match to your index. Considering the penultimate restriction site of the transcripts for the index could help to get some tags back. Also, as the annotation of the ends of plant genes is less reliable than for instance in human (lack of polyA site consensus), considering potentially longer transcripts by virtually extending a bit the 3'UTRs could also help (like done here).

        However, I am afraid I can not help regarding the lack of correlation array/DGE..
        Do you by chance have technical or biological replicates of one of these experiments?
        Does the protocol only retain polyA+ long RNAs? No possible contamination by sRNAs?
        Maybe something wrong with the index table (genome/annotation version..)? Did you double check the position of a few tags "manually"?
        good luck..

        Comment


        • #5
          May be you can find the answer at this papers "H. Alexander Ebhardt etc. al., Meta-analysis of small RNA-sequencing errors reveals ubiquitous post-transcriptional RNA modifications, Nucleic Acids Research, 2009, 1–10, doi:10.1093/nar/gkp093", I am not a biologiest, so I don't know what the key point is in this matter.

          Comment


          • #6
            Originally posted by steven View Post
            Thanks for the info!
            Indeed, a partial digestion can explain a part of the tags that do not match to your index. Considering the penultimate restriction site of the transcripts for the index could help to get some tags back. Also, as the annotation of the ends of plant genes is less reliable than for instance in human (lack of polyA site consensus), considering potentially longer transcripts by virtually extending a bit the 3'UTRs could also help (like done here).

            However, I am afraid I can not help regarding the lack of correlation array/DGE..
            Do you by chance have technical or biological replicates of one of these experiments?
            Does the protocol only retain polyA+ long RNAs? No possible contamination by sRNAs?
            Maybe something wrong with the index table (genome/annotation version..)? Did you double check the position of a few tags "manually"?
            good luck..
            Thanks for your advice.
            Another sample from different treatment on the same flewcell seem to be bad as this.
            The protocol only retains polyA+ long RNAs.
            Last edited by xile; 09-15-2009, 05:34 AM.

            Comment


            • #7
              Originally posted by BENM View Post
              May be you can find the answer at this papers "H. Alexander Ebhardt etc. al., Meta-analysis of small RNA-sequencing errors reveals ubiquitous post-transcriptional RNA modifications, Nucleic Acids Research, 2009, 1–10, doi:10.1093/nar/gkp093", I am not a biologiest, so I don't know what the key point is in this matter.
              Thanks. I think the proportion of RNA modifications is lower.

              Comment


              • #8
                Also, when you annotate a tag to a transcript, do you make sure that it is the 3' most tag? We saw that though 3' most is the major tag, partial digestion causes other tags to be present as well, which need to be all summed up to get the final tag-count for the gene.
                --
                bioinfosm

                Comment


                • #9
                  Originally posted by bioinfosm View Post
                  Also, when you annotate a tag to a transcript, do you make sure that it is the 3' most tag? We saw that though 3' most is the major tag, partial digestion causes other tags to be present as well, which need to be all summed up to get the final tag-count for the gene.
                  I saw that about 20 % of total reads is not the 3' most tag.

                  Comment


                  • #10
                    Try to log2/log10 transform reads and intensity, and do plotting / linear modelling then...

                    Comment


                    • #11
                      Originally posted by xile View Post
                      I saw that about 20 % of total reads is not the 3' most tag.
                      Does including/excluding them make a difference to the correlation plot?
                      --
                      bioinfosm

                      Comment


                      • #12
                        Originally posted by jwaage View Post
                        Try to log2/log10 transform reads and intensity, and do plotting / linear modelling then...
                        Thanks. I redraw the plotClick image for larger version

Name:	log.JPG
Views:	1
Size:	71.0 KB
ID:	303143.and it's looks better. Is it normal? I have no experience for this DGE-tag data.
                        Last edited by xile; 09-16-2009, 06:49 PM.

                        Comment


                        • #13
                          Originally posted by bioinfosm View Post
                          Does including/excluding them make a difference to the correlation plot?
                          Thanks, It's not significant different.

                          Comment


                          • #14
                            Your graph looks decent. What's your correlation coefficient (R2) - our experience (and other litterature) typically shows between 55 and 70 depending on tissue and experimental setup.

                            Comment


                            • #15
                              Originally posted by jwaage View Post
                              Your graph looks decent. What's your correlation coefficient (R2) - our experience (and other litterature) typically shows between 55 and 70 depending on tissue and experimental setup.
                              After log2 transforming, the correlation coefficient is 0.6637376.

                              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
                              7 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, Yesterday, 06:07 PM
                              0 responses
                              7 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, 03-22-2024, 10:03 AM
                              0 responses
                              49 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, 03-21-2024, 07:32 AM
                              0 responses
                              66 views
                              0 likes
                              Last Post seqadmin  
                              Working...
                              X