A High-Throughput PIXUL-Matrix-Based Toolbox to Profile Frozen and Formalin-Fixed Paraffin-Embedded Tissues Multiomes.

Large-scale high-dimensional multiomics studies are essential to unravel molecular complexity in health and disease. We developed an integrated system for tissue sampling (CryoGrid), analytes preparation (PIXUL), and downstream multiomic analysis in a 96-well plate format (Matrix), MultiomicsTracks96, which we used to interrogate matched frozen and formalin-fixed paraffin-embedded (FFPE) mouse organs. Using this system, we generated 8-dimensional omics data sets encompassing 4 molecular layers of intracellular organization: epigenome (H3K27Ac, H3K4m3, RNA polymerase II, and 5mC levels), transcriptome (messenger RNA levels), epitranscriptome (m6A levels), and proteome (protein levels) in brain, heart, kidney, and liver. There was a high correlation between data from matched frozen and FFPE organs. The Segway genome segmentation algorithm applied to epigenomic profiles confirmed known organ-specific superenhancers in both FFPE and frozen samples. Linear regression analysis showed that proteomic profiles, known to be poorly correlated with transcriptomic data, can be more accurately predicted by the full suite of multiomics data, compared with using epigenomic, transcriptomic, or epitranscriptomic measurements individually.

CryoGrid-PIXUL-RNA: high throughput RNA isolation platform for tissue transcript analysis.

BACKGROUND: Disease molecular complexity requires high throughput workflows to map disease pathways through analysis of vast tissue repositories. Great progress has been made in tissue multiomics analytical technologies. To match the high throughput of these advanced analytical platforms, we have previously developed a multipurpose 96-well microplate sonicator, PIXUL, that can be used in multiple workflows to extract analytes from cultured cells and tissue fragments for various downstream molecular assays. And yet, the sample preparation devices, such as PIXUL, along with the downstream multiomics analytical capabilities have not been fully exploited to interrogate tissues because storing and sampling of such biospecimens remain, in comparison, inefficient. RESULTS: To mitigate this tissue interrogation bottleneck, we have developed a low-cost user-friendly system, CryoGrid, to catalog, cryostore and sample tissue fragments. TRIzol is widely used to isolate RNA but it is labor-intensive, hazardous, requires fume-hoods, and is an expensive reagent. Columns are also commonly used to extract RNA but they involve many steps, are prone to human errors, and are also expensive. Both TRIzol and column protocols use test tubes. We developed a microplate PIXUL-based TRIzol-free and column-free RNA isolation protocol that uses a buffer containing proteinase K (PK buffer). We have integrated the CryoGrid system with PIXUL-based PK buffer, TRIzol, and PureLink column methods to isolate RNA for gene-specific qPCR and genome-wide transcript analyses. CryoGrid-PIXUL, when integrated with either PK buffer, TRIzol or PureLink column RNA isolation protocols, yielded similar transcript profiles in frozen organs (brain, heart, kidney and liver) from a mouse model of sepsis. CONCLUSIONS: RNA isolation using the CryoGrid-PIXUL system combined with the 96-well microplate PK buffer method offers an inexpensive user-friendly high throughput workflow to study transcriptional responses in tissues in health and disease as well as in therapeutic interventions.

Single cell epigenetic visualization assay.

Characterization of the epigenetic status of individual cells remains a challenge. Current sequencing approaches have limited coverage, and it is difficult to assign an epigenetic status to the transcription state of individual gene alleles in the same cell. To address these limitations, a targeted microscopy-based epigenetic visualization assay (EVA) was developed for detection and quantification of epigenetic marks at genes of interest in single cells. The assay is based on an in situ biochemical reaction between an antibody-conjugated alkaline phosphatase bound to the epigenetic mark of interest, and a 5'-phosphorylated fluorophore-labeled DNA oligo tethered to a target gene by gene-specific oligonucleotides. When the epigenetic mark is present at the gene, phosphate group removal by the phosphatase protects the oligo from λ-exonuclease activity providing a quantitative fluorescent readout. We applied EVA to measure 5-methylcytosine (5mC) and H3K9Ac levels at different genes and the HIV-1 provirus in human cell lines. To link epigenetic marks to gene transcription, EVA was combined with RNA-FISH. Higher 5mC levels at the silenced compared to transcribed XIST gene alleles in female somatic cells validated this approach and demonstrated that EVA can be used to relate epigenetic marks to the transcription status of individual gene alleles.

