Biological states are controlled by orchestrated transcriptional factors (TFs) within gene regulatory networks. Here we show TFs responsible for the dynamic changes of biological states can be prioritized with temporal PageRank. We further show such TF prioritization can be extended by integrating gene regulatory networks reverse engineered from multi-omics profiles, e.g. gene expression, chromatin accessibility, and chromosome conformation assays, using multiplex PageRank.

Transfer RNAs (tRNAs) are a central component for the biological synthesis of proteins, and they are among the most highly conserved and frequently transcribed genes in all living things. Despite their clear significance for fundamental cellular processes, the forces governing tRNA evolution are poorly understood. We present evidence that transcription-associated mutagenesis and strong purifying selection are key determinants of patterns of sequence variation within and surrounding tRNA genes in humans and diverse model organisms. Remarkably, the mutation rate at broadly expressed cytosolic tRNA loci is likely between 7 and 10 times greater than the nuclear genome average. Furthermore, evolutionary analyses provide strong evidence that tRNA genes, but not their flanking sequences, experience strong purifying selection acting against this elevated mutation rate. We also find a strong correlation between tRNA expression levels and the mutation rates in their immediate flanking regions, suggesting a simple method for estimating individual tRNA gene activity. Collectively, this study illuminates the extreme competing forces in tRNA gene evolution and indicates that mutations at tRNA loci contribute disproportionately to mutational load and have unexplored fitness consequences in human populations.

Background

Ozone exposure induces airway neutrophilia and modifies innate immune monocytic cell-surface phenotypes in healthy individuals. High-dose inhaled corticosteroids can reduce O(3)-induced airway inflammation, but their effect on innate immune activation is unknown.

Objectives

We used a human O(3) inhalation challenge model to examine the effectiveness of clinically relevant doses of inhaled corticosteroids on airway inflammation and markers of innate immune activation in healthy volunteers.

Methods

Seventeen O(3)-responsive subjects [>10% increase in the percentage of polymorphonuclear leukocytes (PMNs) in sputum, PMNs per milligram vs. baseline sputum] received placebo, or either a single therapeutic dose (0.5 mg) or a high dose (2 mg) of inhaled fluticasone proprionate (FP) 1 hr before a 3-hr O(3) challenge (0.25 ppm) on three separate occasions at least 2 weeks apart. Lung function, exhaled nitric oxide, sputum, and systemic biomarkers were assessed 1-5 hr after the O(3) challenge. To determine the effect of FP on cellular function, we assessed sputum cells from seven subjects by flow cytometry for cell-surface marker activation.

Results

FP had no effect on O(3)-induced lung function decline. Compared with placebo, 0.5 mg and 2 mg FP reduced O(3)-induced sputum neutrophilia by 18% and 35%, respectively. A similar effect was observed on the airway-specific serum biomarker Clara cell protein 16 (CCP16). Furthermore, FP pretreatment significantly reduced O(3)-induced modification of CD11b, mCD14, CD64, CD16, HLA-DR, and CD86 on sputum monocytes in a dose-dependent manner.

Conclusions

This study confirmed and extended data demonstrating the protective effect of FP against O(3)-induced airway inflammation and immune cell activation.

Background

Altered kinase localization is gaining appreciation as a mechanism of cardiovascular disease. Previous work suggests GSK-3β (glycogen synthase kinase 3β) localizes to and regulates contractile function of the myofilament. We aimed to discover GSK-3β's in vivo role in regulating myofilament function, the mechanisms involved, and the translational relevance.

Methods

Inducible cardiomyocyte-specific GSK-3β knockout mice and left ventricular myocardium from nonfailing and failing human hearts were studied.

Results

