The pig (Sus scrofa) is an important model for evolutionary, comparative, and translational research; however, current functional genomic resources in pigs remain largely limited to one-dimensional genomic annotation and are therefore insufficient for systematically resolving regulatory region-gene relationships, particularly distal ones. Here, we present the Pig Matrix database, a comprehensive 3D regulatory genomics database for pigs, available at https://pigmatrix.kiz.ac.cn/. Built on a standardized experimental framework, Pig Matrix integrates matched multiomics datasets across tissues, developmental stages, and porcine cell lines, including genomic, transcriptomic, epigenomic, and 3D genome information. In total, it contains 16 library types across 7 omics layers and 7,959 processed files from 1,170 libraries. By integrating epigenomic and 3D genome information, Pig Matrix links cis-regulatory elements (CREs) to putative proximal and distal target genes, thereby facilitating interpretation of noncoding variants and genomic signals. This database provides modules for genes, candidate CREs, 3D genome architecture, genome browsing, and single-cell transcriptomics, together with dedicated evolution and comparative resources and user-oriented Genome Annotation and LiftOver tools. A representative use case illustrates how 3D regulatory annotation extends interpretation beyond linear annotation alone, recovering additional candidate genes in domestication-related signals, notably including the classical domestication gene KIT. Pig Matrix also incorporates xenotransplantation-related resources and may support benchmarking of AI models for regulatory genomics. Together, Pig Matrix provides an integrated platform for regulatory interpretation, evolutionary analysis, and comparative genomics in pigs.

Abstract Hepatocellular carcinoma (HCC) is a deadly cancer that emerges from a continuous progression of liver cells from normal to abnormal, often following infections by hepatitis B/C viruses (HBV/HCV), liver fibrosis, and liver cirrhosis (LC), ultimately culminating in cancer. However, there is currently limited systematic molecular analysis of biomarkers at different stages of HCC progression using multi‐omics approaches. We carried out an innovative pipeline by utilizing targeted proteomics and metabolomics to identify potential biomarkers for early detection of HCC in 316 participants, including healthy adults and patients diagnosed with HBV, HCV, LC, and HCC from three independent cohorts. We first established a detailed database of candidate biomarkers for HCC containing 3059 proteins and 103 metabolites, and identified pivotal candidates implicated in the progressive trajectory of liver cancers. Through our developed DeepPRM, scheduled multiple reaction monitoring (MRM)‐targeted approach, and machine learning‐based computational pipeline, we identified an eight‐biomolecular‐based combination with an accuracy rate of 91.43% for early diagnosis of HCC, and a 12‐biomolecular‐based combination with an accuracy rate of 80.00% for detecting changes in HBV–LC progression. These two biomarker combinations significantly improved accuracy compared to traditional tumor biomarkers. Our extensive analysis provides valuable proteomic and metabolomic data resources that will contribute to a deeper understanding of liver disease progression and enhance the identification of potential therapeutic targets.

Non-histone protein acetylation is involved in key cellular processes both in eukaryotes and prokaryotes. Acetylation in bacteria is used to modify proteins involved in metabolism and allow the bacteria to adapt to their environment. TTE (Thermoanaerobacter tengcongensis) is an anaerobic, thermophilic saccharolytic bacterium that grows at extreme temperature range between 50 and 80 ℃. The annotated TTE proteome contains less than 3000 proteins. We analyzed the proteome and acetylome of TTE using 2DLC-MS/MS (2-dimensional liquid chromatography mass spectrum). We evaluated the ability of mass spectrometry technology to cover a relatively small proteome as much as possible. And we also observed wide spread of acetylation in TTE, which changed under different temperatures. A total of 2082 proteins were identified, which accounts for about 82% of the database. A total of 2050 (~ 98%) proteins were quantified in at least one culture condition and 1818 proteins were quantified in all 4 conditions. The result also consisted 3457 acetylation sites corresponding to 827 distinct proteins, which covered 40% of the proteins identified. Bioinformatics analysis reported that proteins related to replication, recombination, repair, and extracellular structure cell wall biogenesis had more than half members acetylated, while energy production, carbohydrate transport, and metabolism related proteins were least acetylated. Our result suggested that acetylation affects the ATP-related energy metabolism and energy-dependent biosynthesis process. Comparing the enzymes related with lysine acetylation and acetyl-CoA (acetyl-coenzyme A) metabolism, we suggested that the acetylation of TTE took a non-enzymatic mechanism and affected by abundance of acetyl-CoA.

