Highly interacting regions of the human genome are enriched with enhancers and bound by DNA repair proteins

Role of the redox state of the Pirin-bound cofactor on interaction with the master regulators of inflammation and other pathways

Persistent cellular stress induced perpetuation and uncontrolled amplification of inflammatory response results in a shift from tissue repair toward collateral damage, significant alterations of tissue functions, and derangements of homeostasis which in turn can lead to a large number of acute and chronic pathological conditions, such as chronic heart failure, atherosclerosis, myocardial infarction, neurodegenerative diseases, diabetes, rheumatoid arthritis, and cancer. Keeping the vital role of balanced inflammation in maintaining tissue integrity in mind, the way to combating inflammatory diseases may be through identification and characterization of mediators of inflammation that can be targeted without hampering normal body function. Pirin (PIR) is a non-heme iron containing protein having two different conformations depending on the oxidation state of the iron. Through exploration of the Pirin interactome and using molecular docking approaches, we identified that the Fe2+-bound Pirin directly interacts with BCL3, NFKBIA, NFIX and SMAD9 with more resemblance to the native binding pose and higher affinity than the Fe3+-bound form. In addition, Pirin appears to have a function in the regulation of inflammation, the transition between the canonical and non-canonical NF-κB pathways, and the remodeling of the actin cytoskeleton. Moreover, Pirin signaling appears to have a critical role in tumor invasion and metastasis, as well as metabolic and neuro-pathological complications. There are regulatory variants in PIR that can influence expression of not only PIR but also other genes, including VEGFD and ACE2. Disparity exists between South Asian and European populations in the frequencies of variant alleles at some of these regulatory loci that may lead to differential occurrence of Pirin-mediated pathogenic conditions.

C-terminal determinants for RNA binding motif 7 protein stability and RNA recognition

The 7SK ribonucleoprotein (RNP) is a critical regulator of eukaryotic transcription. Recently, RNA binding motif 7 (RBM7) containing an RNA recognition motif (RRM) was reported to associate with 7SK RNA and core 7SK RNP protein components in response to DNA damage. However, little is known about the mode of RBM7-7SK RNA recognition. Here, we found that RRM constructs containing extended C-termini have increased solubility compared to a minimal RRM construct, although these constructs aggregate in a temperature and concentration-dependent manner. Using solution NMR dynamics experiments, we identified additional structural features observed previously in crystal but not in solution structures. To identify potential RBM7-7SK RNA binding sites, we analyzed deposited data from in cellulo crosslinking experiments and found that RBM7 primarily crosslinks to the distal region of 7SK stem-loop 3 (SL3). Electrophoretic mobility shift assays and NMR chemical shift perturbation experiments showed weak binding to 7SK SL3 constructs in vitro. Together, these results provide new insights into RBM7 RRM folding and recognition of 7SK RNA.

Long-distance association of topological boundaries through nuclear condensates

The eukaryotic genome is partitioned into distinct topological domains separated by boundary elements. Emerging data support the concept that several well-established nuclear compartments are ribonucleoprotein condensates assembled through the physical process of phase separation. Here, based on our demonstration that chemical disruption of nuclear condensate assembly weakens the insulation properties of a specific subset (∼20%) of topologically associated domain (TAD) boundaries, we report that the disrupted boundaries are characterized by a high level of transcription and striking spatial clustering. These topological boundary regions tend to be spatially associated, even interchromosomally, segregate with nuclear speckles, and harbor a specific subset of "housekeeping" genes widely expressed in diverse cell types. These observations reveal a previously unappreciated mode of genome organization mediated by conserved boundary elements harboring highly and widely expressed transcription units and associated transcriptional condensates.

3D chromatin remodelling in the germ line modulates genome evolutionary plasticity

Chromosome folding has profound impacts on gene regulation, whose evolutionary consequences are far from being understood. Here we explore the relationship between 3D chromatin remodelling in mouse germ cells and evolutionary changes in genome structure. Using a comprehensive integrative computational analysis, we (i) reconstruct seven ancestral rodent genomes analysing whole-genome sequences of 14 species representatives of the major phylogroups, (ii) detect lineage-specific chromosome rearrangements and (iii) identify the dynamics of the structural and epigenetic properties of evolutionary breakpoint regions (EBRs) throughout mouse spermatogenesis. Our results show that EBRs are devoid of programmed meiotic DNA double-strand breaks (DSBs) and meiotic cohesins in primary spermatocytes, but are associated in post-meiotic cells with sites of DNA damage and functional long-range interaction regions that recapitulate ancestral chromosomal configurations. Overall, we propose a model that integrates evolutionary genome reshuffling with DNA damage response mechanisms and the dynamic spatial genome organisation of germ cells.

The role of genome folding in the heritability and evolvability of structural variations is not well understood. Here the authors investigate the impact of the three-dimensional genome topology of germ cells in the formation and transmission of gross structural genomic changes detected from comparing whole-genome sequences of 14 rodent species.

