Regulon active landscape reveals cell development and functional state changes of human primary osteoblasts in vivo.

BACKGROUND: While transcription factor (TF) regulation is known to play an important role in osteoblast development, differentiation, and bone metabolism, the molecular features of TFs in human osteoblasts at the single-cell resolution level have not yet been characterized. Here, we identified modules (regulons) of co-regulated genes by applying single-cell regulatory network inference and clustering to the single-cell RNA sequencing profiles of human osteoblasts. We also performed cell-specific network (CSN) analysis, reconstructed regulon activity-based osteoblast development trajectories, and validated the functions of important regulons both in vivo and in vitro. RESULTS: We identified four cell clusters: preosteoblast-S1, preosteoblast-S2, intermediate osteoblasts, and mature osteoblasts. CSN analysis results and regulon activity-based osteoblast development trajectories revealed cell development and functional state changes of osteoblasts. CREM and FOSL2 regulons were mainly active in preosteoblast-S1, FOXC2 regulons were mainly active in intermediate osteoblast, and RUNX2 and CREB3L1 regulons were most active in mature osteoblasts. CONCLUSIONS: This is the first study to describe the unique features of human osteoblasts in vivo based on cellular regulon active landscapes. Functional state changes of CREM, FOSL2, FOXC2, RUNX2, and CREB3L1 regulons regarding immunity, cell proliferation, and differentiation identified the important cell stages or subtypes that may be predominantly affected by bone metabolism disorders. These findings may lead to a deeper understanding of the mechanisms underlying bone metabolism and associated diseases.

Divergent clonal selection dominates medulloblastoma at recurrence.

The development of targeted anti-cancer therapies through the study of cancer genomes is intended to increase survival rates and decrease treatment-related toxicity. We treated a transposon-driven, functional genomic mouse model of medulloblastoma with 'humanized' in vivo therapy (microneurosurgical tumour resection followed by multi-fractionated, image-guided radiotherapy). Genetic events in recurrent murine medulloblastoma exhibit a very poor overlap with those in matched murine diagnostic samples (<5%). Whole-genome sequencing of 33 pairs of human diagnostic and post-therapy medulloblastomas demonstrated substantial genetic divergence of the dominant clone after therapy (<12% diagnostic events were retained at recurrence). In both mice and humans, the dominant clone at recurrence arose through clonal selection of a pre-existing minor clone present at diagnosis. Targeted therapy is unlikely to be effective in the absence of the target, therefore our results offer a simple, proximal, and remediable explanation for the failure of prior clinical trials of targeted therapy.

Subgroup-specific structural variation across 1,000 medulloblastoma genomes.

Medulloblastoma, the most common malignant paediatric brain tumour, is currently treated with nonspecific cytotoxic therapies including surgery, whole-brain radiation, and aggressive chemotherapy. As medulloblastoma exhibits marked intertumoural heterogeneity, with at least four distinct molecular variants, previous attempts to identify targets for therapy have been underpowered because of small samples sizes. Here we report somatic copy number aberrations (SCNAs) in 1,087 unique medulloblastomas. SCNAs are common in medulloblastoma, and are predominantly subgroup-enriched. The most common region of focal copy number gain is a tandem duplication of SNCAIP, a gene associated with Parkinson's disease, which is exquisitely restricted to Group 4α. Recurrent translocations of PVT1, including PVT1-MYC and PVT1-NDRG1, that arise through chromothripsis are restricted to Group 3. Numerous targetable SCNAs, including recurrent events targeting TGF-ß signalling in Group 3, and NF-κB signalling in Group 4, suggest future avenues for rational, targeted therapy.

How the mountain pine beetle (Dendroctonus ponderosae) breached the Canadian Rocky Mountains.

