Multi-ancestry GWAS analysis identifies two novel loci associated with diabetic eye disease and highlights APOL1 as a high risk locus in patients with diabetic macular edema.

Diabetic retinopathy (DR) is a common complication of diabetes. Approximately 20% of DR patients have diabetic macular edema (DME) characterized by fluid leakage into the retina. There is a genetic component to DR and DME risk, but few replicable loci. Because not all DR cases have DME, we focused on DME to increase power, and conducted a multi-ancestry GWAS to assess DME risk in a total of 1,502 DME patients and 5,603 non-DME controls in discovery and replication datasets. Two loci reached GWAS significance (p<5x10-8). The strongest association was rs2239785, (K150E) in APOL1. The second finding was rs10402468, which co-localized to PLVAP and ANKLE1 in vascular / endothelium tissues. We conducted multiple sensitivity analyses to establish that the associations were specific to DME status and did not reflect diabetes status or other diabetic complications. Here we report two novel loci for risk of DME which replicated in multiple clinical trial and biobank derived datasets. One of these loci, containing the gene APOL1, is a risk factor in African American DME and DKD patients, indicating that this locus plays a broader role in diabetic complications for multiple ancestries. Trial Registration: NCT00473330, NCT00473382, NCT03622580, NCT03622593, NCT04108156.

Diffusion tensor MRI is sensitive to fibrotic injury in a mouse model of oxalate-induced chronic kidney disease.

Chronic kidney disease (CKD) is characterized by inflammation and fibrosis in the kidney. Renal biopsies and estimated glomerular filtration rate (eGFR) remain the standard of care, but these endpoints have limitations in detecting the stage, progression, and spatial distribution of fibrotic pathology in the kidney. MRI diffusion tensor imaging (DTI) has emerged as a promising noninvasive technology to evaluate renal fibrosis in vivo both in clinical and preclinical studies. However, these imaging studies have not systematically identified fibrosis particularly deeper in the kidney where biopsy sampling is limited, or completed an extensive analysis of whole organ histology, blood biomarkers, and gene expression to evaluate the relative strengths and weaknesses of MRI for evaluating renal fibrosis. In this study, we performed DTI in the sodium oxalate mouse model of CKD. The DTI parameters fractional anisotropy, apparent diffusion coefficient, and axial diffusivity were compared between the control and oxalate groups with region of interest (ROI) analysis to determine changes in the cortex and medulla. In addition, voxel-based analysis (VBA) was implemented to systematically identify local regions of injury over the whole kidney. DTI parameters were found to be significantly different in the medulla by both ROI analysis and VBA, which also spatially matched with collagen III immunohistochemistry (IHC). The DTI parameters in this medullary region exhibited moderate to strong correlations with histology, blood biomarkers, hydroxyproline, and gene expression. Our results thus highlight the sensitivity of DTI to the heterogeneity of renal fibrosis and importance of whole kidney noninvasive imaging.NEW & NOTEWORTHY Chronic kidney disease (CKD) can be characterized by inflammation and fibrosis of the kidney. Although standard of care methods have been limited in scope, safety, and spatial distribution, MRI diffusion tensor imaging (DTI) has emerged as a promising noninvasive technology to evaluate renal fibrosis in vivo. In this study, we performed DTI in an oxalate mouse model of CKD to systematically identify local kidney injury. DTI parameters strongly correlated with histology, blood biomarkers, hydroxyproline, and gene expression.

Unraveling Pre-filled Syringe Needle Clogging: Exploring a Fresh Outlook Through Innovative Techniques.

