Hi, Sean I have one more question. Could you explain what the multiple sequences labeled as 'upstream1000, upstream2000, upstream5000' in the galGla3 assembly? Thanks. 2011/5/4 Sean MacEachern <sean.maceach@gmail.com> > Hi Namyoung, > > unfortunately the Chicken genome is not as complete as many of the > other mammalian genomes. There is a new release that is due to come > out later this year, but that most likely will still be missing many > of the micro chromosomes. There is hope that when single molecule > sequencing is reliably running that this may improve things. > > Sean > > On Tue, May 3, 2011 at 1:26 PM, Namyoung Jung <jnamyoung@gmail.com> wrote: > > Hi, > > I'm a graduate student at Hopkns and working with chicken BSgenome > > Has anyone noticed there are some missing chromosome like 21-28 in the > > package? > > If anyone can explian this , it'll be helpful for my project. > > Thanks. > > > > All the best, > > -- > > Namyoung Jung > > PhD candidate, Department of Cellular and Molecular Medicine > > Johns Hopkins School of Medicine, Baltimore, MD, 21205 > > njung1@jhmi.edu > > > > [[alternative HTML version deleted]] > > > > _______________________________________________ > > Bioconductor mailing list > > Bioconductor@r-project.org > > https://stat.ethz.ch/mailman/listinfo/bioconductor > > Search the archives: > http://news.gmane.org/gmane.science.biology.informatics.conductor > > > -- Namyoung Jung PhD candidate, Department of Cellular and Molecular Medicine Johns Hopkins School of Medicine, Baltimore, MD, 21205 njung1@jhmi.edu [[alternative HTML version deleted]]

Hi Namyoung, Sean, On 11-05-16 01:49 PM, Sean MacEachern wrote: > Hi, > > I've never used them, so I may be incorrect, but I believe they are > the corresponding 1000, 2000 or 5000 bases upstream of a particular > feature on the genome (eg a probe or a gene), these may be useful for > identifying potential regulatory regions or designing custom primers > etc.. Correct. Like the chromosome sequences, the "upstream" sequences were downloaded from here http://hgdownload.cse.ucsc.edu/goldenPath/galGal3/bigZips/ at the time the BSgenome.Ggallus.UCSC.galGal3 was made (i.e. a long time ago). There is some information at the above URL about these upstream sequences. Note that UCSC updates those sequences weekly (they are based on RefSeq genes) so the sequences in the BSgenome data package are probably very outdated (I'm talking about the "upstream" sequences here, the chromosome sequences for galGal3 will never change). Including these upstream sequences (when they are available) is something we've been doing for a very long time, as a convenience, for many organisms. But those days if you want something accurate/up-to-date you should consider retrieving the latest RefSeq gene/transcript/exon locations with the makeTranscriptDbFromUCSC() function (from the GenomicFeatures package), compute the genomic ranges corresponding to the upstream regions (in a GRanges object), and then use the getSeq() function from the BSgenome package to extract the upstream sequences. Cheers, H. > > Cheers, > > Sean > > On Mon, May 16, 2011 at 3:22 PM, Namyoung Jung<jnamyoung at="" gmail.com=""> wrote: >> Hi, Sean >> >> I have one more question. >> Could you explain what the multiple sequences labeled as 'upstream1000, >> upstream2000, upstream5000' in the galGla3 assembly? >> Thanks. >> >> >> 2011/5/4 Sean MacEachern<sean.maceach at="" gmail.com=""> >>> >>> Hi Namyoung, >>> >>> unfortunately the Chicken genome is not as complete as many of the >>> other mammalian genomes. There is a new release that is due to come >>> out later this year, but that most likely will still be missing many >>> of the micro chromosomes. There is hope that when single molecule >>> sequencing is reliably running that this may improve things. >>> >>> Sean >>> >>> On Tue, May 3, 2011 at 1:26 PM, Namyoung Jung<jnamyoung at="" gmail.com=""> wrote: >>>> Hi, >>>> I'm a graduate student at Hopkns and working with chicken BSgenome >>>> Has anyone noticed there are some missing chromosome like 21-28 in the >>>> package? >>>> If anyone can explian this , it'll be helpful for my project. >>>> Thanks. >>>> >>>> All the best, >>>> -- >>>> Namyoung Jung >>>> PhD candidate, Department of Cellular and Molecular Medicine >>>> Johns Hopkins School of Medicine, Baltimore, MD, 21205 >>>> njung1 at jhmi.edu >>>> >>>> [[alternative HTML version deleted]] >>>> >>>> _______________________________________________ >>>> Bioconductor mailing list >>>> Bioconductor at r-project.org >>>> https://stat.ethz.ch/mailman/listinfo/bioconductor >>>> Search the archives: >>>> http://news.gmane.org/gmane.science.biology.informatics.conductor >>>> >> >> >> >> -- >> Namyoung Jung >> PhD candidate, Department of Cellular and Molecular Medicine >> Johns Hopkins School of Medicine, Baltimore, MD, 21205 >> njung1 at jhmi.edu >> > > _______________________________________________ > Bioconductor mailing list > Bioconductor at r-project.org > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor -- Hervé Pagès Program in Computational Biology Division of Public Health Sciences Fred Hutchinson Cancer Research Center 1100 Fairview Ave. N, M1-B514 P.O. Box 19024 Seattle, WA 98109-1024 E-mail: hpages at fhcrc.org Phone: (206) 667-5791 Fax: (206) 667-1319