Heatmaps are preferred visualization modes for biologists to display high-dimensional information from high-throughput omics data. Many software including website services and R packages are available to generate various types of customized heatmaps. Here, we describe oppHeatmap (omics pilot platform of heatmap), a new tool constructed to render different kinds of heatmaps through MATLAB. oppHeatmap is available for plotting ordinary heatmaps, hierarchical clustering, TreeMaps, microplates graph, sample correlation (full heatmap and upper and lower triangle parts), gene correlation (between columns or tables), and polar heatmaps. oppHeatmap can support the modification of borders, fonts, and colors to customize the final plots. oppHeatmap can not only read data from Microsoft Excel to generate specific heatmaps but also make Excel heatmaps by coloring each cell in Excel. The graphs can be stored in SVG (supported vector graph) format and modified by other SVG recognition software. oppHeatmap is designed by MATLAB AppDesigner with GUI (graphical user interface) operation. The program for oppHeatmap was available at https://github.com/HangZhouSheep/oppHeatmap .are diagramed by.

BACKGROUND:The pig is an economically important livestock species and is a widely applied large animal model in medical research. Enhancers are critical regulatory elements that have fundamental functions in evolution, development and disease. Genome-wide quantification of functional enhancers in the pig is needed.RESULTS:We performed self-transcribing active regulatory region sequencing (STARR-seq) in the porcine kidney epithelial PK15 and testicular ST cell lines, and reliably identified 2576 functional enhancers. Most of these enhancers were located in repetitive sequences and were enriched within silent and lowly expressed genes. Enhancers poorly overlapped with chromatin accessibility regions and were highly enriched in chromatin with the repressive histone modification H3K9me3, which is different from predicted pig enhancers detected using ChIP-seq for H3K27ac or/and H3K4me1 modified histones. This suggests that most pig enhancers identified with STARR-seq are endogenously repressed at the chromatin level and may function during cell type-specific development or at specific developmental stages. Additionally, the PPP3CA gene is associated with the loin muscle area trait and the QKI gene is associated with alkaline phosphatase activity that may be regulated by distal functional enhancers.CONCLUSIONS:In summary, we generated the first functional enhancer map in PK15 and ST cells for the pig genome and highlight its potential roles in pig breeding.

Cohesin regulates sister chromatid cohesion but also contributes to chromosome folding by promoting the formation of chromatin loops, a process mediated by loop extrusion. Although PDS5 regulates cohesin dynamics on chromatin, the exact function of PDS5 in cohesin-mediated chromatin looping remains unclear. Two paralogs of PDS5 exist in vertebrates, PDS5A and PDS5B. Here we show that PDS5A and PDS5B co-localize with RAD21 and CTCF at loop anchors. Rapid PDS5A or PDS5B degradation in liver cancer cells using an inducible degron system reduces chromatin loops and increases loop size. RAD21 enrichment at loop anchors is decreased upon depletion of PDS5A or PDS5B. PDS5B loss also reduces CTCF signals at loop anchors and has a stronger effect on loop enlargement compared with PDS5A. Co-depletion of PDS5A and PDS5B reduces RAD21 levels at loop anchors although the amount of cohesin on chromatin is increased. Our study provides insight into how PDS5 proteins regulate cohesin-mediated chromatin looping.

