Toward a fine-scale population health monitoring system.

Understanding population health disparities is an essential component of equitable precision health efforts. Epidemiology research often relies on definitions of race and ethnicity, but these population labels may not adequately capture disease burdens and environmental factors impacting specific sub-populations. Here, we propose a framework for repurposing data from electronic health records (EHRs) in concert with genomic data to explore the demographic ties that can impact disease burdens. Using data from a diverse biobank in New York City, we identified 17 communities sharing recent genetic ancestry. We observed 1,177 health outcomes that were statistically associated with a specific group and demonstrated significant differences in the segregation of genetic variants contributing to Mendelian diseases. We also demonstrated that fine-scale population structure can impact the prediction of complex disease risk within groups. This work reinforces the utility of linking genomic data to EHRs and provides a framework toward fine-scale monitoring of population health.

Trained Immunity-Promoting Nanobiologic Therapy Suppresses Tumor Growth and Potentiates Checkpoint Inhibition.

Trained immunity, a functional state of myeloid cells, has been proposed as a compelling immune-oncological target. Its efficient induction requires direct engagement of myeloid progenitors in the bone marrow. For this purpose, we developed a bone marrow-avid nanobiologic platform designed specifically to induce trained immunity. We established the potent anti-tumor capabilities of our lead candidate MTP10-HDL in a B16F10 mouse melanoma model. These anti-tumor effects result from trained immunity-induced myelopoiesis caused by epigenetic rewiring of multipotent progenitors in the bone marrow, which overcomes the immunosuppressive tumor microenvironment. Furthermore, MTP10-HDL nanotherapy potentiates checkpoint inhibition in this melanoma model refractory to anti-PD-1 and anti-CTLA-4 therapy. Finally, we determined MTP10-HDL's favorable biodistribution and safety profile in non-human primates. In conclusion, we show that rationally designed nanobiologics can promote trained immunity and elicit a durable anti-tumor response either as a monotherapy or in combination with checkpoint inhibitor drugs.

Nuclear RNA catabolism controls endogenous retroviruses, gene expression asymmetry, and dedifferentiation.

Endogenous retroviruses (ERVs) are remnants of ancient parasitic infections and comprise sizable portions of most genomes. Although epigenetic mechanisms silence most ERVs by generating a repressive environment that prevents their expression (heterochromatin), little is known about mechanisms silencing ERVs residing in open regions of the genome (euchromatin). This is particularly important during embryonic development, where induction and repression of distinct classes of ERVs occur in short temporal windows. Here, we demonstrate that transcription-associated RNA degradation by the nuclear RNA exosome and Integrator is a regulatory mechanism that controls the productive transcription of most genes and many ERVs involved in preimplantation development. Disrupting nuclear RNA catabolism promotes dedifferentiation to a totipotent-like state characterized by defects in RNAPII elongation and decreased expression of long genes (gene-length asymmetry). Our results indicate that RNA catabolism is a core regulatory module of gene networks that safeguards RNAPII activity, ERV expression, cell identity, and developmental potency.

Clonally expanded CD8 T cells characterize amyotrophic lateral sclerosis-4.

Amyotrophic lateral sclerosis (ALS) is a heterogenous neurodegenerative disorder that affects motor neurons and voluntary muscle control1. ALS heterogeneity includes the age of manifestation, the rate of progression and the anatomical sites of symptom onset. Disease-causing mutations in specific genes have been identified and define different subtypes of ALS1. Although several ALS-associated genes have been shown to affect immune functions2, whether specific immune features account for ALS heterogeneity is poorly understood. Amyotrophic lateral sclerosis-4 (ALS4) is characterized by juvenile onset and slow progression3. Patients with ALS4 show motor difficulties by the time that they are in their thirties, and most of them require devices to assist with walking by their fifties. ALS4 is caused by mutations in the senataxin gene (SETX). Here, using Setx knock-in mice that carry the ALS4-causative L389S mutation, we describe an immunological signature that consists of clonally expanded, terminally differentiated effector memory (TEMRA) CD8 T cells in the central nervous system and the blood of knock-in mice. Increased frequencies of antigen-specific CD8 T cells in knock-in mice mirror the progression of motor neuron disease and correlate with anti-glioma immunity. Furthermore, bone marrow transplantation experiments indicate that the immune system has a key role in ALS4 neurodegeneration. In patients with ALS4, clonally expanded TEMRA CD8 T cells circulate in the peripheral blood. Our results provide evidence of an antigen-specific CD8 T cell response in ALS4, which could be used to unravel disease mechanisms and as a potential biomarker of disease state.

