Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Retrieving promoter sequences using gene symbol

    Hi all,

    I've been working on this for a few days and don't seem to be getting anywhere.

    I have a list of gene symbols (example list below) that I need to use to retrieve promoter sequences for a promoter analysis. Basically I want to use the gene symbol to identify a promoter (for which I expect there may be several promoters) and then use that location to retrieve 1000 nucleotides upstream and maybe 200-500 nt downstream of the promoter.

    The two main strategies I have tried were:

    1. Download extracted promoter sequences from UCSC download site. Convert gene symbol to refseq ID. Match my gene list to promoters in the pre-compiled fasta of promoter sequences.
    PROBLEM: At least some of my refseq IDs don't seem to be found in this precompiled promoter sequence dataset.

    2. Use my GTF annotation file to select promoter coordinates from my gene symbols.
    PROBLEM: My UCSC GTF files don't appear to contain 5'UTR or whole transcript intervals (only exon and intron intervals). My Annotation file does have the refseq NM_00xxxxx ID though, so I could retrieve those, but where do I find transcript intervals from that? And I only want the primary promoter for each transcript.

    If it is helpful, I can program in python - I just need specific help with the direction.

    Thanks for the help guys. I really appreciate it.

    Paul

    Appologies if this is a repost - I've seen MANY similar posts, but nothing that I've found particularly helpful (that didn't lead to a dead end).

    Example list of gene symbols for which I need promoter sequences:
    Snrpd2
    Snrpe
    Snrpg
    Snrpn
    Snx11
    Socs1
    Sod1
    Sox11
    Sox12
    Sox4
    Sphk1
    Spin2c

  • #2
    Paul: Have you tried BioMart from Ensembl? You can find some help/video's on this page.

    Comment


    • #3
      Using R...

      Ref_annotations is your gff file you have to import using the function import.gff2 (with asRangedData=FALSE)
      Ref_genome is your genome imported using read.DNAStringSet

      The following code should give you the starting base of the first annotated exon of each gene

      Code:
      B <- Ref_annotations[which(seqnames(Ref_annotations) %in% names(Ref_genome))]
      C <- B[which(strand(B) == "+")]
      f <- as.factor(elementMetadata(C)$gene_name)
      rg <- split(C,f)
      rh <- unlist(range(rg))
      end(rh) <- start(rh)
      start(rh) <- start(rh)
      names(rh) <- levels(f)
      D <- rh
      C <- B[which(strand(B) == "-")]
      f <- as.factor(elementMetadata(C)$gene_name)
      rg <- split(C,f)
      rh <- unlist(range(rg))
      start(rh) <- end(rh)
      end(rh) <- end(rh)
      names(rh) <- levels(f)
      E <- rh
      F <- sort(c(D, E))
      Then you can export F as a bed file (function export.bed)

      Hope it helps...
      Last edited by SylvainL; 07-21-2016, 06:52 AM.

      Comment


      • #4
        Hi GenoMax and SylvainL,

        Thanks so much for your suggestions and time! They were both very helpful.

        For anyone later who comes across this post - I strongly urge you to familiarize yourself with biomaRt. Its a powerful tool for extracting all kinds of useful information.

        Comment

        Latest Articles

        Collapse

        • seqadmin
          Genetic Variation in Immunogenetics and Antibody Diversity
          by seqadmin



          The field of immunogenetics explores how genetic variations influence immune responses and susceptibility to disease. In a recent SEQanswers webinar, Oscar Rodriguez, Ph.D., Postdoctoral Researcher at the University of Louisville, and Ruben Martínez Barricarte, Ph.D., Assistant Professor of Medicine at Vanderbilt University, shared recent advancements in immunogenetics. This article discusses their research on genetic variation in antibody loci, antibody production processes,...
          11-06-2024, 07:24 PM
        • seqadmin
          Choosing Between NGS and qPCR
          by seqadmin



          Next-generation sequencing (NGS) and quantitative polymerase chain reaction (qPCR) are essential techniques for investigating the genome, transcriptome, and epigenome. In many cases, choosing the appropriate technique is straightforward, but in others, it can be more challenging to determine the most effective option. A simple distinction is that smaller, more focused projects are typically better suited for qPCR, while larger, more complex datasets benefit from NGS. However,...
          10-18-2024, 07:11 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, 11-01-2024, 06:09 AM
        0 responses
        30 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 10-30-2024, 05:31 AM
        0 responses
        21 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 10-24-2024, 06:58 AM
        0 responses
        26 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 10-23-2024, 08:43 AM
        0 responses
        57 views
        0 likes
        Last Post seqadmin  
        Working...
        X