Aberrant DNA methylation is an epigenetic hallmark of malignant tumours. The DNA methylation level is regulated by not only DNA methyltransferases (DNMTs) but also Ten-Eleven Translocation (TET) family proteins. However, the exact role of TET genes in breast cancer remains controversial. Here, we uncover that the ERα-positive breast cancer patients with high TET2 mRNA expression had better overall survival rates. Consistently, knockout of TET2 promotes the tumorigenesis of ERα-positive MCF7 breast cancer cells. Mechanistically, TET2 loss leads to aberrant DNA methylation (gain of 5mC) at a large proportion of enhancers, accompanied by significant reduction in H3K4me1 and H3K27ac enrichment. By analysing the epigenetically reprogrammed enhancers, we identify oestrogen responsive element (ERE) as one of the enriched motifs of transcriptional factors. Importantly, TET2 loss impairs 17beta-oestradiol (E2)-induced transcription of the epigenetically reprogrammed EREs-associated genes through attenuating the binding of ERα. Taken together, these findings shed light on our understanding of the epigenetic mechanisms underlying the enhancer reprogramming during breast cancer pathogenesis.

Myeloid-derived suppressor cells (MDSCs) are pathologically activated neutrophils and monocytes with potent immunosuppressive activity that regulate immune responses in the tumor microenvironment. We identified a novel long noncoding RNA (lncRNA), named as lnc57Rik, in the MDSCs that controls their immunosuppressive functions. Lnc57Rik was induced in in vitro and in vivo inflammatory settings and upregulated the genes related to MDSC-mediated immunosuppression, including Arg-1, NOS2, NOX2, and COX2. Furthermore, Lnc57Rik can not only bind with the C/EBPβ isoform liver-enriched activator protein to activate C/EBPβ but also with the methyltransferase WD repeat-containing protein 5 that enables the enrichment of histone H3 trimethylated lysine 4 marks on the promoter regions of Arg-1, NOS2, NOX2, and COX2, eventually resulting in their transcriptional activation. Furthermore, the conserved human lnc57Rik has a similar function as murine lnc57Rik. Taken together, upregulation of lnc57Rik in the tumor microenvironment promotes the immunosuppressive function of MDSCs.

Chromatin dynamics as well as genetic evolution underlies the adaptability of tumour cells to environmental cues. Three subtypes of tumour cells have been identified in glioblastoma, one of the commonest malignant brain tumours in adults. During tumour progression or under therapeutic pressure, the non-mesenchymal subtypes may progress to the mesenchymal subtype, leading to unfavourable prognosis. However, the molecular mechanisms for this transition remain poorly understood. Here taking a TNFα-induced cellular model, we profile the chromatin accessibility dynamics during mesenchymal transition. Moreover, we identify the JUN family as one of the key driving transcription factors for the gained chromatin accessibility. Accordingly, inhibition of JUN phosphorylation and therefore its transcription activity successfully impedes TNFα-induced chromatin remodelling and mesenchymal transition. In line with these findings based on experimental models, JUN activity is positively correlated with mesenchymal features in clinical glioblastoma specimens. Together, this study unveils a deregulated transcription regulatory network in glioblastoma progression and hopefully provides a rationale for anti-glioblastoma therapy.

Generation of functional transcripts requires transcriptional initiation at regular start sites, avoiding production of aberrant and potentially hazardous aberrant RNAs. The mechanisms maintaining transcriptional fidelity and the impact of spurious transcripts on cellular physiology and organ function have not been fully elucidated. Here we show that TET3, which successively oxidizes 5-methylcytosine to 5-hydroxymethylcytosine (5hmC) and other derivatives, prevents aberrant intragenic entry of RNA polymerase II pSer5 into highly expressed genes of airway smooth muscle cells, assuring faithful transcriptional initiation at canonical start sites. Loss of TET3-dependent 5hmC production in SMCs results in accumulation of spurious transcripts, which stimulate the endosomal nucleic-acid-sensing TLR7/8 signaling pathway, thereby provoking massive inflammation and airway remodeling resembling human bronchial asthma. Furthermore, we found that 5hmC levels are substantially lower in human asthma airways compared with control samples. Suppression of spurious transcription might be important to prevent chronic inflammation in asthma.