PIXUL-ChIP: integrated high-throughput sample preparation and analytical platform for epigenetic studies.

Chromatin immunoprecipitation (ChIP) is the most widely used approach for identification of genome-associated proteins and their modifications. We have previously introduced a microplate-based ChIP platform, Matrix ChIP, where the entire ChIP procedure is done on the same plate without sample transfers. Compared to conventional ChIP protocols, the Matrix ChIP assay is faster and has increased throughput. However, even with microplate ChIP assays, sample preparation and chromatin fragmentation (which is required to map genomic locations) remains a major bottleneck. We have developed a novel technology (termed 'PIXUL') utilizing an array of ultrasound transducers for simultaneous shearing of samples in standard 96-well microplates. We integrated PIXUL with Matrix ChIP ('PIXUL-ChIP'), that allows for fast, reproducible, low-cost and high-throughput sample preparation and ChIP analysis of 96 samples (cell culture or tissues) in one day. Further, we demonstrated that chromatin prepared using PIXUL can be used in an existing ChIP-seq workflow. Thus, the high-throughput capacity of PIXUL-ChIP provides the means to carry out ChIP-qPCR or ChIP-seq experiments involving dozens of samples. Given the complexity of epigenetic processes, the use of PIXUL-ChIP will advance our understanding of these processes in health and disease, as well as facilitate screening of epigenetic drugs.

Selective Extracellular Signal-Regulated Kinase 1/2 (ERK1/2) Inhibition by the SCH772984 Compound Attenuates In Vitro and In Vivo Inflammatory Responses and Prolongs Survival in Murine Sepsis Models.

Sepsis is the leading cause of death in intensive care units worldwide. Current treatments of sepsis are largely supportive and clinical trials using specific pharmacotherapy for sepsis have failed to improve outcomes. Here, we used the lipopolysaccharide (LPS)-stimulated mouse RAW264.7 cell line and AlphaLisa assay for TNFa as a readout to perform a supervised drug repurposing screen for sepsis treatment with compounds targeting epigenetic enzymes, including kinases. We identified the SCH772984 compound, an extracellular signal-regulated kinase (ERK) 1/2 inhibitor, as an effective blocker of TNFa production in vitro. RNA-Seq of the SCH772984-treated RAW264.7 cells at 1, 4, and 24 h time points of LPS challenge followed by functional annotation of differentially expressed genes highlighted the suppression of cellular pathways related to the immune system. SCH772984 treatment improved survival in the LPS-induced lethal endotoxemia and cecal ligation and puncture (CLP) mouse models of sepsis, and reduced plasma levels of Ccl2/Mcp1. Functional analyses of RNA-seq datasets for kidney, lung, liver, and heart tissues from SCH772984-treated animals collected at 6 h and 12 h post-CLP revealed a significant downregulation of pathways related to the immune response and platelets activation but upregulation of the extracellular matrix organization and retinoic acid signaling pathways. Thus, this study defined transcriptome signatures of SCH772984 action in vitro and in vivo, an agent that has the potential to improve sepsis outcome.

DNA methylation yields epigenetic clues into the diabetic nephropathy of Pima Indians.

Environmental factors drive epigenetic programming. DNA methylation is the best studied modification transmitting epigenetic information. A study by Qiu et al. examined potential epigenetic roots for the decline of renal function in Pima Indians. A genomewide survey of blood leukocytes uncovered differentially methylated DNA sites in regulatory regions of genes associated with chronic kidney disease. This longitudinal study provides the first clues on epigenetic links between environmental factors and a high prevalence of diabetic kidney disease in Pima Indians.

