<p>[Old question, but possibly still relevant to answer.]</p>
<p>This will not necessarily be rare. One or the other of the mates could span the ligation junction itself, in which case your read start and end could be entirely different chromosomes. In fact, <a href="https://journals.plos.org/ploscompbiol/article?id=10.1371/journal.pcbi.1008839" rel="nofollow noreferrer">some tools</a> specifically measure the quality of Hi-C libraries by measuring the number of such ligations directly observable within reads. This will be a much larger number of reads than bin-crossing reads.</p>
<p>The safest thing to do in all such cases (IMO) is to take the read start position as the location of interest for each mate. That way, if a junction does occur, you still have an unambiguous position which is your best guess of where each end of the ligation product actually originates. This gets around both the minor issue of the bin-crossing problem and the major issue of the ligation junction problem. There is probably some imprecision here and you throw out some information, but sorting out the complex possible situations otherwise seems like a lot of hassle for minimal gain.</p>

CC BY-SA 4.0
2023-01-05T17:20:26.557Z
I think what I was referring to here were undigested contiguous reads, like [Figure 1 of the HiCUP method paper](https://f1000research.com/articles/4-1310/v1) describes. In this case, HiCUP marks this as an invalid read pair and filters them out. From my experience, there weren't enough of these occurrences (thankfully) that including or excluding them would have affected the contact matrix that much. At the time I asked this question, the qc3C method you referenced wasn't published yet, but I think that's a good reference to link here, now that it is available.
James Hawley
CC BY-SA 4.0
2022-12-18T05:57:40.153Z
2022-12-18T05:57:40.153Z
Maximilian Press
<p>Chromosome capture experiments, like Hi-C, typically bin the genome into distinct bins. A contact matrix, <span class="math-container">$M$</span>, is given by <span class="math-container">$M_{i,j} = |\{pairs | mate_1 \in bin_i, mate_2 \in bin_j\}|$</span>.</p>

<p>If one of the read mates spans a bin boundary, where should that pair be assigned?</p>

<p>For example, if <span class="math-container">$mate_1$</span> overlaps <span class="math-container">$bin_i$</span> and <span class="math-container">$mate_2$</span> overlaps both <span class="math-container">$bin_j$</span> and <span class="math-container">$bin_{j+1}$</span>, should this read pair be assigned to <span class="math-container">$M_{i,j}$</span> or <span class="math-container">$M_{i,j+1}$</span>? Should it be randomly assigned to either bin with equal probability? Should it be assigned to both?</p>

<p>I know this may be a rare event, but if Hi-C datasets are generated using billions of reads (such as <a href="https://doi.org/10/xqj" rel="nofollow noreferrer">Rao <em>et al.</em>, Cell, 2014</a>), this is bound to happen.</p>

<hi-c><chromosome-capture>
Where to assign a read pair when a mate spans a HiC bin boundary?
# Where to assign a read pair when a mate spans a HiC bin boundary?

(Asked by: username_0 on 2019-03-26)

Chromosome capture experiments, like Hi-C, typically bin the genome into distinct bins. A contact matrix, $M$, is given by $M\_{i,j} = |\{pairs | mate\_1 \in bin\_i, mate\_2 \in bin\_j\}|$.

If one of the read mates spans a bin boundary, where should that pair be assigned?

For example, if $mate\_1$ overlaps $bin\_i$ and $mate\_2$ overlaps both $bin\_j$ and $bin\_{j+1}$, should this read pair be assigned to $M\_{i,j}$ or $M\_{i,j+1}$? Should it be randomly assigned to either bin with equal probability? Should it be assigned to both?

I know this may be a rare event, but if Hi-C datasets are generated using billions of reads (such as [Rao *et al.*, Cell, 2014](https://doi.org/10/xqj)), this is bound to happen.

**Tags:** hi-c, chromosome-capture

---

## 1 answers

### Answer by username_1 on 2022-12-18. Score: 1 | Accepted answer

[Old question, but possibly still relevant to answer.]

This will not necessarily be rare. One or the other of the mates could span the ligation junction itself, in which case your read start and end could be entirely different chromosomes. In fact, [some tools](https://journals.plos.org/ploscompbiol/article?id=10.1371/journal.pcbi.1008839) specifically measure the quality of Hi-C libraries by measuring the number of such ligations directly observable within reads. This will be a much larger number of reads than bin-crossing reads.

The safest thing to do in all such cases (IMO) is to take the read start position as the location of interest for each mate. That way, if a junction does occur, you still have an unambiguous position which is your best guess of where each end of the ligation product actually originates. This gets around both the minor issue of the bin-crossing problem and the major issue of the ligation junction problem. There is probably some imprecision here and you throw out some information, but sorting out the complex possible situations otherwise seems like a lot of hassle for minimal gain.

> **Comments:**
> **username_0**: I think what I was referring to here were undigested contiguous reads, like [Figure 1 of the HiCUP method paper](https://f1000research.com/articles/4-1310/v1) describes. In this case, HiCUP marks this as an invalid read pair and filters them out. From my experience, there weren't enough of these occurrences (thankfully) that including or excluding them would have affected the contact matrix that much. At the time I asked this question, the qc3C method you referenced wasn't published yet, but I think that's a good reference to link here, now that it is available.  

---