Fast Optimal Circular Clustering and Applications on Round Genomes

Round genomes are found in bacteria, plant chloroplasts, and mitochondria. Genetic or epigenetic marks can present biologically interesting clusters along a circular genome. The circular data clustering problem groups N points on a circle into K clusters to minimize the within-cluster sum of squared distances. Repeatedly applying the K-means algorithm takes quadratic time, impractical for large circular datasets. To overcome this issue, we developed a reproducible fast optimal circular clustering (FOCC) algorithm of worst-case O(KN log² N) time. The core is a fast optimal framed clustering algorithm, which we designed by integrating two divide-and-conquer and one bracket dynamic programming strategies. The algorithm is optimal based on a property of monotonic increasing cluster borders over frames on linearized data. On clustering 50,000 circular data points, FOCC outruns brute-force or heuristic circular clustering by three orders of magnitude in time. We produced clusters of CpG sites and genes along three round genomes, exhibiting higher quality than heuristic clustering. More broadly, the presented subquadratic-time algorithms offer the fastest known solution to not only framed and circular clustering, but also angular, periodical, and looped clustering. We implemented these algorithms in the R package 'OptCirClust' (https://CRAN.R-project.org/package=OptCirClust).

Integration of DNA damage responses with dynamic spatial genome organization

The maintenance of genome stability and cellular homeostasis depends on the temporal and spatial coordination of successive events constituting the classical DNA damage response (DDR). Recent findings suggest close integration and coordination of DDR signaling with specific cellular processes. The mechanisms underlying such coordination remain unclear. We review emerging crosstalk between DNA repair factors, chromatin remodeling, replication, transcription, spatial genome organization, cytoskeletal forces, and liquid-liquid phase separation (LLPS) in mediating DNA repair. We present an overarching DNA repair framework within which these dynamic processes intersect in nuclear space over time. Collectively, this interplay ensures the efficient assembly of DNA repair proteins onto shifting genome structures to preserve genome stability and cell survival.

Disparities in COVID-19 severities and casualties across ethnic groups around the globe and patterns of ACE2 and PIR variants

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) mediated Coronavirus disease-19 (COVID-19) has affected millions of individuals around all corners of the globe. Symptoms and severities of infection with this highly contagious virus vary among individuals and there is disparity in the number of COVID-19-related casualties across different ethnic groups. The primary receptor for SARS-CoV-2 entry into the host cells is angiotensin-converting enzyme 2 (ACE2). Certain variants of ACE2 are known to be associated with COVID-19 comorbidities such as hypertension, cardiovascular complications, diabetes, chronic lung disease, etc. In this study, we looked into the geographic distribution of disease-associated variants of ACE2 as well as closely located PIR gene to explore any possible correlation with the disparities in COVID-19 severities and casualties across ethnic groups. Frequencies of the ACE2 variants associated with COVID-19 comorbidities are higher in the European and the admixed American populations. These variants are also present with stronger pairwise linkage disequilibrium (LD) in the European and the admixed American populations. On the other hand, the variants with protective role are more prevalent in the East and the South Asian populations. Strong pairwise LD exists among the activity modifying (modifier) variants of the PIR and ACE2 genes only in the European super-population. Absence of these PIR variants in the South Asian population may contribute to the overall lower COVID-19 case fatality rates (CFR) despite the dense population in this region.

Super hotspots and super coldspots in the repair of UV-induced DNA damage in the human genome

The formation of UV-induced DNA damage and its repair are influenced by many factors that modulate lesion formation and the accessibility of repair machinery. However, it remains unknown which genomic sites are prioritized for immediate repair after UV damage induction, and whether these prioritized sites overlap with hotspots of UV damage. We identified the super hotspots subject to the earliest repair for (6-4) pyrimidine-pyrimidone photoproduct by using the eXcision Repair-sequencing (XR-seq) method. We further identified super coldspots for (6-4) pyrimidine-pyrimidone photoproduct repair and super hotspots for cyclobutane pyrimidine dimer repair by analyzing available XR-seq time-course data. By integrating datasets of XR-seq, Damage-seq, adductSeq, and cyclobutane pyrimidine dimer-seq, we show that neither repair super hotspots nor repair super coldspots overlap hotspots of UV damage. Furthermore, we demonstrate that repair super hotspots are significantly enriched in frequently interacting regions and superenhancers. Finally, we report our discovery of an enrichment of cytosine in repair super hotspots and super coldspots. These findings suggest that local DNA features together with large-scale chromatin features contribute to the orders of magnitude variability in the rates of UV damage repair.

Cell-type-specific 3D epigenomes in the developing human cortex

