Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • Azazel
    Member
    • Oct 2010
    • 52

    #1

    Comparing detected genes. Coords for probesets of Affy GeneChip Mouse 430 2.0 Array.

    Does anyone know where I can download genomic coordinates (chromosome, strand, start, end) for all probesets of this microarray? I was searching the Affy website for a couple of hours, can't find it. The microarray in question is:

    Affymetrix GeneChip Mouse Genome 430 2.0 Array

    Admittedly this is a very microarray-centric question, but the reason I'm asking is this: I want to check if genes detected by the microarray are also detected by HT-Sequencing (not RNA-Seq actually, but similar).

    I.e. I want to know if genes detected by the microarray have also been detected by the HT-Sequencing. The question is just whether the same genes have been detected or not, NOT whether similar expression levels have been detected.

    From the HT-Sequencing, I have coordinates of sequenced tags. I want to map them to some intermediate (RefSeq) and compare with expression levels from the probesets, which is data that I also have. Because directly comparing sequenced tags against probesets is not very informative I think.

    Suggestions for alternative approaches are also very welcome.

    Thanks.
  • ETHANol
    Senior Member
    • Feb 2010
    • 308

    #2
    With the R/Bioconductor package BiomaRt use this script

    > library(biomaRt)
    > ensembl = useMart("ensembl")
    > ensembl = useMart("ensembl", dataset = "hsapiens_gene_ensembl")
    > affyids <- read.delim("mytextfile.txt", header=F,stringsAsFactors=F)[,1]
    > genelocs <- getBM(attributes = c("chromosome_name", "start_position",
    "end_position", "strand"), filters = "affy_hg_u133_plus_2", values = affyids, mart = ensembl)
    #### you can also use "affy_hg_u133_plus_2", "hgnc_symbol" to get affy probe set and gene symbol. "mytextfile.txt" is a list of Affy probe sets and you need to change the filter to your gene chip.


    > data.frame(genelocs)
    > write.table(genelocs, file = "genelocs.bed", append = FALSE, sep = "\t")
    --------------
    Ethan

    Comment

    Latest Articles

    Collapse

    • SEQadmin2
      16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
      by SEQadmin2




      NGS has revolutionized clinical microbiology by offering new ways to identify and surveil pathogens without prior knowledge of the target. One popular approach is 16S rRNA and ITS amplicon sequencing, which targets specific microbial markers. But studies are increasingly incorporating shotgun metagenomics, a technique that sequences all DNA from a sample without targeting specific microbial markers. This approach allows more detailed taxonomic classification, detection...
      Today, 06:46 AM
    • SEQadmin2
      Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
      by SEQadmin2



      Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

      Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
      ...
      10-14-2026, 09:55 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by SEQadmin2, 10-28-2026, 09:57 AM
    0 responses
    23 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 10-28-2026, 09:51 AM
    0 responses
    15 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 10-22-2026, 11:13 AM
    0 responses
    68 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 10-22-2026, 10:48 AM
    0 responses
    48 views
    0 reactions
    Last Post SEQadmin2  
    Working...