The mountain pine beetle (MPB; Dendroctonus ponderosae Hopkins), a major pine forest pest native to western North America, has extended its range north and eastward during an ongoing outbreak. Determining how the MPB has expanded its range to breach putative barriers, whether physical (nonforested prairie and high elevation of the Rocky Mountains) or climatic (extreme continental climate where temperatures can be below -40 °C), may contribute to our general understanding of range changes as well as management of the current epidemic. Here, we use a panel of 1,536 single nucleotide polymorphisms (SNPs) to assess population genetic structure, connectivity, and signals of selection within this MPB range expansion. Biallelic SNPs in MPB from southwestern Canada revealed higher genetic differentiation and lower genetic connectivity than in the northern part of its range. A total of 208 unique SNPs were identified using different outlier detection tests, of which 32 returned annotations for products with putative functions in cholesterol synthesis, actin filament contraction, and membrane transport. We suggest that MPB has been able to spread beyond its previous range by adjusting its cellular and metabolic functions, with genome scale differentiation enabling populations to better withstand cooler climates and facilitate longer dispersal distances. Our study is the first to assess landscape-wide selective adaptation in an insect. We have shown that interrogation of genomic resources can identify shifts in genetic diversity and putative adaptive signals in this forest pest species.

Efficient Determination of the Enantiomeric Purity and Absolute Configuration of Flavanones by Using (S)-3,3'-Dibromo-1,1'-bi-2-naphthol as a Chiral Solvating Agent.

The enantiomeric purity and absolute configuration of flavanones were first determined using (S)-3,3'-dibromo-1,1'-bi-2-naphthol as a chiral solvating agent by means of (1)H NMR spectroscopy. The enantiomeric purity results closely matched those based on chiral HPLC analysis.

[Advances in post-translational modifications of cGAS in innate immunity].

As a member of the nucleotidyltransferase family, cyclic guanosine monophosphate-adenosine monophosphate synthase (cyclic GMP-AMP synthase, or cGAS) is primarily involved in innate immunity as a nucleic acids sensor that activates its downstream pathway and regulates type I interferon synthesis. The regulation of cGAS function is correlated with the bacterial and viral infections, autoimmune diseases, tumors, and other diseases. Besides, post-translational modification is one of the most in-depth and extensive ways of cGAS function adjustment. There are mainly six types of post-translational modifications (PTMs) of cGAS, including phosphorylation, acetylation, ubiquitination, sumoylation, peptide chain cleavage, and glutamylation. This article not only systematically summarizes how PTMs of cGAS regulate the functions of cGAS under different physiological and pathological conditions, but also probes deep into the potential of PTMs as therapeutic targets.

Spatial Dissection of the Distinct Cellular Responses to Normal Aging and Alzheimer's Disease in Human Prefrontal Cortex at Single-Nucleus Resolution.

Aging significantly elevates the risk for Alzheimer's disease (AD), contributing to the accumulation of AD pathologies, such as amyloid-ß (Aß), inflammation, and oxidative stress. The human prefrontal cortex (PFC) is highly vulnerable to the impacts of both aging and AD. Unveiling and understanding the molecular alterations in PFC associated with normal aging (NA) and AD is essential for elucidating the mechanisms of AD progression and developing novel therapeutics for this devastating disease. In this study, for the first time, we employed a cutting-edge spatial transcriptome platform, STOmics® SpaTial Enhanced Resolution Omics-sequencing (Stereo-seq), to generate the first comprehensive, subcellular resolution spatial transcriptome atlas of the human PFC from six AD cases at various neuropathological stages and six age, sex, and ethnicity matched controls. Our analyses revealed distinct transcriptional alterations across six neocortex layers, highlighted the AD-associated disruptions in laminar architecture, and identified changes in layer-to-layer interactions as AD progresses. Further, throughout the progression from NA to various stages of AD, we discovered specific genes that were significantly upregulated in neurons experiencing high stress and in nearby non-neuronal cells, compared to cells distant from the source of stress. Notably, the cell-cell interactions between the neurons under the high stress and adjacent glial cells that promote Aß clearance and neuroprotection were diminished in AD in response to stressors compared to NA. Through cell-type specific gene co-expression analysis, we identified three modules in excitatory and inhibitory neurons associated with neuronal protection, protein dephosphorylation, and negative regulation of Aß plaque formation. These modules negatively correlated with AD progression, indicating a reduced capacity for toxic substance clearance in AD subject samples. Moreover, we have discovered a novel transcription factor, ZNF460, that regulates all three modules, establishing it as a potential new therapeutic target for AD. Overall, utilizing the latest spatial transcriptome platform, our study developed the first transcriptome-wide atlas with subcellular resolution for assessing the molecular alterations in the human PFC due to AD. This atlas sheds light on the potential mechanisms underlying the progression from NA to AD.