Regulation of ribosomal RNA expression across the lifespan is fine-tuned by maternal diet before implantation.

Cells and organisms respond to nutrient deprivation by decreasing global rates of transcription, translation and DNA replication. To what extent such changes can be reversed is largely unknown. We examined the effect of maternal dietary restriction on RNA synthesis in the offspring. Low protein diet fed either throughout gestation or for the preimplantation period alone reduced cellular RNA content across fetal somatic tissues during challenge and increased it beyond controls in fetal and adult tissues after challenge release. Changes in transcription of ribosomal RNA, the major component of cellular RNA, were responsible for this phenotype as evidenced by matching alterations in RNA polymerase I density and DNA methylation at ribosomal DNA loci. Cellular levels of the ribosomal transcription factor Rrn3 mirrored the rRNA expression pattern. In cell culture experiments, Rrn3 overexpression reduced rDNA methylation and increased rRNA expression; the converse occurred after inhibition of Rrn3 activity. These observations define novel mechanism where poor nutrition before implantation irreversibly alters basal rates of rRNA transcription thereafter in a process mediated by rDNA methylation and Rrn3 factor.

Quadruplex formation by both G-rich and C-rich DNA strands of the C9orf72 (GGGGCC)8•(GGCCCC)8 repeat: effect of CpG methylation.

Unusual DNA/RNA structures of the C9orf72 repeat may participate in repeat expansions or pathogenesis of amyotrophic lateral sclerosis and frontotemporal dementia. Expanded repeats are CpG methylated with unknown consequences. Typically, quadruplex structures form by G-rich but not complementary C-rich strands. Using CD, UV and electrophoresis, we characterized the structures formed by (GGGGCC)8 and (GGCCCC)8 strands with and without 5-methylcytosine (5mCpG) or 5-hydroxymethylcytosine (5hmCpG) methylation. All strands formed heterogenous mixtures of structures, with features of quadruplexes (at pH 7.5, in K(+), Na(+) or Li(+)), but no feature typical of i-motifs. C-rich strands formed quadruplexes, likely stabilized by G•C•G•C-tetrads and C•C•C•C-tetrads. Unlike G•G•G•G-tetrads, some G•C•G•C-tetrad conformations do not require the N7-Guanine position, hence C9orf72 quadruplexes still formed when N7-deazaGuanine replace all Guanines. 5mCpG and 5hmCpG increased and decreased the thermal stability of these structures. hnRNPK, through band-shift analysis, bound C-rich but not G-rich strands, with a binding preference of unmethylated > 5hmCpG > 5mCpG, where methylated DNA-protein complexes were retained in the wells, distinct from unmethylated complexes. Our findings suggest that for C-rich sequences interspersed with G-residues, one must consider quadruplex formation and that methylation of quadruplexes may affect epigenetic processes.

Heterogeneity of epigenetic changes at ischemia/reperfusion- and endotoxin-induced acute kidney injury genes.

Aberrant gene expression is a molecular hallmark of acute kidney injury (AKI). As epigenetic processes control gene expression in a cell- and environment-defined manner, understanding the epigenetic pathways that regulate genes altered by AKI may open vital new insights into the complexities of disease pathogenesis and identify possible therapeutic targets. Here we used matrix chromatin immunoprecipitation and integrative analysis to study 20 key permissive and repressive epigenetic histone marks at transcriptionally induced Tnf, Ngal, Kim-1, and Icam-1 genes in mouse models of AKI; unilateral renal ischemia/reperfusion, lipopolysaccharide (LPS), and their synergistically injurious combination. Results revealed unexpected heterogeneity of transcriptional and epigenetic responses. Tnf and Ngal were transcriptionally upregulated in response to both treatments individually, and to combination treatment. Kim-1 was induced by ischemia/reperfusion and Icam-1 by LPS only. Epigenetic alterations at these genes exhibited distinct time-dependent changes that shared some similarities, such as reduction in repressive histone modifications, and also had major ischemia/reperfusion versus endotoxin differences. Thus, diversity of changes at AKI genes in response to different insults indicates involvement of several epigenetic pathways. This could be exploited pharmacologically through rational-drug design to alter the course and improve clinical outcomes of this syndrome.

