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.

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.

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.

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.

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.

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.

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.

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.

Nowadays, the use of cryopreserved fat tissue for soft tissue augmentation is common, except for its unpredictable fat graft absorption, and the toxicity of the cryoprotective agent remains a limitation. In this study, the effects of freezing stored fat tissue without a cryoprotector, and the addition of the stromal vascular fraction (SVF) on the survival of cryopreserved transplants was studied.

Synthetic biology approach has been frequently applied to produce plant rare bioactive compounds in microbial cell factories by fermentation. However, to reach an ideal manufactural efficiency, it is necessary to optimize the microbial cell factories systemically by boosting sufficient carbon flux to the precursor synthesis and tuning the expression level and efficiency of key bioparts related to the synthetic pathway. We previously developed a yeast cell factory to produce ginsenoside Rh2 from glucose. However, the ginsenoside Rh2 yield was too low for commercialization due to the low supply of the ginsenoside aglycone protopanaxadiol (PPD) and poor performance of the key UDP-glycosyltransferase (UGT) (biopart UGTPg45) in the final step of the biosynthetic pathway. In the present study, we constructed a PPD-producing chassis via modular engineering of the mevalonic acid pathway and optimization of P450 expression levels. The new yeast chassis could produce 529.0 mg/L of PPD in shake flasks and 11.02 g/L in 10 L fed-batch fermentation. Based on this high PPD-producing chassis, we established a series of cell factories to produce ginsenoside Rh2, which we optimized by improving the C3-OH glycosylation efficiency. We increased the copy number of UGTPg45, and engineered its promoter to increase expression levels. In addition, we screened for more efficient and compatible UGT bioparts from other plant species and mutants originating from the direct evolution of UGTPg45. Combining all engineered strategies, we built a yeast cell factory with the greatest ginsenoside Rh2 production reported to date, 179.3 mg/L in shake flasks and 2.25 g/L in 10 L fed-batch fermentation. The results set up a successful example for improving yeast cell factories to produce plant rare natural products, especially the glycosylated ones.

Pseudomonas putida SJTE-1 can utilize 17β-estradiol (E2) as its carbon source, while the enzymes for E2 transformation in this strain is still unclear. 17β-hydroxysteroid dehydrogenases (17β-HSD) can catalyze the reduction/oxidation at C17 site of steroid hormone specifically, critical for steroid transformation. Here a novel 3-oxoacyl-(acyl-carrier protein) (ACP) reductase (ANI02794.1) was identified as it could bβ-estradiol, and was proved to be capable of functioning as 17β-HSD. Sequences alignment showed it contained the two consensus regions and the conserved residues of short-chain dehydrogenase/reductase (SDR). Its encoding gene was cloned and over-expressed in Escherichia coli BL21(DE3) strain, and the recombinant protein was purified by the metal-ion affinity chromatography with the yield of 18 mg/L culture. HPLC (High Performance Liquid Chromatography) detection showed this enzyme could convert 17β-estradiol into estrone using NAD+ as cofactor. Its Km value was 0.082 mM and its Vmax value was 0.81 mM/s; its transformation efficiency of 17β-estradiol into estrone was over 96.6% in five minutes. Its optimal temperature was 37 °C and optimal was pH 9.0; the divalent ions had different effects on the enzymatic activity. In conclusion, this 3-oxoacyl-ACP reductase functioned as 17β-HSD in P. putida SJTE-1 and played important role in its estrogen metabolism.

Alzheimer's disease (AD) and Parkinson's disease (PD) are the top two common neurodegenerative diseases in elderly. Recent studies found the α-synuclein have a key role in AD. Although many clinical and pathological features between AD and PD are shared, the genetic association between them remains unclear, especially whether α-synuclein in PD genetically alters AD risk.

Animal feeding is controlled by external sensory cues and internal metabolic states. Does it also depend on enteric neurons that sense mechanical cues to signal fullness of the digestive tract? Here, we identify a group of piezo-expressing neurons innervating the Drosophila crop (the fly equivalent of the stomach) that monitor crop volume to avoid food overconsumption. These neurons reside in the pars intercerebralis (PI), a neuro-secretory center in the brain involved in homeostatic control, and express insulin-like peptides with well-established roles in regulating food intake and metabolism. Piezo knockdown in these neurons of wild-type flies phenocopies the food overconsumption phenotype of piezo-null mutant flies. Conversely, expression of either fly Piezo or mammalian Piezo1 in these neurons of piezo-null mutants suppresses the overconsumption phenotype. Importantly, Piezo+ neurons at the PI are activated directly by crop distension, thus conveying a rapid satiety signal along the "brain-gut axis" to control feeding.

Accurate prediction of lymph-node metastasis in cancers is pivotal for the next targeted clinical interventions that allow favorable prognosis for patients. Different molecular profiles (mRNA and non-coding RNAs) have been widely used to establish classifiers for cancer prediction (e.g., tumor origin, cancerous or non-cancerous state, cancer subtype). However, few studies focus on lymphatic metastasis evaluation using these profiles, and the performance of classifiers based on different profiles has also not been compared. Here, differentially expressed mRNAs, miRNAs, and lncRNAs between lymph-node metastatic and non-metastatic groups were identified as molecular signatures to construct classifiers for lymphatic metastasis prediction in different cancers. With this similar feature selection strategy, support vector machine (SVM) classifiers based on different profiles were systematically compared in their prediction performance. For representative cancers (a total of nine types), these classifiers achieved comparative overall accuracies of 81.00% (67.96-92.19%), 81.97% (70.83-95.24%), and 80.78% (69.61-90.00%) on independent mRNA, miRNA, and lncRNA datasets, with a small set of biomarkers (6, 12, and 4 on average). Therefore, our proposed feature selection strategies are economical and efficient to identify biomarkers that aid in developing competitive classifiers for predicting lymph-node metastasis in cancers. A user-friendly webserver was also deployed to help researchers in metastasis risk determination by submitting their expression profiles of different origins.

Cardiovascular benefits of empagliflozin in patients with type 2 diabetes mellitus (T2DM) have been reported; however, the underlying mechanism remains unknown.