Hidden network generating rules from partially observed complex networks

Complex biological, neuroscience, geoscience and social networks exhibit heterogeneous self-similar higher order topological structures that are usually characterized as being multifractal in nature. However, describing their topological complexity through a compact mathematical description and deciphering their topological governing rules has remained elusive and prevented a comprehensive understanding of networks. To overcome this challenge, we propose a weighted multifractal graph model capable of capturing the underlying generating rules of complex systems and characterizing their node heterogeneity and pairwise interactions. To infer the generating measure with hidden information, we introduce a variational expectation maximization framework. We demonstrate the robustness of the network generator reconstruction as a function of model properties, especially in noisy and partially observed scenarios. The proposed network generator inference framework is able to reproduce network properties, differentiate varying structures in brain networks and chromosomal interactions, and detect topologically associating domain regions in conformation maps of the human genome.

Reactivation of a developmentally silenced embryonic globin gene

The α- and β-globin loci harbor developmentally expressed genes, which are silenced throughout post-natal life. Reactivation of these genes may offer therapeutic approaches for the hemoglobinopathies, the most common single gene disorders. Here, we address mechanisms regulating the embryonically expressed α-like globin, termed ζ-globin. We show that in embryonic erythroid cells, the ζ-gene lies within a ~65 kb sub-TAD (topologically associating domain) of open, acetylated chromatin and interacts with the α-globin super-enhancer. By contrast, in adult erythroid cells, the ζ-gene is packaged within a small (~10 kb) sub-domain of hypoacetylated, facultative heterochromatin within the acetylated sub-TAD and that it no longer interacts with its enhancers. The ζ-gene can be partially re-activated by acetylation and inhibition of histone de-acetylases. In addition to suggesting therapies for severe α-thalassemia, these findings illustrate the general principles by which reactivation of developmental genes may rescue abnormalities arising from mutations in their adult paralogues.

Hierarchical dinucleotide distribution in genome along evolution and its effect on chromatin packing

Dinucleotide densities and their distribution patterns vary significantly among species. Previous studies revealed that CpG is susceptible to methylation, enriched at topologically associating domain boundaries and its distribution along the genome correlates with chromatin compartmentalization. However, the multi-scale organizations of CpG in the linear genome, their role in chromatin organization, and how they change along the evolution are only partially understood. By comparing the CpG distribution at different genomic length scales, we quantify the difference between the CpG distributions of different species and evaluate how the hierarchical uneven CpG distribution appears in evolution. The clustering of species based on the CpG distribution is consistent with the phylogenetic tree. Interestingly, we found the CpG distribution and chromatin structure to be correlated in many different length scales, especially for mammals and avians, consistent with the mosaic CpG distribution in the genomes of these species.

Topologically Associating Domain Boundaries are Commonly Required for Normal Genome Function

Topologically associating domain (TAD) boundaries are thought to partition the genome into distinct regulatory territories. Anecdotal evidence suggests that their disruption may interfere with normal gene expression and cause disease phenotype, but the overall extent to which this occurs remains unknown. Here we show that TAD boundary deletions commonly disrupt normal genome function in vivo. We used CRISPR genome editing in mice to individually delete eight TAD boundaries (11-80kb in size) from the genome in mice. All deletions examined resulted in at least one detectable molecular or organismal phenotype, which included altered chromatin interactions or gene expression, reduced viability, and anatomical phenotypes. For 5 of 8 (62%) loci examined, boundary deletions were associated with increased embryonic lethality or other developmental phenotypes. For example, a TAD boundary deletion near Smad3/Smad6 caused complete embryonic lethality, while a deletion near Tbx5/Lhx5 resulted in a severe lung malformation. Our findings demonstrate the importance of TAD boundary sequences for in vivo genome function and suggest that noncoding deletions affecting TAD boundaries should be carefully considered for potential pathogenicity in clinical genetics screening.

Comparative genome analyses highlight transposon-mediated genome expansion shapes the evolutionary architecture of 3D genomic folding in cotton

Transposable element (TE) amplification has been recognized as a driving force mediating genome size expansion and evolution, but the effect on shaping of 3D genomic architecture remains largely unknown in plants. Here, we report three reference-grade cotton genome assemblies of Gossypium rotundifolium (K2), G. arboreum (A2) and G. raimondii (D5) using Oxford Nanopore sequencing technology. Comparative genome analyses document the details of lineage-specific TE amplification contributing to three-fold change of genome size (K2, 2.44 Gb; A2, 1.62 Gb; D5, 750.19 Mb), and indicate a relatively conserved gene content and synteny relationship among genomes. We found that approximately 17% of syntenic genes exhibit chromatin status switching of A/B compartment, and active TE amplification increases the proportion of A compartment in gene regions in K2 and A2 relative to D5. We also found that only 42% of topologically associating domain (TAD) boundaries were conserved by comparing three genomes, and abundant TE amplification was linked to the organization of lineage-specific TADs. This study sheds light on the regulatory role of transposon-mediated genome expansion in the evolution of higher-order chromatin structure in plants.

DiffGR: Detecting differentially interacting genomic regions from Hi-C contact maps