Beads-free protein immunoprecipitation for a mass spectrometry-based interactome and posttranslational modifications analysis.

BACKGROUND: Protein immunoprecipitation (IP) coupled with MS provides means to interrogate protein complexes and their posttranslational modifications (PTMs). In a typical protein IP assay antibodies are conjugated to protein A/G beads requiring large amounts of antibodies, tube transfers and centrifugations. RESULTS: As an alternative, we present Matrix-IP, beads-free microplate-based platform with surface-immobilized antibodies. Assay utilizes standard 96-well polypropylene PCR plates that are laboratory-fabricated with UV-C light and then protein A/G coated prior to IP reaction. We demonstrate application of Matrix-IP platform in MS analysis of heterogeneous nuclear ribonucleoprotein K (hnRNP K) interactome and PTMs. CONCLUSION: Matrix-IP is time-saving, easy to use high throughput method adaptable for low sample amounts and automation.

Experimental acute lung injury induces multi-organ epigenetic modifications in key angiogenic genes implicated in sepsis-associated endothelial dysfunction.

INTRODUCTION: The Tie2/angiopoietin (Tie2/Ang) and vascular endothelial growth factor receptor-ligand systems (VEGFR/VEGF) are recognized to play important roles in the regulation of microvascular endothelial function. Downregulation of these genes during sepsis has been implicated in the pathogenesis of sepsis-related microvascular leak and multiple organ dysfunction syndrome. Mechanisms responsible for dysregulation of angiogenic genes in sepsis are poorly defined. METHODS: Western blot, reverse transcription-polymerase chain reaction, and multiplex chromatin immunoprecipitation platform (Matrix ChIP) were used to investigate serum albumin leak, changes in gene expression, and associated epigenetic alterations in a murine model of acute lung injury-induced sepsis (ALI-sepsis). RESULTS: Experimental ALI-sepsis induced microvascular leak and downregulation of expression of Angpt1 (Ang1), Tek (Tie2), and Kdr (Vegfr2 or Flk-1) genes in the lung, kidney, and liver. These changes correlate with a decrease in RNA polymerase II density at these genes, and the greatest response was observed in the lung. ALI-sepsis reduced levels of transcription-permissive histone H3 lysine acetylation (H3KAc) at these loci in all examined tissues. Decreases in permissive H3K4m3 and H3Km2 marks were detected only in the lung. In contrast, only minimal alterations in transcription-repressive histone modifications (H3K27m3, H3K9m2, H3K9m3, and H4K20m3) were observed in all tissues. CONCLUSIONS: Our results demonstrate that decreases in transcription-permissive, but not increases in transcription-repressive, histone modifications at Angpt1, Tek, and Kdr are a systemic, rather than a lung-restricted, response, involving key end-organs in experimental ALI-sepsis. Given that ventilator-associated pneumonia is a major cause of sepsis in critically ill patients, elucidation of mechanisms mediating epigenetic alterations during sepsis provides fundamental new insights into the pathogenesis of sepsis-induced microvascular leak and subsequent end-organ injury/dysfunction.

Heterogeneous nuclear ribonucleoprotein (HnRNP) K genome-wide binding survey reveals its role in regulating 3'-end RNA processing and transcription termination at the early growth response 1 (EGR1) gene through XRN2 exonuclease.

