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.

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.

A multi-omic analysis of MCF10A cells provides a resource for integrative assessment of ligand-mediated molecular and phenotypic responses.

The phenotype of a cell and its underlying molecular state is strongly influenced by extracellular signals, including growth factors, hormones, and extracellular matrix proteins. While these signals are normally tightly controlled, their dysregulation leads to phenotypic and molecular states associated with diverse diseases. To develop a detailed understanding of the linkage between molecular and phenotypic changes, we generated a comprehensive dataset that catalogs the transcriptional, proteomic, epigenomic and phenotypic responses of MCF10A mammary epithelial cells after exposure to the ligands EGF, HGF, OSM, IFNG, TGFB and BMP2. Systematic assessment of the molecular and cellular phenotypes induced by these ligands comprise the LINCS Microenvironment (ME) perturbation dataset, which has been curated and made publicly available for community-wide analysis and development of novel computational methods ( synapse.org/LINCS_MCF10A ). In illustrative analyses, we demonstrate how this dataset can be used to discover functionally related molecular features linked to specific cellular phenotypes. Beyond these analyses, this dataset will serve as a resource for the broader scientific community to mine for biological insights, to compare signals carried across distinct molecular modalities, and to develop new computational methods for integrative data analysis.

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α.

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.

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.

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.

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.

Validation of the Orgasm Rating Scale in the Context of Masturbation / Validación de la Orgasm Rating Scale en el Contexto de la Masturbación

Background: The Orgasm Rating Scale (ORS) assesses the subjectiveexperience of orgasm. Its psychometric properties have only been examinedin the context of sexual intercourse. This study aims to validate the ORSin the context of solitary masturbation. Methods: A sample of 1,171 menand 1,424 women aged 18-83 years (M = 40.51, SD = 12.07) completedthe ORS in the solitary masturbation context along with other scales toassess sexual attitudes, solitary sexual desire, propensity for sexual arousal/inhibition, and sexual functioning. Results: A four-dimensional structurewas confirmed, similar to the homologous version for the context of sexualintercourse. Measures obtained from the ORS were sex and age invariant,exhibited adequate internal consistency, discriminated between peoplewith orgasmic difficulties, and were associated with related variables.Conclusions: The ORS is a multidimensional measure that providesreliable, valid measures of the subjective experience of orgasm in thecontext of solitary masturbation.

Antecedentes: la Orgasm Rating Scale (ORS) evalúa la experienciasubjetiva del orgasmo. Sus propiedades psicométricas únicamente seexaminaron en el contexto de las relaciones sexuales. El objetivo de esteestudio es validar la ORS en el contexto de la masturbación en solitario.Método: una muestra formada por 1.171 hombres y 1.424 mujeres de 18a 83 años (M = 40,51; DT = 12,07) completaron la ORS en el contexto dela masturbación en solitario junto con otras escalas para evaluar actitudessexuales, deseo sexual solitario, propensión a la excitación/inhibiciónsexual y funcionamiento sexual. Resultados: se confirma una estructurade cuatro dimensiones, al igual que su versión homóloga para el contextode las relaciones sexuales. Las medidas obtenidas mediante la ORS soninvariantes por sexo y edad, sus dimensiones muestran una adecuadaconsistencia interna, discriminan entre personas con y sin dificultadesorgásmicas y se relacionan con variables afines. Conclusiones: la ORSes una medida multidimensional que aporta medidas fiables y válidas dela experiencia subjetiva del orgasmo en el contexto de la masturbación ensolitario.

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.

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.

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.

Transcriptome analysis of human dermal fibroblasts following red light phototherapy.

Fibrosis occurs when collagen deposition and fibroblast proliferation replace healthy tissue. Red light (RL) may improve skin fibrosis via photobiomodulation, the process by which photosensitive chromophores in cells absorb visible or near-infrared light and undergo photophysical reactions. Our previous research demonstrated that high fluence RL reduces fibroblast proliferation, collagen deposition, and migration. Despite the identification of several cellular mechanisms underpinning RL phototherapy, little is known about the transcriptional changes that lead to anti-fibrotic cellular responses. Herein, RNA sequencing was performed on human dermal fibroblasts treated with RL phototherapy. Pathway enrichment and transcription factor analysis revealed regulation of extracellular matrices, proliferation, and cellular responses to oxygen-containing compounds following RL phototherapy. Specifically, RL phototherapy increased the expression of MMP1, which codes for matrix metalloproteinase-1 (MMP-1) and is responsible for remodeling extracellular collagen. Differential regulation of MMP1 was confirmed with RT-qPCR and ELISA. Additionally, RL upregulated PRSS35, which has not been previously associated with skin activity, but has known anti-fibrotic functions. Our results suggest that RL may benefit patients by altering fibrotic gene expression.

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.

PRDM15 is a key regulator of metabolism critical to sustain B-cell lymphomagenesis.

PRDM (PRDI-BF1 and RIZ homology domain containing) family members are sequence-specific transcriptional regulators involved in cell identity and fate determination, often dysregulated in cancer. The PRDM15 gene is of particular interest, given its low expression in adult tissues and its overexpression in B-cell lymphomas. Despite its well characterized role in stem cell biology and during early development, the role of PRDM15 in cancer remains obscure. Herein, we demonstrate that while PRDM15 is largely dispensable for mouse adult somatic cell homeostasis in vivo, it plays a critical role in B-cell lymphomagenesis. Mechanistically, PRDM15 regulates a transcriptional program that sustains the activity of the PI3K/AKT/mTOR pathway and glycolysis in B-cell lymphomas. Abrogation of PRDM15 induces a metabolic crisis and selective death of lymphoma cells. Collectively, our data demonstrate that PRDM15 fuels the metabolic requirement of B-cell lymphomas and validate it as an attractive and previously unrecognized target in oncology.

Interoperable RNA-Seq analysis in the cloud.

RNA-Sequencing (RNA-Seq) is currently the leading technology for genome-wide transcript quantification. Mapping the raw reads to transcript and gene level counts can be achieved by different aligners. Here we report an in-depth comparison of transcript quantification methods. Our goal is the specific use of cost-efficient RNA-Seq analysis for deployment in a cloud infrastructure composed of interacting microservices. The individual modules cover file transfer into the cloud and APIs to handle the cloud alignment jobs. We next demonstrate how newly generated RNA-Seq data can be placed in the context of thousands of previously published datasets in near real time. With in-depth benchmarks, we identify suitable gene count quantification methods to facilitate cost-effective, accurate, and cloud-based RNA-Seq analysis service. Pseudo-alignment algorithms such as kallisto and Salmon combine high read quality estimation with cost efficient runtime performance. HISAT2 is the fastest of the classical aligners with good alignment quality. This article is part of a Special Issue entitled: Transcriptional Profiles and Regulatory Gene Networks edited by Dr. Federico Manuel Giorgi and Dr. Shaun Mahony.