OBJECTIVE: This study aimed to investigate the movement of liquid in the needle region of staked-in-needle pre-filled syringes using neutron imaging and synchrotron X-ray tomography. The objective was to gain insights into the dynamics of liquid presence and understand the factors contributing to needle clogging. METHODS: Staked-in-needle pre-filled syringes were examined using neutron radiography and synchrotron X-ray phase-contrast computed tomography. Neutron radiography provided a 2D visualization of liquid presence in the needle, while synchrotron X-ray tomography offered high-resolution 3D imaging to study detailed morphological features of the liquid. RESULTS: Neutron radiography revealed liquid presence in the needle region for as-received samples and after temperature and pressure cycling. Pressure cycling had a more pronounced effect on liquid formation. Synchrotron X-ray tomography confirmed the presence of liquid and revealed various morphologies, including droplets of different sizes, liquid segments blocking sections of the needle, and a thin layer covering the needle wall. Liquid presence was also observed between the steel needle and the glass barrel. CONCLUSIONS: The combination of neutron imaging and synchrotron X-ray tomography provided valuable insights into the dynamics of liquid movement in staked-in-needle pre-filled syringes. Temperature and pressure cycling were found to contribute to additional liquid formation, with pressure changes playing a significant role. The detailed morphological analysis enhanced the understanding of microstructural arrangements within the needle. This research contributes to addressing the issue of needle clogging and can guide the development of strategies to improve pre-filled syringe performance.

Quantifying KRAS G12C Covalent Drug Inhibitor Activity in Mouse Tumors Using Mass Spectrometry.

The growing opportunities recognized for covalent drug inhibitors, like KRAS G12C inhibitors, are driving the need for mass spectrometry methods that can quickly and robustly measure therapeutic drug activity in vivo for drug discovery research and development. Effective front-end sample preparation is critical for proteins extracted from tumors but is generally labor intensive and impractical for large sample numbers typical in pharmacodynamic (PD) studies. Herein, we describe an automated and integrated sample preparation method for the measurement of activity levels of KRAS G12C drug inhibitor alkylation from complex tumor samples involving high throughput detergent removal and preconcentration followed by quantitation using mass spectrometry. We introduce a robust assay with an average intra-assay coefficient of variation (CV) of 4% and an interassay CV of 6% obtained from seven studies, enabling us to understand the relationship between KRAS G12C target occupancy and the therapeutic PD effect from mouse tumor samples. Further, the data demonstrated that the drug candidate GDC-6036, a KRAS G12C covalent inhibitor, shows dose-dependent target inhibition (KRAS G12C alkylation) and MAPK pathway inhibition, which correlate with high antitumor potency in the MIA PaCa-2 pancreatic xenograft model.

USP7 small-molecule inhibitors interfere with ubiquitin binding.

The ubiquitin system regulates essential cellular processes in eukaryotes. Ubiquitin is ligated to substrate proteins as monomers or chains and the topology of ubiquitin modifications regulates substrate interactions with specific proteins. Thus ubiquitination directs a variety of substrate fates including proteasomal degradation. Deubiquitinase enzymes cleave ubiquitin from substrates and are implicated in disease; for example, ubiquitin-specific protease-7 (USP7) regulates stability of the p53 tumour suppressor and other proteins critical for tumour cell survival. However, developing selective deubiquitinase inhibitors has been challenging and no co-crystal structures have been solved with small-molecule inhibitors. Here, using nuclear magnetic resonance-based screening and structure-based design, we describe the development of selective USP7 inhibitors GNE-6640 and GNE-6776. These compounds induce tumour cell death and enhance cytotoxicity with chemotherapeutic agents and targeted compounds, including PIM kinase inhibitors. Structural studies reveal that GNE-6640 and GNE-6776 non-covalently target USP7 12 Å distant from the catalytic cysteine. The compounds attenuate ubiquitin binding and thus inhibit USP7 deubiquitinase activity. GNE-6640 and GNE-6776 interact with acidic residues that mediate hydrogen-bond interactions with the ubiquitin Lys48 side chain, suggesting that USP7 preferentially interacts with and cleaves ubiquitin moieties that have free Lys48 side chains. We investigated this idea by engineering di-ubiquitin chains containing differential proximal and distal isotopic labels and measuring USP7 binding by nuclear magnetic resonance. This preferential binding protracted the depolymerization kinetics of Lys48-linked ubiquitin chains relative to Lys63-linked chains. In summary, engineering compounds that inhibit USP7 activity by attenuating ubiquitin binding suggests opportunities for developing other deubiquitinase inhibitors and may be a strategy more broadly applicable to inhibiting proteins that require ubiquitin binding for full functional activity.

