MAP-RSeq: Mayo Analysis Pipeline for RNA sequencing

The impact of hypoxia preconditioning on mesenchymal stem cells performance in hypertensive kidney disease

BACKGROUND

Autologous mesenchymal stem cells (MSCs) have emerged as a therapeutic option for many diseases. Hypertensive kidney disease (HKD) might impair MSCs' reparative ability by altering the biomolecular properties, but the characteristics of this impairment are unclear. In our previous pre-clinical studies, we found hypoxic preconditioning (HPC) enhanced angiogenesis and suppressed senescence gene expression. Thus, we hypothesize that HPC would improve human MSCs by enhancing their functionality and angiogenesis, creating an anti-inflammatory and anti-senescence environment.

METHODS

MSC samples (n = 12 each) were collected from the abdominal fat of healthy kidney donors (HC), hypertensive patients (HTN), and patients with hypertensive kidney disease (HKD). MSCs were harvested and cultured in Normoxic (20% O2) or Hypoxic (1% O2) conditions. MSC functionality was measured by proliferation assays and cytokine released in conditioned media. Senescence was evaluated by senescence-associated beta-galactosidase (SA-beta-gal) activity. Additionally, transcriptome analysis using RNA-sequencing and quantitative PCR (qPCR) were performed.

RESULTS

At baseline, normoxic HTN-MSCs had higher proliferation capacity compared to HC. However, HPC augmented proliferation in HC. HPC did not affect the release of pro-angiogenic protein VEGF, but increased EGF in HC-MSC, and decreased HGF in HC and HKD MSCs. Under HPC, SA-β-gal activity tended to decrease, particularly in HC group. HPC upregulated mostly the pro-angiogenic and inflammatory genes in HC and HKD and a few senescence genes in HKD.

CONCLUSIONS

HPC has a more favorable functional effect on HC- than on HKD-MSC, reflected in increased proliferation and EGF release, and modest decrease in senescence, whereas it has little effect on HTN or HKD MSCs.

Chronic kidney disease and left ventricular diastolic dysfunction (CKD-LVDD) alter cardiac expression of mitochondria-related genes in swine

Cardiovascular disease and heart failure doubles in patients with chronic kidney disease (CKD), but the underlying mechanisms remain obscure. Mitochondria are central to maintaining cellular respiration and modulating cardiomyocyte function. We took advantage of our novel swine model of CKD and left ventricular diastolic dysfunction (CKD-LVDD) to investigate the expression of mitochondria-related genes and potential mechanisms regulating their expression. CKD-LVDD and normal control pigs (n=6/group, 3 males/3 females) were studied for 14 weeks. Renal and cardiac hemodynamics were quantified by multidetector-CT, echocardiography, and pressure-volume loop studies, respectively. Mitochondrial morphology (electron microscopy) and function (Oroboros) were assessed ex vivo. In randomly selected pigs (n=3/group), cardiac mRNA-, MeDIP-, and miRNA-sequencing (seq) were performed to identify mitochondria-related genes and study their pre- and post -transcriptional regulation. CKD-LVDD exhibited cardiac mitochondrial structural abnormalities and elevated mitochondrial H2O2 emission but preserved mitochondrial function. Cardiac mRNA-seq identified 862 mitochondria-related genes, of which 69 were upregulated and 33 downregulated (fold-change ≥2, false discovery rate≤0.05). Functional analysis showed that upregulated genes were primarily implicated in processes associated with oxidative stress, whereas those downregulated mainly participated in respiration and ATP synthesis. Integrated mRNA/miRNA/MeDIP-seq analysis showed that upregulated genes were modulated predominantly by miRNAs, whereas those downregulated were by miRNA and epigenetic mechanisms. CKD-LVDD alters cardiac expression of mitochondria-related genes, associated with mitochondrial structural damage but preserved respiratory function, possibly reflecting intrinsic compensatory mechanisms. Our findings may guide the development of early interventions at stages of cardiac dysfunction in which mitochondrial injury could be prevented, and the development of LVDD ameliorated.

Bridging the Gap: Multi-Omics Profiling of Brain Tissue in Alzheimer's Disease and Older Controls in Multi-Ethnic Populations

INTRODUCTION

Multi-omics studies in Alzheimer's disease (AD) revealed many potential disease pathways and therapeutic targets. Despite their promise of precision medicine, these studies lacked African Americans (AA) and Latin Americans (LA), who are disproportionately affected by AD.

METHODS

To bridge this gap, Accelerating Medicines Partnership in AD (AMP-AD) expanded brain multi-omics profiling to multi-ethnic donors.

RESULTS

We generated multi-omics data and curated and harmonized phenotypic data from AA (n=306), LA (n=326), or AA and LA (n=4) brain donors plus Non-Hispanic White (n=252) and other (n=20) ethnic groups, to establish a foundational dataset enriched for AA and LA participants. This study describes the data available to the research community, including transcriptome from three brain regions, whole genome sequence, and proteome measures.

DISCUSSION

Inclusion of traditionally underrepresented groups in multi-omics studies is essential to discover the full spectrum of precision medicine targets that will be pertinent to all populations affected with AD.

Abundant transcriptomic alterations in the human cerebellum of patients with a C9orf72 repeat expansion

The most prominent genetic cause of both amyotrophic lateral sclerosis (ALS) and frontotemporal lobar degeneration (FTLD) is a repeat expansion in the gene C9orf72. Importantly, the transcriptomic consequences of the C9orf72 repeat expansion remain largely unclear. Here, we used short-read RNA sequencing (RNAseq) to profile the cerebellar transcriptome, detecting alterations in patients with a C9orf72 repeat expansion. We focused on the cerebellum, since key C9orf72-related pathologies are abundant in this neuroanatomical region, yet TDP-43 pathology and neuronal loss are minimal. Consistent with previous work, we showed a reduction in the expression of the C9orf72 gene and an elevation in homeobox genes, when comparing patients with the expansion to both patients without the C9orf72 repeat expansion and control subjects. Interestingly, we identified more than 1000 alternative splicing events, including 4 in genes previously associated with ALS and/or FTLD. We also found an increase of cryptic splicing in C9orf72 patients compared to patients without the expansion and controls. Furthermore, we demonstrated that the expression level of select RNA-binding proteins is associated with cryptic splice junction inclusion. Overall, this study explores the presence of widespread transcriptomic changes in the cerebellum, a region not confounded by severe neurodegeneration, in post-mortem tissue from C9orf72 patients.

Identification of skewed X chromosome inactivation using exome and transcriptome sequencing in patients with suspected rare genetic disease

BACKGROUND

X-chromosome inactivation (XCI) is an epigenetic process that occurs during early development in mammalian females by randomly silencing one of two copies of the X chromosome in each cell. The preferential inactivation of either the maternal or paternal copy of the X chromosome in a majority of cells results in a skewed or non-random pattern of X inactivation and is observed in over 25% of adult females. Identifying skewed X inactivation is of clinical significance in patients with suspected rare genetic diseases due to the possibility of biased expression of disease-causing genes present on the active X chromosome. The current clinical test for the detection of skewed XCI relies on the methylation status of the methylation-sensitive restriction enzyme (Hpall) binding site present in proximity of short tandem polymorphic repeats on the androgen receptor (AR) gene. This approach using one locus results in uninformative or inconclusive data for 10-20% of tests. Further, recent studies have shown inconsistency between methylation of the AR locus and the state of inactivation of the X chromosome. Herein, we develop a method for estimating X inactivation status, using exome and transcriptome sequencing data derived from blood in 227 female samples. We built a reference model for evaluation of XCI in 135 females from the GTEx consortium. We tested and validated the model on 11 female individuals with different types of undiagnosed rare genetic disorders who were clinically tested for X-skew using the AR gene assay and compared results to our outlier-based analysis technique.

RESULTS

In comparison to the AR clinical test for identification of X inactivation, our method was concordant with the AR method in 9 samples, discordant in 1, and provided a measure of X inactivation in 1 sample with uninformative clinical results. We applied this method on an additional 81 females presenting to the clinic with phenotypes consistent with different hereditary disorders without a known genetic diagnosis.

CONCLUSIONS

This study presents the use of transcriptome and exome sequencing data to provide an accurate and complete estimation of X-inactivation and skew status in a cohort of female patients with different types of suspected rare genetic disease.

Early B cell transcriptomic markers of measles-specific humoral immunity following a 3rd dose of MMR vaccine

B cell transcriptomic signatures hold promise for the early prediction of vaccine-induced humoral immunity and vaccine protective efficacy. We performed a longitudinal study in 232 healthy adult participants before/after a 3rd dose of MMR (MMR3) vaccine. We assessed baseline and early transcriptional patterns in purified B cells and their association with measles-specific humoral immunity after MMR vaccination using two analytical methods ("per gene" linear models and joint analysis). Our study identified distinct early transcriptional signatures/genes following MMR3 that were associated with measles-specific neutralizing antibody titer and/or binding antibody titer. The most significant genes included: the interleukin 20 receptor subunit beta/IL20RB gene (a subunit receptor for IL-24, a cytokine involved in the germinal center B cell maturation/response); the phorbol-12-myristate-13-acetate-induced protein 1/PMAIP1, the brain expressed X-linked 2/BEX2 gene and the B cell Fas apoptotic inhibitory molecule/FAIM, involved in the selection of high-affinity B cell clones and apoptosis/regulation of apoptosis; as well as IL16 (encoding the B lymphocyte-derived IL-16 ligand of CD4), involved in the crosstalk between B cells, dendritic cells and helper T cells. Significantly enriched pathways included B cell signaling, apoptosis/regulation of apoptosis, metabolic pathways, cell cycle-related pathways, and pathways associated with viral infections, among others. In conclusion, our study identified genes/pathways linked to antigen-induced B cell proliferation, differentiation, apoptosis, and clonal selection, that are associated with, and impact measles virus-specific humoral immunity after MMR vaccination.

Predictive value of CDC37 gene expression for targeted therapy in metastatic colorectal cancer

BACKGROUND

CDC37 is a key determinant of client kinase recruitment to the HSP90 chaperoning system. We hypothesized that kinase-specific dependency on CDC37 alters the efficacy of targeted therapies for metastatic colorectal cancer (mCRC).

MATERIAL AND METHODS

Two independent mCRC cohorts were analyzed to compare the survival outcomes between CDC37-high and CDC37-low patients (stratified by the median cutoff values): the CALGB/SWOG 80405 trial (226 and 207 patients receiving first-line bevacizumab- and cetuximab-containing chemotherapies, respectively) and Japanese retrospective (50 refractory patients receiving regorafenib) cohorts. A dataset of specimens submitted to a commercial CLIA-certified laboratory was utilized to characterize molecular profiles of CDC37-high (top quartile, N = 5055) and CDC37-low (bottom quartile, N = 5055) CRCs.