Leveraging health systems data to characterize a large effect variant conferring risk for liver disease in Puerto Ricans.

The integration of genomic data into health systems offers opportunities to identify genomic factors underlying the continuum of rare and common disease. We applied a population-scale haplotype association approach based on identity-by-descent (IBD) in a large multi-ethnic biobank to a spectrum of disease outcomes derived from electronic health records (EHRs) and uncovered a risk locus for liver disease. We used genome sequencing and in silico approaches to fine-map the signal to a non-coding variant (c.2784-12T>C) in the gene ABCB4. In vitro analysis confirmed the variant disrupted splicing of the ABCB4 pre-mRNA. Four of five homozygotes had evidence of advanced liver disease, and there was a significant association with liver disease among heterozygotes, suggesting the variant is linked to increased risk of liver disease in an allele dose-dependent manner. Population-level screening revealed the variant to be at a carrier rate of 1.95% in Puerto Rican individuals, likely as the result of a Puerto Rican founder effect. This work demonstrates that integrating EHR and genomic data at a population scale can facilitate strategies for understanding the continuum of genomic risk for common diseases, particularly in populations underrepresented in genomic medicine.

KEA3: improved kinase enrichment analysis via data integration.

Phosphoproteomics and proteomics experiments capture a global snapshot of the cellular signaling network, but these methods do not directly measure kinase state. Kinase Enrichment Analysis 3 (KEA3) is a webserver application that infers overrepresentation of upstream kinases whose putative substrates are in a user-inputted list of proteins. KEA3 can be applied to analyze data from phosphoproteomics and proteomics studies to predict the upstream kinases responsible for observed differential phosphorylations. The KEA3 background database contains measured and predicted kinase-substrate interactions (KSI), kinase-protein interactions (KPI), and interactions supported by co-expression and co-occurrence data. To benchmark the performance of KEA3, we examined whether KEA3 can predict the perturbed kinase from single-kinase perturbation followed by gene expression experiments, and phosphoproteomics data collected from kinase-targeting small molecules. We show that integrating KSIs and KPIs across data sources to produce a composite ranking improves the recovery of the expected kinase. The KEA3 webserver is available at https://maayanlab.cloud/kea3.

Transcriptomic profiling of recessive dystrophic epidermolysis bullosa wounded skin highlights drug repurposing opportunities to improve wound healing.

Chronic wounds present a major disease burden in people with recessive dystrophic epidermolysis bullosa (RDEB), an inherited blistering skin disorder caused by mutations in COL7A1 encoding type VII collagen, the major component of anchoring fibrils at the dermal-epidermal junction. Treatment of RDEB wounds is mostly symptomatic, and there is considerable unmet need in trying to improve and accelerate wound healing. In this study, we defined transcriptomic profiles and gene pathways in RDEB wounds and compared these to intact skin in RDEB and healthy control subjects. We then used a reverse transcriptomics approach to discover drugs or compounds, which might restore RDEB wound profiles towards intact skin. Differential expression analysis identified >2000 differences between RDEB wounds and intact skin, with RDEB wounds displaying aberrant cytokine-cytokine interactions, Toll-like receptor signalling, and JAK-STAT signalling pathways. In-silico prediction for compounds that reverse gene expression signatures highlighted methotrexate as a leading candidate. Overall, this study provides insight into the molecular profiles of RDEB wounds and underscores the possible clinical value of reverse transcriptomics data analysis in RDEB, and the potential of this approach in discovering or repurposing drugs for other diseases.

Geneshot: search engine for ranking genes from arbitrary text queries.

The frequency by which genes are studied correlates with the prior knowledge accumulated about them. This leads to an imbalance in research attention where some genes are highly investigated while others are ignored. Geneshot is a search engine developed to illuminate this gap and to promote attention to the under-studied genome. Through a simple web interface, Geneshot enables researchers to enter arbitrary search terms, to receive ranked lists of genes relevant to the search terms. Returned ranked gene lists contain genes that were previously published in association with the search terms, as well as genes predicted to be associated with the terms based on data integration from multiple sources. The search results are presented with interactive visualizations. To predict gene function, Geneshot utilizes gene-gene similarity matrices from processed RNA-seq data, or from gene-gene co-occurrence data obtained from multiple sources. In addition, Geneshot can be used to analyze the novelty of gene sets and augment gene sets with additional relevant genes. The Geneshot web-server and API are freely and openly available from https://amp.pharm.mssm.edu/geneshot.