MiDAS-Meaningful Immunogenetic Data at Scale.

Human immunogenetic variation in the form of HLA and KIR types has been shown to be strongly associated with a multitude of immune-related phenotypes. However, association studies involving immunogenetic loci most commonly involve simple analyses of classical HLA allelic diversity, resulting in limitations regarding the interpretability and reproducibility of results. We here present MiDAS, a comprehensive R package for immunogenetic data transformation and statistical analysis. MiDAS recodes input data in the form of HLA alleles and KIR types into biologically meaningful variables, allowing HLA amino acid fine mapping, analyses of HLA evolutionary divergence as well as experimentally validated HLA-KIR interactions. Further, MiDAS enables comprehensive statistical association analysis workflows with phenotypes of diverse measurement scales. MiDAS thus closes the gap between the inference of immunogenetic variation and its efficient utilization to make relevant discoveries related to immune and disease biology. It is freely available under a MIT license.

Characterization of Hepatic UDP-Glucuronosyltransferase Enzyme Abundance-Activity Correlations and Population Variability Using a Proteomics Approach and Comparison with Cytochrome P450 Enzymes.

The expression of ten major drug-metabolizing UDP-glucuronosyltransferase (UGT) enzymes in a panel of 130 human hepatic microsomal samples was measured using a liquid chromatography-tandem mass spectrometry-based approach. Simultaneously, ten cytochromes P450 and P450 reductase were also measured, and activity-expression relationships were assessed for comparison. The resulting data sets demonstrated that, with the exception of UGT2B17, 10th to 90th percentiles of UGT expression spanned 3- to 8-fold ranges. These ranges were small relative to ranges of reported mean UGT enzyme expression across different laboratories. We tested correlation of UGT expression with enzymatic activities using selective probe substrates. A high degree of abundance-activity correlation (Spearman's rank correlation coefficient > 0.6) was observed for UGT1As (1A1, 3, 4, 6) and cytochromes P450. In contrast, protein abundance and activity did not correlate strongly for UGT1A9 and UGT2B enzymes (2B4, 7, 10, 15, and 17). Protein abundance was strongly correlated for UGTs 2B7, 2B10, and 2B15. We suggest a number of factors may contribute to these differences including incomplete selectivity of probe substrates, correlated expression of these UGT2B isoforms, and the impact of splice and polymorphic variants on the peptides used in proteomics analysis, and exemplify this in the case of UGT2B10. Extensive correlation analyses identified important criteria for validating the fidelity of proteomics and enzymatic activity approaches for assessing UGT variability, population differences, and ontogenetic changes. SIGNIFICANCE STATEMENT: Protein expression data allow detailed assessment of interindividual variability and enzyme ontogeny. This study has observed that expression and enzyme activity are well correlated for hepatic UGT1A enzymes and cytochromes P450. However, for the UGT2B family, caution is advised when assuming correlation of expression and activity as is often done in physiologically based pharmacokinetic modeling. This can be due to incomplete probe substrate specificities, but may also be related to presence of inactive UGT protein materials and the effect of splicing variations.

Germline genetic polymorphisms influence tumor gene expression and immune cell infiltration.

Cancer immunotherapy has emerged as an effective therapy in a variety of cancers. However, a key challenge in the field is that only a subset of patients who receive immunotherapy exhibit durable response. It has been hypothesized that host genetics influences the inherent immune profiles of patients and may underlie their differential response to immunotherapy. Herein, we systematically determined the association of common germline genetic variants with gene expression and immune cell infiltration of the tumor. We identified 64,094 expression quantitative trait loci (eQTLs) that associated with 18,210 genes (eGenes) across 24 human cancers. Overall, eGenes were enriched for their being involved in immune processes, suggesting that expression of immune genes can be shaped by hereditary genetic variants. We identified the endoplasmic reticulum aminopeptidase 2 (ERAP2) gene as a pan-cancer type eGene whose expression levels stratified overall survival in a subset of patients with bladder cancer receiving anti-PD-L1 (atezolizumab) therapy. Finally, we identified 103 gene signature QTLs (gsQTLs) that were associated with predicted immune cell abundance within the tumor microenvironment. Our findings highlight the impact of germline SNPs on cancer-immune phenotypes and response to therapy; and these analyses provide a resource for integration of germline genetics as a component of personalized cancer immunotherapy.

Characterizing and Quantitating Therapeutic Tethered Multimeric Antibody Degradation Using Affinity Capture Mass Spectrometry.

There are many pharmacokinetic challenges associated with administering protein therapeutics, including biotransformation via clipping, deamidation, isomerization, oxidation, etc. In the case of engineered multivalent tethered antibody formats, proteolysis or deconjugation at the fusion or conjugation site present further issues. Unlike degradations associated with antibody drug conjugates, such biotransformations of tethered antibody formats usually result in degraded products with large mass differences. These large differences can result in processing or mass spectrometry response bias among the resulting product species that can lead to inaccurate stability quantitation. Herein, we describe an assay strategy for characterizing and quantitating degradations accurately for multivalent antibodies by incorporating response bias corrections. For the multivalent tethered antibody molecules selected, an ∼30-80% difference in response, compared to the cleaved product, was observed. To correct for the response bias, selected tethered multivalent antibodies and an IgG antibody (representing the stable intact and the degraded product species, respectively) were spiked in serum at known ratios for analysis. Following affinity capture, we generated calibration curves (five-parameter logistic fit p < 0.05) by plotting the measured ratios of the MS ion responses against the known spiked-in ratios (CVs < 8% for calibration standards). The qualified calibration curve (accuracy within 8% and 2% for measuring degradations of 5% and 15% product, respectively) was then used, through interpolation, to determine stability profiles for the same multivalent tethered antibody formats from both in vitro serum and pharmacokinetic study samples.

Interpreting In Vitro Release Performance from Long-Acting Parenteral Nanosuspensions Using USP-4 Dissolution and Spectroscopic Techniques.

Injectable sustained release dosage forms have emerged as desirable therapeutic routes for patients that require life-long treatments. The prevalence of drug molecules with low aqueous solubility and bioavailability has added momentum toward the development of suspension-based long-acting parenteral (LAP) formulations; the previously undesirable physicochemical properties of Biopharmaceutics Classification System (BCS) Class II/IV compounds are best suited for extended release applications. Effective in vitro release (IVR) testing of crystalline suspensions affirms product quality during early-stage development and provides connections with in vivo performance. However, before in vitro-in vivo correlations (IVIVCs) can be established, it is necessary to evaluate formulation attributes that directly affect IVR properties. In this work, a series of crystalline LAP nanosuspensions were formulated with different stabilizing polymers and applied to a continuous flow-through (USP-4) dissolution method. This technique confirmed the role of salt effects on the stability of polymer-coated nanoparticles through the detection of disparate active pharmaceutical ingredient (API) release profiles. The polymer stabilizers with extended hydrophilic chains exhibited elevated intrapolymer activity from the loss of hydrogen-bond cushioning in dissolution media with heightened ionic strength, confirmed through one-dimensional (1D) 1H NMR and two-dimensional nuclear Overhauser effect spectroscopy (2D NOESY) experiments. Thus, steric repulsion within the affected nanosuspensions was limited and release rates decreased. Additionally, the strength of interaction between hydrophobic polymer components and the API crystalline surface contributed to suspension dissolution properties, confirmed through solution- and solid-state spectroscopic analyses. This study provides a unique perspective on the dynamic interface between the crystalline drug and aqueous microenvironment during dissolution.

Probing in Vitro Release Kinetics of Long-Acting Injectable Nanosuspensions via Flow-NMR Spectroscopy.

Novel treatment routes are emerging for an array of diseases and afflictions. Complex dosage forms, based on active pharmaceutical ingredients (APIs) with previously undesirable physicochemical characteristics, are becoming mainstream and actively pursued in various pipeline initiatives. To fundamentally understand how constituents in these dosage forms interact on a molecular level, analytical methods need to be developed that encompass selectivity and sensitivity requirements previously reserved for a myriad of in vitro techniques. The knowledge of precise chemical interactions between drugs and excipients in a dosage form can streamline formulation development and process screening capabilities through the identification of properties that influence rates and mechanisms of drug release in a cost-effective manner, relative to long-term in vivo studies. Through this work, a noncompendial in vitro release (IVR) method was developed that distinguished the presence of individual components in a complex crystalline nanosuspension environment. Doravirine was formulated as a series of long-acting injectable nanosuspensions with assorted excipients, using low- and high-energy wet media milling methods. IVR behavior of all formulation components were monitored using a robust continuous flow-through (CFT) dissolution setup (USP-4 apparatus) with on-line 1H NMR end-analysis (flow-NMR). Results from this investigation led to a better understanding of formulation parameter influences on nanosuspension stability, surface chemistry, and dissolution behavior. Flow-NMR can be applied to a broad range of dosage forms in which specific molecular interactions from the solution microenvironment require further insight to enhance product development capabilities.

Direct Drug Analysis in Polymeric Implants Using Desorption Electrospray Ionization - Mass Spectrometry Imaging (DESI-MSI).

PURPOSE: Desorption electrospray ionization mass spectrometry imaging (DESI-MSI) coupled with gas-phase ion mobility spectrometry was used to characterize the drug distribution in polymeric implants before and after exposure to accelerated in vitro release (IVR) media. DESI-MSI provides definitive chemical identification and localization of formulation components, including 2D chemical mapping of individual components with essentially no sample preparation. METHODS: Polymeric implants containing 40% (w/w) entecavir and poly(D,L-lactide) (PLA) were prepared and then exposed to either acidified PBS (pH 2.5) or MeOH:H2O (50:50, v/v) medias during a 7-day IVR test using continuous flow-through (CFT) cell dissolution. The amount of drug released from the polymer matrix during the 7-day IVR test was monitored by online-ultraviolet spectroscopy (UV) and HPLC-UV. After that period, intact implants and radial sections of implants were analyzed by DESI-MSI with ion mobility spectrometry. The active ingredient along with impurities and contaminants were used to generate chemical maps before and after exposure to the release medias. RESULTS: Bi-phasic release profiles were observed for implants during IVR release using both medias. During the second phase of release, implants exposed to PBS, pH 2.5, released the entecavir faster than the implants exposed to MeOH:H2O (50:50, v/v). Radial images of the polymer interior show that entecavir is localized along the central core of the implant after exposure to MeOH:H2O (50:50, v/v) and that the drug is more uniformly distributed throughout the implant after exposure to acidified PBS (pH 2.5). CONCLUSIONS: DESI-MSI coupled with ion mobility analysis produced chemical images of the drug distribution on the exterior and interior of cylindrical polymeric implants before and after exposure to various release medias. These results demonstrated the utility of this technique for rapid characterization of drug and impurity/degradant distribution within polymeric implants with direct implications for formulation development as well as analytical method development activities for various solid parenteral and oral dosage forms. These results are especially meaningful since samples were analyzed with essentially no preparative procedures.

Previously reported placebo-response-associated variants do not predict patient outcomes in inflammatory disease Phase III trial placebo arms.

In clinical trials, a placebo response refers to improvement in disease symptoms arising from the psychological effect of receiving a treatment rather than the actual treatment under investigation. Previous research has reported genomic variation associated with the likelihood of observing a placebo response, but these studies have been limited in scope and have not been validated. Here, we analyzed whole-genome sequencing data from 784 patients undergoing placebo treatment in Phase III Asthma or Rheumatoid Arthritis trials to assess the impact of previously reported variation on patient outcomes in the placebo arms and to identify novel variants associated with the placebo response. Contrary to expectations based on previous reports, we did not observe any statistically significant associations between genomic variants and placebo treatment outcome. Our findings suggest that the biological origin of the placebo response is complex and likely to be variable between disease areas.

Dual functionality of ß-tryptase protomers as both proteases and cofactors in the active tetramer.

Human ß-tryptase, a tetrameric trypsin-like serine protease, is an important mediator of the allergic inflammatory responses in asthma. During acute hypersensitivity reactions, mast cells degranulate, releasing active tetramer as a complex with proteoglycans. Extensive efforts have focused on developing therapeutic ß-tryptase inhibitors, but its unique activation mechanism is less well-explored. Tryptase is active only after proteolytic removal of the pro-domain followed by tetramer formation via two distinct symmetry-related interfaces. We show that the cleaved I16G mutant cannot tetramerize, likely due to impaired insertion of its N terminus into its "activation pocket," indicating allosteric linkage at multiple sites on each protomer. We engineered cysteines into each of the two distinct interfaces (Y75C for small or I99C for large) to assess the activity of each tetramer and disulfide-locked dimer. Using size-exclusion chromatography and enzymatic assays, we demonstrate that the two large tetramer interfaces regulate enzymatic activity, elucidating the importance of this protein-protein interaction for allosteric regulation. Notably, the I99C large interface dimer is active, even in the absence of heparin. We show that a monomeric ß-tryptase mutant (I99C*/Y75A/Y37bA, where C* is cysteinylated Cys-99) cannot form a dimer or tetramer, yet it is active but only in the presence of heparin. Thus heparin both stabilizes the tetramer and allosterically conditions the active site. We hypothesize that each ß-tryptase protomer in the tetramer has two distinct roles, acting both as a protease and as a cofactor for its neighboring protomer, to allosterically regulate enzymatic activity, providing a rationale for direct correlation of tetramer stability with proteolytic activity.

Uncovering a Dual Regulatory Role for Caspases During Endoplasmic Reticulum Stress-induced Cell Death.

Many diseases are associated with endoplasmic reticulum (ER) stress, which results from an accumulation of misfolded proteins. This triggers an adaptive response called the "unfolded protein response" (UPR), and prolonged exposure to ER stress leads to cell death. Caspases are reported to play a critical role in ER stress-induced cell death but the underlying mechanisms by which they exert their effect continue to remain elusive. To understand the role caspases play during ER stress, a systems level approach integrating analysis of the transcriptome, proteome, and proteolytic substrate profile was employed. This quantitative analysis revealed transcriptional profiles for most human genes, provided information on protein abundance for 4476 proteins, and identified 445 caspase substrates. Based on these data sets many caspase substrates were shown to be downregulated at the protein level during ER stress suggesting caspase activity inhibits their cellular function. Additionally, RNA sequencing revealed a role for caspases in regulation of ER stress-induced transcriptional pathways and gene set enrichment analysis showed expression of multiple gene targets of essential transcription factors to be upregulated during ER stress upon inhibition of caspases. Furthermore, these transcription factors were degraded in a caspase-dependent manner during ER stress. These results indicate that caspases play a dual role in regulating the cellular response to ER stress through both post-translational and transcriptional regulatory mechanisms. Moreover, this study provides unique insight into progression of the unfolded protein response into cell death, which may help identify therapeutic strategies to treat ER stress-related diseases.

USP Apparatus 4: a Valuable In Vitro Tool to Enable Formulation Development of Long-Acting Parenteral (LAP) Nanosuspension Formulations of Poorly Water-Soluble Compounds.

Long-acting or extended release parenteral dosage forms have attracted extensive attention due to their ability to maintain therapeutic drug concentrations over long periods of time and reduce administration frequency, thus improving patient compliance. It is essential to have an in vitro release (IVR) testing method that can be used to assure product quality during routine production as well as predict and understand the in vivo performance of a formulation. The purpose of this work was to develop a discriminatory in vitro release method to guide formulation and process development of long-acting parenteral (LAP) nanosuspension formulations composed of poorly water-soluble drugs (BCS class II). Injectable nanosuspension formulations were developed to serve as test articles for method development. Several different IVR methods were evaluated for their application to the formulation screening and process development including (1) USP apparatus 2, (2) dialysis and reverse dialysis sac, and (3) continuous flow-through cell (USP apparatus 4). Preliminary data shows the promising results to support the utilization of USP 4 over more widely accepted USP 2 and dialysis methods. A combination of more representative in vivo hydrodynamics and ease of maintaining sink conditions yields the USP 4 flow-through cell method a more suitable in vitro release method for nanosuspension-based LAP formulations of poorly water-soluble compounds, such as compounds A and B.

Identifying and mitigating batch effects in whole genome sequencing data.

BACKGROUND: Large sample sets of whole genome sequencing with deep coverage are being generated, however assembling datasets from different sources inevitably introduces batch effects. These batch effects are not well understood and can be due to changes in the sequencing protocol or bioinformatics tools used to process the data. No systematic algorithms or heuristics exist to detect and filter batch effects or remove associations impacted by batch effects in whole genome sequencing data. RESULTS: We describe key quality metrics, provide a freely available software package to compute them, and demonstrate that identification of batch effects is aided by principal components analysis of these metrics. To mitigate batch effects, we developed new site-specific filters that identified and removed variants that falsely associated with the phenotype due to batch effect. These include filtering based on: a haplotype based genotype correction, a differential genotype quality test, and removing sites with missing genotype rate greater than 30% after setting genotypes with quality scores less than 20 to missing. This method removed 96.1% of unconfirmed genome-wide significant SNP associations and 97.6% of unconfirmed genome-wide significant indel associations. We performed analyses to demonstrate that: 1) These filters impacted variants known to be disease associated as 2 out of 16 confirmed associations in an AMD candidate SNP analysis were filtered, representing a reduction in power of 12.5%, 2) In the absence of batch effects, these filters removed only a small proportion of variants across the genome (type I error rate of 3%), and 3) in an independent dataset, the method removed 90.2% of unconfirmed genome-wide SNP associations and 89.8% of unconfirmed genome-wide indel associations. CONCLUSIONS: Researchers currently do not have effective tools to identify and mitigate batch effects in whole genome sequencing data. We developed and validated methods and filters to address this deficiency.

