Ionizing radiation damage and repair from 3D-genomic perspective.

Ionizing radiation (IR)-induced DNA damage and repair are complex and occur at hierarchical chromatin structures; radiobiology needs to be studied from a 3D-genomic perspective. Differences in IR damage and repair throughout the 3D genome may help to explain differences in radiosensitivity.

A transcriptome-based single-cell biological age model and resource for tissue-specific aging measures.

Accurately measuring biological age is crucial for improving healthcare for the elderly population. However, the complexity of aging biology poses challenges in how to robustly estimate aging and interpret the biological significance of the traits used for estimation. Here we present SCALE, a statistical pipeline that quantifies biological aging in different tissues using explainable features learned from literature and single-cell transcriptomic data. Applying SCALE to the "Mouse Aging Cell Atlas" (Tabula Muris Senis) data, we identified tissue-level transcriptomic aging programs for more than 20 murine tissues and created a multitissue resource of mouse quantitative aging-associated genes. We observe that SCALE correlates well with other age indicators, such as the accumulation of somatic mutations, and can distinguish subtle differences in aging even in cells of the same chronological age. We further compared SCALE with other transcriptomic and methylation "clocks" in data from aging muscle stem cells, Alzheimer's disease, and heterochronic parabiosis. Our results confirm that SCALE is more generalizable and reliable in assessing biological aging in aging-related diseases and rejuvenating interventions. Overall, SCALE represents a valuable advancement in our ability to measure aging accurately, robustly, and interpretably in single cells.

Stratifying TAD boundaries pinpoints focal genomic regions of regulation, damage, and repair.

Advances in chromatin mapping have exposed the complex chromatin hierarchical organization in mammals, including topologically associating domains (TADs) and their substructures, yet the functional implications of this hierarchy in gene regulation and disease progression are not fully elucidated. Our study delves into the phenomenon of shared TAD boundaries, which are pivotal in maintaining the hierarchical chromatin structure and regulating gene activity. By integrating high-resolution Hi-C data, chromatin accessibility, and DNA double-strand breaks (DSBs) data from various cell lines, we systematically explore the complex regulatory landscape at high-level TAD boundaries. Our findings indicate that these boundaries are not only key architectural elements but also vibrant hubs, enriched with functionally crucial genes and complex transcription factor binding site-clustered regions. Moreover, they exhibit a pronounced enrichment of DSBs, suggesting a nuanced interplay between transcriptional regulation and genomic stability. Our research provides novel insights into the intricate relationship between the 3D genome structure, gene regulation, and DNA repair mechanisms, highlighting the role of shared TAD boundaries in maintaining genomic integrity and resilience against perturbations. The implications of our findings extend to understanding the complexities of genomic diseases and open new avenues for therapeutic interventions targeting the structural and functional integrity of TAD boundaries.

The developmental and evolutionary characteristics of transcription factor binding site clustered regions based on an explainable machine learning model.

Gene expression is temporally and spatially regulated by the interaction of transcription factors (TFs) and cis-regulatory elements (CREs). The uneven distribution of TF binding sites across the genome poses challenges in understanding how this distribution evolves to regulate spatio-temporal gene expression and consequent heritable phenotypic variation. In this study, chromatin accessibility profiles and gene expression profiles were collected from several species including mammals (human, mouse, bovine), fish (zebrafish and medaka), and chicken. Transcription factor binding sites clustered regions (TFCRs) at different embryonic stages were characterized to investigate regulatory evolution. The study revealed dynamic changes in TFCR distribution during embryonic development and species evolution. The synchronization between TFCR complexity and gene expression was assessed across species using RegulatoryScore. Additionally, an explainable machine learning model highlighted the importance of the distance between TFCR and promoter in the coordinated regulation of TFCRs on gene expression. Our results revealed the developmental and evolutionary dynamics of TFCRs during embryonic development from fish, chicken to mammals. These data provide valuable resources for exploring the relationship between transcriptional regulation and phenotypic differences during embryonic development.

dbEmbryo multi-omics database for analyses of synergistic regulation in early mammalian embryo development.