The heterogeneous nuclear ribonucleoprotein K (hnRNPK) is a nucleic acid-binding protein that acts as a docking platform integrating signal transduction pathways to nucleic acid-related processes. Given that hnRNPK could be involved in other steps that compose gene expression the definition of its genome-wide occupancy is important to better understand its role in transcription and co-transcriptional processes. Here, we used chromatin immunoprecipitation followed by deep sequencing (ChIP-Seq) to analyze the genome-wide hnRNPK-DNA interaction in colon cancer cell line HCT116. 9.1/3.6 and 7.0/3.4 million tags were sequenced/mapped, then 1809 and 642 hnRNPK binding sites were detected in quiescent and 30-min serum-stimulated cells, respectively. The inspection of sequencing tracks revealed inducible hnRNPK recruitment along a number of immediate early gene loci, including EGR1 and ZFP36, with the highest densities present at the transcription termination sites. Strikingly, hnRNPK knockdown with siRNA resulted in increased pre-RNA levels transcribed downstream of the EGR1 polyadenylation (A) site suggesting altered 3'-end pre-RNA degradation. Further ChIP survey of hnRNPK knockdown uncovered decreased recruitment of the 5'-3' exonuclease XRN2 along EGR1 and downstream of the poly(A) signal without altering RNA polymerase II density at these sites. Immunoprecipitation of hnRNPK and XRN2 from intact and RNase A-treated nuclear extracts followed by shotgun mass spectrometry revealed the presence of hnRNPK and XRN2 in the same complexes along with other spliceosome-related proteins. Our data suggest that hnRNPK may play a role in recruitment of XRN2 to gene loci thus regulating coupling 3'-end pre-mRNA processing to transcription termination.

Acyl-CoA synthetase 1 is induced by Gram-negative bacteria and lipopolysaccharide and is required for phospholipid turnover in stimulated macrophages.

The enzyme acyl-CoA synthetase 1 (ACSL1) is induced by peroxisome proliferator-activated receptor α (PPARα) and PPARγ in insulin target tissues, such as skeletal muscle and adipose tissue, and plays an important role in ß-oxidation in these tissues. In macrophages, however, ACSL1 mediates inflammatory effects without significant effects on ß-oxidation. Thus, the function of ACSL1 varies in different tissues. We therefore investigated the signals and signal transduction pathways resulting in ACSL1 induction in macrophages as well as the consequences of ACSL1 deficiency for phospholipid turnover in LPS-activated macrophages. LPS, Gram-negative bacteria, IFN-γ, and TNFα all induce ACSL1 expression in macrophages, whereas PPAR agonists do not. LPS-induced ACSL1 expression is dependent on Toll-like receptor 4 (TLR4) and its adaptor protein TRIF (Toll-like receptor adaptor molecule 1) but does not require the MyD88 (myeloid differentiation primary response gene 88) arm of TLR4 signaling; nor does it require STAT1 (signal transducer and activator of transcription 1) for maximal induction. Furthermore, ACSL1 deletion attenuates phospholipid turnover in LPS-stimulated macrophages. Thus, the regulation and biological function of ACSL1 in macrophages differ markedly from that in insulin target tissues. These results suggest that ACSL1 may have an important role in the innate immune response. Further, these findings illustrate an interesting paradigm in which the same enzyme, ACSL1, confers distinct biological effects in different cell types, and these disparate functions are paralleled by differences in the pathways that regulate its expression.

Losartan reverses permissive epigenetic changes in renal glomeruli of diabetic db/db mice.

Epigenetic mechanisms such as chromatin histone H3 lysine methylation and acetylation have been implicated in diabetic vascular complications. However, histone modification profiles at pathologic genes associated with diabetic nephropathy in vivo and their regulation by the angiotensin II type 1 receptor (AT1R) are not clear. Here we tested whether treatment of type 2 diabetic db/db mice with the AT1R blocker losartan not only ameliorates diabetic nephropathy, but also reverses epigenetic changes. As expected, the db/db mice had increased blood pressure, mesangial hypertrophy, proteinuria, and glomerular expression of RAGE and PAI-1 vs. control db/+ mice. This was associated with increased RNA polymerase II recruitment and permissive histone marks as well as decreased repressive histone marks at these genes, and altered expression of relevant histone modification enzymes. Increased MCP-1 mRNA levels were not associated with such epigenetic changes, suggesting post-transcriptional regulation. Losartan attenuated key parameters of diabetic nephropathy and gene expression, and reversed some but not all the epigenetic changes in db/db mice. Losartan also attenuated increased H3K9/14Ac at RAGE, PAI-1, and MCP-1 promoters in mesangial cells cultured under diabetic conditions. Our results provide novel information about the chromatin state at key pathologic genes in vivo in diabetic nephropathy mediated in part by AT1R. Thus, combination therapies targeting epigenetic regulators and AT1R could be evaluated for more effective treatment of diabetic nephropathy.