Lyn attenuates sepsis-associated acute kidney injury by inhibition of phospho-STAT3 and apoptosis.

Sepsis-associated acute kidney injury (SA-AKI) is a life-threatening condition associated with high mortality and morbidity. However, the underlying pathogenesis of SA-AKI is still unclear. Lyn belongs to Src family kinases (SFKs), which exert numerous biological functions including modulation in receptor-mediated intracellular signaling and intercellular communication. Previous studies demonstrated that Lyn gene deletion obviously aggravates LPS-induced lung inflammation, but the role and possible mechanism of Lyn in SA-AKI have not been reported yet. Here, we found that Lyn protected against renal tubular injury in cecal ligation and puncture (CLP) induced AKI mouse model by inhibition of signal transducer and activator of transcription 3 (STAT3) phosphorylation and cell apoptosis. Moreover, Lyn agonist MLR-1023 pretreatment improved renal function, inhibited STAT3 phosphorylation and decreased cell apoptosis. Thus, Lyn appears to play a crucial role in orchestrating STAT3-mediated inflammation and cell apoptosis in SA-AKI. Hence, Lyn kinase may be a promising therapeutic target for SA-AKI.

The clinical utility of integrative genomics in childhood cancer extends beyond targetable mutations.

We conducted integrative somatic-germline analyses by deeply sequencing 864 cancer-associated genes, complete genomes and transcriptomes for 300 mostly previously treated children and adolescents/young adults with cancer of poor prognosis or with rare tumors enrolled in the SickKids Cancer Sequencing (KiCS) program. Clinically actionable variants were identified in 56% of patients. Improved diagnostic accuracy led to modified management in a subset. Therapeutically targetable variants (54% of patients) were of unanticipated timing and type, with over 20% derived from the germline. Corroborating mutational signatures (SBS3/BRCAness) in patients with germline homologous recombination defects demonstrates the potential utility of PARP inhibitors. Mutational burden was significantly elevated in 9% of patients. Sequential sampling identified changes in therapeutically targetable drivers in over one-third of patients, suggesting benefit from rebiopsy for genomic analysis at the time of relapse. Comprehensive cancer genomic profiling is useful at multiple points in the care trajectory for children and adolescents/young adults with cancer, supporting its integration into early clinical management.

Endothelial cell polarity and extracellular matrix composition require functional ATP6AP2 during developmental and pathological angiogenesis.

The (Pro)renin receptor ([P]RR), also known as ATP6AP2, is a single-transmembrane protein that is implicated in a multitude of biological processes. However, the exact role of ATP6AP2 during blood vessel development remains largely undefined. Here, we use an inducible endothelial cell-specific (EC-specific) Atp6ap2-KO mouse model to investigate the role of ATP6AP2 during both physiological and pathological angiogenesis in vivo. We observed that postnatal deletion of Atp6ap2 in ECs results in cell migration defects, loss of tip cell polarity, and subsequent impairment of retinal angiogenesis. In vitro, Atp6ap2-deficient ECs similarly displayed reduced cell migration, impaired sprouting, and defective cell polarity. Transcriptional profiling of ECs isolated from Atp6ap2 mutant mice further indicated regulatory roles in angiogenesis, cell migration, and extracellular matrix composition. Mechanistically, we provided evidence that expression of various extracellular matrix components is controlled by ATP6AP2 via the ERK pathway. Furthermore, Atp6ap2-deficient retinas exhibited reduced revascularization in an oxygen-induced retinopathy model. Collectively, our results demonstrate a critical role of ATP6AP2 as a regulator of developmental and pathological angiogenesis.