During early mammalian embryo development, different epigenetic marks undergo reprogramming and play crucial roles in the mediation of gene expression. Currently, several databases provide multi-omics information on early embryos. However, how interconnected epigenetic markers function together to coordinate the expression of the genetic code in a spatiotemporal manner remains difficult to analyze, markedly limiting scientific and clinical research. Here, we present dbEmbryo, an integrated and interactive multi-omics database for human and mouse early embryos. dbEmbryo integrates data on gene expression, DNA methylation, histone modifications, chromatin accessibility, and higher-order chromatin structure profiles for human and mouse early embryos. It incorporates customized analysis tools, such as "multi-omics visualization," "Gene&Peak annotation," "ZGA gene cluster," "cis-regulation," "synergistic regulation," "promoter signal enrichment," and "3D genome." Users can retrieve gene expression and epigenetic profile patterns to analyze synergistic changes across different early embryo developmental stages. We showed the uniqueness of dbEmbryo among extant databases containing data on early embryo development and provided an overview. Using dbEmbryo, we obtained a phase-separated model of transcriptional control during early embryo development. dbEmbryo offers web-based analytical tools and a comprehensive resource for biologists and clinicians to decipher molecular regulatory mechanisms of human and mouse early embryo development.

Spatial density of open chromatin: an effective metric for the functional characterization of topologically associated domains.

Topologically associated domains (TADs) are spatial and functional units of metazoan chromatin structure. Interpretation of the interplay between regulatory factors and chromatin structure within TADs is crucial to understand the spatial and temporal regulation of gene expression. However, a computational metric for the sensitive characterization of TAD regulatory landscape is lacking. Here, we present the spatial density of open chromatin (SDOC) metric as a quantitative measurement of intra-TAD chromatin state and structure. SDOC sensitively reflects epigenetic properties and gene transcriptional activity in TADs. During mouse T-cell development, we found that TADs with decreased SDOC are enriched in repressed developmental genes, and the joint effect of SDOC-decreasing and TAD clustering corresponds to the highest level of gene repression. In addition, we revealed a pervasive preference for TADs with similar SDOC to interact with each other, which may reflect the principle of chromatin organization.

Computational methods for the prediction of chromatin interaction and organization using sequence and epigenomic profiles.

The exploration of three-dimensional chromatin interaction and organization provides insight into mechanisms underlying gene regulation, cell differentiation and disease development. Advances in chromosome conformation capture technologies, such as high-throughput chromosome conformation capture (Hi-C) and chromatin interaction analysis by paired-end tag (ChIA-PET), have enabled the exploration of chromatin interaction and organization. However, high-resolution Hi-C and ChIA-PET data are only available for a limited number of cell lines, and their acquisition is costly, time consuming, laborious and affected by theoretical limitations. Increasing evidence shows that DNA sequence and epigenomic features are informative predictors of regulatory interaction and chromatin architecture. Based on these features, numerous computational methods have been developed for the prediction of chromatin interaction and organization, whereas they are not extensively applied in biomedical study. A systematical study to summarize and evaluate such methods is still needed to facilitate their application. Here, we summarize 48 computational methods for the prediction of chromatin interaction and organization using sequence and epigenomic profiles, categorize them and compare their performance. Besides, we provide a comprehensive guideline for the selection of suitable methods to predict chromatin interaction and organization based on available data and biological question of interest.

Prognosis prediction of icotinib as targeted therapy for advanced EGFR-positive non-small cell lung cancer patients.

