Hypoxia-Inducible Factor 1α Stabilization Restores Epigenetic Control of Nitric Oxide Synthase 1 Expression and Reverses Gastroparesis in Female Diabetic Mice.

BACKGROUND & AIMS: Although depletion of neuronal nitric oxide synthase (NOS1)-expressing neurons contributes to gastroparesis, stimulating nitrergic signaling is not an effective therapy. We investigated whether hypoxia-inducible factor 1α (HIF1A), which is activated by high O2 consumption in central neurons, is a Nos1 transcription factor in enteric neurons and whether stabilizing HIF1A reverses gastroparesis. METHODS: Mice with streptozotocin-induced diabetes, human and mouse tissues, NOS1+ mouse neuroblastoma cells, and isolated nitrergic neurons were studied. Gastric emptying of solids and volumes were determined by breath test and single-photon emission computed tomography, respectively. Gene expression was analyzed by RNA-sequencing, microarrays, immunoblotting, and immunofluorescence. Epigenetic assays included chromatin immunoprecipitation sequencing (13 targets), chromosome conformation capture sequencing, and reporter assays. Mechanistic studies used Cre-mediated recombination, RNA interference, and clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated protein 9 (Cas9)-mediated epigenome editing. RESULTS: HIF1A signaling from physiological intracellular hypoxia was active in mouse and human NOS1+ myenteric neurons but reduced in diabetes. Deleting Hif1a in Nos1-expressing neurons reduced NOS1 protein by 50% to 92% and delayed gastric emptying of solids in female but not male mice. Stabilizing HIF1A with roxadustat (FG-4592), which is approved for human use, restored NOS1 and reversed gastroparesis in female diabetic mice. In nitrergic neurons, HIF1A up-regulated Nos1 transcription by binding and activating proximal and distal cis-regulatory elements, including newly discovered super-enhancers, facilitating RNA polymerase loading and pause-release, and by recruiting cohesin to loop anchors to alter chromosome topology. CONCLUSIONS: Pharmacologic HIF1A stabilization is a novel, translatable approach to restoring nitrergic signaling and treating diabetic gastroparesis. The newly recognized effects of HIF1A on chromosome topology may provide insights into physioxia- and ischemia-related organ function.

Interactive enhancer hubs (iHUBs) mediate transcriptional reprogramming and adaptive resistance in pancreatic cancer.

OBJECTIVE: Pancreatic ductal adenocarcinoma (PDAC) displays a remarkable propensity towards therapy resistance. However, molecular epigenetic and transcriptional mechanisms enabling this are poorly understood. In this study, we aimed to identify novel mechanistic approaches to overcome or prevent resistance in PDAC. DESIGN: We used in vitro and in vivo models of resistant PDAC and integrated epigenomic, transcriptomic, nascent RNA and chromatin topology data. We identified a JunD-driven subgroup of enhancers, called interactive hubs (iHUBs), which mediate transcriptional reprogramming and chemoresistance in PDAC. RESULTS: iHUBs display characteristics typical for active enhancers (H3K27ac enrichment) in both therapy sensitive and resistant states but exhibit increased interactions and production of enhancer RNA (eRNA) in the resistant state. Notably, deletion of individual iHUBs was sufficient to decrease transcription of target genes and sensitise resistant cells to chemotherapy. Overlapping motif analysis and transcriptional profiling identified the activator protein 1 (AP1) transcription factor JunD as a master transcription factor of these enhancers. JunD depletion decreased iHUB interaction frequency and transcription of target genes. Moreover, targeting either eRNA production or signaling pathways upstream of iHUB activation using clinically tested small molecule inhibitors decreased eRNA production and interaction frequency, and restored chemotherapy responsiveness in vitro and in vivo. Representative iHUB target genes were found to be more expressed in patients with poor response to chemotherapy compared with responsive patients. CONCLUSION: Our findings identify an important role for a subgroup of highly connected enhancers (iHUBs) in regulating chemotherapy response and demonstrate targetability in sensitisation to chemotherapy.

Bioinformatics and DNA-extraction strategies to reliably detect genetic variants from FFPE breast tissue samples.

