Tweet
Hi guys. Could someone show me on some example how we identify that certain node(s) in the graph represent one single contig?
-
-
I'm not completely clear on what you are asking but I think you want to know how to decide when to merge distinct nodes in the de Bruijn graph into a single node. This is done when two nodes are unambiguously connected. In other words, if two nodes x and y are connected by an edge and neither x nor y branches, then they can be merged. I've attached a simple diagram of this. In the diagram the connected blue nodes can be merged together, so can the middle red nodes. In this case 19 nodes would be merged into 5 "contigs".
Let me know if that isn't clear or if you are asking something else. -
Thank you very much. That's exactly what I'm interested in. Identifying and acquiring contigs and then using scaffolding to merge them and get the original DNA.Let me know if that isn't clear or if you are asking something else.
Ok. If I'm getting it correctly then the image below represents just 1 single contig, because this graph fragment(although it is a full graph itself) can be unambiguously assembled into just one node, namely the one having the following sequence: TAGTCGAGGCTTTAGATCCGATGAGGCTTTAGAGACAG. Or I'm getting it wrong and there are actually 4 contigs?In this case 19 nodes would be merged into 5 "contigs".Comment
-
Not quite since the graph contains a loop. Label the nodes as follows:Originally posted by bioinf View Post
x = TAGTCGAG
y = GAGGCTTTAGA
z = AGAGACAG
w = AGATCCGAGATGAG
Note that node y branches in both directions (to x/w on its left and to w/z on its right). This branch means that the graph cannot be unambiguously simplified further.
In particular, the path of the assembly that you suggest is:
x -> y -> w -> y -> z
This is ambiguous as the following is also a valid assembly:
x -> y -> w -> y -> w -> y -> z
The second assembly is the same as the first except it travels through the y/w loop twice. In general there is no way to know how many times to travel through the loop so most assemblers will output 4 contigs here.Last edited by jts; 01-10-2011, 05:05 AM.Comment
-
I see. I guess even the mate-pairs can't help in such graph. The only solution in this case is to have the information about the length of the original DNA strand. Then we can deduce the number of times the repetition occured.
What is generally done in such cases? What is the common approach?Comment
-
It varies assembler to assembler. Some will select the most likely number of copies of the repeat based on read pairs spanning the loop and the insert size distribution. Others will just built a scaffold of x,z and leave the sequence inbetween as a run of "N"s.Comment
-
Now everything is clear. Thank you.Comment
-
Simple model of De Bruijn graph
Dear jts or bioinf,
I have a trouble in understanding de Bruijn graph because the concept of graph is vague for me. What does it mean when we say k-mer acts as edge and k-1 mers act as edges? How the use of k-1 mer is used instead of only k mer? Do you have simpler model of assembly using de Bruijn graph? I would appreciate if you can help.
Thanks for your time,
Scientist1Comment
-
Sorry, my first question is:
What does it mean when we say k-mer acts as edge and k-1 mers act as nodes?Comment
-
Amplicon sequencing is a targeted approach that allows researchers to investigate specific regions of the genome. This technique is routinely used in applications such as variant identification, clinical research, and infectious disease surveillance. The amplicon sequencing process begins by designing primers that flank the regions of interest. The DNA sequences are then amplified through PCR (typically multiplex PCR) to produce amplicons complementary to the targets. RNA targets...03-21-2023, 01:49 PM -
Targeted sequencing is an effective way to sequence and analyze specific genomic regions of interest. This method enables researchers to focus their efforts on their desired targets, as opposed to other methods like whole genome sequencing that involve the sequencing of total DNA. Utilizing targeted sequencing is an attractive option for many researchers because it is often faster, more cost-effective, and only generates applicable data. While there are many approaches...03-10-2023, 05:31 AM
Comment