Faculty Bibliography

Unbiased assays such as shotgun proteomics and RNA-seq provide high-resolution molecular characterization of tumors. These assays measure molecules with highly varied distributions, making interpretation and hypothesis testing challenging. Samples with the most extreme measurements for a molecule can reveal the most interesting biological insights yet are often excluded from analysis. Furthermore, rare disease subtypes are, by definition, underrepresented in cancer cohorts. To provide a strategy for identifying molecules aberrantly enriched in small sample cohorts, we present BlackSheep, a package for nonparametric description and differential analysis of genome-wide data, available from Bioconductor (https://www.bioconductor.org/packages/release/bioc/html/blacksheepr.html) and Bioconda (https://bioconda.github.io/recipes/blksheep/README.html). BlackSheep is a complementary tool to other differential expression analysis methods, which is particularly useful when analyzing small subgroups in a larger cohort.

Introduction: Although a substantial body of research documents a relationship between traumatic stress in childhood and the initiation of substance use later in the life course, only limited research has examined potential linkages between adverse childhood experiences (ACEs) and the initiation of non-medical prescription opioid use and other opioid use behaviors. The present study contributes to this growing body of work by investigating the association of childhood trauma with early initiation of a series of opioid use behaviors. Methods: New York City young adults (n = 539) ages 18-29 who reported non-medical use of prescription opioids or heroin use in the past 30 days were recruited using Respondent-Driven Sampling in 2014-16. Ten ACEs were assessed via self-report with the ACE Questionnaire. Associations between number of ACEs and self-reported ages of initiating seven opioid use behaviors (e.g., non-medical prescription opioid use, heroin use, heroin injection) were estimated with multivariable logistic regression. Results: Eighty nine percent of participants reported at least one ACE, and 46% reported four or more ACEs, a well-supported threshold indicating elevated risk for negative health consequences. Every increase of one trauma was associated with a 12-23% increase in odds of early initiation across the seven opioid use behaviors. Findings also document that the mean age at initiation increased with increasing risk severity across the behaviors, contributing to evidence of a trajectory from opioid pill misuse to opioid injection. Discussion: Increasing number of childhood traumas was associated with increased odds of earlier initiation of multiple opioid misuse behaviors. In light of prior research linking earlier initiation of substance use with increased substance use severity, present findings suggest the importance of ACEs as individual-level determinants of increased opioid use severity. Efforts to prevent onset and escalation of opioid use among at-risk youth may benefit from trauma prevention programs and trauma-focused screening and treatment, as well as increased attention to ameliorating upstream socio-structural drivers of childhood trauma.

Advances in mass-spectrometry have generated increasingly large-scale proteomics datasets containing tens of thousands of phosphorylation sites (phosphosites) that require prioritization. We develop a bioinformatics tool called HotPho and systematically discover 3D co-clustering of phosphosites and cancer mutations on protein structures. HotPho identifies 474 such hybrid clusters containing 1255 co-clustering phosphosites, including RET p.S904/Y928, the conserved HRAS/KRAS p.Y96, and IDH1 p.Y139/IDH2 p.Y179 that are adjacent to recurrent mutations on protein structures not found by linear proximity approaches. Hybrid clusters, enriched in histone and kinase domains, frequently include expression-associated mutations experimentally shown as activating and conferring genetic dependency. Approximately 300 co-clustering phosphosites are verified in patient samples of 5 cancer types or previously implicated in cancer, including CTNNB1 p.S29/Y30, EGFR p.S720, MAPK1 p.S142, and PTPN12 p.S275. In summary, systematic 3D clustering analysis highlights nearly 3,000 likely functional mutations and over 1000 cancer phosphosites for downstream investigation and evaluation of potential clinical relevance.

Within the last decade, numerous studies have demonstrated changes in the gut microbiome associated with specific autoimmune diseases. Due to differences in study design, data quality control, analysis and statistical methods, many results of these studies are inconsistent and incomparable. To better understand the relationship between the intestinal microbiome and autoimmunity, we have completed a comprehensive re-analysis of 42 studies focusing on the gut microbiome in 12 autoimmune diseases to identify a microbial signature predictive of multiple sclerosis (MS), inflammatory bowel disease (IBD), rheumatoid arthritis (RA) and general autoimmune disease using both 16S rRNA sequencing data and shotgun metagenomics data. To do this, we used four machine learning algorithms, random forest, eXtreme Gradient Boosting (XGBoost), ridge regression, and support vector machine with radial kernel and recursive feature elimination to rank disease predictive taxa comparing disease vs. healthy participants and pairwise comparisons of each disease. Comparing the performance of these models, we found the two tree-based methods, XGBoost and random forest, most capable of handling sparse multidimensional data, to consistently produce the best results. Through this modeling, we identified a number of taxa consistently identified as dysregulated in a general autoimmune disease model including Odoribacter, Lachnospiraceae Clostridium, and Mogibacteriaceae implicating all as potential factors connecting the gut microbiome to autoimmune response. Further, we computed pairwise comparison models to identify disease specific taxa signatures highlighting a role for Peptostreptococcaceae and Ruminococcaceae Gemmiger in IBD and Akkermansia, Butyricicoccus, and Mogibacteriaceae in MS. We then connected a subset of these taxa with potential metabolic alterations based on metagenomic/metabolomic correlation analysis, identifying 215 metabolites associated with autoimmunity-predictive taxa.

Over-accumulation of oxalate in humans may lead to nephrolithiasis and nephrocalcinosis. Humans lack endogenous oxalate degradation pathways (ODP), but intestinal microbes can degrade oxalate using multiple ODPs and protect against its absorption. The exact oxalate-degrading taxa in the human microbiota and their ODP have not been described. We leverage multi-omics data (>3000 samples from >1000 subjects) to show that the human microbiota primarily uses the type II ODP, rather than type I. Further, among the diverse ODP-encoding microbes, an oxalate autotroph, Oxalobacter formigenes, dominates this function transcriptionally. Patients with Inflammatory Bowel Disease (IBD) frequently suffer from disrupted oxalate homeostasis and calcium oxalate nephrolithiasis. We show that the enteric oxalate level is elevated in IBD patients, with highest levels in Crohn's disease patients with both ileal and colonic involvement consistent with known nephrolithiasis risk. We show that the microbiota ODP expression is reduced in IBD patients, which may contribute to the disrupted oxalate homeostasis. The specific changes in ODP expression by several important taxa suggest that they play distinct roles in IBD-induced nephrolithiasis risk. Lastly, we colonize mice that are maintained in the gnotobiotic facility with O. formigenes, using either a laboratory isolate or an isolate we cultured from human stools, and observed a significant reduction in host fecal and urine oxalate levels, supporting our in silico prediction of the importance of the microbiome, particularly O. formigenes in host oxalate homeostasis.

OBJECTIVE:Hydroxychloroquine (HCQ) is a mainstay of therapy in the treatment of SLE. The effect of HCQ on platelets and vascular health is uncertain. We investigated the relationship between HCQ use and dose with platelet activity, platelet transcriptomics and vascular health in patients with SLE. METHODS:Platelet aggregation, platelet mRNA expression and vascular health (sublingual capillary perfused boundary region (PBR), red blood cell filling (RBCF) and brachial artery reactivity testing) were analysed by HCQ use and dose. RESULTS:Among 132 subjects with SLE (age: 39.7±12.9 years, 97% female), 108 were on HCQ. SLE disease activity was similar between subjects on and off HCQ. Platelet aggregation in response to multiple agonists was significantly lower in patients on HCQ. There were inverse relationships between HCQ dose and gene expression pathways of platelet activity. Gene expression of P-selectin (SELP) was inversely correlated with HCQ dose (r=-0.41, p=0.003), which was validated at the protein level. Subjects on HCQ had improved vascular function correlating with HCQ dose as measured by lower PBR (r=-0.52, p=0.007), higher RBCF (r=0.55, p=0.004) and greater brachial artery reactivity (r=0.43, p=0.056). CONCLUSION/CONCLUSIONS:HCQ use was associated with decreased platelet activation and activation-related transcripts and improved vascular health in SLE.

Atherosclerosis is a disease of chronic inflammation. We investigated the roles of the cytokines IL-4 and IL-13, the classical activators of STAT6, in the resolution of atherosclerosis inflammation. Using Il4^-/-Il13^-/- mice, resolution was impaired, and in control mice, in both progressing and resolving plaques, levels of IL-4 were stably low, and IL-13 was undetectable. This suggested that IL-4 is required for atherosclerosis resolution, but collaborates with other factors. We had observed increased Wnt signaling in macrophages in resolving plaques, and human genetic data from others showed that a loss-of-function Wnt mutation was associated with premature atherosclerosis. We now find an inverse association between activation of Wnt signaling and disease severity in mice and humans. Wnt enhanced the expression of inflammation resolving factors after treatment with plaque-relevant low concentrations of IL-4. Mechanistically, activation of the Wnt pathway following lipid lowering potentiates IL-4 responsiveness in macrophages via a PGE₂/STAT3 axis.

OBJECTIVE:<0.05, |log2foldchange| >0.5) and analyzed using weighted gene co-expression network analysis. Weighted gene co-expression network analysis revealed blood modules enriched for immune activation, secretory granules, and coagulation in patients with PAD. Of these 127 differentially expressed transcripts, 40 were significantly associated with MACLE (log-rank false discovery rate <0.1). MicroRNA-4477b was significantly increased in patients with PAD with subsequent MACLE and in a mouse hindlimb ischemia model. CONCLUSIONS:A whole blood transcript signature identified patients with symptomatic PAD and PAD patients at increased risk of MACLE. A previously uncharacterized transcript microRNA-4477b was overexpressed in prevalent PAD, incident MACLE, and in a mouse hindlimb ischemia model. Our novel transcriptomic signature provides insight into potential mechanisms of patients with severe symptomatic PAD.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the etiological agent of Coronavirus Disease 2019 (COVID-19). There is a dire need for novel effective antivirals to treat COVID-19, as the only approved direct-acting antiviral to date is remdesivir, targeting the viral polymerase complex. A potential alternate target in the viral life cycle is the main SARS-CoV-2 protease 3CL^pro (M^pro). The drug candidate PF-00835231 is the active compound of the first anti-3CL^pro regimen in clinical trials. Here, we perform a comparative analysis of PF-00835231, the pre-clinical 3CL^pro inhibitor GC-376, and the polymerase inhibitor remdesivir, in alveolar basal epithelial cells modified to express ACE2 (A549^+ACE2 cells). We find PF-00835231 with at least similar or higher potency than remdesivir or GC-376. A time-of-drug-addition approach delineates the timing of early SARS-CoV-2 life cycle steps in A549^+ACE2 cells and validates PF-00835231's early time of action. In a model of the human polarized airway epithelium, both PF-00835231 and remdesivir potently inhibit SARS-CoV-2 at low micromolar concentrations. Finally, we show that the efflux transporter P-glycoprotein, which was previously suggested to diminish PF-00835231's efficacy based on experiments in monkey kidney Vero E6 cells, does not negatively impact PF-00835231 efficacy in either A549^+ACE2 cells or human polarized airway epithelial cultures. Thus, our study provides in vitro evidence for the potential of PF-00835231 as an effective SARS-CoV-2 antiviral and addresses concerns that emerged based on prior studies in non-human in vitro models.Importance:The arsenal of SARS-CoV-2 specific antiviral drugs is extremely limited. Only one direct-acting antiviral drug is currently approved, the viral polymerase inhibitor remdesivir, and it has limited efficacy. Thus, there is a substantial need to develop additional antiviral compounds with minimal side effects and alternate viral targets. One such alternate target is its main protease, 3CL^pro (M^pro), an essential component of the SARS-CoV-2 life cycle processing the viral polyprotein into the components of the viral polymerase complex. In this study, we characterize a novel antiviral drug, PF-00835231, which is the active component of the first-in-class 3CL^pro-targeting regimen in clinical trials. Using 3D in vitro models of the human airway epithelium, we demonstrate the antiviral potential of PF-00835231 for inhibition of SARS-CoV-2.

Glioblastoma (GBM) is the most aggressive nervous system cancer. Understanding its molecular pathogenesis is crucial to improving diagnosis and treatment. Integrated analysis of genomic, proteomic, post-translational modification and metabolomic data on 99 treatment-naive GBMs provides insights to GBM biology. We identify key phosphorylation events (e.g., phosphorylated PTPN11 and PLCG1) as potential switches mediating oncogenic pathway activation, as well as potential targets for EGFR-, TP53-, and RB1-altered tumors. Immune subtypes with distinct immune cell types are discovered using bulk omics methodologies, validated by snRNA-seq, and correlated with specific expression and histone acetylation patterns. Histone H2B acetylation in classical-like and immune-low GBM is driven largely by BRDs, CREBBP, and EP300. Integrated metabolomic and proteomic data identify specific lipid distributions across subtypes and distinct global metabolic changes in IDH-mutated tumors. This work highlights biological relationships that could contribute to stratification of GBM patients for more effective treatment.