Acetylation within the globular core domain of histone H3 on lysine 56 (H3K56) has recently been shown to have a critical role in packaging DNA into chromatin following DNA replication and repair in budding yeast. However, the function or occurrence of this specific histone mark has not been studied in multicellular eukaryotes, mainly because the Rtt109 enzyme that is known to mediate acetylation of H3K56 (H3K56ac) is fungal-specific. Here we demonstrate that the histone acetyl transferase CBP (also known as Nejire) in flies and CBP and p300 (Ep300) in humans acetylate H3K56, whereas Drosophila Sir2 and human SIRT1 and SIRT2 deacetylate H3K56ac. The histone chaperones ASF1A in humans and Asf1 in Drosophila are required for acetylation of H3K56 in vivo, whereas the histone chaperone CAF-1 (chromatin assembly factor 1) in humans and Caf1 in Drosophila are required for the incorporation of histones bearing this mark into chromatin. We show that, in response to DNA damage, histones bearing acetylated K56 are assembled into chromatin in Drosophila and human cells, forming foci that colocalize with sites of DNA repair. Furthermore, acetylation of H3K56 is increased in multiple types of cancer, correlating with increased levels of ASF1A in these tumours. Our identification of multiple proteins regulating the levels of H3K56 acetylation in metazoans will allow future studies of this critical and unique histone modification that couples chromatin assembly to DNA synthesis, cell proliferation and cancer.

The histone chaperone Chromatin Assembly Factor 1 (CAF-1) deposits tetrameric (H3/H4)2 histones onto newly-synthesized DNA during DNA replication. To understand the mechanism of the tri-subunit CAF-1 complex in this process, we investigated the protein-protein interactions within the CAF-1-H3/H4 architecture using biophysical and biochemical approaches. Hydrogen/deuterium exchange and chemical cross-linking coupled to mass spectrometry reveal interactions that are essential for CAF-1 function in budding yeast, and importantly indicate that the Cac1 subunit functions as a scaffold within the CAF-1-H3/H4 complex. Cac1 alone not only binds H3/H4 with high affinity, but also promotes histone tetramerization independent of the other subunits. Moreover, we identify a minimal region in the C-terminus of Cac1, including the structured winged helix domain and glutamate/aspartate-rich domain, which is sufficient to induce (H3/H4)2 tetramerization. These findings reveal a key role of Cac1 in histone tetramerization, providing a new model for CAF-1-H3/H4 architecture and function during eukaryotic replication.

The hypoxic tumor microenvironment serves as a niche for maintaining the glioma-initiating cells (GICs) that are critical for glioblastoma (GBM) occurrence and recurrence. Here, we report that hypoxia-induced miR-215 is vital for reprograming GICs to fit the hypoxic microenvironment via suppressing the expression of an epigenetic regulator KDM1B and modulating activities of multiple pathways. Interestingly, biogenesis of miR-215 and several miRNAs is accelerated post-transcriptionally by hypoxia-inducible factors (HIFs) through HIF-Drosha interaction. Moreover, miR-215 expression correlates inversely with KDM1B while correlating positively with HIF1α and GBM progression in patients. These findings reveal a direct role of HIF in regulating miRNA biogenesis and consequently activating the miR-215-KDM1B-mediated signaling required for GIC adaptation to hypoxia.

Adverse drug reactions (ADRs) are responsible for drug candidate failure during clinical trials. It is crucial to investigate biological pathways contributing to ADRs. Here, we applied a large-scale analysis to identify overrepresented ADR-pathway combinations through merging clinical phenotypic data, biological pathway data, and drug-target relations. Evaluation was performed by scientific literature review and defining a pathway-based ADR-ADR similarity measure. The results showed that our method is efficient for finding the associations between ADRs and pathways. To more systematically understand the mechanisms of ADRs, we constructed an ADR-pathway network and an ADR-ADR network. Through network analysis on biology and pharmacology, it was found that frequent ADRs were associated with more pathways than infrequent and rare ADRs. Moreover, environmental information processing pathways contributed most to the observed ADRs. Integrating the system organ class of ADRs, we found that most classes tended to interact with other classes instead of themselves. ADR classes were distributed promiscuously in all the ADR cliques. These results reflected that drug perturbation to a certain pathway can cause changes in multiple organs, rather than in one specific organ. Our work not only provides a global view of the associations between ADRs and pathways, but also is helpful to understand the mechanisms of ADRs.