Abstract The dynamic changes of RNA N6-methyl-adenosine (m6A) during cancer progression contribute to quick adaption to microenvironmental changes. Here, we profiled the cancer cell m6A dynamics in the hypoxic tumor niche and its pathological consequences in glioblastoma multiforme (GBM). The m6A demethylase ALKBH5 was induced in GBM models under hypoxic conditions and was associated with a hypoxic gene signature in GBM patient samples. Depletion or inactivation of ALKBH5 in GBM cells significantly suppressed hypoxia-induced tumor-associated macrophage (TAM) recruitment and immunosuppression in allograft tumors. Expression and secretion of CXCL8/IL8 were significantly suppressed in ALKBH5-deficient tumors. However, ALKBH5 did not regulate CXCL8 m6A directly. Instead, hypoxia-induced ALKBH5 erased m6A deposition from the lncRNA NEAT1, stabilizing the transcript and facilitating NEAT1-mediated paraspeckle assembly, which led to relocation of the transcriptional repressor SFPQ from the CXCL8 promoter to paraspeckles and, ultimately, upregulation of CXCL8/IL8 expression. Accordingly, ectopic expression of CXCL8 in ALKBH5-deficient GBM cells partially restored TAM recruitment and tumor progression. Together, this study links hypoxia-induced epitranscriptomic changes to the emergence of an immunosuppressive microenvironment facilitating tumor evasion. Significance: Hypoxia induces tumor immune microenvironment remodeling through an ALKBH5-mediated epigenetic and epitranscriptomic mechanism, providing potential immunotherapeutic strategies for treating glioblastoma.

BACKGROUND:The pathogenesis of life-threatening cardiopulmonary diseases such as pulmonary hypertension (PH) and chronic obstructive pulmonary disease (COPD) originates from a complex interplay of environmental factors and genetic predispositions that is not fully understood. Likewise, little is known about developmental abnormalities or epigenetic dysregulations that might predispose for PH or COPD in adult individuals.METHODS:To identify pathology-associated epigenetic alteration in diseased lung tissues, we screened a cohort of human patients with PH and COPD for changes of histone modifications by immunofluorescence staining. To analyze the function of H4K20me2/3 in lung pathogenesis, we developed a series of Suv4-20h1 knockout mouse lines targeting cardiopulmonary progenitor cells and different heart and lung cell types, followed by hemodynamic studies and morphometric assessment of tissue samples. Molecular, cellular, and biochemical techniques were applied to analyze the function of Suv4-20h1-dependent epigenetic processes in cardiopulmonary progenitor cells and their derivatives.RESULTS:We discovered a strong reduction of the histone modifications of H4K20me2/3 in human patients with COPD but not patients with PH that depend on the activity of the H4K20 di-methyltransferase SUV4-20H1. Loss of Suv4-20h1 in cardiopulmonary progenitor cells caused a COPD-like/PH phenotype in mice including the formation of perivascular tertiary lymphoid tissue and goblet cell hyperplasia, hyperproliferation of smooth muscle cells/myofibroblasts, impaired alveolarization and maturation defects of the microvasculature leading to massive right ventricular dilatation and premature death. Mechanistically, SUV4-20H1 binds directly to the 5'-upstream regulatory element of the superoxide dismutase 3 (Sod3) gene to repress its expression. Increased levels of the extracellular SOD3 enzyme in Suv4-20h1 mutants increases hydrogen peroxide concentrations, causing vascular defects and impairing alveolarization.CONCLUSIONS:Our findings reveal a pivotal role of the histone modifier SUV4-20H1 in cardiopulmonary codevelopment and uncover the developmental origins of cardiopulmonary diseases. We assume that the study will facilitate the understanding of pathogenic events causing PH and COPD and aid the development of epigenetic drugs for the treatment of cardiopulmonary diseases.