Lineage-specific epigenomic changes during human corticogenesis have remained elusive due to challenges with sample availability and tissue heterogeneity. For example, previous studies used single-cell RNA sequencing to identify at least nine major cell types and up to 26 distinct subtypes in the dorsal cortex alone^1,2. Here, we characterize cell type-specific cis-regulatory chromatin interactions, open chromatin peaks, and transcriptomes for radial glia, intermediate progenitor cells, excitatory neurons, and interneurons isolated from mid-gestational human cortex samples. We show that chromatin interactions underlie multiple aspects of gene regulation, with transposable elements and disease-associated variants enriched at distal interacting regions in a cell type-specific manner. In addition, promoters with significantly increased levels of chromatin interactivity, termed super interactive promoters, are enriched for lineage-specific genes, suggesting that interactions at these loci contribute to the fine-tuning of transcription. Finally, we develop CRISPRview, a novel technique integrating immunostaining, CRISPRi, RNAscope, and image analysis for validating cell type-specific cis-regulatory elements in heterogeneous populations of primary cells. Our study presents the first cell type-specific characterization of 3D epigenomes in the developing human cortex, advancing our understanding of gene regulation and lineage specification during this critical developmental window.

Assessing stationary distributions derived from chromatin contact maps

Background

The spatial configuration of chromosomes is essential to various cellular processes, notably gene regulation, while architecture related alterations, such as translocations and gene fusions, are often cancer drivers. Thus, eliciting chromatin conformation is important, yet challenging due to compaction, dynamics and scale. However, a variety of recent assays, in particular Hi-C, have generated new details of chromatin structure, spawning a number of novel biological findings. Many findings have resulted from analyses on the level of native contact data as generated by the assays. Alternatively, reconstruction based approaches often proceed by first converting contact frequencies into distances, then generating a three dimensional (3D) chromatin configuration that best recapitulates these distances. Subsequent analyses can enrich contact level analyses via superposition of genomic attributes on the reconstruction. But, such advantages depend on the accuracy of the reconstruction which, absent gold standards, is inherently difficult to assess. Attempts at accuracy evaluation have relied on simulation and/or FISH imaging that typically features a handful of low resolution probes. While newly advanced multiplexed FISH imaging offers possibilities for refined 3D reconstruction accuracy evaluation, availability of such data is limited due to assay complexity and the resolution thereof is appreciably lower than the reconstructions being assessed. Accordingly, there is demand for new methods of reconstruction accuracy appraisal.

Results

Here we explore the potential of recently proposed stationary distributions, hereafter StatDns, derived from Hi-C contact matrices, to serve as a basis for reconstruction accuracy assessment. Current usage of such StatDns has focussed on the identification of highly interactive regions (HIRs): computationally defined regions of the genome purportedly involved in numerous long-range intra-chromosomal contacts. Consistent identification of HIRs would be informative with respect to inferred 3D architecture since the corresponding regions of the reconstruction would have an elevated number of k nearest neighbors (kNNs). More generally, we anticipate a monotone decreasing relationship between StatDn values and kNN distances. After initially evaluating the reproducibility of StatDns across replicate Hi-C data sets, we use this implied StatDn - kNN relationship to gauge the utility of StatDns for reconstruction validation, making recourse to both real and simulated examples.

Conclusions

Our analyses demonstrate that, as constructed, StatDns do not provide a suitable measure for assessing the accuracy of 3D genome reconstructions. Whether this is attributable to specific choices surrounding normalization in defining StatDns or to the logic underlying their very formulation remains to be determined.

The Cajal Body Protein WRAP53β Prepares the Scene for Repair of DNA Double-Strand Breaks by Regulating Local Ubiquitination

Proper repair of DNA double-strand breaks is critical for maintaining genome integrity and avoiding disease. Modification of damaged chromatin has profound consequences for the initial signaling and regulation of repair. One such modification involves ubiquitination by E3 ligases RNF8 and RNF168 within minutes after DNA double-strand break formation, altering chromatin structure and recruiting factors such as 53BP1 and BRCA1 for repair via non-homologous end-joining (NHEJ) and homologous recombination (HR), respectively. The WD40 protein WRAP53β plays an essential role in localizing RNF8 to DNA breaks by scaffolding its interaction with the upstream factor MDC1. Loss of WRAP53β impairs ubiquitination at DNA lesions and reduces downstream repair by both NHEJ and HR. Intriguingly, WRAP53β depletion attenuates repair of DNA double-strand breaks more than depletion of RNF8, indicating functions other than RNF8-mediated ubiquitination. WRAP53β plays key roles with respect to the nuclear organelles Cajal bodies, including organizing the genome to promote associated transcription and collecting factors involved in maturation of the spliceosome and telomere elongation within these organelles. It is possible that similar functions may aid also in DNA repair. Here we describe the involvement of WRAP53β in Cajal bodies and DNA double-strand break repair in detail and explore whether and how these processes may be linked. We also discuss the possibility that the overexpression of WRAP53β detected in several cancer types may reflect its normal participation in the DNA damage response rather than oncogenic properties.