ChEA3: transcription factor enrichment analysis by orthogonal omics integration.

Identifying the transcription factors (TFs) responsible for observed changes in gene expression is an important step in understanding gene regulatory networks. ChIP-X Enrichment Analysis 3 (ChEA3) is a transcription factor enrichment analysis tool that ranks TFs associated with user-submitted gene sets. The ChEA3 background database contains a collection of gene set libraries generated from multiple sources including TF-gene co-expression from RNA-seq studies, TF-target associations from ChIP-seq experiments, and TF-gene co-occurrence computed from crowd-submitted gene lists. Enrichment results from these distinct sources are integrated to generate a composite rank that improves the prediction of the correct upstream TF compared to ranks produced by individual libraries. We compare ChEA3 with existing TF prediction tools and show that ChEA3 performs better. By integrating the ChEA3 libraries, we illuminate general transcription factor properties such as whether the TF behaves as an activator or a repressor. The ChEA3 web-server is available from https://amp.pharm.mssm.edu/ChEA3.

eXpression2Kinases (X2K) Web: linking expression signatures to upstream cell signaling networks.

While gene expression data at the mRNA level can be globally and accurately measured, profiling the activity of cell signaling pathways is currently much more difficult. eXpression2Kinases (X2K) computationally predicts involvement of upstream cell signaling pathways, given a signature of differentially expressed genes. X2K first computes enrichment for transcription factors likely to regulate the expression of the differentially expressed genes. The next step of X2K connects these enriched transcription factors through known protein-protein interactions (PPIs) to construct a subnetwork. The final step performs kinase enrichment analysis on the members of the subnetwork. X2K Web is a new implementation of the original eXpression2Kinases algorithm with important enhancements. X2K Web includes many new transcription factor and kinase libraries, and PPI networks. For demonstration, thousands of gene expression signatures induced by kinase inhibitors, applied to six breast cancer cell lines, are provided for fetching directly into X2K Web. The results are displayed as interactive downloadable vector graphic network images and bar graphs. Benchmarking various settings via random permutations enabled the identification of an optimal set of parameters to be used as the default settings in X2K Web. X2K Web is freely available from http://X2K.cloud.

Data Portal for the Library of Integrated Network-based Cellular Signatures (LINCS) program: integrated access to diverse large-scale cellular perturbation response data.