Lysine acetylation (Kac) as an important posttranslational modification of histones is essential for the regulation of gene expression in hepatocellular carcinoma (HCC). However, the atlas of whole acetylated proteins in HCC tissues and the difference in protein acetylation between normal human tissues and HCC tissues are unknown. In this report, we characterized the proteome and acetyl proteome (acetylome) profile of normal, paracancerous, and HCC liver tissues in human clinical samples by quantitative proteomics techniques. We identified 6781 acetylation sites of 2582 proteins and quantified 2492 acetylation sites of 1190 proteins in normal, paracancerous, and HCC liver tissues. Among them, 15 proteins were multiacetylated with more than 10 lysine residues. The histone acetyltransferases p300 and CBP were found to be hyperacetylated in hepatitis B virus pathway. Moreover, we found that 250 Kac sites of 214 proteins were upregulated and 662 Kac sites of 451 proteins were downregulated in HCC compared with normal liver tissues. Additionally, the acetylation levels of lysine 120 in histone H2B (H2BK120ac), lysine 18 in histone H3.3 (H3.3K18ac), and lysine 77 in histone H4 (H4K77ac) were increased in HCC. Interestingly, the higher levels of H2BK120ac, H3.3K18ac, and H4K77ac were significantly associated with worse prognosis, such as poorer survival and higher recurrence in an independent clinical cohort of HCC patients. Overall, this study lays a foundation for understanding the functions of acetylation in HCC and provides potential prognostic factors for the diagnosis and therapy of HCC.

The host immune response is a fundamental mechanism for attenuating cancer progression. Here we report a role for the DNA demethylase and tumor suppressor TET2 in host anti-tumor immunity. Deletion of Tet2 in mice elevates IL-6 levels upon tumor challenge. Elevated IL-6 stimulates immunosuppressive granulocytic myeloid-derived suppressor cells (G-MDSCs), which in turn reduce CD8+ T cells upon tumor challenge. Consequently, systematic knockout of Tet2 in mice leads to accelerated syngeneic tumor growth, which is constrained by anti-PD-1 blockade. Removal of G-MDSCs by the anti-mouse Ly6g antibodies restores CD8+ T-cell numbers in Tet2-/- mice and reboots their anti-tumor activity. Importantly, anti-IL-6 antibody treatment blocks the expansion of G-MDSCs and inhibits syngeneic tumor growth. Collectively, these findings reveal a TET2-mediated IL-6/G-MDSCs/CD8+ T-cell immune response cascade that safeguards host adaptive anti-tumor immunity, offering a cell non-autonomous mechanism of TET2 for tumor suppression.

The author would like to add the below information in this correction. A similar study from Chao Lu group was published online on 5 September 2019 in Nature, entitled "The histone mark H3K36me2 recruits DNMT3A and shapes the intergenic DNA methylation landscape" (Weinberg et al., 2019). Although both the studies reported the preferential recognition of H3K36me2 by DNMT3A PWWP, ours in addition uncovered a stimulation function by such interaction on the activity of DNMT3A. On the disease connections, we used a NSD2 gain-of-function model which led to the discovery of potential therapeutic implication of DNA inhibitors in the related cancers, while the other study only used NSD1 and DNMT3A loss-of-function models.

Alzheimer's disease (AD), a progressive neurodegenerative disorder, is the most common untreatable form of dementia. Identifying molecular biomarkers that allow early detection remains a key challenge in the diagnosis, treatment, and prognostic evaluation of the disease. Here, we report a novel experimental and analytical model characterizing epigenetic alterations during AD onset and progression. We generated the first integrated base-resolution genome-wide maps of the distribution of 5-methyl-cytosine (5mC), 5-hydroxymethyl-cytosine (5hmC), and 5-formyl/carboxy-cytosine (5fC/caC) in normal and AD neurons. We identified 27 AD region-specific and 39 CpG site-specific epigenetic signatures that were independently validated across our familial and sporadic AD models, and in an independent clinical cohort. Thus, our work establishes a new model and strategy to study the epigenetic alterations underlying AD onset and progression and provides a set of highly reliable AD-specific epigenetic signatures that may have early diagnostic and prognostic implications.