Hominoid-specific enzyme GLUD2 promotes growth of IDH1R132H glioma.

Somatic mutation of isocitrate dehydrogenase 1 (IDH1) is now recognized as the most common initiating event for secondary glioblastoma, a brain tumor type arising with high frequency in the frontal lobe. A puzzling feature of IDH1 mutation is the selective manifestation of glioma as the only neoplasm frequently associated with early postzygotic occurrence of this genomic alteration. We report here that IDH1(R132H) exhibits a growth-inhibitory effect that is abrogated in the presence of glutamate dehydrogenase 2 (GLUD2), a hominoid-specific enzyme purportedly optimized to facilitate glutamate turnover in human forebrain. Using murine glioma progenitor cells, we demonstrate that IDH1(R132H) exerts a growth-inhibitory effect that is paralleled by deficiency in metabolic flux from glucose and glutamine to lipids. Examining human gliomas, we find that glutamate dehydrogenase 1 (GLUD1) and GLUD2 are overexpressed in IDH1-mutant tumors and that orthotopic growth of an IDH1-mutant glioma line is inhibited by knockdown of GLUD1/2. Strikingly, introduction of GLUD2 into murine glioma progenitor cells reverses deleterious effects of IDH1 mutation on metabolic flux and tumor growth. Further, we report that glutamate, a substrate of GLUD2 and a neurotransmitter abundant in mammalian neocortex, can support growth of glioma progenitor cells irrespective of IDH1 mutation status. These findings suggest that specialization of human neocortex for high glutamate neurotransmitter flux creates a metabolic niche conducive to growth of IDH1 mutant tumors.