Skinned cardiomyocytes from knockout mice failed to exhibit calcium sensitization with stretch indicating a loss of length-dependent activation (LDA), the mechanism underlying the Frank-Starling Law. Titin acts as a length sensor for LDA, and knockout mice had decreased titin stiffness compared with control mice, explaining the lack of LDA. Knockout mice exhibited no changes in titin isoforms, titin phosphorylation, or other thin filament phosphorylation sites known to affect passive tension or LDA. Mass spectrometry identified several z-disc proteins as myofilament phospho-substrates of GSK-3β. Agreeing with the localization of its targets, GSK-3β that is phosphorylated at Y216 binds to the z-disc. We showed pY216 was necessary and sufficient for z-disc binding using adenoviruses for wild-type, Y216F, and Y216E GSK-3β in neonatal rat ventricular cardiomyocytes. One of GSK-3β's z-disc targets, abLIM-1 (actin-binding LIM protein 1), binds to the z-disc domains of titin that are important for maintaining passive tension. Genetic knockdown of abLIM-1 via siRNA in human engineered heart tissues resulted in enhancement of LDA, indicating abLIM-1 may act as a negative regulator that is modulated by GSK-3β. Last, GSK-3β myofilament localization was reduced in left ventricular myocardium from failing human hearts, which correlated with depressed LDA.

Conclusions

We identified a novel mechanism by which GSK-3β localizes to the myofilament to modulate LDA. Importantly, z-disc GSK-3β levels were reduced in patients with heart failure, indicating z-disc localized GSK-3β is a possible therapeutic target to restore the Frank-Starling mechanism in patients with heart failure.

Genomic sequence interpretation can miss clinically relevant missense variants for several reasons. Rare missense variants are numerous in the exome and difficult to prioritise. Affected genes may also not have existing disease association. To improve variant prioritisation, we leverage population exome data to identify intragenic missense-depleted regions (MDRs) genome-wide that may be important in disease. We then use missense depletion analyses to help prioritise undiagnosed disease exome variants. We demonstrate application of this strategy to identify a novel gene association for human brain malformation. We identified de novo missense variants that affect the GDP/GTP-binding site of ARF1 in three unrelated patients. Corresponding functional analysis suggests ARF1 GDP/GTP-activation is affected by the specific missense mutations associated with heterotopia. These findings expand the genetic pathway underpinning neurologic disease that classically includes FLNA. ARF1 along with ARFGEF2 add further evidence implicating ARF/GEFs in the brain. Using functional ontology, top MDR-containing genes were highly enriched for nucleotide-binding function, suggesting these may be candidates for human disease. Routine consideration of MDR in the interpretation of exome data for rare diseases may help identify strong genetic factors for many severe conditions, infertility/reduction in reproductive capability, and embryonic conditions contributing to preterm loss.

Abnormal activation of the epidermal growth factor receptor (EGFR) due to a deletion of exons 2-7 of EGFR (EGFRvIII) is a common alteration in glioblastoma (GBM). While this alteration can drive gliomagenesis, tumors harboring EGFRvIII are heterogeneous. To investigate the role for EGFRvIII activation in tumor phenotype we used a neural progenitor cell-based murine model of GBM driven by EGFR signaling and generated tumor progenitor cells with high and low EGFRvIII activation, pEGFRHi and pEGFRLo. In vivo, ex vivo, and in vitro studies suggested a direct association between EGFRvIII activity and increased tumor cell proliferation, decreased tumor cell adhesion to the extracellular matrix, and altered progenitor cell phenotype. Time-lapse confocal imaging of tumor cells in brain slice cultures demonstrated blood vessel co-option by tumor cells and highlighted differences in invasive pattern. Inhibition of EGFR signaling in pEGFRHi promoted cell differentiation and increased cell-matrix adhesion. Conversely, increased EGFRvIII activation in pEGFRLo reduced cell-matrix adhesion. Our study using a murine model for GBM driven by a single genetic driver, suggests differences in EGFR activation contribute to tumor heterogeneity and aggressiveness.