Analysis of gliomas DNA methylation: Assessment of pre-analytical variables.

Precision oncology is driven by molecular biomarkers. For glioblastoma multiforme (GBM), the most common malignant adult primary brain tumor, O6-methylguanine-DNA methyltransferase ( MGMT ) gene DNA promoter methylation is an important prognostic and treatment clinical biomarker. Time consuming pre-analytical steps such as biospecimen storage before fixing, sampling, and processing are major sources of errors and batch effects, that are further confounded by intra-tumor heterogeneity of MGMT promoter methylation. To assess the effect of pre-analytical variables on GBM DNA methylation, tissue storage/sampling (CryoGrid), sample preparation multi-sonicator (PIXUL) and 5-methylcytosine (5mC) DNA immunoprecipitation (Matrix MeDIP-qPCR/seq) platforms were used. MGMT promoter CpG methylation was examined in 173 surgical samples from 90 individuals, 50 of these were used for intra-tumor heterogeneity studies. MGMT promoter methylation levels in paired frozen and formalin fixed paraffin embedded (FFPE) samples were very close, confirming suitability of FFPE for MGMT promoter methylation analysis in clinical settings. Matrix MeDIP-qPCR yielded similar results to methylation specific PCR (MS-PCR). Warm ex-vivo ischemia (37°C up to 4hrs) and 3 cycles of repeated sample thawing and freezing did not alter 5mC levels at MGMT promoter, exon and upstream enhancer regions, demonstrating the resistance of DNA methylation to the most common variations in sample processing conditions that might be encountered in research and clinical settings. 20-30% of specimens exhibited intratumor heterogeneity in the MGMT DNA promoter methylation. Collectively these data demonstrate that variations in sample fixation, ischemia duration and temperature, and DNA methylation assay technique do not have significant impact on assessment of MGMT promoter methylation status. However, intratumor methylation heterogeneity underscores the need for histologic verification and value of multiple biopsies at different GBM geographic tumor sites in assessment of MGMT promoter methylation. Matrix-MeDIP-seq analysis revealed that MGMT promoter methylation status clustered with other differentially methylated genomic loci (e.g. HOXA and lncRNAs), that are likewise resilient to variation in above post-resection pre-analytical conditions. These MGMT -associated global DNA methylation patterns offer new opportunities to validate more granular data-based epigenetic GBM clinical biomarkers where the CryoGrid-PIXUL-Matrix toolbox could prove to be useful.

Regulation of Versican Expression in Macrophages is Mediated by Canonical Type I Interferon Signaling via ISGF3.

Growing evidence supports a role for versican as an important component of the inflammatory response, with both pro- and anti-inflammatory roles depending on the specific context of the system or disease under investigation. Our goal is to understand the regulation of macrophage-derived versican and the role it plays in innate immunity. In previous work, we showed that LPS triggers a signaling cascade involving TLR4, the Trif adaptor, type I interferons, and the type I interferon receptor, leading to increased versican expression by macrophages. In the present study, we used a combination of chromatin immunoprecipitation, siRNA, chemical inhibitors, and mouse model approaches to investigate the regulatory events downstream of the type I interferon receptor to better define the mechanism controlling versican expression. Results indicate that transcriptional regulation by canonical type I interferon signaling via the heterotrimeric transcription factor, ISGF3, controls versican expression in macrophages exposed to LPS. This pathway is not dependent on MAPK signaling, which has been shown to regulate versican expression in other cell types. The stability of versican mRNA may also contribute to prolonged versican expression in macrophages. These findings strongly support a role for macrophage-derived versican as a type I interferon-stimulated gene and further our understanding of versican's role in regulating inflammation.