BACKGROUND: Archived formalin fixed paraffin embedded (FFPE) samples are valuable clinical resources to examine clinically relevant morphology features and also to study genetic changes. However, DNA quality and quantity of FFPE samples are often sub-optimal, and resulting NGS-based genetics variant detections are prone to false positives. Evaluations of wet-lab and bioinformatics approaches are needed to optimize variant detection from FFPE samples. RESULTS: As a pilot study, we designed within-subject triplicate samples of DNA derived from paired FFPE and fresh frozen breast tissues to highlight FFPE-specific artifacts. For FFPE samples, we tested two FFPE DNA extraction methods to determine impact of wet-lab procedures on variant calling: QIAGEN QIAamp DNA Mini Kit ("QA"), and QIAGEN GeneRead DNA FFPE Kit ("QGR"). We also used negative-control (NA12891) and positive control samples (Horizon Discovery Reference Standard FFPE). All DNA sample libraries were prepared for NGS according to the QIAseq Human Breast Cancer Targeted DNA Panel protocol and sequenced on the HiSeq 4000. Variant calling and filtering were performed using QIAGEN Gene Globe Data Portal. Detailed variant concordance comparisons and mutational signature analysis were performed to investigate effects of FFPE samples compared to paired fresh frozen samples, along with different DNA extraction methods. In this study, we found that five times or more variants were called with FFPE samples, compared to their paired fresh-frozen tissue samples even after applying molecular barcoding error-correction and default bioinformatics filtering recommended by the vendor. We also found that QGR as an optimized FFPE-DNA extraction approach leads to much fewer discordant variants between paired fresh frozen and FFPE samples. Approximately 92% of the uniquely called FFPE variants were of low allelic frequency range (< 5%), and collectively shared a "C > T|G > A" mutational signature known to be representative of FFPE artifacts resulting from cytosine deamination. Based on control samples and FFPE-frozen replicates, we derived an effective filtering strategy with associated empirical false-discovery estimates. CONCLUSIONS: Through this study, we demonstrated feasibility of calling and filtering genetic variants from FFPE tissue samples using a combined strategy with molecular barcodes, optimized DNA extraction, and bioinformatics methods incorporating genomics context such as mutational signature and variant allelic frequency.

Indel detection from RNA-seq data: tool evaluation and strategies for accurate detection of actionable mutations.

Driver somatic mutations are a hallmark of a tumor that can be used for diagnosis and targeted therapy. Mutations are primarily detected from tumor DNA. As dynamic molecules of gene activities, transcriptome profiling by RNA sequence (RNA-seq) is becoming increasingly popular, which not only measures gene expression but also structural variations such as mutations and fusion transcripts. Although single-nucleotide variants (SNVs) can be easily identified from RNA-seq, intermediate long insertions/deletions (indels > 2 bases and less than sequence reads) cause significant challenges and are ignored by most RNA-seq analysis tools. This study evaluates commonly used RNA-seq analysis programs along with variant and somatic mutation callers in a series of data sets with simulated and known indels. The aim is to develop strategies for accurate indel detection. Our results show that the RNA-seq alignment is the most important step for indel identification and the evaluated programs have a wide range of sensitivity to map sequence reads with indels, from not at all to decently sensitive. The sensitivity is impacted by sequence read lengths. Most variant calling programs rely on hard evidence indels marked in the alignment and the programs with realignment may use soft-clipped reads for indel inferencing. Based on the observations, we have provided practical recommendations for indel detection when different RNA-seq aligners are used and demonstrated the best option with highly reliable results. With careful customization of bioinformatics algorithms, RNA-seq can be reliably used for both SNV and indel mutation detection that can be used for clinical decision-making.

High frequency electrical stimulation promotes expression of extracellular matrix proteins from human astrocytes.

Therapeutic benefits of deep brain stimulation (DBS), a neurosurgical treatment for certain movement disorders and other neurologic conditions, are well documented, but DBS mechanisms remain largely unexplained. DBS is thought to modulate pathological neural activity. However, although astrocytes, the most numerous cell type in the brain, play a significant role in neurotransmission, chemical homeostasis and synaptic plasticity, their role in DBS has not been fully examined. To investigate astrocytic function in DBS, we applied DBS-like high frequency electrical stimulation for 24 h to human astrocytes in vitro and analyzed single cell transcriptome mRNA profile. We found that DBS-like high frequency stimulation negatively impacts astrocyte metabolism and promotes the release of extracellular matrix (matricellular) proteins, including IGFBP3, GREM1, IGFBP5, THBS1, and PAPPA. Our results suggest that astrocytes are involved in the long-term modulation of extra cellular matrix environments and that they may influence persistent cell-to-cell interaction and help maintain neuromodulation over time.