Conformational dynamics control ubiquitin-deubiquitinase interactions and influence in vivo signaling.

Ubiquitin is a highly conserved eukaryotic protein that interacts with a diverse set of partners to act as a cellular signaling hub. Ubiquitin's conformational flexibility has been postulated to underlie its multifaceted recognition. Here we use computational and library-based means to interrogate core mutations that modulate the conformational dynamics of human ubiquitin. These ubiquitin variants exhibit increased affinity for the USP14 deubiquitinase, with concomitantly reduced affinity for other deubiquitinases. Strikingly, the kinetics of conformational motion are dramatically slowed in these variants without a detectable change in either the ground state fold or excited state population. These variants can be ligated into substrate-linked chains in vitro and in vivo but cannot solely support growth in eukaryotic cells. Proteomic analyses reveal nearly identical interaction profiles between WT ubiquitin and the variants but identify a small subset of altered interactions. Taken together, these results show that conformational dynamics are critical for ubiquitin-deubiquitinase interactions and imply that the fine tuning of motion has played a key role in the evolution of ubiquitin as a signaling hub.

Small molecule inhibition of group I p21-activated kinases in breast cancer induces apoptosis and potentiates the activity of microtubule stabilizing agents.

INTRODUCTION: Breast cancer, the most common cause of cancer-related deaths worldwide among women, is a molecularly and clinically heterogeneous disease. Extensive genetic and epigenetic profiling of breast tumors has recently revealed novel putative driver genes, including p21-activated kinase (PAK)1. PAK1 is a serine/threonine kinase downstream of small GTP-binding proteins, Rac1 and Cdc42, and is an integral component of growth factor signaling networks and cellular functions fundamental to tumorigenesis. METHODS: PAK1 dysregulation (copy number gain, mRNA and protein expression) was evaluated in two cohorts of breast cancer tissues (n=980 and 1,108). A novel small molecule inhibitor, FRAX1036, and RNA interference were used to examine PAK1 loss of function and combination with docetaxel in vitro. Mechanism of action for the therapeutic combination, both cellular and molecular, was assessed via time-lapse microscopy and immunoblotting. RESULTS: We demonstrate that focal genomic amplification and overexpression of PAK1 are associated with poor clinical outcome in the luminal subtype of breast cancer (P=1.29×10(-4) and P=0.015, respectively). Given the role for PAK1 in regulating cytoskeletal organization, we hypothesized that combination of PAK1 inhibition with taxane treatment could be combined to further interfere with microtubule dynamics and cell survival. Consistent with this, administration of docetaxel with either a novel small molecule inhibitor of group I PAKs, FRAX1036, or PAK1 small interfering RNA oligonucleotides dramatically altered signaling to cytoskeletal-associated proteins, such as stathmin, and induced microtubule disorganization and cellular apoptosis. Live-cell imaging revealed that the duration of mitotic arrest mediated by docetaxel was significantly reduced in the presence of FRAX1036, and this was associated with increased kinetics of apoptosis. CONCLUSIONS: Taken together, these findings further support PAK1 as a potential target in breast cancer and suggest combination with taxanes as a viable strategy to increase anti-tumor efficacy.