Multi-modal molecular determinants of clinically relevant osteoporosis subtypes.

Due to a rapidly aging global population, osteoporosis and the associated risk of bone fractures have become a wide-spread public health problem. However, osteoporosis is very heterogeneous, and the existing standard diagnostic measure is not sufficient to accurately identify all patients at risk of osteoporotic fractures and to guide therapy. Here, we constructed the first prospective multi-omics atlas of the largest osteoporosis cohort to date (longitudinal data from 366 participants at three time points), and also implemented an explainable data-intensive analysis framework (DLSF: Deep Latent Space Fusion) for an omnigenic model based on a multi-modal approach that can capture the multi-modal molecular signatures (M3S) as explicit functional representations of hidden genotypes. Accordingly, through DLSF, we identified two subtypes of the osteoporosis population in Chinese individuals with corresponding molecular phenotypes, i.e., clinical intervention relevant subtypes (CISs), in which bone mineral density benefits response to calcium supplements in 2-year follow-up samples. Many snpGenes associated with these molecular phenotypes reveal diverse candidate biological mechanisms underlying osteoporosis, with xQTL preferences of osteoporosis and its subtypes indicating an omnigenic effect on different biological domains. Finally, these two subtypes were found to have different relevance to prior fracture and different fracture risk according to 4-year follow-up data. Thus, in clinical application, M3S could help us further develop improved diagnostic and treatment strategies for osteoporosis and identify a new composite index for fracture prediction, which were remarkably validated in an independent cohort (166 participants).

Oxidative stress-triggered Wnt signaling perturbation characterizes the tipping point of lung adeno-to-squamous transdifferentiation.

Lkb1 deficiency confers the Kras-mutant lung cancer with strong plasticity and the potential for adeno-to-squamous transdifferentiation (AST). However, it remains largely unknown how Lkb1 deficiency dynamically regulates AST. Using the classical AST mouse model (Kras LSL-G12D/+;Lkb1flox/flox, KL), we here comprehensively analyze the temporal transcriptomic dynamics of lung tumors at different stages by dynamic network biomarker (DNB) and identify the tipping point at which the Wnt signaling is abruptly suppressed by the excessive accumulation of reactive oxygen species (ROS) through its downstream effector FOXO3A. Bidirectional genetic perturbation of the Wnt pathway using two different Ctnnb1 conditional knockout mouse strains confirms its essential role in the negative regulation of AST. Importantly, pharmacological activation of the Wnt pathway before but not after the tipping point inhibits squamous transdifferentiation, highlighting the irreversibility of AST after crossing the tipping point. Through comparative transcriptomic analyses of mouse and human tumors, we find that the lineage-specific transcription factors (TFs) of adenocarcinoma and squamous cell carcinoma form a "Yin-Yang" counteracting network. Interestingly, inactivation of the Wnt pathway preferentially suppresses the adenomatous lineage TF network and thus disrupts the "Yin-Yang" homeostasis to lean towards the squamous lineage, whereas ectopic expression of NKX2-1, an adenomatous lineage TF, significantly dampens such phenotypic transition accelerated by the Wnt pathway inactivation. The negative correlation between the Wnt pathway and AST is further observed in a large cohort of human lung adenosquamous carcinoma. Collectively, our study identifies the tipping point of AST and highlights an essential role of the ROS-Wnt axis in dynamically orchestrating the homeostasis between adeno- and squamous-specific TF networks at the AST tipping point.

Dynamic network biomarker analysis and system pharmacology methods to explore the therapeutic effects and targets of Xiaoyaosan against liver cirrhosis.