Tumor necrosis factor overcomes immune evasion in p53-mutant medulloblastoma.

Many immunotherapies act by enhancing the ability of cytotoxic T cells to kill tumor cells. Killing depends on T cell recognition of antigens presented by class I major histocompatibility complex (MHC-I) proteins on tumor cells. In this study, we showed that medulloblastomas lacking the p53 tumor suppressor do not express surface MHC-I and are therefore resistant to immune rejection. Mechanistically, this is because p53 regulates expression of the peptide transporter Tap1 and the aminopeptidase Erap1, which are required for MHC-I trafficking to the cell surface. In vitro, tumor necrosis factor (TNF) or lymphotoxin-ß receptor agonist can rescue expression of Erap1, Tap1 and MHC-I on p53-mutant tumor cells. In vivo, low doses of TNF prolong survival and synergize with immune checkpoint inhibitors to promote tumor rejection. These studies identified p53 as a key regulator of immune evasion and suggest that TNF could be used to enhance sensitivity of tumors to immunotherapy.

Structural evolution of LiN n + (n = 2, 4, 6, 8, and 10) clusters: mass spectrometry and theoretical calculations.

Mixed nitrogen-lithium cluster cations LiN n + were generated by laser vaporization and analyzed by time-of-flight mass spectrometry. It is found that LiN8 + has the highest ion abundance among the LiN n + ions in the mass spectrum. Density functional calculations were conducted to search for the stable structures of the Li-N clusters. The theoretical results show that the most stable isomers of LiN n + clusters are in the form of Li+(N2) n/2, and the order of their calculated binding energies is consistent with that of Li-N2 bond lengths. The most stable structures of LiN n + evolve from one-dimensional linear type (C ∞v, n = 2; D ∞h, n = 4), to two-dimensional branch type (D 3h, n = 6), then to three-dimensional tetrahedral (T d, n = 8) and square pyramid (C 4v, n = 10) types. Further natural bond orbital analyses show that electrons are transferred from the lone pair on Nα of every N2 unit to the empty orbitals of lithium atom in LiN2-8 +, while in LiN10 +, electrons are transferred from the bonding orbital of the Li-Nα bonds to the antibonding orbital of the other Li-Nα bonds. In both cases, the N2 units become dipoles and strongly interact with Li+. The average second-order perturbation stabilization energy for LiN8 + is the highest among the observed LiN n + clusters. For neutral LiN2-8 clusters, the most stable isomers were also formed by a Li atom and n/2 number of N2 units, while that of LiN10 is in the form of Li+(N2)3(η1-N4).

[Effects of recruitment maneuver on elderly patients after major operations].

OBJECTIVE: To study the effect of recruitment maneuver (RM) in preventing atelectasis and lung injury in elderly patients after major operations, and to evaluate the safety of RMs. METHODS: Forty elderly patients after major operations were admitted to intensive care unit (ICU) of Beijing Air Force General Hospital from February 2007 to February 2008 were randomized into RM group and control group. The patients were still under the effect of anesthesia and muscle relaxation when admitted. All of them were under invasive blood pressure monitoring, which was continued for over 6 hours. RM was conducted by regulating inspired oxygen concentration (FiO2) to 0.60, respiratory rate 20/min, tidal volume (VT) 5 ml/kg, with 25 cm H2O (1 cm H2O=0.098 kPa) of continuous positive end-expiratory pressure (PEEP) for 30 seconds, and then the previous ventilator setting was resumed. The above modality was repeated once after 1 hour. Heart rate (HR), central venous pressure (CVP), mean arterial pressure (MAP), platform airway pressure (Pplat), percutaneous oxygen saturation (SpO2) were measured before and after the RM. Arterial blood gas analysis was done before and after RMs. The presence of pulmonary atelectasis or pulmonary infection was looked for after RMs. RESULTS: (1) There were significant changes in HR, CVP and MAP during RM (all P<0.05), circulation function was not affected. (2) There were no significant changes in HR, MAP, CVP and SpO2 before and after RMs (all P>0.05). Pplat was significantly reduced after RMs (P<0.05). (3) The incidence of pulmonary atelectasis or pulmonary infection was significantly lower in RM group (both P<0.05). (4) Oxygenation index (PaO2/FiO2) in RM group was significant increased (P<0.05) . CONCLUSION: RM is safe when used in elderly patients. It can significantly improve oxygenation in elderly patients.