The protein anti-silencing function 1 (Asf1) chaperones histones H3/H4 for assembly into nucleosomes every cell cycle as well as during DNA transcription and repair. Asf1 interacts directly with H4 through the C-terminal tail of H4, which itself interacts with the docking domain of H2A in the nucleosome. The structure of this region of the H4 C-terminus differs greatly in these two contexts.

The modulator of retrovirus infection (MRI or CYREN) is a 30-kDa protein with a conserved N-terminal Ku-binding motif (KBM) and a C-terminal XLF-like motif (XLM). We show that MRI is intrinsically disordered and interacts with many DNA damage response (DDR) proteins, including the kinases ataxia telangiectasia mutated (ATM) and DNA-PKcs and the classical non-homologous end joining (cNHEJ) factors Ku70, Ku80, XRCC4, XLF, PAXX, and XRCC4. MRI forms large multimeric complexes that depend on its N and C termini and localizes to DNA double-strand breaks (DSBs), where it promotes the retention of DDR factors. Mice deficient in MRI and XLF exhibit embryonic lethality at a stage similar to those deficient in the core cNHEJ factors XRCC4 or DNA ligase IV. Moreover, MRI is required for cNHEJ-mediated DSB repair in XLF-deficient lymphocytes. We propose that MRI is an adaptor that, through multivalent interactions, increases the avidity of DDR factors to DSB-associated chromatin to promote cNHEJ.

Here, we use ChAP-MS (chromatin affinity purification with mass spectrometry), for the affinity purification of a sequence-specific single-copy endogenous chromosomal locus containing a DNA double-strand break (DSB). We found multiple new histone post-translational modifications enriched on chromatin bearing a DSB from budding yeast. One of these, methylation of histone H3 on lysine 125, has not previously been reported. Among over 100 novel proteins enriched at a DSB were the phosphatase Sit4, the RNA pol II degradation factor Def1, the mRNA export protein Yra1 and the HECT E3 ligase Tom1. Each of these proteins was required for resistance to radiomimetics, and many were required for resistance to heat, which we show here to cause a defect in DSB repair in yeast. Yra1 and Def1 were required for DSB repair per se, while Sit4 was required for rapid inactivation of the DNA damage checkpoint after DSB repair. Thus, our unbiased proteomics approach has led to the unexpected discovery of novel roles for these and other proteins in the DNA damage response.

More than 150 ginsenosides have been isolated and identified from Panax plants. Ginsenosides with different glycosylation degrees have demonstrated different chemical properties and bioactivity. In this study, we systematically cloned and characterized 46 UGT94 family UDP-glycosyltransferases (UGT94s) from a mixed Panax ginseng/callus cDNA sample with high amino acid identity. These UGT94s were found to catalyze sugar chain elongation at C3-O-Glc and/or C20-O-Glc of protopanaxadiol (PPD)-type, C20-O-Glc or C6-O-Glc of protopanaxatriol (PPT)-type or both C3-O-Glc of PPD-type and C6-O-Glc of PPT-type or C20-O-Glc of PPD-type and PPT-type ginsenosides with different efficiencies. We also cloned 26 and 51 UGT94s from individual P. ginseng and P. notoginseng plants, respectively; our characterization results suggest that there is a group of UGT94s with high amino acid identity but diverse functions or catalyzing activities even within individual plants. These UGT94s were classified into three clades of the phylogenetic tree and consistent with their catalytic function. Based on these UGT94s, we elucidated the biosynthetic pathway of a group of ginsenosides. Our present results reveal a series of UGTs involved in second sugar chain elongation of saponins in Panax plants, and provide a scientific basis for understanding the diverse evolution mechanisms of UGT94s among plants.