Recent advances in Hi-C techniques have allowed us to map genome-wide chromatin interactions and uncover higher-order chromatin structures, thereby shedding light on the principles of genome architecture and functions. However, statistical methods for detecting changes in large-scale chromatin organization such as topologically-associating domain (TAD) are still lacking. We proposed a new statistical method, DiffGR, for detecting differentially interacting genomic regions at the TAD level between Hi-C contact maps. We utilized the stratum-adjusted correlation coefficient (SCC) to measure similarity of local TAD regions. We then developed a nonparametric approach to identify statistically significant changes of genomic interacting regions. Through simulation studies, we demonstrated that DiffGR can robustly and effectively discover differential genomic regions under various conditions. Furthermore, we successfully revealed cell type-specific changes in genomic interacting regions using real Hi-C datasets. DiffGR is publicly available at https://github.com/wmalab/DiffGR.

Tight DNA-protein complexes isolated from barley seedlings are rich in potential guanine quadruplex sequences

Background The concept of chromatin domains attached to the nuclear matrix is being revisited, with nucleus described as a set of topologically associating domains. The significance of the tightly bound to DNA proteins (TBP), a protein group that remains attached to DNA after its deproteinization should be also revisited, as the existence of these interactions is in good agreement with the concept of the topologically associating domain. The work aimed to characterize the DNA component of TBP isolated from barley seedlings. Methods The tight DNA-protein complexes from the first leaves, coleoptiles, and roots of barley seedlings were isolated by purification with chromatography on nitrocellulose or exhaustive digestion of DNA with DNase I. Cloning and transformation were performed using pMOSBBlue Blunt Ended Cloning Kit. Inserts were amplified by PCR, and sequencing was performed on the MegaBace 1000 Sequencing System. The BLAST search was performed using sequence databases at NCBI, CR-EST, and TREP and Ensembl Plants databases. Comparison to MAR/SAR sequences was performed using http://smartdb.bioinf.med.uni-goettingen.de/cgi-bin/SMARtDB/smar.cgi database. The prediction of G quadruplexes (GQ) was performed with the aid of R-studio library pqsfinder. CD spectra were recorded on a Chirascan CS/3D spectrometer. Results Although the barley genome is AT-rich (43% of GC pairs), most DNA fragments associated with TBP were GC-rich (up to 70% in some fractions). Both fractionation procedures yielded a high proportion of CT-motif sequences presented predominantly by the 16-bp CC(TCTCCC)2 TC fragment present in clones derived from the TBP-bound DNA and absent in free DNA. BLAST analysis revealed alignment with different barley repeats. Some clones, however, aligned with both nuclear and chloroplast structural genes. Alignments with MAR/SAR motifs were very few. The analysis produced by the pqsfinder program revealed numerous potential quadruplex-forming sites in the TBP-bound sequences. A set of oligonucleotides containing sites of possible GQs were designed and ordered. Three of them represented the minus strand of the CT-repeat. Two were derived from sequences of two clones of nitrocellulose retained fraction from leaves and contained GC-rich motifs different from the CT motif. Circular dichroism spectroscopy revealed profound changes in spectra when oligonucleotides were incubated with 100 mM KCl. There was either an increase of positive band in the area of 260 nm or the formation of a positive band at 290 nm. In the former case, changes are typical for parallel G-quadruplexes and, in the latter, 3 + 1 structures. Discussion The G-quadruplexes anchor proteins are probably involved in the maintenance of the topologically associated domain structure.

Constitutively bound CTCF sites maintain 3D chromatin architecture and long-range epigenetically regulated domains

The architectural protein CTCF is a mediator of chromatin conformation, but how CTCF binding to DNA is orchestrated to maintain long-range gene expression is poorly understood. Here we perform RNAi knockdown to reduce CTCF levels and reveal a shared subset of CTCF-bound sites are robustly resistant to protein depletion. The ‘persistent’ CTCF sites are enriched at domain boundaries and chromatin loops constitutive to all cell types. CRISPR-Cas9 deletion of 2 persistent CTCF sites at the boundary between a long-range epigenetically active (LREA) and silenced (LRES) region, within the Kallikrein (KLK) locus, results in concordant activation of all 8 KLK genes within the LRES region. CTCF genome-wide depletion results in alteration in Topologically Associating Domain (TAD) structure, including the merging of TADs, whereas TAD boundaries are not altered where persistent sites are maintained. We propose that the subset of essential CTCF sites are involved in cell-type constitutive, higher order chromatin architecture.

NSD2 overexpression drives clustered chromatin and transcriptional changes in a subset of insulated domains

CTCF and cohesin play a key role in organizing chromatin into topologically associating domain (TAD) structures. Disruption of a single CTCF binding site is sufficient to change chromosomal interactions leading to alterations in chromatin modifications and gene regulation. However, the extent to which alterations in chromatin modifications can disrupt 3D chromosome organization leading to transcriptional changes is unknown. In multiple myeloma, a 4;14 translocation induces overexpression of the histone methyltransferase, NSD2, resulting in expansion of H3K36me2 and shrinkage of antagonistic H3K27me3 domains. Using isogenic cell lines producing high and low levels of NSD2, here we find oncogene activation is linked to alterations in H3K27ac and CTCF within H3K36me2 enriched chromatin. A logistic regression model reveals that differentially expressed genes are significantly enriched within the same insulated domain as altered H3K27ac and CTCF peaks. These results identify a bidirectional relationship between 2D chromatin and 3D genome organization in gene regulation.