ETHNOPHARMACOLOGICAL RELEVANCE: Xiaoyaosan is a traditional Chinese herbal formula that has long been used to treat liver cirrhosis, liver failure, and hepatocarcinoma (HCC). However, little is known about its mechanism of action and targets in treating chronic liver disease. AIM OF THE STUDY: This study aimed to detect the critical transition of HCC progression and to explore the regulatory mechanism and targets of Xiaoyaosan treating liver cirrhosis (cirrhosis) using integrative medicinal research involving system biology and pharmacology. MATERIALS AND METHODS: We recruited chronic liver disease participants to obtain gene expression data and applied the dynamic network biomarker (DNB) method to identify molecular markers and the critical transition. We combined network pharmacology and DNB analysis to locate the potential DNBs (targets). Then we validated the DNBs in the liver cirrhosis rat models using Xiaoyaosan treatment. The expression of genes encoding the four DNBs, including Cebpa, Csf1, Egfr, and Il7r, were further validated in rat liver tissue using Western blot analysis. RESULTS: We found EGFR, CEBPA, Csf1, Ccnb1, Rrmm2, C3, Il7r, Ccna2, and Peg10 overlap in the DNB list and Xiaoyaosan-Target-Disease (XTD) network constructed using network pharmacology databases. We investigated the diagnostic ability of each member in the DNB cluster and found EGFR, CEBPA, CSF1, and IL7R had high diagnostic abilities with AUC >0.7 and P-value < 0.05. We validated these findings in rats and found that liver function improved significantly and fibrotic changes were relieved in the Xiaoyaosan treatment group. The expression levels of CSF1 and IL7R in the Xiaoyaosan group were significantly lower than those in the cirrhosis model group. In contrast, CEBPA expression in the Xiaoyaosan group was significantly higher than that in the cirrhosis model group. The expression of EGFR in the Xiaoyaosan group was slightly decreased than in the model group but not significantly. CONCLUSION: Using the DNB method and network pharmacology approach, this study revealed that CEBPA, IL7R, EGFR, and CSF1 expression was remarkably altered in chronic liver disease and thus, may play an important role in driving the progression of cirrhosis. Therefore, CEBPA, IL7R, EGFR, and CSF1 may be important targets of Xiaoyaosan in treating cirrhosis and can be considered for developing novel therapeutics.

Dynamical network biomarkers: Theory and applications.

This paper reviews theory of DNB (Dynamical Network Biomarkers) and its applications including both modern medicine and traditional medicine. We show that omics data such as gene/protein expression profiles can be effectively used to detect pre-disease states before critical transitions from healthy states to disease states by using the DNB theory. The DNB theory with big biological data is expected to lead to ultra-early precision and preventive medicine.

Corrigendum to "Suppression of Dynamical Network Biomarker Signals at the Predisease State (Mibyou) before Metabolic Syndrome in Mice by a Traditional Japanese Medicine (Kampo Formula) Bofutsushosan".

[This corrects the article DOI: 10.1155/2020/9129134.].

The Effects of Graded Levels of Calorie Restriction: XVI. Metabolomic Changes in the Cerebellum Indicate Activation of Hypothalamocerebellar Connections Driven by Hunger Responses.

Calorie restriction (CR) remains the most robust intervention to extend life span and improve healthspan. Though the cerebellum is more commonly associated with motor control, it has strong links with the hypothalamus and is thought to be associated with nutritional regulation and adiposity. Using a global mass spectrometry-based metabolomics approach, we identified 756 metabolites that were significantly differentially expressed in the cerebellar region of the brain of C57BL/6J mice, fed graded levels of CR (10, 20, 30, and 40 CR) compared to mice fed ad libitum for 12 hours a day. Pathway enrichment indicated changes in the pathways of adenosine and guanine (which are precursors of DNA production), aromatic amino acids (tyrosine, phenylalanine, and tryptophan) and the sulfur-containing amino acid methionine. We also saw increases in the tricarboxylic acid cycle (TCA) cycle, electron donor, and dopamine and histamine pathways. In particular, changes in l-histidine and homocarnosine correlated positively with the level of CR and food anticipatory activity and negatively with insulin and body temperature. Several metabolic and pathway changes acted against changes seen in age-associated neurodegenerative disorders, including increases in the TCA cycle and reduced l-proline. Carnitine metabolites contributed to discrimination between CR groups, which corroborates previous work in the liver and plasma. These results indicate the conservation of certain aspects of metabolism across tissues with CR. Moreover, this is the first study to indicate CR alters the cerebellar metabolome, and does so in a graded fashion, after only a short period of restriction.

Suppression of Dynamical Network Biomarker Signals at the Predisease State (Mibyou) before Metabolic Syndrome in Mice by a Traditional Japanese Medicine (Kampo Formula) Bofutsushosan.

Due to the increasing incidence of metabolic syndrome, the development of new therapeutic strategies is urgently required. One promising approach is to focus on the predisease state (so-called Mibyou in traditional Japanese medicine) before metabolic syndrome as a preemptive medical target. We recently succeeded in detecting a predisease state before metabolic syndrome using a mathematical theory called the dynamical network biomarker (DNB) theory. The detected predisease state was characterized by 147 DNB genes among a total of 24,217 genes in TSOD (Tsumura-Suzuki Obese Diabetes) mice, a well-accepted model of metabolic syndrome, at 5 weeks of age. The timing of the predisease state was much earlier than the onset of metabolic syndrome in TSOD mice reported to be at approximately 8-12 weeks of age. In the present study, we investigated whether the predisease state in TSOD mice can be inhibited by the oral administration of a Kampo formula, bofutsushosan (BTS), which is usually used to treat obese patients with metabolic syndrome in Japan, from 3 to 7 weeks of age. We found the comprehensive suppression of the early warning signals of the DNB genes by BTS at 5 weeks of age and later. Specifically, the standard deviations of 134 genes among the 147 DNB genes decreased at 5 weeks of age as compared to the nontreatment control group, and 80 of them showed more than 50% reduction. In addition, at 7 weeks of age, the body weight and blood glucose level were significantly lower in the BTS-treated group than in the nontreatment control group. The results of our study suggest a novel mechanism of BTS; it suppressed fluctuations of the DNB genes at the predisease state before metabolic syndrome and thus prevented the subsequent transition to metabolic syndrome. In conclusion, this study demonstrated the preventive and preemptive effects of a Kampo formula on Mibyou before metabolic syndrome for the first time based on scientific evaluation.