Opposing Effects of CREBBP Mutations Govern the Phenotype of Rubinstein-Taybi Syndrome and Adult SHH Medulloblastoma.

Recurrent mutations in chromatin modifiers are specifically prevalent in adolescent or adult patients with Sonic hedgehog-associated medulloblastoma (SHH MB). Here, we report that mutations in the acetyltransferase CREBBP have opposing effects during the development of the cerebellum, the primary site of origin of SHH MB. Our data reveal that loss of Crebbp in cerebellar granule neuron progenitors (GNPs) during embryonic development of mice compromises GNP development, in part by downregulation of brain-derived neurotrophic factor (Bdnf). Interestingly, concomitant cerebellar hypoplasia was also observed in patients with Rubinstein-Taybi syndrome, a congenital disorder caused by germline mutations of CREBBP. By contrast, loss of Crebbp in GNPs during postnatal development synergizes with oncogenic activation of SHH signaling to drive MB growth, thereby explaining the enrichment of somatic CREBBP mutations in SHH MB of adult patients. Together, our data provide insights into time-sensitive consequences of CREBBP mutations and corresponding associations with human diseases.

Pyruvate Kinase Inhibits Proliferation during Postnatal Cerebellar Neurogenesis and Suppresses Medulloblastoma Formation.

Aerobic glycolysis supports proliferation through unresolved mechanisms. We have previously shown that aerobic glycolysis is required for the regulated proliferation of cerebellar granule neuron progenitors (CGNP) and for the growth of CGNP-derived medulloblastoma. Blocking the initiation of glycolysis via deletion of hexokinase-2 (Hk2) disrupts CGNP proliferation and restricts medulloblastoma growth. Here, we assessed whether disrupting pyruvate kinase-M (Pkm), an enzyme that acts in the terminal steps of glycolysis, would alter CGNP metabolism, proliferation, and tumorigenesis. We observed a dichotomous pattern of PKM expression, in which postmitotic neurons throughout the brain expressed the constitutively active PKM1 isoform, while neural progenitors and medulloblastomas exclusively expressed the less active PKM2. Isoform-specific Pkm2 deletion in CGNPs blocked all Pkm expression. Pkm2-deleted CGNPs showed reduced lactate production and increased SHH-driven proliferation. 13C-flux analysis showed that Pkm2 deletion reduced the flow of glucose carbons into lactate and glutamate without markedly increasing glucose-to-ribose flux. Pkm2 deletion accelerated tumor formation in medulloblastoma-prone ND2:SmoA1 mice, indicating the disrupting PKM releases CGNPs from a tumor-suppressive effect. These findings show that distal and proximal disruptions of glycolysis have opposite effects on proliferation, and that efforts to block the oncogenic effect of aerobic glycolysis must target reactions upstream of PKM. Cancer Res; 77(12); 3217-30. ©2017 AACR.

Spatial heterogeneity in medulloblastoma.