RVboost: RNA-seq variants prioritization using a boosting method.

MOTIVATION: RNA-seq has become the method of choice to quantify genes and exons, discover novel transcripts and detect fusion genes. However, reliable variant identification from RNA-seq data remains challenging because of the complexities of the transcriptome, the challenges of accurately mapping exon boundary spanning reads and the bias introduced during the sequencing library preparation. METHOD: We developed RVboost, a novel method specific for RNA variant prioritization. RVboost uses several attributes unique in the process of RNA library preparation, sequencing and RNA-seq data analyses. It uses a boosting method to train a model of 'good quality' variants using common variants from HapMap, and prioritizes and calls the RNA variants based on the trained model. We packaged RVboost in a comprehensive workflow, which integrates tools of variant calling, annotation and filtering. RESULTS: RVboost consistently outperforms the variant quality score recalibration from the Genome Analysis Tool Kit and the RNA-seq variant-calling pipeline SNPiR in 12 RNA-seq samples using ground-truth variants from paired exome sequencing data. Several RNA-seq-specific attributes were identified as critical to differentiate true and false variants, including the distance of the variant positions to exon boundaries, and the percent of the reads supporting the variant in the first six base pairs. The latter identifies false variants introduced by the random hexamer priming during the library construction. AVAILABILITY AND IMPLEMENTATION: The RVboost package is implemented to readily run in Mac or Linux environments. The software and user manual are available at http://bioinformaticstools.mayo.edu/research/rvboost/.

PatternCNV: a versatile tool for detecting copy number changes from exome sequencing data.

MOTIVATION: Exome sequencing (exome-seq) data, which are typically used for calling exonic mutations, have also been utilized in detecting DNA copy number variations (CNVs). Despite the existence of several CNV detection tools, there is still a great need for a sensitive and an accurate CNV-calling algorithm with built-in QC steps, and does not require a paired reference for each sample. RESULTS: We developed a novel method named PatternCNV, which (i) accounts for the read coverage variations between exons while leveraging the consistencies of this variability across different samples; (ii) reduces alignment BAM files to WIG format and therefore greatly accelerates computation; (iii) incorporates multiple QC measures designed to identify outlier samples and batch effects; and (iv) provides a variety of visualization options including chromosome, gene and exon-level views of CNVs, along with a tabular summarization of the exon-level CNVs. Compared with other CNV-calling algorithms using data from a lymphoma exome-seq study, PatternCNV has higher sensitivity and specificity. AVAILABILITY AND IMPLEMENTATION: The software for PatternCNV is implemented using Perl and R, and can be used in Mac or Linux environments. Software and user manual are available at http://bioinformaticstools.mayo.edu/research/patterncnv/, and R package at https://github.com/topsoil/patternCNV/.

HiChIP: a high-throughput pipeline for integrative analysis of ChIP-Seq data.

BACKGROUND: Chromatin immunoprecipitation (ChIP) followed by next-generation sequencing (ChIP-Seq) has been widely used to identify genomic loci of transcription factor (TF) binding and histone modifications. ChIP-Seq data analysis involves multiple steps from read mapping and peak calling to data integration and interpretation. It remains challenging and time-consuming to process large amounts of ChIP-Seq data derived from different antibodies or experimental designs using the same approach. To address this challenge, there is a need for a comprehensive analysis pipeline with flexible settings to accelerate the utilization of this powerful technology in epigenetics research. RESULTS: We have developed a highly integrative pipeline, termed HiChIP for systematic analysis of ChIP-Seq data. HiChIP incorporates several open source software packages selected based on internal assessments and published comparisons. It also includes a set of tools developed in-house. This workflow enables the analysis of both paired-end and single-end ChIP-Seq reads, with or without replicates for the characterization and annotation of both punctate and diffuse binding sites. The main functionality of HiChIP includes: (a) read quality checking; (b) read mapping and filtering; (c) peak calling and peak consistency analysis; and (d) result visualization. In addition, this pipeline contains modules for generating binding profiles over selected genomic features, de novo motif finding from transcription factor (TF) binding sites and functional annotation of peak associated genes. CONCLUSIONS: HiChIP is a comprehensive analysis pipeline that can be configured to analyze ChIP-Seq data derived from varying antibodies and experiment designs. Using public ChIP-Seq data we demonstrate that HiChIP is a fast and reliable pipeline for processing large amounts of ChIP-Seq data.

