Explore
For Enterprise
Join for Free
Log In

艺术与人文
商务
计算机科学
数据科学
信息技术
健康
数学和逻辑
个人发展
物理科学与工程
社会科学
语言学习
学位和证书

浏览 Coursera 的全部课程

浏览
搜索
企业版
登录
免费加入

GenomicRanges - GRanges

Loading...

使用Bioconductor分析基因组科学数据

约翰霍普金斯大学

4（180 个评分） | 10K 名学生已注册

课程 6（共 8 门，基因组数据科学专项课程）

查看授课大纲

Learn to use tools from the Bioconductor project to perform analysis of genomic data. This is the fifth course in the Genomic Big Data Specialization from Johns Hopkins University.

您将学习的技能

Bioinformatics, Bioconductor, Genomics, R Programming

审阅

4（180 个评分）

5 stars
84 ratings
4 stars
51 ratings
3 stars
21 ratings
2 stars
9 ratings
1 star
15 ratings

AB

Apr 12, 2018

Amazing course, motivated me to use bio conductor for descriptive seq studies of DNA.

CZ

Feb 10, 2016

I learned a lot from this class and the instructor is good.

从本节课中

Week One

The class will cover how to install and use Bioconductor software. We will discuss common data structures, including ExpressionSets, SummarizedExperiment and GRanges used across several types of analyses.

What is Bioconductor7:17

Installing Bioconductor3:39

The Bioconductor Website9:55

Useful Online Resources5:08

R Base Types18:11

GRanges - Overview4:41

IRanges - Basic Usage12:11

GenomicRanges - GRanges8:44

GenomicRanges - Basic GRanges Usage8:14

GenomicRanges - seqinfo6:13

AnnotationHub8:22

Usecase: AnnotationHub and GRanges, part 112:55

Usecase: AnnotationHub and GRanges, part 213:45

教学方

Kasper Daniel Hansen, PhD
Assistant Professor, Biostatistics and Genetic Medicine

In this session we shall discuss a quite important and

very useful GRanges class, from the packets, called genomic ranges.

Let's load.

A GRange is a GRange constructor and

essentially it's very similar to an IRanges with some additional

having to do with chromosome and strength.

Gromosones in GRanges are called seqnames.

Say yes, we have a strand "+",

"-", "+" And

then the ranges.

We'll say I ranges start.

So it prints and we can see we have three ranges there.

On two different strands and on a single chromosome.

Strands in G ranges can have three different values.

They can be plus or minus which are the forward and the reverse strand.

And then they can be star which indicates

either that the strand is unknown or that it's there.

An entity that is present on both strands.

A lot of the things that we've learned about, I ranges carried directly to T

ranges we just have the additional booking keeping of the chromosome and the strand.

But there are a couple of new operations that has to do with

the fact that things now have a direction because of the strain.

We can do things, such as get

the flanking sequence on the g range,

and you can see here, that.

The flanking sequence is relative

to the direction of transcription.

Something happens in the positive strand where we go to the left.

It was a negative strand to go to the right.

There are other things such as Permodas and

a couple of other straight forward things there.

Permoda is a default permoda thing.

2,200 bases interval where your 200 bases are downstream and

2,000 bases upstream of the transaction start site.

And g-ranges since it has chromosomes operate within a universe

that is given by something called seqinfo or sequence info.

Right now we haven't given it a lot of information.

It knows that there's a chromosome 1, it doesn't know how long it is,

it doesn't know whether it's circular and what genome it is.

So can we for example give it a length, which is quite useful.

So we're going to say chromosome 1 has a length of 10.

And let's see it there.

There.

We can see the different chromosome names, given by sequence levels.

And now that we have an end of the chromosome,

we can talk two operations that are relative to the entire chromosome.

An example is gaps.

Gaps is a function that basically gives us all the stuff on the chromosome that is

not covered by a range and the g ranges.

So here we get one surprising thing.

The first thing is you can see that the star strand

is treated differently from the past in the minor strand.

There were no ranges originally on the star strand so

therefore when we check the gap, we get the entire star strand back.

Furthermore, you can see that all of the ranges end at 10,

because that's the end of the chromosome so we know, yeah, that's where it stops.

Let's try to expand this with a new chromosome.

It seems straightforward that you should be able to do something like seqnames

of gr, to be equal to chromosome one, chromosome two,

chromosome one.

But, we get an error, and

the reason we get an error is that it has recorded that there's a single chromosome

in this organism, or on this hypothetical organism.

So we need to start off by either reconstructing the object, or

we need to tell it that, hey, there's actually two different levels.

And now when it knows that the seqlevels can take two different values.

We can assign a new vector of C values.

Now we have two different chromosomes here.

Often we sought Granges and asserting order seems to be sensible here,

but it's actually relative to the order of a seqlevel.

So here we have seqlevels that are chromosome one and chromosome two,

let's reverse that, so we're essentially saying that chromosome 2 comes before

chromosome 1 and then we're going to solve two main [INAUDIBLE] and

now the seq chromosome 2 is the first range that comes out.

Things with chromosome one and chromosome ten and sort of stand out computer

sorting chromosome ten would happen before chromosome two.

We could also assign a genome to these different things so let's try to do that.

Let's say that the genome of this g range is some string, well let's say hq19.

And it will print it.

We can see it all in 19, and actually turns out if you look at the seek info

slot that every single chromosome can have its own genome.

That seems very esoteric, but that happens because say we seek from saying organisms.

Or say you're spiking in external DNA into our sequencing experiment.

This external DNA sometimes comes as our quality control.

Or spiking controls.

These external pieces of DNA may come from different organisms.

So it's a little less charitable as possible that each sequence comes from

its own genome.

The whole thing is they will all come from the same genome.

So why is this nice?

Well, this is very nice to label

your GRanges with genome because it gives it a certain kind of protection.

We are all used to dealing with data from different genome builds and

that's a, I feel probably a frequent source of errors in computational.

But here there's some building stuff.

So let's take, let's make a copy of the T ranges and

say that this new copy is actually on HG18.

Now let's do something like a find overlaps,

that we haven't really discussed fully but we'll do that in a different session.

But it works in the same way as the Iranges.

And if we do that we get basically an error out saying that, hey we're trying

to do a final overlap between two GRanges and they are incompatible G notes.

That's a really,

really nice safety feature that I would encourage people to make use of.

探索我们的目录

免费加入并获得个性化推荐、更新和优惠。

Coursera 致力于普及全世界最好的教育，它与全球一流大学和机构合作提供在线课程。

© 2019 Coursera Inc. 保留所有权利。

通过 App Store 下载

通过 Google Play 获取

Coursera
关于
管理团队
工作机会
目录
证书
MasterTrack™ 证书
学位
企业版
政府版

社区
Learners
Partners
Developers
Beta Testers
Translators

连接
博客
Facebook
领英
Twitter
Instagram
YouTube
技术博客

更多
条款
隐私
帮助
内容访问
媒体
联系我们
目录
附属公司