RESULTS

In the bevacizumab-treated group, CDC37-high patients showed significantly better progression-free survival (PFS) (median 13.3 vs 9.6 months, hazard ratio [HR] 0.59, 95% confidence interval [CI] 0.44-0.79, p < 0.01) than CDC37-low patients. In the cetuximab-treated group, CDC37-high and CDC37-low patients had similar outcomes. In the regorafenib-treated group, CDC37-high patients showed significantly better overall survival (median 11.3 vs 6.0 months, HR 0.24, 95% CI 0.11-0.54, p < 0.01) and PFS (median 3.5 vs 1.9 months, HR 0.51, 95% CI 0.28-0.94, p = 0.03). Comprehensive molecular profiling revealed that CDC37-high CRCs were associated with higher VEGFA, FLT1, and KDR expressions and activated hypoxia signature.

CONCLUSIONS

CDC37-high mCRC patients derived more benefit from anti-VEGF therapies, including bevacizumab and regorafenib, but not from cetuximab. Molecular profiles suggested that such tumors were dependent on angiogenesis-relating pathways.

BEERS2: RNA-Seq simulation through high fidelity in silico modeling

Simulation of RNA-seq reads is critical in the assessment, comparison, benchmarking and development of bioinformatics tools. Yet the field of RNA-seq simulators has progressed little in the last decade. To address this need we have developed BEERS2, which combines a flexible and highly configurable design with detailed simulation of the entire library preparation and sequencing pipeline. BEERS2 takes input transcripts (typically fully length messenger RNA transcripts with polyA tails) from either customizable input or from CAMPAREE simulated RNA samples. It produces realistic reads of these transcripts as FASTQ, SAM or BAM formats with the SAM or BAM formats containing the true alignment to the reference genome. It also produces true transcript-level quantification values. BEERS2 combines a flexible and highly configurable design with detailed simulation of the entire library preparation and sequencing pipeline and is designed to include the effects of polyA selection and RiboZero for ribosomal depletion, hexamer priming sequence biases, GC-content biases in polymerase chain reaction (PCR) amplification, barcode read errors and errors during PCR amplification. These characteristics combine to make BEERS2 the most complete simulation of RNA-seq to date. Finally, we demonstrate the use of BEERS2 by measuring the effect of several settings on the popular Salmon pseudoalignment algorithm.

Follicular lymphoma B cells exhibit heterogeneous transcriptional states with associated somatic alterations and tumor microenvironments

Follicular lymphoma (FL) is an indolent non-Hodgkin lymphoma of germinal center origin, which presents with significant biologic and clinical heterogeneity. Using RNA-seq on B cells sorted from 87 FL biopsies, combined with machine-learning approaches, we identify 3 transcriptional states that divide the biological ontology of FL B cells into inflamed, proliferative, and chromatin-modifying states, with relationship to prior GC B cell phenotypes. When integrated with whole-exome sequencing and immune profiling, we find that each state was associated with a combination of mutations in chromatin modifiers, copy-number alterations to TNFAIP3, and T follicular helper cells (Tfh) cell interactions, or primarily by a microenvironment rich in activated T cells. Altogether, these data define FL B cell transcriptional states across a large cohort of patients, contribute to our understanding of FL heterogeneity at the tumor cell level, and provide a foundation for guiding therapeutic intervention.

EHMT2 Inactivation in Pancreatic Epithelial Cells Shapes the Transcriptional Landscape and Inflammation Response of the Whole Pancreas

The Euchromatic Histone Methyl Transferase Protein 2 (EHMT2), also known as G9a, deposits transcriptionally repressive chromatin marks that play pivotal roles in the maturation and homeostasis of multiple organs. Recently, we have shown that EHMT2 inactivation alters growth and immune gene expression networks, antagonizing KRAS-mediated pancreatic cancer initiation and promotion. Here, we elucidate the essential role of EHMT2 in maintaining a transcriptional landscape that protects organs from inflammation. Comparative RNA-seq studies between normal postnatal and young adult pancreatic tissue from EHMT2 conditional knockout animals ( EHMT2 fl/fl ) targeted to the exocrine pancreatic epithelial cells ( Pdx1-Cre and P48 Cre/+ ), reveal alterations in gene expression networks in the whole organ related to injury-inflammation-repair, suggesting an increased predisposition to damage. Thus, we induced an inflammation repair response in the EHMT2 fl/fl pancreas and used a data science-based approach to integrate RNA-seq-derived pathways and networks, deconvolution digital cytology, and spatial transcriptomics. We also analyzed the tissue response to damage at the morphological, biochemical, and molecular pathology levels. The EHMT2 fl/fl pancreas displays an enhanced injury-inflammation-repair response, offering insights into fundamental molecular and cellular mechanisms involved in this process. More importantly, these data show that conditional EHMT2 inactivation in exocrine cells reprograms the local environment to recruit mesenchymal and immunological cells needed to mount an increased inflammatory response. Mechanistically, this response is an enhanced injury-inflammation-repair reaction with a small contribution of specific EHMT2-regulated transcripts. Thus, this new knowledge extends the mechanisms underlying the role of the EHMT2-mediated pathway in suppressing pancreatic cancer initiation and modulating inflammatory pancreatic diseases.

Gene expression biomarkers differentiate overall survival of colorectal cancer upon targeted therapies

While monoclonal antibody-based targeted therapies have substantially improved progression-free survival in cancer patients, the variability in individual responses poses a significant challenge in patient care. Therefore, identifying cancer subtypes and their associated biomarkers is required for assigning effective treatment. In this study, we integrated genotype and pre-treatment tissue RNA-seq data and identified biomarkers causally associated with the overall survival (OS) of colorectal cancer (CRC) patients treated with either cetuximab or bevacizumab. We performed enrichment analysis for specific consensus molecular subtypes (CMS) of colorectal cancer and evaluated differential expression of identified genes using paired tumor and normal tissue from an external cohort. In addition, we replicated the causal effect of these genes on OS using validation cohort and assessed their association with the Cancer Genome Atlas Program data as an external cohort. One of the replicated findings was WDR62, whose overexpression shortened OS of patients treated with cetuximab. Enrichment of its over expression in CMS1 and low expression in CMS4 suggests that patients with CMS4 subtype may drive greater benefit from cetuximab. In summary, this study highlights the importance of integrating different omics data for identifying promising biomarkers specific to a treatment or a cancer subtype.

ABCA7 deficiency causes neuronal dysregulation by altering mitochondrial lipid metabolism

ABCA7 loss-of-function variants are associated with increased risk of Alzheimer's disease (AD). Using ABCA7 knockout human iPSC models generated with CRISPR/Cas9, we investigated the impacts of ABCA7 deficiency on neuronal metabolism and function. Lipidomics revealed that mitochondria-related phospholipids, such as phosphatidylglycerol and cardiolipin were reduced in the ABCA7-deficient iPSC-derived cortical organoids. Consistently, ABCA7 deficiency-induced alterations of mitochondrial morphology accompanied by reduced ATP synthase activity and exacerbated oxidative damage in the organoids. Furthermore, ABCA7-deficient iPSC-derived neurons showed compromised mitochondrial respiration and excess ROS generation, as well as enlarged mitochondrial morphology compared to the isogenic controls. ABCA7 deficiency also decreased spontaneous synaptic firing and network formation in iPSC-derived neurons, in which the effects were rescued by supplementation with phosphatidylglycerol or NAD+ precursor, nicotinamide mononucleotide. Importantly, effects of ABCA7 deficiency on mitochondria morphology and synapses were recapitulated in synaptosomes isolated from the brain of neuron-specific Abca7 knockout mice. Together, our results provide evidence that ABCA7 loss-of-function contributes to AD risk by modulating mitochondria lipid metabolism.

A small molecule MST1/2 inhibitor accelerates murine liver regeneration with improved survival in models of steatohepatitis

Dysfunctional liver regeneration following surgical resection remains a major cause of postoperative mortality and has no therapeutic options. Without targeted therapies, the current treatment paradigm relies on supportive therapy until homeostasis can be achieved. Pharmacologic acceleration of regeneration represents an alternative therapeutic avenue. Therefore, we aimed to generate a small molecule inhibitor that could accelerate liver regeneration with an emphasis on diseased models, which represent a significant portion of patients who require surgical resection and are often not studied. Utilizing a clinically approved small molecule inhibitor as a parent compound, standard medicinal chemistry approaches were utilized to generate a small molecule inhibitor targeting serine/threonine kinase 4/3 (MST1/2) with reduced off-target effects. This compound, mCLC846, was then applied to preclinical models of murine partial hepatectomy, which included models of diet-induced metabolic dysfunction-associated steatohepatitis (MASH). mCLC846 demonstrated on target inhibition of MST1/2 and reduced epidermal growth factor receptor inhibition. The inhibitory effects resulted in restored pancreatic beta-cell function and survival under diabetogenic conditions. Liver-specific cell-line exposure resulted in Yes-associated protein activation. Oral delivery of mCLC846 perioperatively resulted in accelerated murine liver regeneration and improved survival in diet-induced MASH models. Bulk transcriptional analysis of regenerating liver remnants suggested that mCLC846 enhanced the normal regenerative pathways and induced them following liver resection. Overall, pharmacological acceleration of liver regeneration with mCLC846 was feasible, had an acceptable therapeutic index, and provided a survival benefit in models of diet-induced MASH.

Topical Ruxolitinib in the Treatment of Necrobiosis Lipoidica: A Prospective, Open-Label Study

Necrobiosis lipoidica (NL) is a rare granulomatous disease. There are few effective treatments for NL. We sought to investigate the efficacy and safety of the Jak1/2 inhibitor, ruxolitnib, in the treatment of NL and identify the biomarkers associated with the disease and treatment response. We conducted an open-label, phase 2 study of ruxolitinib in 12 patients with NL. We performed transcriptomic analysis of tissue samples before and after treatment. At week 12, the mean NL lesion score decreased by 58.2% (SD = 28.7%, P = .003). Transcriptomic analysis demonstrated enrichment of type I and type II IFN pathways in baseline disease. Weighted gene coexpression network analysis demonstrated post-treatment changes in IFN pathways with key hub genes IFNG and signal transducer and activator of transcription 1 gene STAT1. Limitations include small sample size and a study group limited to patients with <10% body surface area. In conclusion, ruxolitinib is an effective treatment for NL and targets the key pathogenic mediators of the disease.

AdipoRon reduces TGFβ1-mediated collagen deposition in vitro and alleviates knee stiffness in vivo