Diabetes is a complex metabolic syndrome that is characterized by prolonged high blood glucose levels and frequently associated with life-threatening complications1,2. Epidemiological studies have suggested that diabetes is also linked to an increased risk of cancer3–5. High glucose levels may be a prevailing factor that contributes to the link between diabetes and cancer, but little is known about the molecular basis of this link and how the high glucose state may drive genetic and/or epigenetic alterations that result in a cancer phenotype. Here we show that hyperglycaemic conditions have an adverse effect on the DNA 5-hydroxymethylome. We identify the tumour suppressor TET2 as a substrate of the AMP-activated kinase (AMPK), which phosphorylates TET2 at serine 99, thereby stabilizing the tumour suppressor. Increased glucose levels impede AMPK-mediated phosphorylation at serine 99, which results in the destabilization of TET2 followed by dysregulation of both 5-hydroxymethylcytosine (5hmC) and the tumour suppressive function of TET2 in vitro and in vivo. Treatment with the anti-diabetic drug metformin protects AMPK-mediated phosphorylation of serine 99, thereby increasing TET2 stability and 5hmC levels. These findings define a novel ‘phospho-switch’ that regulates TET2 stability and a regulatory pathway that links glucose and AMPK to TET2 and 5hmC, which connects diabetes to cancer. Our data also unravel an epigenetic pathway by which metformin mediates tumour suppression. Thus, this study presents a new model for how a pernicious environment can directly reprogram the epigenome towards an oncogenic state, offering a potential strategy for cancer prevention and treatment. Modulation of DNA 5-hydroxymethylcytosine by glucose reveals an AMPK–TET2–5hmC axis that links diabetes to cancer.

Immune response by T cells is essential for a healthy body against cancer, infection, and pathophysiological alteration. The activation and expansion of T cells can be inhibited by dasatinib, a tyrosine inhibitor, thus improving the outcome of diseases, such as autoimmune disease, graft-versus-host disease, and transplant rejection. The underlying mechanism of inhibition by dasatinib is elusive. Here, we designed and synthesized a CRISPR/Cas9 screening library that includes 6,149 genes. Using the library, we performed dasatinib CRISPR/cas9 screening in Jurkat cell, a T lymphocyte cell. We firstly identified survival essential genes for Jurkat cells. Comparing with other CRISPR/Cas9 screenings, we obtained Jurkat cell specific essential genes. By comparing dasatinib treatment to control, we identified a set of dasatinib targets, which includes known targets: CSK, LCK, ZAP70, and previously unknown targets: ZFP36L2, LRPPRC, CFLAR, PD-1, CD45 et al. Visualizing these target genes on T cell receptor signaling pathway, we found several genes could be inhibited by dasatinib. Furthermore, we introduced a framework, 9-square, to classify genes and found a group of genes that are associated with dasatinib resistance, possibly linking the side effects of dasatinib. These data reveal a set of dasatinib targets and demonstrate the molecular potential functions of dasatinib. Identification of dasatinib targets will broaden our understanding to its molecular mechanism, and thus benefits to clinical outcome.

N6-methyladenosine (m6A) is an abundant modification in eukaryotic mRNA, regulating mRNA dynamics by influencing mRNA stability, splicing, export, and translation. However, the precise m6A regulating machinery still remains incompletely understood. Here we demonstrate that ZC3H13, a zinc-finger protein, plays an important role in modulating RNA m6A methylation in the nucleus. We show that knockdown of Zc3h13 in mouse embryonic stem cell significantly decreases global m6A level on mRNA. Upon Zc3h13 knockdown, a great majority of WTAP, Virilizer, and Hakai translocate to the cytoplasm, suggesting that Zc3h13 is required for nuclear localization of the Zc3h13-WTAP-Virilizer-Hakai complex, which is important for RNA m6A methylation. Finally, Zc3h13 depletion, as does WTAP, Virilizer, or Hakai, impairs self-renewal and triggers mESC differentiation. Taken together, our findings demonstrate that Zc3h13 plays a critical role in anchoring WTAP, Virilizer, and Hakai in the nucleus to facilitate m6A methylation and to regulate mESC self-renewal.