CAP-miRSeq: a comprehensive analysis pipeline for microRNA sequencing data.

BACKGROUND: miRNAs play a key role in normal physiology and various diseases. miRNA profiling through next generation sequencing (miRNA-seq) has become the main platform for biological research and biomarker discovery. However, analyzing miRNA sequencing data is challenging as it needs significant amount of computational resources and bioinformatics expertise. Several web based analytical tools have been developed but they are limited to processing one or a pair of samples at time and are not suitable for a large scale study. Lack of flexibility and reliability of these web applications are also common issues. RESULTS: We developed a Comprehensive Analysis Pipeline for microRNA Sequencing data (CAP-miRSeq) that integrates read pre-processing, alignment, mature/precursor/novel miRNA detection and quantification, data visualization, variant detection in miRNA coding region, and more flexible differential expression analysis between experimental conditions. According to computational infrastructure, users can install the package locally or deploy it in Amazon Cloud to run samples sequentially or in parallel for a large number of samples for speedy analyses. In either case, summary and expression reports for all samples are generated for easier quality assessment and downstream analyses. Using well characterized data, we demonstrated the pipeline's superior performances, flexibility, and practical use in research and biomarker discovery. CONCLUSIONS: CAP-miRSeq is a powerful and flexible tool for users to process and analyze miRNA-seq data scalable from a few to hundreds of samples. The results are presented in the convenient way for investigators or analysts to conduct further investigation and discovery.

Lactic acid induces transcriptional repression of macrophage inflammatory response via histone acetylation.

Lactic acid has emerged as an important modulator of immune cell function. It can be produced by both gut microbiota and the host metabolism at homeostasis and during disease states. The production of lactic acid in the gut microenvironment is vital for tissue homeostasis. In the present study, we examined how lactic acid integrates cellular metabolism to shape the epigenome of macrophages during pro-inflammatory response. We found that lactic acid serves as a primary fuel source to promote histone H3K27 acetylation, which allows the expression of immunosuppressive gene program including Nr4a1. Consequently, macrophage pro-inflammatory function was transcriptionally repressed. Furthermore, the histone acetylation induced by lactic acid promotes a form of long-term immunosuppression ("trained immunosuppression"). Pre-exposure to lactic acid induces lipopolysaccharide tolerance. These findings thus indicate that lactic acid sensing and its effect on chromatin remodeling in macrophages represent a key homeostatic mechanism that can provide a tolerogenic tissue microenvironment.

Intestinal epigenomic alterations are associated with a dysregulated nutrient absorption phenotype in obesity.

Obesity is an epidemic with myriad health effects, but little is understood regarding individual obese phenotypes and how they may respond to therapy. Epigenetic changes associated with obesity have been detected in blood, liver, pancreas, and adipose tissues. Previous work found that dietary glucose hyperabsorption occurs in some obese subjects, but detailed transcriptional or epigenomic features of the intestine associated with this phenotype are unknown. This study evaluated differentially expressed genes and relative chromatin accessibility in intestinal organoids established from donors classified as lean, obese, or obese hyperabsorptive by body mass index and glucose transport assays. Transcriptomic analysis indicated that obese hyperabsorptive subjects have significantly upregulated dietary nutrient absorption proteins and downregulated type I interferon targets. Chromatin accessibility and transcription factor footprinting suggested that enhanced binding of HNF4G promotes the obese hyperabsorption phenotype. Quantitative PCR assessment in a larger subject cohort suggested that intestinal epithelial expression of CUBN, GIP, and SLC2A5 have high correlation with hyperabsorption. The obese hyperabsorption phenotype is characterized by transcriptional changes that support increased nutrient uptake and may be driven by differentially accessible chromatin. Recognizing unique intestinal phenotypes in obesity provides new perspective in considering therapeutic targets and options to manage the disease.