Arthrofibrosis, which causes joint motion restrictions, is a common complication following total knee arthroplasty (TKA). Key features associated with arthrofibrosis include myofibroblast activation, knee stiffness, and excessive scar tissue formation. We previously demonstrated that adiponectin levels are suppressed within the knee tissues of patients affected by arthrofibrosis and showed that AdipoRon, an adiponectin receptor agonist, exhibited anti-fibrotic properties in human mesenchymal stem cells. In this study, the therapeutic potential of AdipoRon was evaluated on TGFβ1-mediated myofibroblast differentiation of primary human knee fibroblasts and in a mouse model of knee stiffness. Picrosirius red staining revealed that AdipoRon reduced TGFβ1-induced collagen deposition in primary knee fibroblasts derived from patients undergoing primary TKA and revision TKA for arthrofibrosis. AdipoRon also reduced mRNA and protein levels of ACTA2, a key myofibroblast marker. RNA-seq analysis corroborated the anti-myofibrogenic effects of AdipoRon. In our knee stiffness mouse model, 6 weeks of knee immobilization, to induce a knee contracture, in conjunction with daily vehicle (DMSO) or AdipoRon (1, 5, and 25 mg/kg) via intraperitoneal injections were well tolerated based on animal behavior and weight measurements. Biomechanical testing demonstrated that passive extension angles (PEAs) of experimental knees were similar between vehicle and AdipoRon treatment groups in mice evaluated immediately following immobilization. Interestingly, relative to vehicle-treated mice, 5 mg/kg AdipoRon therapy improved the PEA of the experimental knees in mice that underwent 4 weeks of knee remobilization following the immobilization and therapy. Together, these studies revealed that AdipoRon may be an effective therapeutic modality for arthrofibrosis.

miRNA family miR-29 inhibits PINK1-PRKN dependent mitophagy via ATG9A

Loss-of-function mutations in the genes encoding PINK1 and PRKN result in early-onset Parkinson disease (EOPD). Together the encoded enzymes direct a neuroprotective pathway that ensures the elimination of damaged mitochondria via autophagy. We performed a genome-wide high content imaging miRNA screen for inhibitors of the PINK1-PRKN pathway and identified all three members of the miRNA family 29 (miR-29). Using RNAseq we identified target genes and found that siRNA against ATG9A phenocopied the effects of miR-29 and inhibited the initiation of PINK1-PRKN mitophagy. Furthermore, we discovered two rare, potentially deleterious, missense variants (p.R631W and p.S828L) in our EOPD cohort and tested them experimentally in cells. While expression of wild-type ATG9A was able to rescue the effects of miR-29a, the EOPD-associated variants behaved like loss-of-function mutations. Together, our study validates miR-29 and its target gene ATG9A as novel regulators of mitophagy initiation. It further serves as proof-of-concept of finding novel, potentially disease-causing EOPD-linked variants specifically in mitophagy regulating genes. The nomination of genetic variants and biological pathways is important for the stratification and treatment of patients that suffer from devastating diseases, such as EOPD.

Cryptic splicing of stathmin-2 and UNC13A mRNAs is a pathological hallmark of TDP-43-associated Alzheimer's disease

Nuclear clearance and cytoplasmic accumulations of the RNA-binding protein TDP-43 are pathological hallmarks in almost all patients with amyotrophic lateral sclerosis (ALS) and up to 50% of patients with frontotemporal dementia (FTD) and Alzheimer's disease. In Alzheimer's disease, TDP-43 pathology is predominantly observed in the limbic system and correlates with cognitive decline and reduced hippocampal volume. Disruption of nuclear TDP-43 function leads to abnormal RNA splicing and incorporation of erroneous cryptic exons in numerous transcripts including Stathmin-2 (STMN2, also known as SCG10) and UNC13A, recently reported in tissues from patients with ALS and FTD. Here, we identify both STMN2 and UNC13A cryptic exons in Alzheimer's disease patients, that correlate with TDP-43 pathology burden, but not with amyloid-β or tau deposits. We also demonstrate that processing of the STMN2 pre-mRNA is more sensitive to TDP-43 loss of function than UNC13A. In addition, full-length RNAs encoding STMN2 and UNC13A are suppressed in large RNA-seq datasets generated from Alzheimer's disease post-mortem brain tissue. Collectively, these results open exciting new avenues to use STMN2 and UNC13A as potential therapeutic targets in a broad range of neurodegenerative conditions with TDP-43 proteinopathy including Alzheimer's disease.

Malignant female genital tract smooth muscle tumors with adipocytic differentiation: A morphologic, immunohistochemical, MDM2 fluorescence in situ hybridization and molecular genetic study of 6 lipoleiomyosarcomas

Leiomyosarcoma with adipocytic differentiation or lipoleiomyosarcoma is an uncommon sarcoma of the female genital tract with only a few individual reports in the literature. We therefore performed a morphologic, immunohistochemical, MDM2 gene amplification and RNA and DNA sequencing analysis of a series of gynecologic lipoleiomyosarcoma to better define the clinicopathologic spectrum. Six tumors from 6 patients were identified and classified as spindled lipoleiomyosarcoma (n = 2), mixed spindled and myxoid lipoleiomyosarcoma (n = 1), epithelioid lipoleiomyosarcoma with focal myxoid features (n = 1) and mixed spindled and epithelioid lipoleiomyosarcoma (n = 2). Patient age ranged from 41 to 64 years (mean: 49; median: 50). Primary location included uterine corpus (3), uterine corpus/cervix (2) and broad ligament (1). Tumor size ranged from 4.5 to 22 cm (mean: 11.2; median: 9.8). Four patients had metastasis at presentation or subsequently developed recurrent or distant disease. Patient status was known for 5: 2 dead of disease, 2 alive with disease and 1 alive without evidence of disease. Immunohistochemical expression of smooth muscle markers, ER, PR and WT-1 showed patterns similar to non-adipocytic gynecologic leiomyosarcomas. MDM2 amplification fluorescence in situ hybridization performed on 2 tumors was negative in 1 and equivocal in 1. Sequencing studies performed on 3 tumors found TP53 mutations in 3, with 1 tumor also having an ATRX alteration. No gene fusions were identified. Although lipoleiomyosarcomas have a diverse morphologic spectrum, our findings suggest the smooth muscle component shares morphologic and immunohistochemical features with female genital tract non-adipocytic leiomyosarcomas. Lipoleiomyosarcomas also have genetic alterations associated with non-adipocytic gynecologic leiomyosarcomas.

Direct early growth response-1 knockdown decreases melanoma viability independent of mitogen-activated extracellular signal-related kinase inhibition

To investigate downstream molecular changes caused by mitogen-activated protein kinase (MEK) inhibitor treatment and further explore the impact of direct knockdown of early growth response-1 ( EGR1 ) in melanoma cell culture. RNA-sequencing (RNA-Seq) was performed to determine gene expression changes with MEK inhibitor treatment. Treatment with MEK inhibitor (trametinib) was then assessed in two cutaneous (MEL888, MEL624) and one conjunctival (YUARGE 13-3064) melanoma cell line. Direct knockdown of EGR1 was accomplished using lentiviral vectors containing shRNA. Cell viability was measured using PrestoBlueHS Cell Viability Reagent. Total RNA and protein were assessed by qPCR and SimpleWestern. RNA-Seq demonstrated a profound reduction in EGR1 with MEK inhibitor treatment, prompting further study of melanoma cell lines. Following trametinib treatment of melanoma cells, viability was reduced in both cutaneous (MEL888 26%, P  < 0.01; MEL624 27%, P  < 0.001) and conjunctival (YUARGE 13-3064 33%, P  < 0.01) melanoma compared with DMSO control, with confirmed EGR1 knockdown to 0.04-, 0.01-, and 0.16-fold DMSO-treated levels (all P  < 0.05) in MEL888, MEL624, and YUARGE 13-3064, respectively. Targeted EGR1 knockdown using shRNA reduced viability in both cutaneous (MEL624 78%, P  = 0.05) and conjunctival melanoma (YUARGE-13-3064 67%, P  = 0.02). RNA-Sequencing in MEK inhibitor-treated cells identified EGR1 as a candidate effector molecule of interest. In a malignant melanoma cell population, MEK inhibition reduced viability in both cutaneous and conjunctival melanoma with a profound downstream reduction in EGR1 expression. Targeted knockdown of EGR1 reduced both cutaneous and conjunctival melanoma cell viability independent of MEK inhibition, suggesting a key role for EGR1 in melanoma pathobiology.

Hypomethylation and Overexpression of Th17-Associated Genes is a Hallmark of Intestinal CD4+ Lymphocytes in Crohn's Disease

BACKGROUND

The development of Crohn's disease [CD] involves immune cell signalling pathways regulated by epigenetic modifications. Aberrant DNA methylation has been identified in peripheral blood and bulk intestinal tissue from CD patients. However, the DNA methylome of disease-associated intestinal CD4+ lymphocytes has not been evaluated.

MATERIALS AND METHODS

Genome-wide DNA methylation sequencing was performed from terminal ileum CD4+ cells from 21 CD patients and 12 age- and sex-matched controls. Data were analysed for differentially methylated CpGs [DMCs] and methylated regions [DMRs]. Integration was performed with RNA-sequencing data to evaluate the functional impact of DNA methylation changes on gene expression. DMRs were overlapped with regions of differentially open chromatin [by ATAC-seq] and CCCTC-binding factor [CTCF] binding sites [by ChIP-seq] between peripherally derived Th17 and Treg cells.

RESULTS

CD4+ cells in CD patients had significantly increased DNA methylation compared to those from the controls. A total of 119 051 DMCs and 8113 DMRs were detected. While hypermethylated genes were mostly related to cell metabolism and homeostasis, hypomethylated genes were significantly enriched within the Th17 signalling pathway. The differentially enriched ATAC regions in Th17 cells [compared to Tregs] were hypomethylated in CD patients, suggesting heightened Th17 activity. There was significant overlap between hypomethylated DNA regions and CTCF-associated binding sites.

CONCLUSIONS

The methylome of CD patients shows an overall dominant hypermethylation yet hypomethylation is more concentrated in proinflammatory pathways, including Th17 differentiation. Hypomethylation of Th17-related genes associated with areas of open chromatin and CTCF binding sites constitutes a hallmark of CD-associated intestinal CD4+ cells.

Transcriptomic changes during the replicative senescence of human articular chondrocytes

