AnnBuilder questions...
3
0
Entering edit mode
@dipl-ing-johannes-rainer-846
Last seen 10.3 years ago
hi, i wanted to build an annotation package for the Affymetrix rat2302 chip. because i have problems to access ftp servers through the http_proxy i am behind (it works perfectly with http, but not through the ftp protocol) i downloaded the LocusLink and UniGene database files and created my UG and LL objects with the option "fromWeb=FALSE". i supposed that it would work, but i got an error "could not open connection to ftp://ftp.ncbi.nih.gov/repository/UniGene/Rn.info". i supposed the function would not try to establish connections to the internet when i submit the attribute fromWeb=FALSE and define a local url beginning with file:// . maybe this is a bug, it may be better if one can submit also the url of a local info file to get the information about the built that was used... i wrote a wrapper function that overwrites the getUGBuilt() function to read a local info file, but it took me some time to understand why the call everytime failed another question: locuslink has been substituted by the Gene database... are there plans that AnnBuilder will annotate using the new database? thanks, jo
Annotation annotate AnnBuilder Annotation annotate AnnBuilder • 2.1k views
ADD COMMENT
0
Entering edit mode
Seth Falcon ★ 7.4k
@seth-falcon-992
Last seen 10.3 years ago
Hi Jo, "Dipl.-Ing. Johannes Rainer" <johannes.rainer@tugraz.at> writes: > i wanted to build an annotation package for the Affymetrix rat2302 > chip. because i have problems to access ftp servers through the > http_proxy i am behind (it works perfectly with http, but not through > the ftp protocol) i downloaded the LocusLink and UniGene database > files and created my UG and LL objects with the option > "fromWeb=FALSE". i supposed that it would work, but i got an error > "could not open connection to > ftp://ftp.ncbi.nih.gov/repository/UniGene/Rn.info". i supposed the > function would not try to establish connections to the internet when i > submit the attribute fromWeb=FALSE and define a local url beginning > with file:// . maybe this is a bug, We will look into this. > another question: locuslink has been substituted by the Gene > database... are there plans that AnnBuilder will annotate using the > new database? Yes, we are actively working on this and expect to have something very soon. Best, + seth
ADD COMMENT
0
Entering edit mode
John Zhang ★ 2.9k
@john-zhang-6
Last seen 10.3 years ago
>i wanted to build an annotation package for the Affymetrix rat2302 >chip. because i have problems to access ftp servers through the >http_proxy i am behind (it works perfectly with http, but not through >the ftp protocol) i downloaded the LocusLink and UniGene database files >and created my UG and LL objects with the option "fromWeb=FALSE". i >supposed that it would work, but i got an error "could not open >connection to ftp://ftp.ncbi.nih.gov/repository/UniGene/Rn.info". i >supposed the function would not try to establish connections to the >internet when i submit the attribute fromWeb=FALSE and define a local >url beginning with file:// . maybe this is a bug, it may be better if >one can submit also the url of a local info file to get the information >about the built that was used... i wrote a wrapper function that >overwrites the getUGBuilt() function to read a local info file, but it >took me some time to understand why the call everytime failed The code was trying to get the build information. Each of the objects (e. g. LL, UG, ...) has to be provided with a value for built to completely turn off the connection). We will try to figure out a better way of handling this. Thanks. > >another question: locuslink has been substituted by the Gene >database... are there plans that AnnBuilder will annotate using the new >database? > >thanks, jo > >_______________________________________________ >Bioconductor mailing list >Bioconductor@stat.math.ethz.ch >https://stat.ethz.ch/mailman/listinfo/bioconductor Jianhua Zhang Department of Medical Oncology Dana-Farber Cancer Institute 44 Binney Street Boston, MA 02115-6084
ADD COMMENT
0
Entering edit mode
@dipl-ing-johannes-rainer-846
Last seen 10.3 years ago
this sounds nice! i will try it and report the outcome to you, thanks, jo Quoting Steffen Durinck <sdurinck@ebi.ac.uk>: > Hi, > > If you want, you could also try out the biomaRt package see > developmental packages. biomaRt retrieves its annotation via direct > mysql queries to Ensembl. I don't know if that goes through your > firewall... > rat2302 is present in Ensembl so you could query that for annotation > such as gene information, go, sequence,snp,... Ensembl also has > entrez-gene id's (in the package they are still called locuslink but > in the database they are already entrez-gene). > This package is still in an early phase so if you decide to use it > and have any comments, those would be very helpful for directing > further development. > > Below an example for rat230_2.. > > Cheers, > Steffen > > >library(biomaRt) > > mart<-martConnect() > - Connected to ensembl_mart_29, ensembl_mart_29, ensembl_mart_29 and > ensembl_mart_29 - > > > getGene(id="1370752_a_at",array="rat230_2",mart=mart) > An object of class "martTable" > Slot "id": > [1] "1370752_a_at" > > Slot "table": > $symbol > [1] "P53_RAT" > > $description > [1] "CELLULAR TUMOR ANTIGEN P53 (TUMOR SUPPRESSOR P53). > [Source:SWISSPROT;Acc:P10361]" > > $band > [1] "q24" > > $chromosome > [1] "10" > > $start > [1] 56420181 > > $end > [1] 56424324 > > > getGO(id="1370752_a_at",array="rat230_2",mart=mart) > An object of class "martTable" > Slot "id": > [1] "1370752_a_at" "1370752_a_at" "1370752_a_at" "1370752_a_at" > "1370752_a_at" > [6] "1370752_a_at" "1370752_a_at" > > Slot "table": > $GOID > [1] "GO:0003700" "GO:0005634" "GO:0006350" "GO:0006355" "GO:0006915" > [6] "GO:0045786" NA > > $description > [1] "transcription factor activity" > [2] "nucleus" > [3] "transcription" > [4] "regulation of transcription, DNA-dependent" > [5] "apoptosis" > [6] "negative regulation of cell cycle" > [7] NA > > $evidence > [1] "IEA" "IEA" "IEA" "IEA" "IEA" "IEA" NA > > $martID > [1] "ENSRNOG00000010756" "ENSRNOG00000010756" "ENSRNOG00000010756" > [4] "ENSRNOG00000010756" "ENSRNOG00000010756" "ENSRNOG00000010756" > [7] "ENSRNOG00000010756" > > > > > > Dipl.-Ing. Johannes Rainer wrote: > >> hi, >> >> i wanted to build an annotation package for the Affymetrix rat2302 >> chip. because i have problems to access ftp servers through the >> http_proxy i am behind (it works perfectly with http, but not >> through the ftp protocol) i downloaded the LocusLink and UniGene >> database files and created my UG and LL objects with the option >> "fromWeb=FALSE". i supposed that it would work, but i got an error >> "could not open connection to >> ftp://ftp.ncbi.nih.gov/repository/UniGene/Rn.info". i supposed the >> function would not try to establish connections to the internet when >> i submit the attribute fromWeb=FALSE and define a local url >> beginning with file:// . maybe this is a bug, it may be better if >> one can submit also the url of a local info file to get the >> information about the built that was used... i wrote a wrapper >> function that overwrites the getUGBuilt() function to read a local >> info file, but it took me some time to understand why the call >> everytime failed >> >> another question: locuslink has been substituted by the Gene >> database... are there plans that AnnBuilder will annotate using the >> new database? >> >> thanks, jo >> >> _______________________________________________ >> Bioconductor mailing list >> Bioconductor@stat.math.ethz.ch >> https://stat.ethz.ch/mailman/listinfo/bioconductor >> > >
ADD COMMENT
0
Entering edit mode
Hi all, My goal is to discover regions with multiple adjacent genes coordinately regulated. Are there any packages that would allow one to test for genomic loci that are overrepresented in a subset of genes versus the remaining genes on the chip? I have tried using the geneplotter package to visualize the genomic distribution of differentially expressed genes, however, I have found it difficult to identify candidate loci from these plots because neighboring genes are indistinguishable from each other on the output graphic. Thanks, Jeff. __________________________________ Jeffrey Rasmussen Research Consultant, Bioinformatics Department of Immunology University of Washington
ADD REPLY
0
Entering edit mode
You might look at GLAD, DNAcopy, and aCGH as possibilities, but this isn't exactly what you want, I suppose given that you have a set of differentially expressed genes already. Another tact is to look at something like WebGestalt (online tool). It allows you to give a set of genes as input and can plot where these genes are on the genome. Finally, you could use something like Fisher exact test by region or chromosome arm to determine the "significance" of clustering by region but this would require some coding on your part. Sean On Mar 29, 2005, at 4:07 PM, jeffrey rasmussen wrote: > Hi all, > > My goal is to discover regions with multiple adjacent genes > coordinately regulated. Are there any packages that would allow one to > test for genomic loci that are overrepresented in a subset of genes > versus the remaining genes on the chip? I have tried using the > geneplotter package to visualize the genomic distribution of > differentially expressed genes, however, I have found it difficult to > identify candidate loci from these plots because neighboring genes are > indistinguishable from each other on the output graphic. > > Thanks, > Jeff. > > __________________________________ > Jeffrey Rasmussen > Research Consultant, Bioinformatics > Department of Immunology > University of Washington > > _______________________________________________ > Bioconductor mailing list > Bioconductor@stat.math.ethz.ch > https://stat.ethz.ch/mailman/listinfo/bioconductor
ADD REPLY
0
Entering edit mode
It depends on what you mean by that. Two relatively easy things to do, are to look at chromosomes (possibly 1/2 of chromosomes) this information is kept in the xxxCHR environment for each data package. The other relatively easy calculation is to use xxxMAP files, which group probes by cytochrome band. These would then let you construct two way tables, one direction is expressed/not (if that is really the way you want to do things - let me suggest a Bayesian approach - EBarrays or limma) and the other direction in the table is the category (either chromosome or cytochrome band). You could then use virtually any two way table test (but do be careful with Fisher's exact test, as you may end up with a very large and hence long running computation). Doing something more (like looking for amplicons) is a bit harder - but possible, best wishes Robert On Mar 29, 2005, at 1:07 PM, jeffrey rasmussen wrote: > Hi all, > > My goal is to discover regions with multiple adjacent genes > coordinately regulated. Are there any packages that would allow one to > test for genomic loci that are overrepresented in a subset of genes > versus the remaining genes on the chip? I have tried using the > geneplotter package to visualize the genomic distribution of > differentially expressed genes, however, I have found it difficult to > identify candidate loci from these plots because neighboring genes are > indistinguishable from each other on the output graphic. > > Thanks, > Jeff. > > __________________________________ > Jeffrey Rasmussen > Research Consultant, Bioinformatics > Department of Immunology > University of Washington > > _______________________________________________ > Bioconductor mailing list > Bioconductor@stat.math.ethz.ch > https://stat.ethz.ch/mailman/listinfo/bioconductor > > +--------------------------------------------------------------------- -- ----------------+ | Robert Gentleman phone: (206) 667-7700 | | Head, Program in Computational Biology fax: (206) 667-1319 | | Division of Public Health Sciences office: M2-B865 | | Fred Hutchinson Cancer Research Center | | email: rgentlem@fhcrc.org | +--------------------------------------------------------------------- -- ----------------+
ADD REPLY
0
Entering edit mode
We have recently put up a package called "ideogram" that assists in plotting data in a genomic context. The package may be of help in your analysis and can be found at http://www.vai.org/vari/supplemental/kyle.furge/download.html. The package also includes a function named (for better or worse) "ideograb" to semi-interactively select sets of adjacent genes for further analysis. While the previously described two-table testing methods may be more appropriate for your particular analysis, in addition to aCGH, GLAD, etc. you also may want to consider the CLAC package. Also, be aware that if a cluster of differentially expressed genes is biased towards one sample group, this often indicts that a common underlying cytogenetic event has taken place (i.e. DNA copy number gain / loss). The dependence of gene expression on DNA copy number can either confound or confirm your regional discriminate analysis depending on your experimental design. -kyle > From: jeffrey rasmussen <rasmuss@u.washington.edu> > Date: Tue, 29 Mar 2005 13:07:49 -0800 (PST) > To: bioconductor@stat.math.ethz.ch > Subject: [BioC] Detecting overrepresented genomic loci > > Hi all, > > My goal is to discover regions with multiple adjacent genes coordinately > regulated. Are there any packages that would allow one to test for genomic > loci that are overrepresented in a subset of genes versus the remaining > genes on the chip? I have tried using the geneplotter package to visualize > the genomic distribution of differentially expressed genes, however, I > have found it difficult to identify candidate loci from these plots > because neighboring genes are indistinguishable from each other on the > output graphic. > > Thanks, > Jeff. > > __________________________________ > Jeffrey Rasmussen > Research Consultant, Bioinformatics > Department of Immunology > University of Washington > > _______________________________________________ > Bioconductor mailing list > Bioconductor@stat.math.ethz.ch > https://stat.ethz.ch/mailman/listinfo/bioconductor > This email message, including any attachments, is for the so...{{dropped}}
ADD REPLY

Login before adding your answer.

Traffic: 399 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6