Clinical trials on icotinib, a first-generation epidermal growth factor receptor tyrosine kinase inhibitor (EGFR-TKI), have shown promising results as targeted therapy for non-small cell lung cancer (NSCLC). This study aimed to establish an effective scoring system to predict the one-year progression-free survival (PFS) of advanced NSCLC patients with EGFR mutations treated with icotinib as targeted therapy. A total of 208 consecutive patients with advanced EGFR-positive NSCLC treated with icotinib were enrolled in this study. Baseline characteristics were collected within 30 days before icotinib treatment. PFS was taken as the primary endpoint and the response rate as the secondary endpoint. Least absolute shrinkage and selection operator (LASSO) regression analysis and Cox proportional hazards regression analysis were used to select the optimal predictors. We evaluated the scoring system using a five-fold cross-validation. PFS events occurred in 175 patients, with a median PFS of 9.9 months (interquartile range, 6.8-14.5). The objective response rate (ORR) was 36.1%, and the disease control rate (DCR) was 67.3%. The final ABC-Score consisted of three predictors: age, bone metastases and carbohydrate antigen 19-9 (CA19-9). Upon comparison of all three factors, the combined ABC-score (area under the curve (AUC)= 0.660) showed a better predictive accuracy than age (AUC = 0.573), bone metastases (AUC = 0.615), and CA19-9 (AUC = 0.608) individually. A five-fold cross-validation showed good discrimination with AUC = 0.623. The ABC-score developed in this study was significantly effective as a prognostic tool for icotinib in advanced NSCLC patients with EGFR mutations.

New insights on human essential genes based on integrated analysis and the construction of the HEGIAP web-based platform.