Osteoarthritis (OA) is a degenerative joint disease and a leading cause of disability worldwide. Aging is a major risk factor for OA, but the specific mechanisms underlying this connection remain unclear. Although chondrocytes rarely divide in adult articular cartilage, they undergo replicative senescence in vitro which provides an opportunity to study changes related to aging under controlled laboratory conditions. In this pilot study, we performed bulk RNA sequencing on early- and late-passage human articular chondrocytes to identify transcriptomic changes associated with cellular aging. Chondrocytes were isolated from the articular cartilage of three donors, two with OA (age 70-80 years) and one with healthy cartilage (age 26 years). Chondrocytes were serially passaged until replicative senescence and RNA extracted from early- and late-passage cells. Principal component analysis of all genes showed clear separation between early- and late-passage chondrocytes, indicating substantial age-related differences in gene expression. Differentially expressed genes (DEGs) analysis confirmed distinct transcriptomic profiles between early- and late-passage chondrocytes. Hierarchical clustering revealed contrasting expression patterns between the two isolates from osteoarthritic samples and the healthy sample. Focused analysis of DEGs on transcripts associated with turnover of the extra-cellular matrix and the senescence-associated secretory phenotype (SASP) showed consistent downregulation of Col2A1 and ACAN, and upregulation of MMP19, ADAMTS4, and ADAMTS8 in late passage chondrocytes across all samples. SASP components including IL-1α, IL-1β, IL-6, IL-7, p16INK4A (CDKN2A) and CCL2 demonstrated significant upregulation in late passage chondrocytes originally isolated from OA samples. Pathway analysis between sexes with OA revealed shared pathways such as extracellular matrix (ECM) organization, collagen formation, skeletal and muscle development, and nervous system development. Sex-specific differences were observed, with males showing distinctions in ECM organization, regulation of the cell cycle process as well as neuron differentiation. In contrast, females exhibited unique variations in the regulation of the cell cycle process, DNA metabolic process, and the PID-PLK1 pathway.

Impact of biological sex and sex hormones on molecular signatures of skeletal muscle at rest and in response to distinct exercise training modes

Substantial divergence in cardio-metabolic risk, muscle size, and performance exists between men and women. Considering the pivotal role of skeletal muscle in human physiology, we investigated and found, based on RNA sequencing (RNA-seq), that differences in the muscle transcriptome between men and women are largely related to testosterone and estradiol and much less related to genes located on the Y chromosome. We demonstrate inherent unique, sex-dependent differences in muscle transcriptional responses to aerobic, resistance, and combined exercise training in young and older cohorts. The hormonal changes with age likely explain age-related differential expression of transcripts. Furthermore, in primary human myotubes we demonstrate the profound but distinct effects of testosterone and estradiol on amino acid incorporation to multiple individual proteins with specific functions. These results clearly highlight the potential of designing exercise programs tailored specifically to men and women and have implications for people who change gender by altering their hormone profile.

Cross species systems biology discovers glial DDR2, STOM, and KANK2 as therapeutic targets in progressive supranuclear palsy