Spatial heterogeneity of transcriptional and genetic markers between physically isolated biopsies of a single tumor poses major barriers to the identification of biomarkers and the development of targeted therapies that will be effective against the entire tumor. We analyzed the spatial heterogeneity of multiregional biopsies from 35 patients, using a combination of transcriptomic and genomic profiles. Medulloblastomas (MBs), but not high-grade gliomas (HGGs), demonstrated spatially homogeneous transcriptomes, which allowed for accurate subgrouping of tumors from a single biopsy. Conversely, somatic mutations that affect genes suitable for targeted therapeutics demonstrated high levels of spatial heterogeneity in MB, malignant glioma, and renal cell carcinoma (RCC). Actionable targets found in a single MB biopsy were seldom clonal across the entire tumor, which brings the efficacy of monotherapies against a single target into question. Clinical trials of targeted therapies for MB should first ensure the spatially ubiquitous nature of the target mutation.

Repositioning organohalogen drugs: a case study for identification of potent B-Raf V600E inhibitors via docking and bioassay.

Drug repositioning has been attracting increasingly attention for its advantages of reducing costs and risks. Statistics showed that around one quarter of the marketed drugs are organohalogens. However, no study has been reported, to the best of our knowledge, to aim at efficiently repositioning organohalogen drugs, which may be attributed to the lack of accurate halogen bonding scoring function. Here, we present a study to show that two organohalogen drugs were successfully repositioned as potent B-Raf V600E inhibitors via molecular docking with halogen bonding scoring function, namely D(3)DOCKxb developed in our lab, and bioassay. After virtual screening by D(3)DOCKxb against the database CMC (Comprehensive Medicinal Chemistry), 3 organohalogen drugs that were predicted to form strong halogen bonding with B-Raf V600E were purchased and tested with ELISA-based assay. In the end, 2 of them, rafoxanide and closantel, were identified as potent inhibitors with IC50 values of 0.07 µM and 1.90 µM, respectively, which are comparable to that of vemurafenib (IC50: 0.17 µM), a marketed drug targeting B-Raf V600E. Single point mutagenesis experiments confirmed the conformations predicted by D(3)DOCKxb. And comparison experiment revealed that halogen bonding scoring function is essential for repositioning those drugs with heavy halogen atoms in their molecular structures.

Medulloblastoma-associated DDX3 variant selectively alters the translational response to stress.

DDX3X encodes a DEAD-box family RNA helicase (DDX3) commonly mutated in medulloblastoma, a highly aggressive cerebellar tumor affecting both children and adults. Despite being implicated in several facets of RNA metabolism, the nature and scope of DDX3's interactions with RNA remain unclear. Here, we show DDX3 collaborates extensively with the translation initiation machinery through direct binding to 5'UTRs of nearly all coding RNAs, specific sites on the 18S rRNA, and multiple components of the translation initiation complex. Impairment of translation initiation is also evident in primary medulloblastomas harboring mutations in DDX3X, further highlighting DDX3's role in this process. Arsenite-induced stress shifts DDX3 binding from the 5'UTR into the coding region of mRNAs concomitant with a general reduction of translation, and both the shift of DDX3 on mRNA and decreased translation are blunted by expression of a catalytically-impaired, medulloblastoma-associated DDX3R534H variant. Furthermore, despite the global repression of translation induced by arsenite, translation is preserved on select genes involved in chromatin organization in DDX3R534H-expressing cells. Thus, DDX3 interacts extensively with RNA and ribosomal machinery to help remodel the translation landscape in response to stress, while cancer-related DDX3 variants adapt this response to selectively preserve translation.

JAGuaR: junction alignments to genome for RNA-seq reads.

JAGuaR is an alignment protocol for RNA-seq reads that uses an extended reference to increase alignment sensitivity. It uses BWA to align reads to the genome and reference transcript models (including annotated exon-exon junctions) specifically allowing for the possibility of a single read spanning multiple exons. Reads aligned to the transcript models are then re-mapped on to genomic coordinates, transforming alignments that span multiple exons into large-gapped alignments on the genome. While JAGuaR does not detect novel junctions, we demonstrate how JAGuaR generates fast and accurate transcriptome alignments, which allows for both sensitive and specific SNV calling.