High incidence and geographic distribution of cleft palate cases in Finland are associated with a regulatory variant in IRF6.

In Finland the frequency of isolated cleft palate (CP) is higher than that of isolated cleft lip with or without cleft palate (CL/P). This trend contrasts to that in other European countries but its genetic underpinnings are unknown. We performed a genome-wide association study for orofacial clefts, which include CL/P and CP, in the Finnish population. We identified rs570516915, a single nucleotide polymorphism that is highly enriched in Finns and Estonians, as being strongly associated with CP ( P = 5.25 × 10 -34 , OR = 8.65, 95% CI 6.11-12.25), but not with CL/P ( P = 7.2 × 10 -5 ), with genome-wide significance. The risk allele frequency of rs570516915 parallels the regional variation of CP prevalence in Finland, and the association was replicated in independent cohorts of CP cases from Finland ( P = 8.82 × 10 -28 ) and Estonia ( P = 1.25 × 10 -5 ). The risk allele of rs570516915 disrupts a conserved binding site for the transcription factor IRF6 within a previously characterized enhancer upstream of the IRF6 gene. Through reporter assay experiments we found that the risk allele of rs570516915 diminishes the enhancer activity. Oral epithelial cells derived from CRISPR-Cas9 edited induced pluripotent stem cells demonstrate that the CP-associated allele of rs570516915 concomitantly decreases the binding of IRF6 and the expression level of IRF6 , suggesting impaired IRF6 autoregulation as a molecular mechanism underlying the risk for CP.

Epigenetic changes in renal genes dysregulated in mouse and rat models of type 1 diabetes.

Epigenetic processes are increasingly being recognized as factors in the pathophysiology of diabetes complications, but few chromatin studies have been done in diabetic nephropathy (DN). We hypothesized that changes in mRNA expression of DN-related genes are associated with epigenetic alterations and aberrant expression of histone-modifying enzymes. RT-PCR and a matrix-chromatin immunoprecipitation platform were used to examine renal mRNA expression, RNA polymerase II (Pol II) recruitment, and epigenetic marks at DN-related genes in the mouse (OVE26) and streptozotocin-induced rat models of type 1 diabetes. Diabetes induced renal expression of Cox2, S100A4/FSP-1, and vimentin genes in both the mouse and the rat models of DN. Mcp-1 and laminin γ1 (Lamc1) expression were increased in diabetic mice but not in rats. Comparison of mRNA and Pol II levels suggested that the diabetes-induced expression of these transcripts is mediated by transcriptional and posttranscriptional processes. Decreases in histone H3 lysine 27 tri-methylation (H3K27m3, silencing mark) and increases in H3 lysine 4 di-methylation (H3K4m2, activating mark) levels were the most consistent epigenetic alterations in the tested genes. In agreement with these results, immunoblot analysis showed increased protein abundance of renal H3K27m2/3 demethylase KDM6A, but no changes in cognate methyltransferase Ezh2 in kidneys of the OVE26 mice compared with controls. In diabetic rats, Ezh2 expression was higher without changes in KDM6A, demonstrating that mechanisms of DN-induced H3K27m3 loss could be species specific. In summary, we show that altered mRNA expression of some DN-related genes is associated with changes in Pol II recruitment and a corresponding decrease in repressive H3K27m3 at the selected loci, and at least in mice with equivalent changes in renal expression of cognate histone-modifying enzymes. This pattern could contribute to diabetes-mediated transitions in chromatin that facilitate transcriptional changes in the diabetic kidney.

MultiomicsTracks96: A high throughput PIXUL-Matrix-based toolbox to profile frozen and FFPE tissues multiomes.

