<p>Something else that could be happening is that contigs that are being collapsed into "heterozygous" groups. This would be a particular problem when a genome has a substantial amount of repeated sequence. Digging deep into the supplementary information of the <a href="https://academic.oup.com/gigascience/article-lookup/doi/10.1186/2047-217X-1-18" rel="nofollow noreferrer">SOAPdenovo2 paper</a>, I've found the following information:</p>

<blockquote>
  <p>In SOAPdenovo2, heterozygous contig pairs are recognized by utilizing the information of contig depth and the locality of contig. The recognized heterozygous contig pairs should obey the following rules: 1) the similarity between contigs should be high enough, for example, ≥ 95%; 2) the depth of both contigs should be near half of the average depth or all contigs, complying Poisson distribution; 3) the two contigs should be located adjacently in a scaffold and have no relationship to each other inferred by paired-end reads information. The normal contigs neighboring the heterozygous regions, if they exist, could be connected to both of the heterozygous contig pairs (H1 and H2). Only the contig with relatively higher depth in a heterozygous contig pair were kept for scaffolding. The method reduces the influence of genome heterozygosity on final scaffold length. All heterozygous contig pairs were outputted to a file to facilitate further analysis. <strong>However, the trade-off of this method is that it might incorrectly remove paralogous contigs</strong>. This problem could be relieved by a gap-filling procedure while the removed copy of paralogous contigs would be represented by gaps during scaffolding.</p>
</blockquote>

CC BY-SA 3.0
2017-06-05T10:30:04.253Z
Nice. However, this can't be a complete answer, because of "2) the depth of both contigs should be near half of the average depth or all contigs, complying Poisson distribution;". Because in the OP's case the reduction of size is from 6.3G to 2.7G, which means that more than a half of contigs would be coverage outliers.
Kamil S Jaron
CC BY-SA 3.0
2017-06-05T10:36:13.833Z
That would be true for a diploid genome that has a low amount of repeats, but if there's substantial gene duplication (as happens with very large genomes), then the collapsed sequence will be a lot smaller. However, I'm not excluding other factors, just pointing out that this is one of the possible explanations.
gringer
CC BY-SA 3.0
2017-06-05T10:59:40.503Z
I just wanted to point out, that there got to be something else as well given those numbers...
Kamil S Jaron
CC BY-SA 3.0
2017-06-05T21:25:13.250Z
Maybe also the filtering of small contigs as @manuel says.
Kamil S Jaron
CC BY-SA 3.0
2017-06-06T05:54:56.240Z
Nice. I see there are at least a few factors we agreed: low coverage of contigs, small contigs, also possible high gene duplication levels. I also have 3 mate-pair libraries with 10X each, will use them for scaffolding and gap closer. Hopefully this will tell us a bit more on this issue. i will keep you update.
user8095614
CC BY-SA 3.0
2017-06-05T02:34:11.050Z
2017-06-05T10:36:44.423Z
2017-06-05T10:36:44.423Z
gringer
<p>I use SOAPdenovo2 to assemble a large genome (4.8G) using ~20X paired-end reads. The  total length of contig sizes is 6.3G while total length of scaffolds is 2.7G. Note that this is a haploid genome, so there is no issue of heterozygosity for scaffolding. </p>

<p>I am wondering what happened during the scaffolding procedure? How comes that sum of scaffolds is so much smaller than sum of contigs?</p>

CC BY-SA 3.0
2017-06-04T14:34:54.073Z
Scaffolding probably throws out small contigs. Without being an expert in soapdenovo2, I would first bwasw the contigs to he scaffolds to see which sequence looks "lost".
Manuel
<assembly><scaffold>
What causes the difference in total length of assembled contigs and scaffolds in SOAPdenovo2?
# What causes the difference in total length of assembled contigs and scaffolds in SOAPdenovo2?

(Asked by: username_0 on 2017-06-04)

I use SOAPdenovo2 to assemble a large genome (4.8G) using ~20X paired-end reads. The total length of contig sizes is 6.3G while total length of scaffolds is 2.7G. Note that this is a haploid genome, so there is no issue of heterozygosity for scaffolding.

I am wondering what happened during the scaffolding procedure? How comes that sum of scaffolds is so much smaller than sum of contigs?

**Tags:** assembly, scaffold
> **Comments:**
> **username_1**: Scaffolding probably throws out small contigs. Without being an expert in soapdenovo2, I would first bwasw the contigs to he scaffolds to see which sequence looks "lost".  

---

## 1 answers

### Answer by username_2 on 2017-06-05. Score: 1

Something else that could be happening is that contigs that are being collapsed into "heterozygous" groups. This would be a particular problem when a genome has a substantial amount of repeated sequence. Digging deep into the supplementary information of the [SOAPdenovo2 paper](https://academic.oup.com/gigascience/article-lookup/doi/10.1186/2047-217X-1-18), I've found the following information:

> In SOAPdenovo2, heterozygous contig pairs are recognized by utilizing the information of contig depth and the locality of contig. The recognized heterozygous contig pairs should obey the following rules: 1) the similarity between contigs should be high enough, for example, ≥ 95%; 2) the depth of both contigs should be near half of the average depth or all contigs, complying Poisson distribution; 3) the two contigs should be located adjacently in a scaffold and have no relationship to each other inferred by paired-end reads information. The normal contigs neighboring the heterozygous regions, if they exist, could be connected to both of the heterozygous contig pairs (H1 and H2). Only the contig with relatively higher depth in a heterozygous contig pair were kept for scaffolding. The method reduces the influence of genome heterozygosity on final scaffold length. All heterozygous contig pairs were outputted to a file to facilitate further analysis. **However, the trade-off of this method is that it might incorrectly remove paralogous contigs**. This problem could be relieved by a gap-filling procedure while the removed copy of paralogous contigs would be represented by gaps during scaffolding.

> **Comments:**
> **username_3**: Nice. However, this can't be a complete answer, because of "2) the depth of both contigs should be near half of the average depth or all contigs, complying Poisson distribution;". Because in the OP's case the reduction of size is from 6.3G to 2.7G, which means that more than a half of contigs would be coverage outliers.  
> **username_2**: That would be true for a diploid genome that has a low amount of repeats, but if there's substantial gene duplication (as happens with very large genomes), then the collapsed sequence will be a lot smaller. However, I'm not excluding other factors, just pointing out that this is one of the possible explanations.  
> **username_3**: I just wanted to point out, that there got to be something else as well given those numbers...  
> **username_3**: Maybe also the filtering of small contigs as @manuel says.  
> **username_0**: Nice. I see there are at least a few factors we agreed: low coverage of contigs, small contigs, also possible high gene duplication levels. I also have 3 mate-pair libraries with 10X each, will use them for scaffolding and gap closer. Hopefully this will tell us a bit more on this issue. i will keep you update.  

---

