DNA Folding in Chromatin at the Supra-nucleosome Level

基因组功能的变化与染色质空间组织的变化相结合。在过去的一年里，我们在核小体水平上获得了关于乳腺癌染色质重组的详细的全基因组信息。真核生物的DNA紧密堆积在细胞核中，然而它的序列却被许多调控所必需的蛋白质因子有效地识别。为了阐明这种识别的结构机制，我们需要有关于第二级DNA组织或30纳米纤维的详细信息。为了解决这个问题，我们早先计算了两个起始染色质纤维的所有可能的构型，其中DNA连接子L=10-70bp(核小体重复长度，nr1=157-217bp)。结果，我们观察到了两类不同的构象(即拓扑异构体)，它们具有不同的DNA拓扑结构。纤维的最佳几何形状取决于连接物的长度：连接物L=10n和10n 5bp的纤维每核小体增量的DNA连接数(Lk)分别为-1.5和-1.0。换句话说，DNA超卷曲的程度与染色质中的核小体间距直接相关。这使得我们能够解决长期存在的被称为联系数悖论的矛盾。我们假设上述染色质纤维的拓扑多态可能在转录过程中发挥作用，众所周知，转录过程会产生不同水平的DNA超卷曲，从RNA聚合酶的上游和下游。对酵母基因中NRL分布的全基因组分析证实了这一假设。我们还发现，两种纤维拓扑异构体(L=10n和10n=5bp)不仅在平衡构型和平均DNA连接数上存在差异，而且在动力学上也存在差异。特别是，新型的10n 5拓扑异构体的特征是可塑性增加，这使得染色质更容易被转录因子(TF)访问。此外，基因组DNA是通过部分解开核小体来获得的。有趣的是，染色质结构、拓扑和基因调控之间的这种结构-功能关系是否也存在于人类中。为此，我们深入研究了乳腺癌(BRC)细胞中的核小体重新定位(与英国埃塞克斯大学的V.Teif合作)。我们在同一BRC患者的肿瘤和正常组织中生成了高分辨率的核小体图谱，并将这些图谱与相应的无细胞DNA(CfDNA)进行了比较。肿瘤组织的特征是(1)关键调控区的单核小体重新定位，(2)全基因组部分核小体解离增加，(3)NRL减少5-10个碱基。这些效应受核苷酸含量和DNA序列重复的存在的调节，并与DNA甲基化的差异以及稳定染色体的连接子组蛋白变异体H1.4和H1X的结合有关。肿瘤组织中NRL缩短的观察打开了一些关于癌症核小体到细胞水平重组的重要问题，这些问题需要在未来进行更多的分析。在实践方面，我们首次分析了来自同一患者的成对正常组织和肿瘤组织中的核小体定位，并为先前在cfDNA中观察到的影响提供了可能的解释，并在临床上进行了经验应用。此外，我们的发现为液体活检的NRL相关分析开辟了一个新的场所。综上所述，这项研究使我们能够为基于cfDNA的核小体分析为患者诊断、监测和分层奠定坚实的理论基础。描述这些结果的论文已提交给癌症研究；另请参阅BioRxiv预印本：https://doi.org/10.1101/2023.04.17.537031.未来的研究需要将这一方法应用于更大的患者队列和更多的癌症类型。因此，我们与W.M.Linehan(NCI泌尿外科肿瘤学分会)合作，比较了几种类型的肾癌(透明细胞癌、肉瘤样癌和嫌色肾癌)的核小体情况。初步结果表明，肿瘤组织中存在显著的核小体重排(与健康组织相比)。重要的是，我们观察到肿瘤细胞的侵袭性水平与核小体位置的变化之间存在很强的相关性，在数值上表现为核小体重复长度的减少(正在进行中)。

Changes in genome functioning are coupled with the changes in spatial organization of chromatin. In the past year, we obtained a detailed genome-wide information about chromatin reorganization in breast cancer at the nucleosome level. Eukaryotic DNA is tightly packed in the nucleus, nevertheless its sequence is effectively recognized by numerous protein factors essential for regulation. To elucidate structural mechanisms of this recognition one needs to have a detailed information about the second level of DNA organization, or 30-nm fibers. To tackle this problem, earlier we computed all possible configurations of the two-start chromatin fibers with DNA linkers L = 10 - 70 bp (nucleosome repeat length, NRL = 157 - 217 bp). As a result, we observed two different families of conformations (i.e., topoisomers) characterized by different DNA topologies. The optimal geometry of a fiber depends on the linker length: the fibers with linkers L = 10n and 10n+5 bp have DNA linking numbers per nucleosome delta(Lk) = -1.5 and -1.0, respectively. In other words, the level of DNA supercoiling is directly related to the nucleosome spacing in chromatin. This allowed us to resolve the long-standing discrepancy known as the linking-number paradox. We hypothesize that topological polymorphism of chromatin fibers described above may play a role in the process of transcription, which is known to generate different levels of DNA supercoiling upstream and downstream from RNA polymerase. A genome-wide analysis of the NRL distribution in yeast genes confirmed this assumption. We also found that the two fiber topoisomers (with L = 10n and 10n+5 bp) differ not only in their equilibrium configuration and average DNA linking number, but also in their dynamics. In particular, the novel 10n+5 topoisomer is characterized by an increased plasticity, which makes chromatin more accessible to transcription factors (TFs). In addition, genomic DNA is made accessible by partial unwrapping of nucleosomes. It is interesting to see if this structure-function relationship between chromatin structure, topology and gene regulation also exists in humans. To this aim, we thoroughly investigated nucleosome repositioning in breast cancer (BRC) cells (in collaboration with V. Teif, Essex University, UK). We generated high-resolution nucleosome maps in paired tumor and normal tissues from the same BRC patients and compared these with the corresponding cell-free DNA (cfDNA). Tumor tissues were characterized by (1) single-nucleosome repositioning at key regulatory regions, (2) genome-wide increase in partial nucleosome unwrapping and (3) decrease in NRL by 5-10 bp. These effects were modulated by the nucleotide content and the presence of DNA sequence repeats and linked to differential DNA methylation and binding of linker histone variants H1.4 and H1X which stabilize chromatosomes. The observation of NRL shortening in tumor tissues opens a number of important questions about the nucleosome-to-cell level reorganization in cancer that require additional analyses in the future. On the practical side, we analyzed for the first time nucleosome positioning in paired normal and tumor tissues from the same patients and provided a possible explanation for the effects previously observed in cfDNA and used empirically in the clinic. In addition, our findings opened a new venue for NRL-related analyses of liquid biopsies. Taken together, this study allows us to develop solid theoretical foundations for cfDNA-based nucleosomics analysis for patient diagnostics, monitoring and stratification. The paper describing these results has been submitted to Cancer Research; see also bioRxiv preprint: https://doi.org/10.1101/2023.04.17.537031. Future studies are needed to apply this methodology to larger patient cohorts and more cancer types. Therefore, we established a collaboration with W.M. Linehan (Urologic Oncology Branch, NCI) to compare nucleosome landscapes in several types of renal cell carcinoma (clear cell, sarcomatoid and chromophobe renal cell carcinoma). The initial results suggest that there are significant rearrangements of nucleosomes in the tumor tissues (compared to the healthy ones). Importantly, we observe a strong correlation between the level of aggressiveness of the tumor cells and the changes in nucleosome positioning, characterized numerically as a decrease in nucleosome repeat length (in progress).