The Effects of Graded Levels of Calorie Restriction XV: Phase Space Attractors Reveal Distinct Behavioral Phenotypes.

Calorie restriction (CR) has a positive impact on health and life span. Previous work, however, does not reveal the whole underlying mechanism of behavioral phenotypes under CR. We propose a new approach based on phase space reconstruction (PSR) to analyze the behavioral responses of mice to graded CR. This involved reconstructing high-dimensional attractors which topologically represent the intrinsic dynamics of mice based on low-dimensional time series of movement counts observed during the 90-day time course of restriction. PSR together with correlation dimensions (CD), Kolmogorov entropy (KE), and multifractal spectra builds a map from internal attractors to the phenotype of mice and reveals the mice with increasing CR levels undergo significant changes from a normal to a new state. Features of the attractors (CD and KE) were significantly associated with gene expression profiles in the hypothalamus of the same individuals.

Chronic hepatitis B: dynamic change in Traditional Chinese Medicine syndrome by dynamic network biomarkers.

BACKGROUND: In traditional Chinese medicine (TCM) clinical practice, TCM syndromes help to understand human homeostasis and guide individualized treatment. However, the TCM syndrome changes with disease progression, of which the scientific basis and mechanism remain unclear. METHODS: To demonstrate the underlying mechanism of dynamic changes in the TCM syndrome, we applied a dynamic network biomarker (DNB) algorithm to obtain the DNBs of changes in the TCM syndrome, based on the transcriptomic data of patients with chronic hepatitis B and typical TCM syndromes, including healthy controls and patients with liver-gallbladder dampness-heat syndrome (LGDHS), liver-depression spleen-deficiency syndrome (LDSDS), and liver-kidney yin-deficiency syndrome (LKYDS). The DNB model exploits collective fluctuations and correlations of the observed genes, then diagnoses the critical state. RESULTS: Our results showed that the DNBs of TCM syndromes were comprised of 52 genes and the tipping point occurred at the LDSDS stage. Meanwhile, there were numerous differentially expressed genes between LGDHS and LKYDS, which highlighted the drastic changes before and after the tipping point, implying the 52 DNBs could serve as early-warning signals of the upcoming change in the TCM syndrome. Next, we validated DNBs by cytokine profiling and isobaric tags for relative and absolute quantitation (iTRAQ). The results showed that PLG (plasminogen) and coagulation factor XII (F12) were significantly expressed during the progression of TCM syndrome from LGDHS to LKYDS. CONCLUSIONS: This study provides a scientific understanding of changes in the TCM syndrome. During this process, the cytokine system was involved all the time. The DNBs PLG and F12 were confirmed to significantly change during TCM-syndrome progression and indicated a potential value of DNBs in auxiliary diagnosis of TCM syndrome in CHB.Trial registration Identifier: NCT03189992. Registered on June 4, 2017. Retrospectively registered (http://www.clinicaltrials.gov).

The effects of graded levels of calorie restriction: VII. Topological rearrangement of hypothalamic aging networks.

Connectivity in a gene-gene network declines with age, typically within gene clusters. We explored the effect of short-term (3 months) graded calorie restriction (CR) (up to 40 %) on network structure of aging-associated genes in the murine hypothalamus by using conditional mutual information. The networks showed a topological rearrangement when exposed to graded CR with a higher relative within cluster connectivity at 40CR. We observed changes in gene centrality concordant with changes in CR level, with Ppargc1a, and Ppt1 having increased centrality and Etfdh, Traf3 and Abcc1 decreased centrality as CR increased. This change in gene centrality in a graded manner with CR, occurred in the absence of parallel changes in gene expression levels. This study emphasizes the importance of augmenting traditional differential gene expression analyses to better understand structural changes in the transcriptome. Overall our results suggested that CR induced changes in centrality of biological relevant genes that play an important role in preventing the age-associated loss of network integrity irrespective of their gene expression levels.

The effects of graded levels of calorie restriction: VI. Impact of short-term graded calorie restriction on transcriptomic responses of the hypothalamic hunger and circadian signaling pathways.

Food intake and circadian rhythms are regulated by hypothalamic neuropeptides and circulating hormones, which could mediate the anti-ageing effect of calorie restriction (CR). We tested whether these two signaling pathways mediate CR by quantifying hypothalamic transcripts of male C57BL/6 mice exposed to graded levels of CR (10 % to 40 %) for 3 months. We found that the graded CR manipulation resulted in upregulation of core circadian rhythm genes, which correlated negatively with circulating levels of leptin, insulin-like growth factor 1 (IGF-1), insulin, and tumor necrosis factor alpha (TNF-α). In addition, key components in the hunger signaling pathway were expressed in a manner reflecting elevated hunger at greater levels of restriction, and which also correlated negatively with circulating levels of insulin, TNF-α, leptin and IGF-1. Lastly, phenotypes, such as food anticipatory activity and body temperature, were associated with expression levels of both hunger genes and core clock genes. Our results suggest modulation of the hunger and circadian signaling pathways in response to altered levels of circulating hormones, that are themselves downstream of morphological changes resulting from CR treatment, may be important elements in the response to CR, driving some of the key phenotypic outcomes.

Analysing the Effect of Mutation on Protein Function and Discovering Potential Inhibitors of CDK4: Molecular Modelling and Dynamics Studies.

The cyclin-dependent kinase 4 (CDK4)-cyclin D1 complex plays a crucial role in the transition from the G1 phase to S phase of the cell cycle. Among the CDKs, CDK4 is one of the genes most frequently affected by somatic genetic variations that are associated with various forms of cancer. Thus, because the abnormal function of the CDK4-cyclin D1 protein complex might play a vital role in causing cancer, CDK4 can be considered a genetically validated therapeutic target. In this study, we used a systematic, integrated computational approach to identify deleterious nsSNPs and predict their effects on protein-protein (CDK4-cyclin D1) and protein-ligand (CDK4-flavopiridol) interactions. This analysis resulted in the identification of possible inhibitors of mutant CDK4 proteins that bind the conformations induced by deleterious nsSNPs. Using computational prediction methods, we identified five nsSNPs as highly deleterious: R24C, Y180H, A205T, R210P, and R246C. From molecular docking and molecular dynamic studies, we observed that these deleterious nsSNPs affected CDK4-cyclin D1 and CDK4-flavopiridol interactions. Furthermore, in a virtual screening approach, the drug 5_7_DIHYDROXY_ 2_ (3_4_5_TRI HYDROXYPHENYL) _4H_CHROMEN_ 4_ONE displayed good binding affinity for proteins with the mutations R24C or R246C, the drug diosmin displayed good binding affinity for the protein with the mutation Y180H, and the drug rutin displayed good binding affinity for proteins with the mutations A205T and R210P. Overall, this computational investigation of the CDK4 gene highlights the link between genetic variation and biological phenomena in human cancer and aids in the discovery of molecularly targeted therapies for personalized treatment.

Comparative genomics reveals new candidate genes involved in selenium metabolism in prokaryotes.

Selenium (Se) is an important micronutrient that mainly occurs in proteins in the form of selenocysteine and in tRNAs in the form of selenouridine. In the past 20 years, several genes involved in Se utilization have been characterized in both prokaryotes and eukaryotes. However, Se homeostasis and the associated regulatory network are not fully understood. In this study, we conducted comparative genomics and phylogenetic analyses to examine the occurrence of all known Se utilization traits in prokaryotes. Our results revealed a highly mosaic pattern of species that use Se (in different forms) in spite that most organisms do not use this element. Further investigation of genomic context of known Se-related genes in different organisms suggested novel candidate genes that may participate in Se metabolism in bacteria and/or archaea. Among them, a membrane protein, YedE, which contains ten transmembrane domains and shows distant similarity to a sulfur transporter, is exclusively found in Se-utilizing organisms, suggesting that it may be involved in Se transport. A LysR-like transcription factor subfamily might be important for the regulation of Sec biosynthesis and/or other Se-related genes. In addition, a small protein family DUF3343 is widespread in Se-utilizing organisms, which probably serves as an important chaperone for Se trafficking within the cells. Finally, we proposed a simple model of Se homeostasis based on our findings. Our study reveals new candidate genes involved in Se metabolism in prokaryotes and should be useful for a further understanding of the complex metabolism and the roles of Se in biology.