Obesity-driven mitochondrial dysfunction in human adipose tissue-derived mesenchymal stem/stromal cells involves epigenetic changes.

Obesity exacerbates tissue degeneration and compromises the integrity and reparative potential of mesenchymal stem/stromal cells (MSCs), but the underlying mechanisms have not been sufficiently elucidated. Mitochondria modulate the viability, plasticity, proliferative capacity, and differentiation potential of MSCs. We hypothesized that alterations in the 5-hydroxymethylcytosine (5hmC) profile of mitochondria-related genes may mediate obesity-driven dysfunction of human adipose-derived MSCs. MSCs were harvested from abdominal subcutaneous fat of obese and age/sex-matched non-obese subjects (n = 5 each). The 5hmC profile and expression of nuclear-encoded mitochondrial genes were examined by hydroxymethylated DNA immunoprecipitation sequencing (h MeDIP-seq) and mRNA-seq, respectively. MSC mitochondrial structure (electron microscopy) and function, metabolomics, proliferation, and neurogenic differentiation were evaluated in vitro, before and after epigenetic modulation. hMeDIP-seq identified 99 peaks of hyper-hydroxymethylation and 150 peaks of hypo-hydroxymethylation in nuclear-encoded mitochondrial genes from Obese- versus Non-obese-MSCs. Integrated hMeDIP-seq/mRNA-seq analysis identified a select group of overlapping (altered levels of both 5hmC and mRNA) nuclear-encoded mitochondrial genes involved in ATP production, redox activity, cell proliferation, migration, fatty acid metabolism, and neuronal development. Furthermore, Obese-MSCs exhibited decreased mitochondrial matrix density, membrane potential, and levels of fatty acid metabolites, increased superoxide production, and impaired neuronal differentiation, which improved with epigenetic modulation. Obesity elicits epigenetic changes in mitochondria-related genes in human adipose-derived MSCs, accompanied by structural and functional changes in their mitochondria and impaired fatty acid metabolism and neurogenic differentiation capacity. These observations may assist in developing novel therapies to preserve the potential of MSCs for tissue repair and regeneration in obese individuals.

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.

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.

Performance comparisons of methylation and structural variants from low-input whole-genome methylation sequencing.

Aim: Whole-genome methylation sequencing carries both DNA methylation and structural variant information (single nucleotide variant [SNV]; copy number variant [CNV]); however, limited data is available on the reliability of obtaining this information simultaneously from low-input DNA using various library preparation and sequencing protocols. Methods: A HapMap NA12878 sample was sequenced with three protocols (EM-sequencing, QIA-sequencing and Swift-sequencing) and their performance was compared on CpG methylation measurement and SNV and CNV detection. Results: At low DNA input (10-25 ng), EM-sequencing was superior in almost all metrics except CNV detection where all protocols were similar. EM-sequencing captured the highest number of CpGs and true SNVs. Conclusion: EM-sequencing is suitable to detect methylation, SNVs and CNVs from single sequencing with low-input DNA.

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.

Epigenomics may begin to explain in vitro differential response to hypomethylating agents in MMR-D hypermethylated endometrial cancer.

This work examines differences in chromatin accessibility, methylation, and response to DNA hypomethylating agents between mismatch repair-deficient and non-mismatch repair-deficient endometrial cancer. Next-generation sequencing of a stage 1B, grade 2 endometrioid endometrial cancer tumor revealed microsatellite instability and a variant of unknown significance in POLE along with global and MLH1 hypermethylation. Inhibition of viability by decitabine in the study and comparison tumors was minimal, as shown by an inhibitory effect of 0 and 17.9, respectively. Conversely, the inhibitory effect of azacitidine on the study tumor was more pronounced, at 72.8 versus 41.2. In vitro, mismatch repair-deficient endometrial cancer with MLH1 hypermethylation respond better to DNA methyltransferase inhibition by azacytidine (DNA/RNA inhibition), than to decitabine (DNA-only inhibition). Additional large studies are needed to substantiate our findings.

Human transcriptomic response to periprosthetic joint infection.