Although our knowledge of human diseases has increased dramatically, the molecular basis, phenotypic traits, and therapeutic targets of most diseases still remain unclear. An increasing number of studies have observed that similar diseases often are caused by similar molecules, can be diagnosed by similar markers or phenotypes, or can be cured by similar drugs. Thus, the identification of diseases similar to known ones has attracted considerable attention worldwide. To this end, the associations between diseases at the molecular, phenotypic, and taxonomic levels were used to measure the pairwise similarity in diseases. The corresponding performance assessment strategies for these methods involving the terms "category-based," "simulated-patient-based," and "benchmark-data-based" were thus further emphasized. Then, frequently used methods were evaluated using a benchmark-data-based strategy. To facilitate the assessment of disease similarity scores, researchers have designed dozens of tools that implement these methods for calculating disease similarity. Currently, disease similarity has been advantageous in predicting noncoding RNA (ncRNA) function and therapeutic drugs for diseases. In this article, we review disease similarity methods, evaluation strategies, tools, and their applications in the biomedical community. We further evaluate the performance of these methods and discuss the current limitations and future trends for calculating disease similarity.

Osteoporosis is the most common bone metabolic disease, characterized by bone mass loss and bone microstructure changes due to unbalanced bone conversion, which increases bone fragility and fracture risk. Glucocorticoids are clinically used to treat a variety of diseases, including inflammation, cancer and autoimmune diseases. However, excess glucocorticoids can cause osteoporosis. Herein we performed an integrated analysis of two glucocorticoid-related microarray datasets. The WGCNA analysis identified 3 and 4 glucocorticoid-related gene modules, respectively. Differential expression analysis revealed 1047 and 844 differentially expressed genes in the two datasets. After integrating differentially expressed glucocorticoid-related genes, we found that most of the robust differentially expressed genes were up-regulated. Through protein-protein interaction analysis, we obtained 158 glucocorticoid-related candidate genes. Enrichment analysis showed that these genes are significantly enriched in the osteoporosis related pathways. Our results provided new insights into glucocorticoid-induced osteoporosis and potential candidate markers of osteoporosis.

Osteoporosis is a skeletal disorder characterized by a systemic impairment of bone mineral density (BMD). Genome-wide association studies (GWAS) have identified hundreds of susceptibility loci for osteoporosis and BMD. However, the vast majority of susceptibility loci are located in non-coding regions of the genome and provide limited information about the genetic mechanisms of osteoporosis. Herein we performed a comprehensive functional analysis to investigate the genetic and epigenetic mechanisms of osteoporosis and BMD. BMD and osteoporosis are found to share many common susceptibility loci, and the corresponding susceptibility genes are significantly enriched in bone-related biological pathways. The regulatory element enrichment analysis indicated that BMD and osteoporosis susceptibility loci are significantly enriched in 5'UTR and DNase I hypersensitive sites (DHSs) of peripheral blood immune cells. By integrating GWAS and expression Quantitative Trait Locus (eQTL) data, we found that 15 protein-coding genes are regulated by the osteoporosis and BMD susceptibility loci. Our analysis provides new clues for a better understanding of the pathogenic mechanisms and offers potential therapeutic targets for osteoporosis.

To explore the relationship between ultrasound (US) features and Ki-67 labeling index (LI) of soft tissue sarcoma (STS).

A dysfunctional immune response in coronavirus disease 2019 (COVID-19) patients is a recurrent theme impacting symptoms and mortality, yet a detailed understanding of pertinent immune cells is not complete. We applied single-cell RNA sequencing to 284 samples from 196 COVID-19 patients and controls and created a comprehensive immune landscape with 1.46 million cells. The large dataset enabled us to identify that different peripheral immune subtype changes are associated with distinct clinical features, including age, sex, severity, and disease stages of COVID-19. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA was found in diverse epithelial and immune cell types, accompanied by dramatic transcriptomic changes within virus-positive cells. Systemic upregulation of S100A8/A9, mainly by megakaryocytes and monocytes in the peripheral blood, may contribute to the cytokine storms frequently observed in severe patients. Our data provide a rich resource for understanding the pathogenesis of and developing effective therapeutic strategies for COVID-19.