Pleomorphic xanthoastrocytoma (PXA) is an astrocytic neoplasm that is typically well circumscribed and can have a relatively favorable prognosis. Tumor progression to anaplastic PXA (WHO grade III), however, is associated with a more aggressive biologic behavior and worse prognosis. The factors that drive anaplastic progression are largely unknown. We performed comprehensive genomic profiling on a set of 23 PXAs from 19 patients, including 15 with anaplastic PXA. Four patients had tumor tissue from multiple recurrences, including two with anaplastic progression. We find that PXAs are genetically defined by the combination of CDKN2A biallelic inactivation and RAF alterations that were present in all 19 cases, most commonly as CDKN2A homozygous deletion and BRAF p.V600E mutation but also occasionally BRAF or RAF1 fusions or other rearrangements. The third most commonly altered gene in anaplastic PXA was TERT, with 47% (7/15) harboring TERT alterations, either gene amplification (n = 2) or promoter hotspot mutation (n = 5). In tumor pairs analyzed before and after anaplastic progression, two had increased copy number alterations and one had TERT promoter mutation at recurrence. Less commonly altered genes included TP53, BCOR, BCORL1, ARID1A, ATRX, PTEN, and BCL6. All PXA in this cohort were IDH and histone H3 wildtype, and did not contain alterations in EGFR. Genetic profiling performed on six regions from the same tumor identified intratumoral genomic heterogeneity, likely reflecting clonal evolution during tumor progression. Overall, anaplastic PXA is characterized by the combination of CDKN2A biallelic inactivation and oncogenic RAF kinase signaling as well as a relatively small number of additional genetic alterations, with the most common being TERT amplification or promoter mutation. These data define a distinct molecular profile for PXA and suggest additional genetic alterations, including TERT, may be associated with anaplastic progression.

Amplification of the epidermal growth factor receptor gene (EGFR) represents one of the most commonly observed genetic lesions in glioblastoma (GBM); however, therapies targeting this signaling pathway have failed clinically. Here, using human tumors, primary patient-derived xenografts (PDX), and a murine model for GBM, we demonstrate that EGFR inhibition leads to increased invasion of tumor cells. Further, EGFR inhibitor-treated GBM demonstrates altered oxidative stress, with increased lipid peroxidation, and generation of toxic lipid peroxidation products. A tumor cell subpopulation with elevated aldehyde dehydrogenase (ALDH) levels was determined to comprise a significant proportion of the invasive cells observed in EGFR inhibitor-treated GBM. Our analysis of the ALDH1A1 protein in newly diagnosed GBM revealed detectable ALDH1A1 expression in 69% (35/51) of the cases, but in relatively low percentages of tumor cells. Analysis of paired human GBM before and after EGFR inhibitor therapy showed an increase in ALDH1A1 expression in EGFR-amplified tumors (P < 0.05, n = 13 tumor pairs), and in murine GBM ALDH1A1-high clones were more resistant to EGFR inhibition than ALDH1A1-low clones. Our data identify ALDH levels as a biomarker of GBM cells with high invasive potential, altered oxidative stress, and resistance to EGFR inhibition, and reveal a therapeutic target whose inhibition should limit GBM invasion.

Amplification of the epidermal growth factor receptor gene (EGFR) represents one of the most commonly observed genetic lesions in glioblastoma (GBM); however, therapies targeting this signaling pathway have failed clinically. Here, using human tumors, primary patient-derived xenografts (PDX), and a murine model for GBM, we demonstrate that EGFR inhibition leads to increased invasion of tumor cells. Further, EGFR inhibitor-treated GBM demonstrates altered oxidative stress, with increased lipid peroxidation, and generation of toxic lipid peroxidation products. A tumor cell subpopulation with elevated aldehyde dehydrogenase (ALDH) levels was determined to comprise a significant proportion of the invasive cells observed in EGFR inhibitor-treated GBM. Our analysis of the ALDH1A1 protein in newly diagnosed GBM revealed detectable ALDH1A1 expression in 69% (35/51) of the cases, but in relatively low percentages of tumor cells. Analysis of paired human GBM before and after EGFR inhibitor therapy showed an increase in ALDH1A1 expression in EGFR-amplified tumors (P < 0.05, n = 13 tumor pairs), and in murine GBM ALDH1A1-high clones were more resistant to EGFR inhibition than ALDH1A1-low clones. Our data identify ALDH levels as a biomarker of GBM cells with high invasive potential, altered oxidative stress, and resistance to EGFR inhibition, and reveal a therapeutic target whose inhibition should limit GBM invasion.