Essential genes are those whose loss of function compromises organism viability or results in profound loss of fitness. Recent gene-editing technologies have provided new opportunities to characterize essential genes. Here, we present an integrated analysis that comprehensively and systematically elucidates the genetic and regulatory characteristics of human essential genes. First, we found that essential genes act as 'hubs' in protein-protein interaction networks, chromatin structure and epigenetic modification. Second, essential genes represent conserved biological processes across species, although gene essentiality changes differently among species. Third, essential genes are important for cell development due to their discriminate transcription activity in embryo development and oncogenesis. In addition, we developed an interactive web server, the Human Essential Genes Interactive Analysis Platform (http://sysomics.com/HEGIAP/), which integrates abundant analytical tools to enable global, multidimensional interpretation of gene essentiality. Our study provides new insights that improve the understanding of human essential genes.

A survey and evaluation of Web-based tools/databases for variant analysis of TCGA data.

The Cancer Genome Atlas (TCGA) is a publicly funded project that aims to catalog and discover major cancer-causing genomic alterations with the goal of creating a comprehensive 'atlas' of cancer genomic profiles. The availability of this genome-wide information provides an unprecedented opportunity to expand our knowledge of tumourigenesis. Computational analytics and mining are frequently used as effective tools for exploring this byzantine series of biological and biomedical data. However, some of the more advanced computational tools are often difficult to understand or use, thereby limiting their application by scientists who do not have a strong computational background. Hence, it is of great importance to build user-friendly interfaces that allow both computational scientists and life scientists without a computational background to gain greater biological and medical insights. To that end, this survey was designed to systematically present available Web-based tools and facilitate the use TCGA data for cancer research.

DeepHiC: A generative adversarial network for enhancing Hi-C data resolution.

Hi-C is commonly used to study three-dimensional genome organization. However, due to the high sequencing cost and technical constraints, the resolution of most Hi-C datasets is coarse, resulting in a loss of information and biological interpretability. Here we develop DeepHiC, a generative adversarial network, to predict high-resolution Hi-C contact maps from low-coverage sequencing data. We demonstrated that DeepHiC is capable of reproducing high-resolution Hi-C data from as few as 1% downsampled reads. Empowered by adversarial training, our method can restore fine-grained details similar to those in high-resolution Hi-C matrices, boosting accuracy in chromatin loops identification and TADs detection, and outperforms the state-of-the-art methods in accuracy of prediction. Finally, application of DeepHiC to Hi-C data on mouse embryonic development can facilitate chromatin loop detection. We develop a web-based tool (DeepHiC, http://sysomics.com/deephic) that allows researchers to enhance their own Hi-C data with just a few clicks.

Stable H3K4me3 is associated with transcription initiation during early embryo development.

MOTIVATION: During development of the mammalian embryo, histone modification H3K4me3 plays an important role in regulating gene expression and exhibits extensive reprograming on the parental genomes. In addition to these dramatic epigenetic changes, certain unchanging regulatory elements are also essential for embryonic development. RESULTS: Using large-scale H3K4me3 chromatin immunoprecipitation sequencing data, we identified a form of H3K4me3 that was present during all eight stages of the mouse embryo before implantation. This 'stable H3K4me3' was highly accessible and much longer than normal H3K4me3. Moreover, most of the stable H3K4me3 was in the promoter region and was enriched in higher chromatin architecture. Using in-depth analysis, we demonstrated that stable H3K4me3 was related to higher gene expression levels and transcriptional initiation during embryonic development. Furthermore, stable H3K4me3 was much more active in blood tumor cells than in normal blood cells, suggesting a potential mechanism of cancer progression. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Exploration of prognosis-related microRNA and transcription factor co-regulatory networks across cancer types.

The study of cancer prognosis serves as an important part of cancer research. Large-scale cancer studies have identified numerous genes and microRNAs (miRNAs) associated with prognosis. These informative genes and miRNAs represent potential biomarkers to predict survival and to elucidate the molecular mechanism of tumour progression. MiRNAs and transcription factors (TFs) can work cooperatively as essential mediators of gene expression, and their dysregulation affects cancer prognosis. A panoramic view of cancer prognosis at the system level, considering the co-regulation roles of miRNA and TF, remains elusive. Here, we establish 12 prognosis-related miRNA-TF co-regulatory networks. The characteristics of prognostic target genes and their regulators in the network are depicted. Although the target genes and co-regulatory patterns exhibit cancer-specific properties, some miRNAs and TFs are highly conserved across cancers. We illustrate and interpret the roles of these conserved regulators by building a model associated with cancer hallmarks, functional enrichment analysis, network community detection, and exhaustive literature research. The elaborated system-level prognostic miRNA-TF co-regulation landscape, including the highlighted roles of conserved regulators, provides a novel and powerful insights into further biological and medical discoveries.

3DSNP: a database for linking human noncoding SNPs to their three-dimensional interacting genes.

The vast noncoding portion of the human genome harbors a rich array of functional elements and disease-causing regulatory variants. Recent high-throughput chromosome conformation capture studies have outlined the principles of these elements interacting and regulating the expression of distal target genes through three-dimensional (3D) chromatin looping. Here we present 3DSNP, an integrated database for annotating human noncoding variants by exploring their roles in the distal interactions between genes and regulatory elements. 3DSNP integrates 3D chromatin interactions, local chromatin signatures in different cell types and linkage disequilibrium (LD) information from the 1000 Genomes Project. 3DSNP provides informative visualization tools to display the integrated local and 3D chromatin signatures and the genetic associations among variants. Data from different functional categories are integrated in a scoring system that quantitatively measures the functionality of SNPs to help select important variants from a large pool. 3DSNP is a valuable resource for the annotation of human noncoding genome sequence and investigating the impact of noncoding variants on clinical phenotypes. The 3DSNP database is available at http://biotech.bmi.ac.cn/3dsnp/.

Exploring spatially adjacent TFBS-clustered regions with Hi-C data.

MOTIVATION: Transcription factor binding sites (TFBSs) are clustered in the human genome, forming the TFBS-clustered regions that regulate gene transcription, which requires dynamic chromatin configurations between promoters and distal regulatory elements. Here, we propose a regulatory model called spatially adjacent TFBS-clustered regions (SATs), in which TFBS-clustered regions are connected by spatial proximity as identified by high-resolution Hi-C data. RESULTS: TFBS-clustered regions forming SATs appeared less frequently in gene promoters than did isolated TFBS-clustered regions, whereas SATs as a whole appeared more frequently. These observations indicate that multiple distal TFBS-clustered regions combined to form SATs to regulate genes. Further examination confirmed that a substantial portion of genes regulated by SATs were located between the paired TFBS-clustered regions instead of the downstream. We reconstructed the chromosomal conformation of the H1 human embryonic stem cell line using the ShRec3D algorithm and proposed the SAT regulatory model. CONTACT: ylu.phd@gmail.com or boxc@bmi.ac.cn. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

[Investigation of Evolution Rules of Phlegm and Blood Stasis Syndrome in Hyperlipidemia and Ath- erosclerosis by NMR-based Metabolic Profiling and Metabonomic Approaches].

OBJECTIVE: To explore evolution rules of phlegm and blood stasis syndrome ( PBSS) in hyperlipidemia and atherosclerosis (AS) using NMR-based metabolic profiling and metabonomic approaches based on formulas corresponding to syndrome. METHODS: Totally 150 SD rats were divided into the normal group, the model group, the Erchen Decoction (ED) group, the Xuefu Zhuyu Decoction (XZD) group, the Lipitor group, 30 in each group. The hyperlipidemia and AS rat model was duplicated by suturing carotid artery, injecting vitamin D3, and feeding with high fat diet. ED and XZD were used as drug probes. Blood samples were withdrawn at week 2, 4, and 8 after modeling. Blood lipids, blood rheology, histopathology and metabolomics were detected and analyzed. Results Results of blood lipids and pathology showed hyperlipidemia and early AS rat models were successfully established. At week 2 after modeling, levels of total cholesterol (TC) and low-density lipoprotein cholesterol (LDL-C) significantly increased, which reached the peak at week 4 and maintained at higher levels at week 8. ED exerted obvious effect in improving TC and LDL-C levels of early models, while XZD could greatly improve levels of TC and LDL-C of late models. Rheological results showed at week 2, there was no significant difference in whole blood viscosity, plasma viscosity, or hematocrit between the model group and the normal group (P > 0.05). At week 4 partial hemorheological indicators (such as plasma viscosity) were abnormal. Till week 8 whole blood viscosity, plasma viscosity, and hematocrit were significantly abnormal (P <0. 05, P < 0.01). As time went by, whole blood viscosity, plasma viscosity, and hematocrit showed gradual increasing tendency in the ED group, while they showed gradual decreasing tendency in the XZD group. Results of metabonomics showed significant difference in spectra of metabolites between the normal group and the model group. As modeling time was prolonged, contents of acetyl glucoprotein and glucose in the model group increased in late stage, which was in. line with results of blood lipids and hemorheology. ED showed more obvious effect in early and mid-term modeling (at week 2 and 4), and increased contents of partial metabolites (such as choline, phosphatidyl choline, glycerophosphocholine), but these changes in the XZD group were consistent with those of the model group. In late modeling (at week 8) XZD showed more obvious effect in improving contents of lactic acid, acetyl glycoprotein, LDL, creatine, choline, and glucose. CONCLUSIONS: ED and XZD not only showed regulatory effects on lipid disorders, but also could improve dysbolism of Chos. In formulas corresponding to syndrome, damp-phlegm was main pathogenesis of hyperlipidema and AS in early and mid stages. Blood stasis syndrome began to occur along with it progressed. Phlegm can result in blood stasis and intermingles with stasis. Phlegm turbidity runs through the whole process.

Dual-Layer Spectral Detector Computed Tomography Quantitative Parameters: A Potential Tool for Lymph Node Activity Determination in Lymphoma Patients.

Dual-energy CT has shown promising results in determining tumor characteristics and treatment effectiveness through spectral data by assessing normalized iodine concentration (nIC), normalized effective atomic number (nZeff), normalized electron density (nED), and extracellular volume (ECV). This study explores the value of quantitative parameters in contrast-enhanced dual-layer spectral detector CT (SDCT) as a potential tool for detecting lymph node activity in lymphoma patients. A retrospective analysis of 55 lymphoma patients with 289 lymph nodes, assessed through 18FDG-PET/CT and the Deauville five-point scale, revealed significantly higher values of nIC, nZeff, nED, and ECV in active lymph nodes compared to inactive ones (p < 0.001). Generalized linear mixed models showed statistically significant fixed-effect parameters for nIC, nZeff, and ECV (p < 0.05). The area under the receiver operating characteristic curve (AUROC) values of nIC, nZeff, and ECV reached 0.822, 0.845, and 0.811 for diagnosing lymph node activity. In conclusion, the use of g nIC, nZeff, and ECV as alternative imaging biomarkers to PET/CT for identifying lymph node activity in lymphoma holds potential as a reliable diagnostic tool that can guide treatment decisions.

Comparative study on genomic and epigenomic profiles of retinoblastoma or tuberous sclerosis complex via nanopore sequencing and a joint screening framework.

Recurrence and extraocular metastasis in advanced intraocular retinoblastoma (RB) are still major obstacles for successful treatment of Chinese children. Tuberous sclerosis complex (TSC) is a very rare, multisystemic genetic disorder characterized by hamartomatous growth. In this study, we aimed to compare genomic and epigenomic profiles with human RB or TSC using recently developed nanopore sequencing, and to identify disease-associated variations or genes. Peripheral blood samples were collected from either RB or RB/TSC patients plus their normal siblings, followed by nanopore sequencing and identification of disease-specific structural variations (SVs) and differentially methylated regions (DMRs) by a systematic biology strategy named as multiomics-based joint screening framework. In total, 316 RB- and 1295 TSC-unique SVs were identified, as well as 1072 RB- and 1114 TSC-associated DMRs, respectively. We eventually identified 6 key genes for RB for further functional validation. Knockdown of CDK19 with specific siRNAs significantly inhibited Y79 cellular proliferation and increased sensitivity to carboplatin, whereas downregulation of AHNAK2 promoted the cell growth as well as drug resistance. Those two genes might serve as potential diagnostic markers or therapeutic targets of RB. The systematic biology strategy combined with functional validation might be an effective approach for rare pediatric malignances with limited samples and challenging collection process.

Ablation alone is noninferior to radiotherapy plus ablation in the patients with early-stage hepatocellular carcinoma: a population-based study.

Recently, the efficacy of two low-invasive treatments, ablation, and radiotherapy, has been fully compared for the patients with the early-stage hepatocellular carcinoma (HCC). However, the comparison between radiotherapy plus ablation and ablation alone has been less frequently reported. Data from the Surveillance, Epidemiology, and End Results (SEER) database were searched for early-stage HCC patients treated with ablation plus radiotherapy or ablation alone. The outcome measures were overall survival (OS) and cancer-specific survival (CSS). The propensity score matching (PSM) was used to reduce selection bias. We included 240 and 6619 patients in the radiotherapy plus ablation group and ablation group before the PSM. After PSM, 240 pairs of patients were included. The median OS (mOS) and median CSS (mCSS) of patients receiving ablation alone were longer than that of receiving radiotherapy plus ablation (mOS: 47 vs. 34 months, P = 0.019; mCSS: 77 vs. 40 months, P = 0.018, after PSM) before and after PSM. The multivariate analysis indicated that radiotherapy plus ablation independent risk factor for OS and CSS before PSM, but the significance disappeared after PSM. The detailed subgroup analyses indicated ablation alone brought more benefit in very early-stage HCC and older patients. In addition, we found different types of radiotherapy might lead to different outcomes when combined with ablation. In conclusion, ablation alone is noninferior to radiotherapy plus ablation in patients with early-stage HCC. The additional radiation prior to ablation may bring survival benefits in the patients with higher tumor stage. However, due to the risk of selection bias in that study, the results should be interpreted cautiously.

SEPDB: a database of secreted proteins.

Detecting changes in the dynamics of secreted proteins in serum has been a challenge for proteomics. Enter secreted protein database (SEPDB), an integrated secretory proteomics database offering human, mouse and rat secretory proteomics datasets collected from serum, exosomes and cell culture media. SEPDB compiles secreted protein information from secreted protein database, UniProt and Human Protein Atlas databases to annotate secreted proteomics data based on protein subcellular localization and disease markers. SEPDB integrates the latest predictive modeling techniques to measure deviations in the distribution of signal peptide structures of secreted proteins, extends signal peptide sequence prediction by excluding transmembrane structural domain proteins and updates the validation analysis pipeline for secreted proteins. To establish tissue-specific profiles, we have also created secreted proteomics datasets associated with different human tissues. In addition, we provide information on heterogeneous receptor network organizational relationships, reflective of the complex functional information inherent in the molecular structures of secreted proteins that serve as ligands. Users can take advantage of the Refreshed Search, Analyze, Browse and Download functions of SEPDB, which is available online at https://sysomics.com/SEPDB/. Database URL:  https://sysomics.com/SEPDB/.