Periprosthetic joint infection (PJI), a devastating complication of total joint replacement, is of incompletely understood pathogenesis and may sometimes be challenging to clinically distinguish from other causes of arthroplasty failure. We characterized human gene expression in 93 specimens derived from surfaces of resected arthroplasties, comparing transcriptomes of subjects with infection- versus non-infection-associated arthroplasty failure. Differential gene expression analysis confirmed 28 previously reported potential biomarkers of PJI, including bactericidal/permeability increasing protein (BPI), cathelicidin antimicrobial peptide (CAMP), C-C-motif chemokine ligand 3 (CCL3), 4(CCL4) and C-X-C-motif chemokine ligand 2 (CXCL2), colony stimulating factor 2 receptor beta (CSF2RB), colony stimulating factor 3 (CSF3), alpha-defensin (DEFA4), Fc fragment of IgG receptor 1B (CD64B), intercellular adhesion molecule 1 (ICAM1), interferon gamma (IFNG), interleukin 13 receptor subunit alpha 2 (IL13RA2), interleukin 17D (IL17D), interleukin 1 (IL1A, IL1B, IL1RN), interleukin 2 receptors (IL2RA, IL2RG), interleukin 5 receptor (IL5RA), interleukin 6 (IL6), interleukin 8 (IL8), lipopolysaccharide binding protein (LBP), lipocalin (LCN2), lactate dehydrogenase C (LDHC), lactotransferrin (LTF), matrix metallopeptidase 3 (MMP3), peptidase inhibitor 3 (PI3), and vascular endothelial growth factor A (VEGFA), and identified three novel molecules of potential diagnostic use for detection of PJI, namely C-C-motif chemokine ligand CCL20, coagulation factor VII (F7), and B cell receptor FCRL4. Comparative analysis of infections caused by staphylococci versus bacteria other than staphylococci and Staphylococcus aureus versus Staphylococcus epidermidis showed elevated expression of interleukin 13 (IL13), IL17D, and MMP3 in staphylococcal infections, and of IL1B, IL8, and platelet factor PF4V1 in S. aureus compared to S. epidermidis infections. Pathway analysis of over-represented genes suggested activation of host immune response and cellular maintenance and repair functions in response to invasion of infectious agents. The data presented provides new potential targets for diagnosis of PJI and for differentiation of PJI caused by different infectious agents.

Dysfunctional ERG signaling drives pulmonary vascular aging and persistent fibrosis.

Vascular dysfunction is a hallmark of chronic diseases in elderly. The contribution of the vasculature to lung repair and fibrosis is not fully understood. Here, we performed an epigenetic and transcriptional analysis of lung endothelial cells (ECs) from young and aged mice during the resolution or progression of bleomycin-induced lung fibrosis. We identified the transcription factor ETS-related gene (ERG) as putative orchestrator of lung capillary homeostasis and repair, and whose function is dysregulated in aging. ERG dysregulation is associated with reduced chromatin accessibility and maladaptive transcriptional responses to injury. Loss of endothelial ERG enhances paracrine fibroblast activation in vitro, and impairs lung fibrosis resolution in young mice in vivo. scRNA-seq of ERG deficient mouse lungs reveales transcriptional and fibrogenic abnormalities resembling those associated with aging and human lung fibrosis, including reduced number of general capillary (gCap) ECs. Our findings demonstrate that lung endothelial chromatin remodeling deteriorates with aging leading to abnormal transcription, vascular dysrepair, and persistent fibrosis following injury.

In vivo lentiviral vector gene therapy to cure hereditary tyrosinemia type 1 and prevent development of precancerous and cancerous lesions.

Conventional therapy for hereditary tyrosinemia type-1 (HT1) with 2-(2-nitro-4-trifluoromethylbenzoyl)-1,3-cyclohexanedione (NTBC) delays and in some cases fails to prevent disease progression to liver fibrosis, liver failure, and activation of tumorigenic pathways. Here we demonstrate cure of HT1 by direct, in vivo administration of a therapeutic lentiviral vector targeting the expression of a human fumarylacetoacetate hydrolase (FAH) transgene in the porcine model of HT1. This therapy is well tolerated and provides stable long-term expression of FAH in pigs with HT1. Genomic integration displays a benign profile, with subsequent fibrosis and tumorigenicity gene expression patterns similar to wild-type animals as compared to NTBC-treated or diseased untreated animals. Indeed, the phenotypic and genomic data following in vivo lentiviral vector administration demonstrate comparative superiority over other therapies including ex vivo cell therapy and therefore support clinical application of this approach.