Question: Fusion with strand orientation
0
gravatar for mangfu100
4.8 years ago by
mangfu100710
Korea, Republic Of
mangfu100710 wrote:

Hi

I have been reading "FACTERA : a practical method for the discovery of genomic rearrangements at breakpoint resolution".

While reading this paper, there is a explanation about its process. 

http://bioinformatics.oxfordjournals.org/content/early/2014/09/01/bioinformatics.btu549/F1.expansion.html

 

anyway, In this paper, soft-clipping with orientation concept is hard to understand. So I am appreciated if you look at it.

Four orientation of R1 and R2 are possible.(Fig.1D).However, only cases 1a and 2a shown in Figure !d can generate valid fusions, as their reads have soft-clipped sequences facing opposite directions. Thus , before k-mer comparison (Fig.1C), the reverse complement of R1 is taken for cases 1b and 2b, respectively, converting them into cases 1a and 2a.

 

In this paragraph, I can confirm the 1a 2a has the opposite directions in soft-clipped reads, but why only fusion is to be opposite direction in soft-clipped. I can't understand.

Can anyone help me understand this further more with picture?

sequence next-gen alignment • 2.8k views
ADD COMMENTlink modified 4.3 years ago by Biostar ♦♦ 20 • written 4.8 years ago by mangfu100710
2
gravatar for Istvan Albert
4.8 years ago by
Istvan Albert ♦♦ 81k
University Park, USA
Istvan Albert ♦♦ 81k wrote:

First I think you are right to be confused - the way the rationale is written should be much improved. 

I believe that the answer to your question is that the situations in 1a and 2a refer to genes on the same strand getting fused. Whereas cases 1b and 2b refer to a situation where genes on different strand are fused. Thus when mapping the different paired ends reads will have appear to have different orientations.

ADD COMMENTlink written 4.8 years ago by Istvan Albert ♦♦ 81k

Thank you for your advice.

Anyway How did you know that 1a and 2a go on the same strand while the others are on the different strand?

ADD REPLYlink written 4.8 years ago by mangfu100710
1

You check their strand on genome. The genes in fusion transcript are in the same strand based RNA-seq reads aligned to transcriptome, yet the initial positioning and orientation of partner genes on genome could vary due to the fact that there is a multitude of mechanisms (large-scale mutations, e.g. deletions, inversions, translocations) that could lead to fusion.

ADD REPLYlink modified 4.8 years ago • written 4.8 years ago by mikhail.shugay3.3k
1

As mikhail.shugay also stated I suspected this from other types of genome reorganization problems. For example genomic inversions will also show up as reads in the wrong say head-to-head orientation. 

ADD REPLYlink written 4.8 years ago by Istvan Albert ♦♦ 81k

As I look at it again, I didn't understand.

why R1 and R2 in case 1a is on the same orientation?

I think that paired end reads is RF orientation and ,therefore, R1 should be reverse strand and R2 also be a forward strand. so they are not on the same orientation. 

ADD REPLYlink written 4.8 years ago by mangfu100710
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 553 users visited in the last hour