The Library of Integrated Network-based Cellular Signatures (LINCS) program is a national consortium funded by the NIH to generate a diverse and extensive reference library of cell-based perturbation-response signatures, along with novel data analytics tools to improve our understanding of human diseases at the systems level. In contrast to other large-scale data generation efforts, LINCS Data and Signature Generation Centers (DSGCs) employ a wide range of assay technologies cataloging diverse cellular responses. Integration of, and unified access to LINCS data has therefore been particularly challenging. The Big Data to Knowledge (BD2K) LINCS Data Coordination and Integration Center (DCIC) has developed data standards specifications, data processing pipelines, and a suite of end-user software tools to integrate and annotate LINCS-generated data, to make LINCS signatures searchable and usable for different types of users. Here, we describe the LINCS Data Portal (LDP) (http://lincsportal.ccs.miami.edu/), a unified web interface to access datasets generated by the LINCS DSGCs, and its underlying database, LINCS Data Registry (LDR). LINCS data served on the LDP contains extensive metadata and curated annotations. We highlight the features of the LDP user interface that is designed to enable search, browsing, exploration, download and analysis of LINCS data and related curated content.

A PRMT5-ZNF326 axis mediates innate immune activation upon replication stress.

DNA replication stress (RS) is a widespread phenomenon in carcinogenesis, causing genomic instability and extensive chromatin alterations. DNA damage leads to activation of innate immune signaling, but little is known about transcriptional regulators mediating such signaling upon RS. Using a chemical screen, we identified protein arginine methyltransferase 5 (PRMT5) as a key mediator of RS-dependent induction of interferon-stimulated genes (ISGs). This response is also associated with reactivation of endogenous retroviruses (ERVs). Using quantitative mass spectrometry, we identify proteins with PRMT5-dependent symmetric dimethylarginine (SDMA) modification induced upon RS. Among these, we show that PRMT5 targets and modulates the activity of ZNF326, a zinc finger protein essential for ISG response. Our data demonstrate a role for PRMT5-mediated SDMA in the context of RS-induced transcriptional induction, affecting physiological homeostasis and cancer therapy.

Isoform-resolved transcriptome of the human preimplantation embryo.

Human preimplantation development involves extensive remodeling of RNA expression and splicing. However, its transcriptome has been compiled using short-read sequencing data, which fails to capture most full-length mRNAs. Here, we generate an isoform-resolved transcriptome of early human development by performing long- and short-read RNA sequencing on 73 embryos spanning the zygote to blastocyst stages. We identify 110,212 unannotated isoforms transcribed from known genes, including highly conserved protein-coding loci and key developmental regulators. We further identify 17,964 isoforms from 5,239 unannotated genes, which are largely non-coding, primate-specific, and highly associated with transposable elements. These isoforms are widely supported by the integration of published multi-omics datasets, including single-cell 8CLC and blastoid studies. Alternative splicing and gene co-expression network analyses further reveal that embryonic genome activation is associated with splicing disruption and transient upregulation of gene modules. Together, these findings show that the human embryo transcriptome is far more complex than currently known, and will act as a valuable resource to empower future studies exploring development.

WNTinib is a multi-kinase inhibitor with specificity against ß-catenin mutant hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) is a leading cause of cancer-related deaths worldwide. ß-Catenin (CTNNB1)-mutated HCC represents 30% of cases of the disease with no precision therapeutics available. Using chemical libraries derived from clinical multi-kinase inhibitor (KI) scaffolds, we screened HCC organoids to identify WNTinib, a KI with exquisite selectivity in CTNNB1-mutated human and murine models, including patient samples. Multiomic and target engagement analyses, combined with rescue experiments and in vitro and in vivo efficacy studies, revealed that WNTinib is superior to clinical KIs and inhibits KIT/mitogen-activated protein kinase (MAPK) signaling at multiple nodes. Moreover, we demonstrate that reduced engagement on BRAF and p38α kinases by WNTinib relative to several multi-KIs is necessary to avoid compensatory feedback signaling-providing a durable and selective transcriptional repression of mutant ß-catenin/Wnt targets through nuclear translocation of the EZH2 transcriptional repressor. Our studies uncover a previously unknown mechanism to harness the KIT/MAPK/EZH2 pathway to potently and selectively antagonize CTNNB1-mutant HCC with an unprecedented wide therapeutic index.

Validation of the Orgasm Rating Scale in the Context of Masturbation.

BACKGROUND: The Orgasm Rating Scale (ORS) assesses the subjective experience of orgasm. Its psychometric properties have only been examined in the context of sexual intercourse. This study aims to validate the ORS in the context of solitary masturbation. METHODS: A sample of 1,171 men and 1,424 women aged 18-83 years ( M = 40.51, SD = 12.07) completed the ORS in the solitary masturbation context along with other scales to assess sexual attitudes, solitary sexual desire, propensity for sexual arousal/inhibition, and sexual functioning. RESULTS: A four-dimensional structure was confirmed, similar to the homologous version for the context of sexual intercourse. Measures obtained from the ORS were sex and age invariant, exhibited adequate internal consistency, discriminated between people with orgasmic difficulties, and were associated with related variables. CONCLUSIONS: The ORS is a multidimensional measure that provides reliable, valid measures of the subjective experience of orgasm in the context of solitary masturbation.

Age alters the oncogenic trajectory toward luminal mammary tumors that activate unfolded proteins responses.

A major limitation in the use of mouse models in breast cancer research is that most mice develop estrogen receptor-alpha (ERα)-negative mammary tumors, while in humans, the majority of breast cancers are ERα-positive. Therefore, developing mouse models that best mimic the disease in humans is of fundamental need. Here, using an inducible MMTV-rtTA/TetO-NeuNT mouse model, we show that despite being driven by the same oncogene, mammary tumors in young mice are ERα-negative, while they are ERα-positive in aged mice. To further elucidate the mechanisms for this observation, we performed RNAseq analysis and identified genes that are uniquely expressed in aged female-derived mammary tumors. We found these genes to be involved in the activation of the ERα axis of the mitochondrial UPR and the ERα-mediated regulation of XBP-1s, a gene involved in the endoplasmic reticulum UPR. Collectively, our results indicate that aging alters the oncogenic trajectory towards the ERα-positive subtype of breast cancers, and that mammary tumors in aged mice are characterized by the upregulation of multiple UPR stress responses regulated by the ERα.