Background: The multiome is an integrated assembly of distinct classes of molecules and molecular properties, or "omes," measured in the same biospecimen. Freezing and formalin-fixed paraffin-embedding (FFPE) are two common ways to store tissues, and these practices have generated vast biospecimen repositories. However, these biospecimens have been underutilized for multi-omic analysis due to the low throughput of current analytical technologies that impede large-scale studies. Methods: Tissue sampling, preparation, and downstream analysis were integrated into a 96-well format multi-omics workflow, MultiomicsTracks96. Frozen mouse organs were sampled using the CryoGrid system, and matched FFPE samples were processed using a microtome. The 96-well format sonicator, PIXUL, was adapted to extract DNA, RNA, chromatin, and protein from tissues. The 96-well format analytical platform, Matrix, was used for chromatin immunoprecipitation (ChIP), methylated DNA immunoprecipitation (MeDIP), methylated RNA immunoprecipitation (MeRIP), and RNA reverse transcription (RT) assays followed by qPCR and sequencing. LC-MS/MS was used for protein analysis. The Segway genome segmentation algorithm was used to identify functional genomic regions, and linear regressors based on the multi-omics data were trained to predict protein expression. Results: MultiomicsTracks96 was used to generate 8-dimensional datasets including RNA-seq measurements of mRNA expression; MeRIP-seq measurements of m6A and m5C; ChIP-seq measurements of H3K27Ac, H3K4m3, and Pol II; MeDIP-seq measurements of 5mC; and LC-MS/MS measurements of proteins. We observed high correlation between data from matched frozen and FFPE organs. The Segway genome segmentation algorithm applied to epigenomic profiles (ChIP-seq: H3K27Ac, H3K4m3, Pol II; MeDIP-seq: 5mC) was able to recapitulate and predict organ-specific super-enhancers in both FFPE and frozen samples. Linear regression analysis showed that proteomic expression profiles can be more accurately predicted by the full suite of multi-omics data, compared to using epigenomic, transcriptomic, or epitranscriptomic measurements individually. Conclusions: The MultiomicsTracks96 workflow is well suited for high dimensional multi-omics studies - for instance, multiorgan animal models of disease, drug toxicities, environmental exposure, and aging as well as large-scale clinical investigations involving the use of biospecimens from existing tissue repositories.

Direct recruitment of ERK cascade components to inducible genes is regulated by heterogeneous nuclear ribonucleoprotein (hnRNP) K.

Components of the ERK cascade are recruited to genes, but it remains unknown how they are regulated at these sites. The RNA-binding protein heterogeneous nuclear ribonucleoprotein (hnRNP) K interacts with kinases and is found along genes including the mitogen-inducible early response gene EGR-1. Here, we used chromatin immunoprecipitations to study co-recruitment of hnRNP K and ERK cascade activity along the EGR-1 gene. These measurements revealed that the spatiotemporal binding patterns of ERK cascade transducers (GRB2, SOS, B-Raf, MEK, and ERK) at the EGR-1 locus resemble both hnRNP K and RNA polymerase II (Pol II). Inhibition of EGR-1 transcription with either serum-responsive factor knockdown or 5,6-dichloro-1-ß-D-ribofuranosylbenzimidazole altered recruitment of all of the above ERK cascade components along this locus that mirrored the changes in Pol II and hnRNP K profiles. siRNA knockdown of hnRNP K decreased the levels of active MEK and ERK at the EGR-1, changes associated with decreased levels of elongating pre-mRNA and less efficient splicing. The hnRNP K dependence and pattern of ERK cascade activation at the c-MYC locus were different from at EGR-1. Ribonucleoprotein immunoprecipitations revealed that hnRNP K was associated with the EGR-1 but not c-MYC mRNAs. These data suggest a model where Pol II transcription-driven recruitment of hnRNP K along the EGR-1 locus compartmentalizes activation of the ERK cascade at these genes, events that regulate synthesis of mature mRNA.