Single-cell RNA sequencing (scRNA-seq) detects whole transcriptome signals for large amounts of individual cells and is powerful for determining cell-to-cell differences and investigating the functional characteristics of various cell types. scRNA-seq datasets are usually sparse and highly noisy. Many steps in the scRNA-seq analysis workflow, including reasonable gene selection, cell clustering and annotation, as well as discovering the underlying biological mechanisms from such datasets, are difficult. In this study, we proposed an scRNA-seq analysis method based on the latent Dirichlet allocation (LDA) model. The LDA model estimates a series of latent variables, i.e. putative functions (PFs), from the input raw cell-gene data. Thus, we incorporated the 'cell-function-gene' three-layer framework into scRNA-seq analysis, as this framework is capable of discovering latent and complex gene expression patterns via a built-in model approach and obtaining biologically meaningful results through a data-driven functional interpretation process. We compared our method with four classic methods on seven benchmark scRNA-seq datasets. The LDA-based method performed best in the cell clustering test in terms of both accuracy and purity. By analysing three complex public datasets, we demonstrated that our method could distinguish cell types with multiple levels of functional specialization, and precisely reconstruct cell development trajectories. Moreover, the LDA-based method accurately identified the representative PFs and the representative genes for the cell types/cell stages, enabling data-driven cell cluster annotation and functional interpretation. According to the literature, most of the previously reported marker/functionally relevant genes were recognized.

Gain-of-function mutations in the catalytic site of EZH2 (Enhancer of Zeste Homologue 2), is observed in about 22% of diffuse large B-cell lymphoma (DLBCL) cases. Here we show that selective inhibition of histone deacetylase 1,2 (HDAC1,2) activity using a small molecule inhibitor causes cytotoxic or cytostatic effects in EZH2 gain-of-function mutant (EZH2GOF) DLBCL cells. Our results show that blocking the activity of HDAC1,2 increases global H3K27ac without causing a concomitant global decrease in H3K27me3 levels. Our data shows that inhibition of HDAC1,2 is sufficient to decrease H3K27me3 present at DSBs, decrease DSB repair and activate the DNA damage response in these cells. In addition to increased H3K27me3, we found that the EZH2GOF DLBCL cells overexpress another chemotherapy resistance factor - B-lymphoma and BAL-associated protein (BBAP). BBAP monoubiquitinates histone H4K91, a residue that is also subjected to acetylation. Our results show that selective inhibition of HDAC1,2 increases H4K91ac, decreases BBAP-mediated H4K91 monoubiquitination, impairs BBAP-dependent DSB repair and sensitizes the refractory EZH2GOF DLBCL cells to treatment with doxorubicin, a chemotherapy agent. Hence, selective HDAC1,2 inhibition provides a novel DNA repair mechanism-based therapeutic approach as it can overcome both EZH2- and BBAP-mediated DSB repair in the EZH2GOF DLBCL cells.

Much of our understanding of the function of histone post-translational modifications in metazoans is inferred from their genomic localization and / or extrapolated from yeast studies. For example, acetylation of histone H3 lysine 56 (H3 K56Ac) is assumed to be important for transcriptional regulation in metazoan cells based on its occurrence at promoters and its function in yeast. Here we directly assess the function of H3 K56Ac during chromatin disassembly from gene regulatory regions during transcriptional induction in human cells by using mutations that either mimic or prevent H3 K56Ac. Although there is rapid histone H3 disassembly during induction of some estrogen receptor responsive genes, depletion of the histone chaperone ASF1A/B, which is required for H3 K56 acetylation, has no effect on chromatin disassembly at these regions. During the course of this work, we found that all the commercially available antibodies to H3 K56Ac are non-specific in human cells and in Drosophila. We used H3-YFP fusions to show that the H3 K56Q mutation can promote chromatin disassembly from regulatory regions of some estrogen responsive genes in the context of transcriptional induction. However, neither the H3 K56R nor K56Q mutation significantly altered chromatin disassembly dynamics by FRAP analysis. These results indicate that unlike the situation in yeast, human cells do not use H3 K56Ac to promote chromatin disassembly from regulatory regions or from the genome in general. Furthermore, our work highlights the need for rigorous characterization of the specificity of antibodies to histone post-translational modifications in vivo.