Development of Potent Cellular and Humoral Immune Responses in Long-Term Hemodialysis Patients After 1273-mRNA SARS-CoV-2 Vaccination.

Long-term hemodialysis (HD) patients are considered vulnerable and at high-risk of developing severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) infection due to their immunocompromised condition. Since COVID-19 associated mortality rates are higher in HD patients, vaccination is critical to protect them. The response towards vaccination against COVID-19 in HD patients is still uncertain and, in particular the cellular immune response is not fully understood. We monitored the humoral and cellular immune responses by analysis of the serological responses and Spike-specific cellular immunity in COVID-19-recovered and naïve HD patients in a longitudinal study shortly after vaccination to determine the protective effects of 1273-mRNA vaccination against SARS-CoV-2 in these high-risk patients. In naïve HD patients, the cellular immune response measured by IL-2 and IFN-É£ secretion needed a second vaccine dose to significantly increase, with a similar pattern for the humoral response. In contrast, COVID-19 recovered HD patients developed a potent and rapid cellular and humoral immune response after the first vaccine dose. Interestingly, when comparing COVID-19 recovered healthy volunteers (HV), previously vaccinated with BNT162b2 vaccine to HD patients vaccinated with 1273-mRNA, these exhibited a more robust immune response that is maintained longitudinally. Our results indicate that HD patients develop strong cellular and humoral immune responses to 1273-mRNA vaccination and argue in favor of personalized immune monitoring studies in HD patients, especially if COVID-19 pre-exposed, to adapt COVID-19 vaccination protocols for this immunocompromised population.

PPARdelta activation induces metabolic and contractile maturation of human pluripotent stem cell-derived cardiomyocytes.

Pluripotent stem-cell-derived cardiomyocytes (PSC-CMs) provide an unprecedented opportunity to study human heart development and disease, but they are functionally and structurally immature. Here, we induce efficient human PSC-CM (hPSC-CM) maturation through metabolic-pathway modulations. Specifically, we find that peroxisome-proliferator-associated receptor (PPAR) signaling regulates glycolysis and fatty acid oxidation (FAO) in an isoform-specific manner. While PPARalpha (PPARa) is the most active isoform in hPSC-CMs, PPARdelta (PPARd) activation efficiently upregulates the gene regulatory networks underlying FAO, increases mitochondrial and peroxisome content, enhances mitochondrial cristae formation, and augments FAO flux. PPARd activation further increases binucleation, enhances myofibril organization, and improves contractility. Transient lactate exposure, which is frequently used for hPSC-CM purification, induces an independent cardiac maturation program but, when combined with PPARd activation, still enhances oxidative metabolism. In summary, we investigate multiple metabolic modifications in hPSC-CMs and identify a role for PPARd signaling in inducing the metabolic switch from glycolysis to FAO in hPSC-CMs.

Rapid, scalable assessment of SARS-CoV-2 cellular immunity by whole-blood PCR.

Fast, high-throughput methods for measuring the level and duration of protective immune responses to SARS-CoV-2 are needed to anticipate the risk of breakthrough infections. Here we report the development of two quantitative PCR assays for SARS-CoV-2-specific T cell activation. The assays are rapid, internally normalized and probe-based: qTACT requires RNA extraction and dqTACT avoids sample preparation steps. Both assays rely on the quantification of CXCL10 messenger RNA, a chemokine whose expression is strongly correlated with activation of antigen-specific T cells. On restimulation of whole-blood cells with SARS-CoV-2 viral antigens, viral-specific T cells secrete IFN-γ, which stimulates monocytes to produce CXCL10. CXCL10 mRNA can thus serve as a proxy to quantify cellular immunity. Our assays may allow large-scale monitoring of the magnitude and duration of functional T cell immunity to SARS-CoV-2, thus helping to prioritize revaccination strategies in vulnerable populations.