Progressive supranuclear palsy (PSP) is a neurodegenerative parkinsonian disorder characterized by cell-type-specific tau lesions in neurons and glia. Prior work uncovered transcriptome changes in human PSP brains, although their cell-specificity is unknown. Further, systematic data integration and experimental validation platforms to prioritize brain transcriptional perturbations as therapeutic targets in PSP are currently lacking. In this study, we combine bulk tissue (n = 408) and single nucleus RNAseq (n = 34) data from PSP and control brains with transcriptome data from a mouse tauopathy and experimental validations in Drosophila tau models for systematic discovery of high-confidence expression changes in PSP with therapeutic potential. We discover, replicate, and annotate thousands of differentially expressed genes in PSP, many of which reside in glia-enriched co-expression modules and cells. We prioritize DDR2, STOM, and KANK2 as promising therapeutic targets in PSP with striking cross-species validations. We share our findings and data via our interactive application tool PSP RNAseq Atlas ( https://rtools.mayo.edu/PSP_RNAseq_Atlas/ ). Our findings reveal robust glial transcriptome changes in PSP, provide a cross-species systems biology approach, and a tool for therapeutic target discoveries in PSP with potential application in other neurodegenerative diseases.

Mild and moderate chronic hypercapnia elicit distinct transcriptomic responses of immune function in cardiorespiratory nuclei

Chronic hypercapnia (CH) is a hallmark of respiratory-related diseases, and the level of hypercapnia can acutely or progressively become more severe. Previously, we have shown time-dependent adaptations in steady-state physiology during mild (arterial Pco2 ∼55 mmHg) and moderate (∼60 mmHg) CH in adult goats, including transient (mild CH) or sustained (moderate CH) suppression of acute chemosensitivity suggesting limitations in adaptive respiratory control mechanisms as the level of CH increases. Changes in specific markers of glutamate receptor plasticity, interleukin-1ß, and serotonergic modulation within key nodes of cardiorespiratory control do not fully account for the physiological adaptations to CH. Here, we used an unbiased approach (bulk tissue RNA sequencing) to test the hypothesis that mild or moderate CH elicits distinct gene expression profiles in important brain stem regions of cardiorespiratory control, which may explain the contrasting responses to CH. Gene expression profiles from the brain regions validated the accuracy of tissue biopsy methodology. Differential gene expression analyses revealed greater effects of CH on brain stem sites compared with the medial prefrontal cortex. Mild CH elicited an upregulation of predominantly immune-related genes and predicted activation of immune-related pathways and functions. In contrast, moderate CH broadly led to downregulation of genes and predicted inactivation of cellular pathways related to the immune response and vascular function. These data suggest that mild CH leads to a steady-state activation of neuroinflammatory pathways within the brain stem, whereas moderate CH drives the opposite response. Transcriptional shifts in immune-related functions may underlie the cardiorespiratory network's capability to respond to acute, more severe hypercapnia when in a state of progressively increased CH.NEW & NOTEWORTHY Mild chronic hypercapnia (CH) broadly upregulated immune-related genes and a predicted activation of biological pathways related to immune cell activity and the overall immune response. In contrast, moderate CH primarily downregulated genes related to major histocompatibility complex signaling and vasculature function that led to a predicted inactivation of pathways involving the immune response and vascular endothelial function. The severity-dependent effect on immune responses suggests that neuroinflammation has an important role in CH and may be important in the maintenance of proper ventilatory responses to acute and chronic hypercapnia.

BMAL1 modulates senescence programming via AP-1

Cellular senescence and circadian dysregulation are biological hallmarks of aging. Whether they are coordinately regulated has not been thoroughly studied. We hypothesize that BMAL1, a pioneer transcription factor and master regulator of the molecular circadian clock, plays a role in the senescence program. Here, we demonstrate BMAL1 is significantly upregulated in senescent cells and has altered rhythmicity compared to non-senescent cells. Through BMAL1-ChIP-seq, we show that BMAL1 is uniquely localized to genomic motifs associated with AP-1 in senescent cells. Integration of BMAL1-ChIP-seq data with RNA-seq data revealed that BMAL1 presence at AP-1 motifs is associated with active transcription. Finally, we showed that BMAL1 contributes to AP-1 transcriptional control of key features of the senescence program, including altered regulation of cell survival pathways, and confers resistance to drug-induced apoptosis. Overall, these results highlight a previously unappreciated role of the core circadian clock component BMAL1 on the molecular phenotype of senescent cells.

Host immunity and KLF 11 deficiency together promote fibrosis in a mouse model of endometriosis

BACKGROUND

Endometriosis is a debilitating disease typically characterized by prolific fibrotic scarring. Earlier we reported downregulation of two transcription factors belonging TGF-βR signaling pathway Sp/Krüppel-like factor 11 (KLF11) and 10 (KLF10) in human endometriosis lesions. Here we investigated the role of these nuclear factors and immunity in the scaring fibrosis associated with endometriosis.

METHODS

We used a well characterized experimental mouse model of endometriosis. WT, KLF10 or KLF11 deficient mice were compared. The lesions were evaluated histologically, fibrosis was quantified with Masons' Trichome staining, immune-infiltrates were quantified by immunohistochemistry, peritoneal adhesions were score, gene expression was evaluated by bulk RNA sequencing.

RESULTS

Intense fibrotic reactions and large changes in gene expression were detected in KLF11 deficient implants associated with squamous metaplasia of the ectopic endometrium, as compared to KLF10 deficient or WT implants. Fibrosis was mitigated with pharmacologic agents that blocked histone acetylation or TGF-βR signaling or with genetic deficiency for SMAD3. The lesions were richly infiltrated with T-cells, regulatory T-cells, and innate immune cells. Fibrosis was exacerbated when implants expressed ectopic genes implicating autoimmunity as a major factor contributing to the scaring fibrosis.

CONCLUSIONS

Our findings identify KLF11 and TGF-βR signaling as cell intrinsic mechanisms and autoimmune responses as cell extrinsic mechanisms of scaring fibrosis in ectopic endometrium lesions.

GENERAL SIGNIFICANCE

Immunological factors associated with inflammation and tissue repair drive scaring fibrosis in experimental endometriosis, providing the rationale for immune therapy of endometriosis.

Tissue Eosinophilia is Superior to an Analysis by Polyp Status for the Chronic Rhinosinusitis Transcriptome: An RNA Study

OBJECTIVE

RNA sequencing (transcriptomics) is used to study biological pathways. However, the yield of data depends on comparing well-characterized cohorts. We compared tissue eosinophilia versus nasal polyp (NP) status as the metric to characterize transcriptomic mechanisms at play in eosinophilic and non-eosinophilic chronic rhinosinusitis (CRS) versus controls.

METHODS

RNA sequencing was conducted on sinonasal tissue samples of CRS and controls. Analyses were conducted based on polyp status [with nasal polyps (CRSwNP) and without nasal polyps (CRSsNP)] as well as tissue eosinophil levels per high power field (eos/hpf)[non-eosinophilic (<10 eos/hpf, neCRS) or eosinophilic (≥10 eos/hpf, eCRS)]. The yield of differentially expressed genes (DEGs) and biological pathways through Ingenuity Pathway Analysis (IPA) were compared.

RESULTS

CRS tissue differed from controls by 736 statistically significant DEGs. Both NP status and tissue eosinophilia were effective in differentiating CRS from controls and into two distinct subgroups. Statistically significant DEGs identified when comparing CRS by NP status were 60, whereas 110 DEGs were identified using eosinophil cutoff ≥10 and <10 eos/hpf. Additionally, heatmaps showed greater homogeneity within each CRS subgroup when analyzed by tissue eosinophilia versus NP status. On IPA, the IL-17 signaling pathway was significantly different only by tissue eosinophilia status, not NP status, being higher in CRS <10 eos/hpf.

CONCLUSION

Tissue eosinophilia is superior to an analysis by NP status for the study of CRS transcriptome by RNA sequencing in identifying DEGs. Classification of CRS samples by eosinophil counts agnostic of NP status may offer advantageous insights into CRS pathogenetic mechanisms.

LEVEL OF EVIDENCE

3 Laryngoscope, 133:2480-2489, 2023.

Lipoblastoma-Like Tumor and Fibrosarcoma-Like Lipomatous Neoplasm Represent the Same Entity: A Clinicopathologic and Molecular Genetic Study of 23 Cases Occurring in Both Men and Women at Diverse Locations

Lipoblastoma-like tumor (LLT) is a benign soft tissue tumor demonstrating mixed morphologic features of lipoblastoma, myxoid liposarcoma, and spindle cell lipoma but lacking genetic alterations associated with those tumors. LLT was originally thought to be specific to the vulva but has since been reported in the paratesticular region. The morphologic features of LLT overlap with those of "fibrosarcoma-like lipomatous neoplasm" (FLLN), a rare, indolent adipocytic neoplasm considered by some to form part of the spectrum of atypical spindle cell and pleomorphic lipomatous tumor. We compared the morphologic, immunohistochemical, and genetic features of 23 tumors previously classified as LLT (n = 17) and FLLN (n = 6). The 23 tumors occurred in 13 women and 10 men (mean age, 42 years; range, 17 to 80 years). Eighteen (78%) cases arose in the inguinogenital region, whereas 5 tumors (22%) involved noninguinogenital soft tissue, including the flank (n = 1), shoulder (n = 1), foot (n = 1), forearm (n = 1), and chest wall (n = 1). Microscopically, the tumors were lobulated and septated, with variably collagenized fibromyxoid stroma, prominent thin-walled vessels, scattered univacuolated or bivacuolated lipoblasts, and a minor component of mature adipose tissue. Using immunohistochemistry, 5 tumors (42%) showed complete RB1 loss, with partial loss in 7 cases (58%). RNA sequencing, chromosomal microarray, and DNA next-generation sequencing study results were negative for significant alterations. There were no clinical, morphologic, immunohistochemical, or molecular genetic differences between cases previously classified as LLT or FLLN. Clinical follow-up (11 patients [48%]; range, 2-276 months; mean, 48.2 months) showed all patients were alive without disease, and only one patient had experienced a single local recurrence. We conclude that LLT and FLLN represent the same entity, for which "LLT" seems most appropriate. LLT may occur in either sex and any superficial soft tissue location. Careful morphologic study and appropriate ancillary testing should allow for the distinction of LLT from its potential mimics.

The microRNomes of Chinese hamster ovary (CHO) cells and their extracellular vesicles, and how they respond to osmotic and ammonia stress

A new area of focus in Chinese hamster ovary (CHO) biotechnology is the role of small (exosomes) and large (microvesicles or microparticles) extracellular vesicles (EVs). CHO cells in culture exchange large quantities of proteins and RNA through these EVs, yet the content and role of these EVs remain elusive. MicroRNAs (miRs or miRNA) are central to adaptive responses to stress and more broadly to changes in culture conditions. Given that EVs are highly enriched in miRs, and that EVs release large quantities of miRs both in vivo and in vitro, EVs and their miR content likely play an important role in adaptive responses. Here we report the miRNA landscape of CHO cells and their EVs under normal culture conditions and under ammonia and osmotic stress. We show that both cells and EVs are highly enriched in five miRs (among over 600 miRs) that make up about half of their total miR content, and that these highly enriched miRs differ significantly between normal and stress culture conditions. Notable is the high enrichment in miR-92a and miR-23a under normal culture conditions, in contrast to the high enrichment in let-7 family miRs (let-7c, let-7b, and let-7a) under both stress conditions. The latter suggests a preserved stress-responsive function of the let-7 miR family, one of the most highly preserved miR families across species, where among other functions, let-7 miRs regulate core oncogenes, which, depending on the biological context, may tip the balance between cell cycle arrest and apoptosis. While the expected-based on their profound enrichment-important role of these highly enriched miRs remains to be dissected, our data and analysis constitute an important resource for exploring the role of miRs in cell adaptation as well as for synthetic applications.

Transcriptomic and immunophenotypic characterization of two cases of adamantinoma-like Ewing sarcoma of the thyroid gland

INTRODUCTION

Adamantinoma-like Ewing sarcoma (ALES) is a rare aggressive malignancy occasionally diagnosed in the thyroid gland. ALES shows basaloid cytomorphology, expresses keratins, p63, p40, frequently CD99, and harbours the t(11;22) EWSR1::FLI1 translocation. There is debate on whether ALES resembles more sarcoma or carcinoma.

METHODS

We performed RNA sequencing from two ALES cases and compared findings with skeletal Ewing's sarcomas and nonneoplastic thyroid tissue. ALES was investigated by in situ hybridization (ISH) for high-risk human papillomavirus (HPV) DNA and immunohistochemistry for the following antigens: keratin 7, keratin 20, keratin 5, keratins (AE1/AE3 and CAM5.2), CD45, CD20, CD5, CD99, chromogranin, synaptophysin, calcitonin, thyroglobulin, PAX8, TTF1, S100, p40, p63, p16, NUT, desmin, ER, FLI1, INI1, and myogenin.

RESULTS

An uncommon EWSR1::FLI transcript with retained EWSR1 exon 8 was detected in both ALES cases. Regulators of EWSR1::FLI1 splicing (HNRNPH1, SUPT6H, SF3B1) necessary for production of a functional fusion oncoprotein, as well as 53 genes (including TNNT1, NKX2.2) activated downstream to the EWSR1::FLI1 cascade, were overexpressed. Eighty-six genes were uniquely overexpressed in ALES, most of which were related to squamous differentiation. Immunohistochemically, ALES strongly expressed keratins 5, AE1/AE3 and CAM5.2, p63, p40, p16, and focally CD99. INI1 was retained. The remaining immunostains and HPV DNA ISH were negative.

CONCLUSION

Comparative transcriptomic profiling reveals overlapping features of ALES with skeletal Ewing's sarcoma and an epithelial carcinoma, as evidenced by immunohistochemical expression of keratin 5, p63, p40, CD99, the transcriptome profile, and detection of EWSR1::FLI1 fusion transcript by RNA sequencing.

Renal ischemia alters the transcriptomic and epigenetic profile of inflammatory genes in swine scattered tubular-like cells

BACKGROUND

Scattered tubular-like cells (STCs) are differentiated renal tubular cells that during recovery from ischemic injury dedifferentiate to repair other injured renal cells. Renal artery stenosis (RAS), often associated with chronic inflammatory injury, compromises the integrity and function of STCs, but the underlying mechanisms remain unknown. We hypothesized that RAS alters the transcriptomic and epigenetic profile of inflammatory genes in swine STCs.

METHODS

STCs were harvested from pig kidneys after 10 weeks of RAS or sham (n=6 each). STC mRNA profiles of inflammatory genes were analyzed using high-throughput mRNA-sequencing (seq) and their DNA methylation (5mC) and hydroxymethylation (5hmC) profiles by DNA immunoprecipitation and next-generation sequencing (MeDIP-seq) (n=3 each), followed by an integrated (mRNA-seq/MeDIP-seq) analysis. STC protein expression of candidate differentially expressed (DE) genes and common proinflammatory proteins were subsequently assessed in vitro before and after epigenetic (Bobcat339) modulation.

RESULTS

mRNA-seq identified 57 inflammatory genes up-regulated in RAS-STCs versus Normal-STCs (>1.4 or <0.7-fold, P<0.05), of which 14% exhibited lower 5mC and 5% higher 5hmC levels in RAS-STCs versus Normal-STCs, respectively. Inflammatory gene and protein expression was higher in RAS-STCs compared with Normal-STCs but normalized after epigenetic modulation.

CONCLUSIONS

These observations highlight a novel modulatory mechanism of this renal endogenous repair system and support development of epigenetic or anti-inflammatory therapies to preserve the reparative capacity of STCs in individuals with RAS.

Cell-specific crosstalk proteomics reveals cathepsin B signaling as a driver of glioblastoma malignancy near the subventricular zone

Glioblastoma (GBM) is the most prevalent and aggressive malignant primary brain tumor. GBM proximal to the lateral ventricles (LVs) is more aggressive, potentially due to subventricular zone (SVZ) contact. Despite this, crosstalk between GBM and neural stem/progenitor cells (NSC/NPCs) is not well understood. Using cell-specific proteomics, we show that LV-proximal GBM prevents neuronal maturation of NSCs through induction of senescence. Additionally, GBM brain tumor initiating cells (BTICs) increase expression of CTSB upon interaction with NPCs. Lentiviral knockdown and recombinant protein experiments reveal both cell-intrinsic and soluble CTSB promote malignancy-associated phenotypes in BTICs. Soluble CTSB stalls neuronal maturation in NPCs while promoting senescence, providing a link between LV-tumor proximity and neurogenesis disruption. Finally, we show LV-proximal CTSB upregulation in patients, showing the relevance of this crosstalk in human GBM biology. These results demonstrate the value of proteomic analysis in tumor microenvironment research and provide direction for new therapeutic strategies in GBM.

Highlights

Periventricular GBM is more malignant and disrupts neurogenesis in a rodent model.Cell-specific proteomics elucidates tumor-promoting crosstalk between GBM and NPCs.NPCs induce upregulated CTSB expression in GBM, promoting tumor progression.GBM stalls neurogenesis and promotes NPC senescence via CTSB.

Virus-specific and shared gene expression signatures in immune cells after vaccination in response to influenza and vaccinia stimulation

Background

In the vaccine era, individuals receive multiple vaccines in their lifetime. Host gene expression in response to antigenic stimulation is usually virus-specific; however, identifying shared pathways of host response across a wide spectrum of vaccine pathogens can shed light on the molecular mechanisms/components which can be targeted for the development of broad/universal therapeutics and vaccines.

Method

We isolated PBMCs, monocytes, B cells, and CD8+ T cells from the peripheral blood of healthy donors, who received both seasonal influenza vaccine (within <1 year) and smallpox vaccine (within 1 - 4 years). Each of the purified cell populations was stimulated with either influenza virus or vaccinia virus. Differentially expressed genes (DEGs) relative to unstimulated controls were identified for each in vitro viral infection, as well as for both viral infections (shared DEGs). Pathway enrichment analysis was performed to associate identified DEGs with KEGG/biological pathways.

Results

We identified 2,906, 3,888, 681, and 446 DEGs in PBMCs, monocytes, B cells, and CD8+ T cells, respectively, in response to influenza stimulation. Meanwhile, 97, 120, 20, and 10 DEGs were identified as gene signatures in PBMCs, monocytes, B cells, and CD8+ T cells, respectively, upon vaccinia stimulation. The majority of DEGs identified in PBMCs were also found in monocytes after either viral stimulation. Of the virus-specific DEGs, 55, 63, and 9 DEGs occurred in common in PBMCs, monocytes, and B cells, respectively, while no DEGs were shared in infected CD8+ T cells after influenza and vaccinia. Gene set enrichment analysis demonstrated that these shared DEGs were over-represented in innate signaling pathways, including cytokine-cytokine receptor interaction, viral protein interaction with cytokine and cytokine receptor, Toll-like receptor signaling, RIG-I-like receptor signaling pathways, cytosolic DNA-sensing pathways, and natural killer cell mediated cytotoxicity.

Conclusion

Our results provide insights into virus-host interactions in different immune cells, as well as host defense mechanisms against viral stimulation. Our data also highlights the role of monocytes as a major cell population driving gene expression in ex vivo PBMCs in response to viral stimulation. The immune response signaling pathways identified in this study may provide specific targets for the development of novel virus-specific therapeutics and improved vaccines for vaccinia and influenza. Although influenza and vaccinia viruses have been selected in this study as pathogen models, this approach could be applicable to other pathogens.

Dynamics and survival associations of T cell receptor clusters in patients with pleural mesothelioma treated with immunotherapy

BACKGROUND

Immune checkpoint inhibitors (ICIs) are now a first-line treatment option for patients with pleural mesothelioma with the recent approval of ipilimumab and nivolumab. Mesothelioma has a low tumor mutation burden and no robust predictors of survival with ICI. Since ICIs enable adaptive antitumor immune responses, we investigated T-cell receptor (TCR) associations with survival in participants from two clinical trials treated with ICI.

METHODS

We included patients with pleural mesothelioma who were treated with nivolumab (NivoMes, NCT02497508) or nivolumab and ipilimumab (INITIATE, NCT03048474) after first-line therapy. TCR sequencing was performed with the ImmunoSEQ assay in 49 and 39 pretreatment and post-treatment patient peripheral blood mononuclear cell (PBMC) samples. These data were integrated with TCR sequences found in bulk RNAseq data by TRUST4 program in 45 and 35 pretreatment and post-treatment tumor biopsy samples and TCR sequences from over 600 healthy controls. The TCR sequences were clustered into groups of shared antigen specificity using GIANA. Associations of TCR clusters with overall survival were determined by cox proportional hazard analysis.

RESULTS

We identified 4.2 million and 12 thousand complementarity-determining region 3 (CDR3) sequences from PBMCs and tumors, respectively, in patients treated with ICI. These CDR3 sequences were integrated with 2.1 million publically available CDR3 sequences from healthy controls and clustered. ICI-enhanced T-cell infiltration and expanded T cell diversity in tumors. Cases with TCR clones in the top tertile in the pretreatment tissue or in circulation had significantly better survival than the bottom two tertiles (p<0.04). Furthermore, a high number of shared TCR clones between pretreatment tissue and in circulation was associated with improved survival (p=0.01). To potentially select antitumor clusters, we filtered for clusters that were (1) not found in healthy controls, (2) recurrent in multiple patients with mesothelioma, and (3) more prevalent in post-treatment than pretreatment samples. The detection of two-specific TCR clusters provided significant survival benefit compared with detection of 1 cluster (HR<0.001, p=0.026) or the detection of no TCR clusters (HR=0.10, p=0.002). These two clusters were not found in bulk tissue RNA-seq data and have not been reported in public CDR3 databases.

CONCLUSIONS

We identified two unique TCR clusters that were associated with survival on treatment with ICI in patients with pleural mesothelioma. These clusters may enable approaches for antigen discovery and inform future targets for design of adoptive T cell therapies.

Obesity and dyslipidemia are associated with partially reversible modifications to DNA hydroxymethylation of apoptosis- and senescence-related genes in swine adipose-derived mesenchymal stem/stromal cells

BACKGROUND

Obesity dysregulates key biological processes underlying the functional homeostasis, fate decisions, and reparative potential of mesenchymal stem/stromal cells (MSCs). Mechanisms directing obesity-induced phenotypic alterations in MSCs remain unclear, but emerging drivers include dynamic modification of epigenetic marks, like 5-hydroxymethylcytosine (5hmC). We hypothesized that obesity and cardiovascular risk factors induce functionally relevant, locus-specific changes in 5hmC of swine adipose-derived MSCs and evaluated their reversibility using an epigenetic modulator, vitamin-C.

METHODS

Female domestic pigs were fed a 16-week Lean or Obese diet (n = 6 each). MSCs were harvested from subcutaneous adipose tissue, and 5hmC profiles were examined through hydroxymethylated DNA immunoprecipitation sequencing (hMeDIP-seq) followed by an integrative (hMeDIP and mRNA sequencing) gene set enrichment analysis. For clinical context, we compared 5hmC profiles of adipose tissue-derived human MSCs harvested from patients with obesity and healthy controls.

RESULTS

hMeDIP-seq revealed 467 hyper- (fold change ≥ 1.4; p-value ≤ 0.05) and 591 hypo- (fold change ≤ 0.7; p-value ≤ 0.05) hydroxymethylated loci in swine Obese- versus Lean-MSCs. Integrative hMeDIP-seq/mRNA-seq analysis identified overlapping dysregulated gene sets and discrete differentially hydroxymethylated loci with functions related to apoptosis, cell proliferation, and senescence. These 5hmC changes were associated with increased senescence in cultured MSCs (p16/CDKN2A immunoreactivity, senescence-associated β-galactosidase [SA-β-Gal] staining), were partly reversed in swine Obese-MSCs treated with vitamin-C, and shared common pathways with 5hmC changes in human Obese-MSCs.

CONCLUSIONS

Obesity and dyslipidemia are associated with dysregulated DNA hydroxymethylation of apoptosis- and senescence-related genes in swine and human MSCs, potentially affecting cell vitality and regenerative functions. Vitamin-C may mediate reprogramming of this altered epigenomic landscape, providing a potential strategy to improve the success of autologous MSC transplantation in obese patients.

Transcriptome-Powered Pluripotent Stem Cell Differentiation for Regenerative Medicine

Pluripotent stem cells are endless sources for in vitro engineering human tissues for regenerative medicine. Extensive studies have demonstrated that transcription factors are the key to stem cell lineage commitment and differentiation efficacy. As the transcription factor profile varies depending on the cell type, global transcriptome analysis through RNA sequencing (RNAseq) has been a powerful tool for measuring and characterizing the success of stem cell differentiation. RNAseq has been utilized to comprehend how gene expression changes as cells differentiate and provide a guide to inducing cellular differentiation based on promoting the expression of specific genes. It has also been utilized to determine the specific cell type. This review highlights RNAseq techniques, tools for RNAseq data interpretation, RNAseq data analytic methods and their utilities, and transcriptomics-enabled human stem cell differentiation. In addition, the review outlines the potential benefits of the transcriptomics-aided discovery of intrinsic factors influencing stem cell lineage commitment, transcriptomics applied to disease physiology studies using patients' induced pluripotent stem cell (iPSC)-derived cells for regenerative medicine, and the future outlook on the technology and its implementation.

BEERS2: RNA-Seq simulation through high fidelity in silico modeling

Simulation of RNA-seq reads is critical in the assessment, comparison, benchmarking, and development of bioinformatics tools. Yet the field of RNA-seq simulators has progressed little in the last decade. To address this need we have developed BEERS2, which combines a flexible and highly configurable design with detailed simulation of the entire library preparation and sequencing pipeline. BEERS2 takes input transcripts (typically fully-length mRNA transcripts with polyA tails) from either customizable input or from CAMPAREE simulated RNA samples. It produces realistic reads of these transcripts as FASTQ, SAM, or BAM formats with the SAM or BAM formats containing the true alignment to the reference genome. It also produces true transcript-level quantification values. BEERS2 combines a flexible and highly configurable design with detailed simulation of the entire library preparation and sequencing pipeline and is designed to include the effects of polyA selection and RiboZero for ribosomal depletion, hexamer priming sequence biases, GC-content biases in PCR amplification, barcode read errors, and errors during PCR amplification. These characteristics combine to make BEERS2 the most complete simulation of RNA-seq to date. Finally, we demonstrate the use of BEERS2 by measuring the effect of several settings on the popular Salmon pseudoalignment algorithm.

Transcriptomic landscapes of effective and failed liver regeneration in humans

Background & Aims

Although extensive experimental evidence on the process of liver regeneration exists, in humans, validation is largely missing. However, liver regeneration is critically affected by underlying liver disease. Within this project, we aimed to systematically assess early transcriptional changes during liver regeneration in humans and further assess how these processes differ in people with dysfunctional liver regeneration.

Methods

Blood samples of 154 patients and intraoperative tissue samples of 46 patients undergoing liver resection were collected and classified with regard to dysfunctional postoperative liver regeneration. Of those, a matched cohort of 21 patients were used for RNA sequencing. Samples were assessed for circulating cytokines, gene expression dynamics, intrahepatic neutrophil accumulation, and spatial transcriptomics.

Results

Individuals with dysfunctional liver regeneration demonstrated an aggravated transcriptional inflammatory response with higher intracellular adhesion molecule-1 induction. Increased induction of this critical leukocyte adhesion molecule was associated with increased intrahepatic neutrophil accumulation and activation upon induction of liver regeneration in individuals with dysfunctional liver regeneration. Comparing baseline gene expression profiles in individuals with and without dysfunctional liver regeneration, we found that dual-specificity phosphatase 4 (DUSP4) expression, a known critical regulator of intracellular adhesion molecule-1 expression in endothelial cells, was markedly reduced in patients with dysfunctional liver regeneration. Mimicking clinical risk factors for dysfunctional liver regeneration, we found liver sinusoidal endothelial cells of two liver disease models to have significantly reduced baseline levels of DUSP4.

Conclusions

Exploring the landscape of early transcriptional changes of human liver regeneration, we observed that people with dysfunctional regeneration experience overwhelming intrahepatic inflammation. Subclinical liver disease might account for DUSP4 reduction in liver sinusoidal endothelial cells, which ultimately primes the liver for an aggravated inflammatory response.

Impact and implications

Using a unique human biorepository, focused on liver regeneration (LR), we explored the landscape of circulating and tissue-level alterations associated with both functional and dysfunctional LR. In contrast to experimental animal models, people with dysfunctional LR demonstrated an aggravated transcriptional inflammatory response, higher intracellular adhesion molecule-1 (ICAM-1) induction, intrahepatic neutrophil accumulation and activation upon induction of LR. Although inflammatory responses appear rapidly after liver resection, people with dysfunctional LR have exaggerated inflammatory responses that appear to be related to decreased levels of LSEC DUSP4, challenging existing concepts of post-resectional LR.

Mitochondrial DNA variation in Alzheimer's disease reveals a unique microprotein called SHMOOSE

Mitochondrial DNA variants have previously associated with disease, but the underlying mechanisms have been largely elusive. Here, we report that mitochondrial SNP rs2853499 associated with Alzheimer's disease (AD), neuroimaging, and transcriptomics. We mapped rs2853499 to a novel mitochondrial small open reading frame called SHMOOSE with microprotein encoding potential. Indeed, we detected two unique SHMOOSE-derived peptide fragments in mitochondria by using mass spectrometry-the first unique mass spectrometry-based detection of a mitochondrial-encoded microprotein to date. Furthermore, cerebrospinal fluid (CSF) SHMOOSE levels in humans correlated with age, CSF tau, and brain white matter volume. We followed up on these genetic and biochemical findings by carrying out a series of functional experiments. SHMOOSE acted on the brain following intracerebroventricular administration, differentiated mitochondrial gene expression in multiple models, localized to mitochondria, bound the inner mitochondrial membrane protein mitofilin, and boosted mitochondrial oxygen consumption. Altogether, SHMOOSE has vast implications for the fields of neurobiology, Alzheimer's disease, and microproteins.

Diabetic kidney disease induces transcriptome alterations associated with angiogenesis activity in human mesenchymal stromal cells

BACKGROUND

Therapeutic interventions that optimize angiogenic activities may reduce rates of end-stage kidney disease, critical limb ischemia, and lower extremity amputations in individuals with diabetic kidney disease (DKD). Infusion of autologous mesenchymal stromal cells (MSC) is a promising novel therapy to rejuvenate vascular integrity. However, DKD-related factors, including hyperglycemia and uremia, might alter MSC angiogenic repair capacity in an autologous treatment approach.

METHODS

To explore the angiogenic activity of MSC in DKD, the transcriptome of adipose tissue-derived MSC obtained from DKD subjects was compared to age-matched controls without diabetes or kidney impairment. Next-generation RNA sequencing (RNA-seq) was performed on MSC (DKD n = 29; Controls n = 9) to identify differentially expressed (DE; adjusted p < 0.05, |log₂fold change|> 1) messenger RNA (mRNA) and microRNA (miRNA) involved in angiogenesis (GeneCards). Paracrine-mediated angiogenic repair capacity of MSC conditioned medium (MSCcm) was assessed in vitro using human umbilical vein endothelial cells incubated in high glucose and indoxyl sulfate for a hyperglycemic, uremic state.

RESULTS

RNA-seq analyses revealed 133 DE mRNAs (77 upregulated and 56 down-regulated) and 208 DE miRNAs (119 up- and 89 down-regulated) in DKD-MSC versus Control-MSC. Interestingly, miRNA let-7a-5p, which regulates angiogenesis and participates in DKD pathogenesis, interacted with 5 angiogenesis-associated mRNAs (transgelin/TAGLN, thrombospondin 1/THBS1, lysyl oxidase-like 4/LOXL4, collagen 4A1/COL4A1 and collagen 8A1/COL8A1). DKD-MSCcm incubation with injured endothelial cells improved tube formation capacity, enhanced migration, reduced adhesion molecules E-selectin, vascular cell adhesion molecule 1 and intercellular adhesion molecule 1 mRNA expression in endothelial cells. Moreover, angiogenic repair effects did not differ between treatment groups (DKD-MSCcm vs. Control-MSCcm).

CONCLUSIONS

MSC from individuals with DKD show angiogenic transcriptome alterations compared to age-matched controls. However, angiogenic repair potential may be preserved, supporting autologous MSC interventions to treat conditions requiring enhanced angiogenic activities such as DKD, diabetic foot ulcers, and critical limb ischemia.

Novel Genetic Determinants of Dental Maturation in Children

Dental occlusion requires harmonious development of teeth, jaws, and other elements of the craniofacial complex, which are regulated by environmental and genetic factors. We performed the first genome-wide association study (GWAS) on dental development (DD) using the Demirjian radiographic method. Radiographic assessments from participants of the Generation R Study (primary study population, N1 = 2,793; mean age of 9.8 y) were correlated with ~30 million genetic variants while adjusting for age, sex, and genomic principal components (proxy for population stratification). Variants associated with DD at genome-wide significant level (P < 5 × 10-8) mapped to 16q12.2 (IRX5) (lead variant rs3922616, B = 0.16; P = 2.2 × 10-8). We used Fisher's combined probability tests weighted by sample size to perform a meta-analysis (N = 14,805) combining radiographic DD at a mean age of 9.8 y from Generation R with data from a previous GWAS (N2 = 12,012) on number of teeth (NT) in infants used as proxy of DD at a mean age of 9.8 y (including the ALSPAC and NFBC1966). This GWAS meta-analysis revealed 3 novel loci mapping to 7p15.3 (IGF2BP3: P = 3.2 × 10-8), 14q13.3 (PAX9: P = 1.9 × 10-8), and 16q12.2 (IRX5: P = 1.2 × 10-9) and validated 8 previously reported NT loci. A polygenic allele score constructed from these 11 loci was associated with radiographic DD in an independent Generation R set of children (N = 703; B = 0.05, P = 0.004). Furthermore, profiling of the identified genes across an atlas of murine and human stem cells observed expression in the cells involved in the formation of bone and/or dental tissues (>0.3 frequency per kilobase of transcript per million mapped reads), likely reflecting functional specialization. Our findings provide biological insight into the polygenic architecture of the pediatric dental maturation process.

Zika virus alters osteogenic lineage progression of human mesenchymal stromal cells

Arboviruses target bone forming osteoblasts and perturb bone remodeling via paracrine factors. We previously reported that Zika virus (ZIKV) infection of early-stage human mesenchymal stromal cells (MSCs) inhibited the osteogenic lineage commitment of MSCs. To understand the physiological interplay between bone development and ZIKV pathogenesis, we employed a primary in vitro model to examine the biological responses of MSCs to ZIKV infection at different stages of osteogenesis. Precommitted MSCs were infected at the late stage of osteogenic stimulation (Day 7) with ZIKV (multiplicity of infection of 5). We observe that MSCs infected at the late stage of differentiation are highly susceptible to ZIKV infection similar to previous observations with early stage infected MSCs (Day 0). However, in contrast to ZIKV infection at the early stage of differentiation, infection at a later stage significantly elevates the key osteogenic markers and calcium content. Comparative RNA sequencing (RNA-seq) of early and late stage infected MSCs reveals that ZIKV infection alters the mRNA transcriptome during osteogenic induction of MSCs (1251 genes). ZIKV infection provokes a robust antiviral response at both stages of osteogenic differentiation as reflected by the upregulation of interferon responsive genes (n > 140). ZIKV infection enhances the expression of immune-related genes in early stage MSCs while increasing cell cycle genes in late stage MSCs. Remarkably, ZIKA infection in early stage MSCs also activates lipid metabolism-related pathways. In conclusion, ZIKV infection has differentiation stage-dependent effects on MSCs and this mechanistic understanding may permit the development of new therapeutic or preventative measures for bone-related effects of ZIKV infection.

G9a Modulates Lipid Metabolism in CD4 T Cells to Regulate Intestinal Inflammation

BACKGROUND & AIMS

Although T-cell intrinsic expression of G9a has been associated with murine intestinal inflammation, mechanistic insight into the role of this methyltransferase in human T-cell differentiation is ill defined, and manipulation of G9a function for therapeutic use against inflammatory disorders is unexplored.

METHODS

Human naive T cells were isolated from peripheral blood and differentiated in vitro in the presence of a G9a inhibitor (UNC0642) before being characterized via the transcriptome (RNA sequencing), chromatin accessibility (assay for transposase-accessible chromatin by sequencing), protein expression (cytometry by time of flight, flow cytometry), metabolism (mitochondrial stress test, ultrahigh performance liquid chromatography-tandem mas spectroscopy) and function (T-cell suppression assay). The in vivo role of G9a was assessed using 3 murine models.

RESULTS

We discovered that pharmacologic inhibition of G9a enzymatic function in human CD4 T cells led to spontaneous generation of FOXP3+ T cells (G9a-inibitors-T regulatory cells [Tregs]) in vitro that faithfully reproduce human Tregs, functionally and phenotypically. Mechanistically, G9a inhibition altered the transcriptional regulation of genes involved in lipid biosynthesis in T cells, resulting in increased intracellular cholesterol. Metabolomic profiling of G9a-inibitors-Tregs confirmed elevated lipid pathways that support Treg development through oxidative phosphorylation and enhanced lipid membrane composition. Pharmacologic G9a inhibition promoted Treg expansion in vivo upon antigen (gliadin) stimulation and ameliorated acute trinitrobenzene sulfonic acid-induced colitis secondary to tissue-specific Treg development. Finally, Tregs lacking G9a expression (G9a-knockout Tregs) remain functional chronically and can rescue T-cell transfer-induced colitis.

CONCLUSION

G9a inhibition promotes cholesterol metabolism in T cells, favoring a metabolic profile that facilitates Treg development in vitro and in vivo. Our data support the potential use of G9a inhibitors in the treatment of immune-mediated conditions including inflammatory bowel disease.

Trem2 H157Y increases soluble TREM2 production and reduces amyloid pathology

BACKGROUND

The rare p.H157Y variant of TREM2 (Triggering Receptor Expressed on Myeloid Cells 2) was found to increase Alzheimer's disease (AD) risk. This mutation is located at the cleavage site of TREM2 extracellular domain. Ectopic expression of TREM2-H157Y in HEK293 cells resulted in increased TREM2 shedding. However, the physiological outcomes of the TREM2 H157Y mutation remain unknown in the absence and presence of AD related pathologies.

METHODS

We generated a novel Trem2 H157Y knock-in mouse model through CRISPR/Cas9 technology and investigated the effects of Trem2 H157Y on TREM2 proteolytic processing, synaptic function, and AD-related amyloid pathologies by conducting biochemical assays, targeted mass spectrometry analysis of TREM2, hippocampal electrophysiology, immunofluorescent staining, in vivo micro-dialysis, and cortical bulk RNA sequencing.

RESULTS

Consistent with previous in vitro findings, Trem2 H157Y increases TREM2 shedding with elevated soluble TREM2 levels in the brain and serum. Moreover, Trem2 H157Y enhances synaptic plasticity without affecting microglial density and morphology, or TREM2 signaling. In the presence of amyloid pathology, Trem2 H157Y accelerates amyloid-β (Aβ) clearance and reduces amyloid burden, dystrophic neurites, and gliosis in two independent founder lines. Targeted mass spectrometry analysis of TREM2 revealed higher ratios of soluble to full-length TREM2-H157Y compared to wild-type TREM2, indicating that the H157Y mutation promotes TREM2 shedding in the presence of Aβ. TREM2 signaling was further found reduced in Trem2 H157Y homozygous mice. Transcriptomic profiling revealed that Trem2 H157Y downregulates neuroinflammation-related genes and an immune module correlated with the amyloid pathology.

CONCLUSION

Taken together, our findings suggest beneficial effects of the Trem2 H157Y mutation in synaptic function and in mitigating amyloid pathology. Considering the genetic association of TREM2 p.H157Y with AD risk, we speculate TREM2 H157Y in humans might increase AD risk through an amyloid-independent pathway, such as its effects on tauopathy and neurodegeneration which merit further investigation.

AT-101 Enhances the Antitumor Activity of Lenalidomide in Patients with Multiple Myeloma

Bcl-2 and Mcl-1 proteins play a role in multiple myeloma (MM) cell survival, for which targeted inhibitors are being developed. AT-101 is an oral drug, which disrupts Bcl-2 and Mcl-1 function, impedes mitochondrial bioenergetic processes and induces apoptosis in MM cells. When combined with lenalidomide and dexamethasone (Rd), AT-101 significantly reduced tumor burden in an in vivo xenograft model of MM. These data provided rationale for a phase I/II study to establish the effective dose of AT-101 in combination with Rd (ARd regimen) in relapsed/refractory MM. A total of 10 patients were enrolled, most with high-risk cytogenetics (80%) and prior stem cell transplant (70%). Three patients were lenalidomide-refractory, 2 were bortezomib-refractory and 3 were daratumumab-refractory. The ARd combination was well tolerated with most common grade 3/4 adverse events being cytopenia's. The overall response rate was 40% and clinical benefit rate was 90%. The median progression free survival was 14.9 months (95% CI 7.1-NE). Patients responsive to ARd showed a decrease in Bcl-2:Bim or Mcl-1:Noxa protein complexes, increased CD8+ T and NK cells and depletion of T and B-regulatory cells. The ARd regimen demonstrated an acceptable safety profile and promising efficacy in patients with relapsed/refractory MM prompting further investigation in additional patients.

Sonicate Fluid Cellularity Predicted by Transcriptomic Deconvolution Differentiates Infectious from Non-Infectious Arthroplasty Failure

BACKGROUND

Although cellularity is traditionally assessed morphologically, deep sequencing approaches being used for microorganism detection may be able to provide information about cellularity. We hypothesized that cellularity predicted using CIBERSORTx (Stanford University), a transcriptomic-based cellular deconvolution tool, would differentiate between infectious and non-infectious arthroplasty failure.

METHODS

CIBERSORTx-derived cellularity profiles of 93 sonicate fluid samples, including 53 from subjects who underwent failed arthroplasties due to periprosthetic joint infection (PJI) (abbreviated for the purpose of this study as PJIF) and 40 from subjects who had undergone non-infectious arthroplasty failure (abbreviated NIAF) that had been subjected to bulk RNA sequencing were evaluated.

RESULTS

Samples from PJIF and NIAF subjects were differentially clustered by principal component analysis based on the cellularity profile. Twelve of the 22 individual predicted cellular fractions were differentially expressed in the PJIF cases compared with the NIAF cases, including increased predicted neutrophils (mean and standard error, 9.73% ± 1.06% and 0.81% ± 0.60%), activated mast cells (17.12% ± 1.51% and 4.11% ± 0.44%), and eosinophils (1.96% ± 0.37% and 0.42% ± 0.21%), and decreased predicted M0 macrophages (21.33% ± 1.51% and 39.75% ± 2.45%), M2 macrophages (3.56% ± 0.52% and 8.70% ± 1.08%), and regulatory T cells (1.57% ± 0.23% and 3.20% ± 0.34%). The predicted total granulocyte fraction was elevated in the PJIF cases (32.97% ± 2.13% and 11.76% ± 1.61%), and the samples from the NIAF cases had elevated predicted total macrophage and monocyte (34.71% ± 1.71% and 55.34% ± 2.37%) and total B cell fractions (5.89% ± 0.30% and 8.62% ± 0.86%). Receiver operating characteristic curve analysis identified predicted total granulocytes, neutrophils, and activated mast cells as highly able to differentiate between the PJIF cases and the NIAF cases. Within the PJIF cases, the total granulocyte, total macrophage and monocyte, M0 macrophage, and M2 macrophage fractions were differentially expressed in Staphylococcus aureus compared with Staphylococcus epidermidis -associated samples. Within the NIAF cases, the predicted total B cell, naïve B cell, plasma cell, and M2 macrophage fractions were differentially expressed among different causes of failure.

CONCLUSIONS

CIBERSORTx can predict the cellularity of sonicate fluid using transcriptomic data, allowing for the evaluation of the underlying immune response during the PJIF and NIAF cases, without a need to phenotypically assess cell composition.

Exploring Gene Expression Profiles in Primary Central Nervous System Vasculitis

OBJECTIVE

This study was undertaken to explore the gene expression profile of primary central nervous system vasculitis (PCNSV).

METHODS

Brain specimens of 4 patients with granulomatous vasculitis (GV), 5 with lymphocytic vasculitis (LV), 4 with amyloid β-related angiitis (ABRA), and 4 normal controls were studied. RNA-sequencing was performed using the Illumina Hiseq-4,000 platform and the Illumina TruSeq Total-RNA library. Student t test and false discovery rate tests were performed for each of the differentially expressed transcripts. Ingenuity Pathway Analysis was used for the pathway expression analysis. CIBERSORT was used to estimate the abundances of different immune cell subsets in the tissues based on gene expression data.

RESULTS

Transcripts differentially expressed between PCNSV and normal brain indicated that endosomal, mitochondrial, and ribosome dysfunction, alterations in protein synthesis, and noncoding RNAs might be involved in PCNSV. Pathway analysis revealed the activation of dendritic cell maturation and antigen processing as well as neuroinflammation in PCNSV versus normal brain, whereas oxidative phosphorylation was inhibited. CIBERSORT estimation of immune cell subsets suggested that activated NK cells, M1 macrophages, memory B cells, and follicular helper T cells were likely to be more prevalent in PCNSV samples. Naïve CD4 T cells and monocytes were mainly estimated to be present in GV and ABRA. Plasma cell and γδ T-cell signatures were mainly found in LV and normal brain. GV showed higher levels of genes associated with macrophage activities and T cells. ABRA showed higher levels of long noncoding RNAs and miR-616. LV showed higher levels of genes encoding immunoglobulins.

INTERPRETATION

RNA sequencing confirmed PCNSV heterogeneity. ANN NEUROL 2023;93:120-130.

Pepsinogen/Proton Pump Co-Expression in Barrett's Esophageal Cells Induces Cancer-Associated Changes

EDUCATIONAL OBJECTIVE

At the conclusion of this presentation, participants should better understand the carcinogenic potential of pepsin and proton pump expression in Barrett's esophagus.

OBJECTIVE

Barrett's esophagus (BE) is a well-known risk factor for esophageal adenocarcinoma (EAC). Gastric H⁺ /K⁺ ATPase proton pump and pepsin expression has been demonstrated in some cases of BE; however, the contribution of local pepsin and proton pump expression to carcinogenesis is unknown. In this study, RNA sequencing was used to examine global transcriptomic changes in a BE cell line ectopically expressing pepsinogen and/or gastric H⁺ /K⁺ ATPase proton pumps.

STUDY DESIGN

In vitro translational.

METHODS

BAR-T, a human BE cell line devoid of expression of pepsinogen or proton pumps, was transduced by lentivirus-encoding pepsinogen (PGA5) and/or gastric proton pump subunits (ATP4A, ATP4B). Changes relative to the parental line were assessed by RNA sequencing.

RESULTS

Top canonical pathways associated with protein-coding genes differentially expressed in pepsinogen and/or proton pump expressing BAR-T cells included those involved in the tumor microenvironment and epithelial-mesenchymal transition. Top upstream regulators of coding transcripts included TGFB1 and ERBB2, which are associated with the pathogenesis and prognosis of BE and EAC. Top upstream regulators of noncoding transcripts included p300-CBP, I-BET-151, and CD93, which have previously described associations with EAC or carcinogenesis. The top associated disease of both coding and noncoding transcripts was cancer.

CONCLUSIONS

These data support the carcinogenic potential of pepsin and proton pump expression in BE and reveal molecular pathways affected by their expression. Further study is warranted to investigate the role of these pathways in carcinogenesis associated with BE.

LEVEL OF EVIDENCE

NA Laryngoscope, 133:59-69, 2023.

Cardiac epigenetic changes in VEGF signaling genes associate with myocardial microvascular rarefaction in experimental chronic kidney disease

Chronic kidney disease (CKD) is common in patients with heart failure and often results in left ventricular diastolic dysfunction (LVDD). However, the mechanisms responsible for cardiac damage in CKD-LVDD remain to be elucidated. Epigenetic alterations may impose long-lasting effects on cellular transcription and function, but their exact role in CKD-LVDD is unknown. We investigate whether changes in cardiac site-specific DNA methylation profiles might be implicated in cardiac abnormalities in CKD-LVDD. CKD-LVDD and normal control pigs (n = 6 each) were studied for 14 wk. Renal and cardiac hemodynamics were quantified by multidetector CT and echocardiography. In randomly selected pigs (n = 3/group), cardiac site-specific 5-methylcytosine (5mC) immunoprecipitation (MeDIP)- and mRNA-sequencing (seq) were performed, followed by integrated (MeDiP-seq/mRNA-seq analysis), and confirmatory ex vivo studies. MeDIP-seq analysis revealed 261 genes with higher (fold change > 1.4; P < 0.05) and 162 genes with lower (fold change < 0.7; P < 0.05) 5mC levels in CKD-LVDD versus normal pigs, which were primarily implicated in vascular endothelial growth factor (VEGF)-related signaling and angiogenesis. Integrated MeDiP-seq/mRNA-seq analysis identified a select group of VEGF-related genes in which 5mC levels were higher, but mRNA expression was lower in CKD-LVDD versus normal pigs. Cardiac VEGF signaling gene and VEGF protein expression were blunted in CKD-LVDD compared with controls and were associated with decreased subendocardial microvascular density. Cardiac epigenetic changes in VEGF-related genes are associated with impaired angiogenesis and cardiac microvascular rarefaction in swine CKD-LVDD. These observations may assist in developing novel therapies to ameliorate cardiac damage in CKD-LVDD.NEW & NOTEWORTHY Chronic kidney disease (CKD) often leads to left ventricular diastolic dysfunction (LVDD) and heart failure. Using a novel translational swine model of CKD-LVDD, we characterize the cardiac epigenetic landscape, identifying site-specific 5-methylcytosine changes in vascular endothelial growth factor (VEGF)-related genes associated with impaired angiogenesis and cardiac microvascular rarefaction. These observations shed light on the mechanisms of cardiac microvascular damage in CKD-LVDD and may assist in developing novel therapies for these patients.