染色體相互作用密度與拓撲域相關分析

1、    染色體相互作用密度與拓撲域相關分析    許希倫摘要：生物信息學是一門交叉科學，利用計算機方法來揭示大量復雜生物數據所包含的生物學意義。染色體相關結構域是染色體上相互作用密集的一種重要結構，但目前缺少相互作用與拓撲域相關性分析。針對此問題，使用計算機分析方法，設計了相互作用密度指標，可以較好地表示相互作用地聚集程度，并分析了相互作用密度和拓撲域的特點和相關性。關鍵詞：生物信息學;染色體相互作用;染色體拓撲結構域;相互作用密度;相關性分析：tp3-05 ：a：1009-3044（2020）03-0011-021 概述在過去，人們主要依賴顯微技術來研

2、究染色體的空間組織模式。受顯微技術的限制，觀察的精度有限。隨著染色質構象捕獲技術的發(fā)展，可以通過hi-c實驗1測序得到染色體在三維空間中相互接近的片段。通過對hi-c數據的分析，可以得到染色體上兩個區(qū)域間的相互作用強度，并計算得到一種染色體結構拓撲相關結構域，簡稱拓撲域2-4。拓撲域是染色體中重要的結構單元，它廣泛存在于多種物種間，是一種保守的結構。這種結構在同一生物的不同細胞系間十分穩(wěn)定，具有很高的相似性。對拓撲域的研究將有助于我們了解染色體發(fā)揮功能的機制，同時也有助于揭露染色體結構相關疾病的發(fā)生原因，為治療相關疾病提供幫助5-7。本文中我們將分析拓撲域與相互作用的性質，并提出相互作用密度的

3、指標，來表示拓撲域中相互作用的密集程度。2 實驗分析本文中，我們將分析相互作用的性質，拓撲域的性質和相互作用密度與拓撲域的相關程度。我們使用兩種常用的細胞系imr90（人胚肺成纖維細胞），并使用di算法和icfinder算法得到的拓撲域來進行分析。為了構建染色體上的相互作用矩陣，我們首先將染色體按照40kb的窗口大小分割成一個個等長且連續(xù)的區(qū)間，可依次記為b1，b2，bn（假設分成n個區(qū)間）。之后根據hi-c數據構建的相互作用矩陣定義為mnxn其中第i行，第j列元素mi，j表示區(qū)間bi與bj的相互作用強度。2.1 相互作用分析本小節(jié)我們比較了imr90細胞中相互作用強度與距離的關系，其結果如圖

4、1所示。其橫坐標表示位點之間的距離l，實線表示兩個位點在當前距離下相互作用的平均值，距離以40kb為一個單位。虛線表示其標準差?？梢钥吹剑S著兩個位點距離的增加，之間的相互作用強度下降明顯。并且其方差在不同距離下都較大，說明即使兩個位點距離相同，其相互作用強度也有很大差異。由此我們可以得出結論，兩個位點之間的相互作用受距離影響明顯，隨著兩個位點距離變遠，其相互作用的強度會下降。2.2 染色體拓撲域分析在本節(jié)中，我們采用兩個經典di算法4和icfinder算法5得到的拓撲域，來探究拓撲域的性質。本節(jié)我們使用imr90細胞系的1號染色體為例，使用兩種不同的算法，得到其拓撲域檢測結果，并進行對比。表

5、1統(tǒng)計了兩個檢測結果的信息，可以看到，兩個結果間的拓撲域平均大小均在20個單位。2.3 相互作用密度與拓撲域分析染色體拓撲域是染色體上相互作用密集的區(qū)域，也即在染色體上的兩個位置間的相互作用相較于域外或者跨域邊界這兩種情況，在拓撲域內其相互作用會明顯較大。為了讓不同距離的兩個位置間的相互作用具有可比性，并且考慮其全局信息，我們對每個相互作用的兩個位置的距離，根據前面介紹的平均值，對其轉換得到一個規(guī)范化的相互作用矩陣。如果ni，j大于1，說明這兩個位置間的相互作用大于全局的均值。根據規(guī)范化后的矩陣，我們提出了相互作用密度的概念。相互作用密度主要用于描述一個區(qū)域內部相互作用相較于全局的密集程度，值

6、越大，說明其密集程度越大。為了驗證其有效性，我們使用imr90細胞系的1號染色體，使用di和icfinder算法得到的拓撲域，我們統(tǒng)計了其結果的基本信息，如表2所示?？梢钥吹?，拓撲域的密度均值大于1，且icfinder拓撲域中密度大于1的拓撲域比例達到74%。我們計算其拓撲域的相互作用密度的分布，得到的分布圖如圖2所示?？梢钥吹剑浯蟛糠滞負溆虻拿芏染性诖笥?的部分。由此我們可以得出結論，我們所定義的相互作用密度可以有效地體現(xiàn)拓撲域內相互作用的密集程度，密度越大，表示其密集程度越高。接下來我們繼續(xù)探究拓撲域長度對相互作用密度的影響，我們統(tǒng)計了兩個方法得到的拓撲域中，不同長度拓撲域的平均密度

7、。其結果如圖3所示，隨著拓撲域長度的增大，其相互作用密度也隨之增大。由此可以得出結論，在拓撲域內，長間隔的兩個區(qū)間的相互作用強度較短間隔相比，其與全局平均水平差異更大。3 結論在本文中，我們分析了相互作用矩陣和已有算法檢測得到的拓撲域的特點，并基于我們的發(fā)現(xiàn)，提出了相互作用密度的指標，該指標考慮到了距離的因素，可用于描述染色體上相互作用的聚集程度。通過分析該指標與已有算法檢測得到的拓撲域的關系，我們發(fā)現(xiàn)該指標有以下特點：1）已有算法檢測的拓撲域其相互作用密度大部分均大于1;2）隨著拓撲域的增大，其相互作用密度也隨之增大。參考文獻：1 lieberman-aiden e，van berkum n

8、 l，williams l，et al.comprehensive mapping of long-range interactions revealsfolding principles of the human genomej. science， 2009， 326（5950）：289-293.2張文力，高通量測序數據分析現(xiàn)狀與挑戰(zhàn)j.集成技術，2012 （3）：20-24.3呂紅強，郝樂樂，劉源，等，基于生物信息學的hi-c研究現(xiàn)狀與發(fā)展趨勢（三維基因組專刊稿件）j.遺傳，2019： 0-0.4 dixon j r，selvaraj s，yue f，et al.topological d

9、omains inmammalian genomes identified by analysis of chromatin inter-actionsj. nature， 2012， 485（7398）：376.5 de laat w. duboule d.topology of mammalian developmen-tal enhancers and their regulatory landscapesj. nature. 2013.502（7472）：499-506.6 pombo a，dillon n.three-dimensional genome architecture：players and mechanismsj. nature reviews molecular cell biol-ogy， 2015， 16（4）：245-257.7 matharu n，ahituv n.minor loops in major folds： enhancerproruoter looping， chroruatin restructuring， and their associa-tion with transcriptional regulation and diseasej. plos genet-ics， 2015， 11（12）：e1005640.8 hadd