Phosphorylation of histone H3 threonine 118 (H3 T118ph) weakens histone DNA-contacts, disrupting the nucleosome structure. We show that Aurora-A mediated H3 T118ph occurs at pericentromeres and chromosome arms during prophase and is lost upon chromosome alignment. Expression of H3 T118E or H3 T118I (a SIN mutation that bypasses the need for the ATP-dependent nucleosome remodeler SWI/SNF) leads to mitotic problems including defects in spindle attachment, delayed cytokinesis, reduced chromatin packaging, cohesion loss, cohesin and condensin I loss in human cells. In agreement, overexpression of Aurora-A leads to increased H3 T118ph levels, causing cohesion loss, and reduced levels of cohesin and condensin I on chromatin. Normal levels of H3 T118ph are important because it is required for development in fruit flies. We propose that H3 T118ph alters the chromatin structure during specific phases of mitosis to promote timely condensin I and cohesin disassociation, which is essential for effective chromosome segregation.

Nucleosome loss from a promoter region has recently been described as a potential mechanism for transcriptional regulation. We investigated whether H3/H4 histone chaperones mediate the loss of nucleosomes from the promoter of the yeast PHO5 gene during transcriptional activation. We found that antisilencing function 1 (Asf1p) mediates nucleosome disassembly from the PHO5 promoter in vivo. We show that nucleosome disassembly also occurs at a second promoter, that of the PHO8 gene, during activation, and we demonstrate that this is also mediated by Asf1p. Furthermore, we show that nucleosome disassembly is essential for PHO5 and PHO8 activation. Contrary to the current dogma, we demonstrate that nucleosome disassembly is not required to enable binding of the Pho4p activator to its PHO5 UASp2 site in vivo. Finally, we show that nucleosomes are reassembled over the PHO5 promoter during repression. As such, nucleosome disassembly and reassembly are important mechanisms for transcriptional activation and repression, respectively.

Alzheimer's disease (AD) is a common neurodegenerative disease. APOE is the strong genetic risk factor of AD. The existing genome-wide association studies have identified many single nucleotide polymorphisms (SNPs) with minor effects on AD risk and the polygenic risk score (PRS) is presented to combine the effect of these SNPs. On the other hand, the volumes of various brain regions in AD patients have significant changes compared to that in normal individuals. Ch4 brain region containing at least 90% cholinergic neurons is the most extensive and conspicuous in the basal forebrain. Here, we investigated the relationship between the combined effect of AD-associated SNPs and Ch4 volume using the PRS approach. Our results showed that Ch4 volume in AD patients is significantly different from that in normal control subjects (p-value < 2.2 × 10-16). AD PRS, is not associated with the Ch4 volume in AD patients, excluding the APOE region (p-value = 0.264) and including the APOE region (p-value = 0.213). However, AD best-fit PRS, excluding the APOE region, is associated with Ch4 volume in normal control subjects (p-value = 0.015). AD PRS based on 8070 SNPs could explain 3.35% variance of Ch4 volume. In addition, the p-value of AD PRS model in normal control subjects, including the APOE region, is 0.006. AD PRS based on 8079 SNPs could explain 4.23% variance of Ch4 volume. In conclusion, PRS based on AD-associated SNPs is significantly related to Ch4 volume in normal subjects but not in patients.

Acetyl-CoA carboxylase (ACC) plays an important role in the metabolism of various cancer cells, but its role in head and neck squamous cell carcinoma (HNSCC) is uncertain. Therefore, in the present study, we explored the role of ACC2 in HNSCC.