Temporal mapping of CEBPA and CEBPB binding during liver regeneration reveals dynamic occupancy and specific regulatory codes for homeostatic and cell cycle gene batteries

CCAAT/Enhancer-Binding Proteins α and β Regulate Ovulation and Gene Expression via Dose- and Stage-Dependent Mechanisms

The preovulatory luteinizing hormone (LH) surge orchestrates complex cellular and molecular events leading to ovulation. CCAAT/enhancer-binding proteins α and β (C/EBPα/β) are transcription factors acutely induced by the LH surge and crucial for ovulation and granulosa cell luteinization. However, biological processes (BPs) and their regulatory mechanisms downstream of C/EBPα/β in the preovulatory ovary are not completely understood. To address this knowledge gap, we generated Cebpa/bfl/fl;Pgr-Cre mutants and compared them with Cebpa/bfl/fl;Cyp19a1-Cre mutant female mice: Cebpa/bfl/fl;Cyp19a1-Cre mutants have undetectable levels of C/EBPα/β throughout the preovulatory stages and do not ovulate, aligning with previous reports; and Cebpa/bfl/fl;Pgr-Cre mutants present gradual depletion of C/EBPα/β during the late preovulatory stage and a reduced ovulation rate. Comparison of these two models indicates that sustained expression of C/EBPα/β throughout the preovulatory stages is necessary for successful ovulation and provides a unique opportunity to interrogate gene regulatory mechanisms by C/EBPα/β during different preovulatory time windows and the effect of dysregulating C/EBPα/β on ovulation-associated BPs. Our study revealed that C/EBPα/β regulate gene expression and distinct biological functions such as vascular remodeling via dose- and preovulatory stage-dependent mechanisms. These findings shed new light on the intricate mechanisms of gene regulation by C/EBPα/β downstream of the LH surge.

Continuous Activation of C/EBPβ Transcription Factor Prevents Fibrosis Resolution After Alcohol Cessation

BACKGROUND & AIMS

Abstinence is an important therapeutic intervention for patients with alcohol-associated liver disease (ALD). However, fibrosis improvement after cessation is not uniform and some patients do not improve.

METHODS

Mice were fed high-fat diet with 20% alcohol in the drinking water for 20 weeks (ALD) followed by 4 weeks of chow diet with plain water (resolution). scATAC-seq dataset was analyzed using Signac R package. Cebpb floxed mice received AAV-TBG-Cre or AAV-control at the time of alcohol cessation. Hepatocyte-macrophage and endothelial cell-hepatocytes crosstalk was investigated using a Transwell coculture system. To test the role of angiopoietin mice were treated with recombinant angiopoietin-1, 1 week after alcohol cessation.

RESULTS

We analyzed differentially accessible regions in hepatocytes from control, ALD, or 4 weeks post alcohol cessation mice and identified transcription factors activated in ALD that remained activated after alcohol withdrawal. The top hit was CCAAT enhancer binding protein beta (C/EBPβ). We found that hepatocyte-specific Cebpb knockout at the time of alcohol cessation promoted fibrosis resolution. The resolution was mediated by altered hepatocyte-macrophage crosstalk. C/EBPβ suppressed the expression of CYP3A family of enzymes in hepatocytes and downstream macrophage collagen degradation ability. Cebpb knockout in hepatocytes promoted a proresolving phenotype in liver macrophages. We further identified upstream events leading to persistent C/EBPβ activation. C/EBPβ was induced by alcohol-mediated endothelial changes during ALD development and resolution. Restoring endothelial cell function with angiopoietin-1 supplementation reduced C/EBPβ and promoted fibrosis resolution.

CONCLUSIONS

Taken together, alcohol-induced C/EBPβ activation is a key driver of poor disease resolution in ALD and a promising target for patients who fail to recover after alcohol abstinence.

A single-nucleus and spatial transcriptomic atlas of the COVID-19 liver reveals topological, functional, and regenerative organ disruption in patients

BACKGROUND

The molecular underpinnings of organ dysfunction in severe COVID-19 and its potential long-term sequelae are under intense investigation. To shed light on these in the context of liver function, we perform single-nucleus RNA-seq and spatial transcriptomic profiling of livers from 17 COVID-19 decedents.

RESULTS

We identify hepatocytes positive for SARS-CoV-2 RNA with an expression phenotype resembling infected lung epithelial cells, and a central role in a pro-fibrotic TGFβ signaling cell-cell communications network. Integrated analysis and comparisons with healthy controls reveal extensive changes in the cellular composition and expression states in COVID-19 liver, providing the underpinning of hepatocellular injury, ductular reaction, pathologic vascular expansion, and fibrogenesis characteristic of COVID-19 cholangiopathy. We also observe Kupffer cell proliferation and erythrocyte progenitors for the first time in a human liver single-cell atlas. Despite the absence of a clinical acute liver injury phenotype, endothelial cell composition is dramatically impacted in COVID-19, concomitantly with extensive alterations and profibrogenic activation of reactive cholangiocytes and mesenchymal cells.

CONCLUSIONS

Our atlas provides novel insights into liver physiology and pathology in COVID-19 and forms a foundational resource for its investigation and understanding.

C/EBPβ transcription factor promotes alcohol-induced liver fibrosis in males via HDL remodeling

BACKGROUND

Alcohol-associated liver disease (ALD) is the main cause of alcohol-associated mortality. However, the mechanism of ALD development is poorly understood. Epigenetic changes are thought to play an important role in ALD. We aimed to define the epigenetic changes induced by alcohol and predict drivers of these changes.

METHODS

Mice were fed high-fat diet with or without 20% of alcohol in the drinking water for 20 weeks (WDA model). scATAC-seq data set was analyzed using Signac R package. To test the role of C/EBPβ, Cebpb-floxed mice were treated with AAV8-TBG-Cre or AAV8-control.

RESULTS

We analyzed differentially accessible regions in livers from control and alcohol-fed mice and found that activity of C/EBPβ transcription factor was associated with alcohol-induced epigenetic changes in hepatocytes. C/EBPβ protein levels were significantly upregulated in multiple models of ALD and human ALD samples. Using hepatocyte-specific Cebpb knockout mice we found that Cebpb loss protected male mice from alcohol-induced fibrosis development. We found no protection in female mice, suggesting that this mechanism is specific to male ALD. In vitro studies suggested that the protective effect of Cebpb loss was mediated by altered hepatocyte-macrophage cross talk. Cebpb knockout in hepatocytes reduced a profibrotic and promoted a pro-resolving phenotype in macrophages, thus modulating ALD development. We further identified the mediators of the cross talk. Cebpb knockout altered the expression of several HDL protein components, increasing APOA1 and apolipoprotein M and reducing apolipoprotein E and SAA levels in male mice. HDL secreted by Cebpb knockout hepatocytes was sufficient to confer anti-inflammatory and antifibrotic changes to macrophages.

CONCLUSIONS

Taken together, alcohol-induced C/EBPβ activation is a key driver of ALD fibrosis in males via C/EBPβ-dependent HDL remodeling.

Transcriptional cascades during fasting amplify gluconeogenesis and instigate a secondary wave of ketogenic gene transcription

BACKGROUND AND AIMS

During fasting, bodily homeostasis is maintained due to hepatic production of glucose (gluconeogenesis) and ketone bodies (ketogenesis). The main hormones governing hepatic fuel production are glucagon and glucocorticoids that initiate transcriptional programs aimed at supporting gluconeogenesis and ketogenesis.

METHODS

Using primary mouse hepatocytes as an ex vivo model, we employed transcriptomic analysis (RNA-seq), genome-wide profiling of enhancer dynamics (ChIP-seq), perturbation experiments (inhibitors, shRNA), hepatic glucose production measurements and computational analyses.

RESULTS

We found that in addition to the known metabolic genes transcriptionally induced by glucagon and glucocorticoids, these hormones induce a set of genes encoding transcription factors (TFs) thereby initiating transcriptional cascades. Upon activation by glucocorticoids, the glucocorticoid receptor (GR) induced the genes encoding two TFs: CCAAT/enhancer-binding protein beta (C/EBPβ) and peroxisome proliferator-activated receptor alpha (PPARα). We found that the GR-C/EBPβ cascade mainly serves as a secondary amplifier of primary hormone-induced gene programs. C/EBPβ augmented gluconeogenic gene expression and hepatic glucose production. Conversely, the GR-PPARα cascade initiated a secondary transcriptional wave of genes supporting ketogenesis. The cascade led to synergistic induction of ketogenic genes which is dependent on protein synthesis. Genome-wide analysis of enhancer dynamics revealed numerous enhancers activated by the GR-PPARα cascade. These enhancers were proximal to ketogenic genes, enriched for the PPARα response element and showed increased PPARα binding.

CONCLUSION

This study reveals abundant transcriptional cascades occurring during fasting. These cascades serve two separated purposes: the amplification of the gluconeogenic transcriptional program and the induction of a gene program aimed at enhancing ketogenesis.

A spatiotemporal atlas of mouse liver homeostasis and regeneration

The mechanism by which mammalian liver cell responses are coordinated during tissue homeostasis and perturbation is poorly understood, representing a major obstacle in our understanding of many diseases. This knowledge gap is caused by the difficulty involved with studying multiple cell types in different states and locations, particularly when these are transient. We have combined Stereo-seq (spatiotemporal enhanced resolution omics-sequencing) with single-cell transcriptomic profiling of 473,290 cells to generate a high-definition spatiotemporal atlas of mouse liver homeostasis and regeneration at the whole-lobe scale. Our integrative study dissects in detail the molecular gradients controlling liver cell function, systematically defining how gene networks are dynamically modulated through intercellular communication to promote regeneration. Among other important regulators, we identified the transcriptional cofactor TBL1XR1 as a rheostat linking inflammation to Wnt/β-catenin signaling for facilitating hepatocyte proliferation. Our data and analytical pipelines lay the foundation for future high-definition tissue-scale atlases of organ physiology and malfunction.

Germline cis variant determines epigenetic regulation of the anti-cancer drug metabolism gene dihydropyrimidine dehydrogenase (DPYD)

Enhancers are critical for regulating tissue-specific gene expression, and genetic variants within enhancer regions have been suggested to contribute to various cancer-related processes, including therapeutic resistance. However, the precise mechanisms remain elusive. Using a well-defined drug-gene pair, we identified an enhancer region for dihydropyrimidine dehydrogenase (DPD, DPYD gene) expression that is relevant to the metabolism of the anti-cancer drug 5-fluorouracil (5-FU). Using reporter systems, CRISPR genome-edited cell models, and human liver specimens, we demonstrated in vitro and vivo that genotype status for the common germline variant (rs4294451; 27% global minor allele frequency) located within this novel enhancer controls DPYD transcription and alters resistance to 5-FU. The variant genotype increases recruitment of the transcription factor CEBPB to the enhancer and alters the level of direct interactions between the enhancer and DPYD promoter. Our data provide insight into the regulatory mechanisms controlling sensitivity and resistance to 5-FU.

Germline cis variant determines epigenetic regulation of the anti-cancer drug metabolism gene dihydropyrimidine dehydrogenase (DPYD)

Enhancers are critical for regulating tissue-specific gene expression, and genetic variants within enhancer regions have been suggested to contribute to various cancer-related processes, including therapeutic resistance. However, the precise mechanisms remain elusive. Using a well-defined drug-gene pair, we identified an enhancer region for dihydropyrimidine dehydrogenase (DPD, DPYD gene) expression that is relevant to the metabolism of the anti-cancer drug 5-fluorouracil (5-FU). Using reporter systems, CRISPR genome edited cell models, and human liver specimens, we demonstrated in vitro and vivo that genotype status for the common germline variant (rs4294451; 27% global minor allele frequency) located within this novel enhancer controls DPYD transcription and alters resistance to 5-FU. The variant genotype increases recruitment of the transcription factor CEBPB to the enhancer and alters the level of direct interactions between the enhancer and DPYD promoter. Our data provide insight into the regulatory mechanisms controlling sensitivity and resistance to 5-FU.

CEBPA Overexpression Enhances β-Cell Proliferation and Survival

A commonality between type 1 and type 2 diabetes is the decline in functional β-cell mass. The transcription factor Nkx6.1 regulates β-cell development and is integral for proper β-cell function. We have previously demonstrated that Nkx6.1 depends on c-Fos mediated upregulation and the nuclear hormone receptors Nr4a1 and Nr4a3 to increase β-cell insulin secretion, survival, and replication. Here, we demonstrate that Nkx6.1 overexpression results in upregulation of the bZip transcription factor CEBPA and that CEBPA expression is independent of c-Fos regulation. In turn, CEBPA overexpression is sufficient to enhance INS-1 832/13 β-cell and primary rat islet proliferation. CEBPA overexpression also increases the survival of β-cells treated with thapsigargin. We demonstrate that increased survival in response to ER stress corresponds with changes in expression of various genes involved in the unfolded protein response, including decreased Ire1a expression. These data show that CEBPA is sufficient to enhance functional β-cell mass by increasing β-cell proliferation and modulating the unfolded protein response.

Peak Scores Significantly Depend on the Relationships between Contextual Signals in ChIP-Seq Peaks

Chromatin immunoprecipitation followed by massively parallel DNA sequencing (ChIP-seq) is a central genome-wide method for in vivo analyses of DNA-protein interactions in various cellular conditions. Numerous studies have demonstrated the complex contextual organization of ChIP-seq peak sequences and the presence of binding sites for transcription factors in them. We assessed the dependence of the ChIP-seq peak score on the presence of different contextual signals in the peak sequences by analyzing these sequences from several ChIP-seq experiments using our fully enumerative GPU-based de novo motif discovery method, Argo_CUDA. Analysis revealed sets of significant IUPAC motifs corresponding to the binding sites of the target and partner transcription factors. For these ChIP-seq experiments, multiple regression models were constructed, demonstrating a significant dependence of the peak scores on the presence in the peak sequences of not only highly significant target motifs but also less significant motifs corresponding to the binding sites of the partner transcription factors. A significant correlation was shown between the presence of the target motifs FOXA2 and the partner motifs HNF4G, which found experimental confirmation in the scientific literature, demonstrating the important contribution of the partner transcription factors to the binding of the target transcription factor to DNA and, consequently, their important contribution to the peak score.

Sepsis-induced endothelial dysfunction drives acute-on-chronic liver failure through Angiopoietin-2-HGF-C/EBPβ pathway

BACKGROUND AND AIMS

Acute-on-chronic liver failure (ACLF) is an acute liver and multisystem failure in patients with previously stable cirrhosis. A common cause of ACLF is sepsis secondary to bacterial infection. Sepsis-associated ACLF involves a loss of differentiated liver function in the absence of direct liver injury, and its mechanism is unknown. We aimed to study the mechanism of sepsis-associated ACLF using a novel mouse model.

APPROACH AND RESULTS

Sepsis-associated ACLF was induced by cecal ligation and puncture procedure (CLP) in mice treated with thioacetamide (TAA). The combination of TAA and CLP resulted in a significant decrease in liver synthetic function and high mortality. These changes were associated with reduced metabolic gene expression and increased CCAAT enhancer binding protein beta (C/EBPβ) transcriptional activity. We found that C/EBPβ binding to its target gene promoters was increased. In humans, C/EBPβ chromatin binding was similarly increased in the ACLF group compared with control cirrhosis. Hepatocyte-specific Cebpb knockout mice had reduced mortality and increased gene expression of hepatocyte differentiation markers in TAA/CLP mice, suggesting that C/EBPβ promotes liver failure in these mice. C/EBPβ activation was associated with endothelial dysfunction, characterized by reduced Angiopoietin-1/Angiopoietin-2 ratio and increased endothelial production of HGF. Angiopoietin-1 supplementation or Hgf knockdown reduced hepatocyte C/EBPβ accumulation, restored liver function, and reduced mortality, suggesting that endothelial dysfunction induced by sepsis drives ACLF through HGF-C/EBPβ pathway.

CONCLUSIONS

The transcription factor C/EBPβ is activated in both mouse and human ACLF and is a potential therapeutic target to prevent liver failure in patients with sepsis and cirrhosis.

Dynamic de novo heterochromatin assembly and disassembly at replication forks ensures fork stability

Chromatin is dynamically reorganized when DNA replication forks are challenged. However, the process of epigenetic reorganization and its implication for fork stability is poorly understood. Here we discover a checkpoint-regulated cascade of chromatin signalling that activates the histone methyltransferase EHMT2/G9a to catalyse heterochromatin assembly at stressed replication forks. Using biochemical and single molecule chromatin fibre approaches, we show that G9a together with SUV39h1 induces chromatin compaction by accumulating the repressive modifications, H3K9me1/me2/me3, in the vicinity of stressed replication forks. This closed conformation is also favoured by the G9a-dependent exclusion of the H3K9-demethylase JMJD1A/KDM3A, which facilitates heterochromatin disassembly upon fork restart. Untimely heterochromatin disassembly from stressed forks by KDM3A enables PRIMPOL access, triggering single-stranded DNA gap formation and sensitizing cells towards chemotherapeutic drugs. These findings may help in explaining chemotherapy resistance and poor prognosis observed in patients with cancer displaying elevated levels of G9a/H3K9me3.

The rhythmic coupling of Egr-1 and Cidea regulates age-related metabolic dysfunction in the liver of male mice

The liver lipid metabolism of older individuals canbecome impaired and the circadian rhythm of genes involved in lipid metabolism is also disturbed. Although the link between metabolism and circadian rhythms is already recognized, how these processes are decoupled in liver during aging is still largely unknown. Here, we show that the circadian rhythm for the transcription factor Egr-1 expression is shifted forward with age in male mice. Egr-1 deletion accelerates liver age-related metabolic dysfunction, which associates with increased triglyceride accumulation, disruption of the opposite rhythmic coupling of Egr-1 and Cidea (Cell Death Inducing DFFA Like Effector A) at the transcriptional level and large lipid droplet formation. Importantly, adjustment of the central clock with light via a 4-hour forward shift in 6-month-old mice, leads to recovery the rhythm shift of Egr-1 during aging and largely ameliorated liver metabolic dysfunction. All our collected data suggest that liver Egr-1 might integrate the central and peripheral rhythms and regulate metabolic homeostasis in the liver.

C/EBPβ Coupled with E2F2 Promoted the Proliferation of hESC-Derived Hepatocytes through Direct Binding to the Promoter Regions of Cell-Cycle-Related Genes

Human embryonic stem cells (hESCs) hold the potential to solve the problem of the shortage of functional hepatocytes in clinical applications and drug development. However, a large number of usable hepatocytes derived from hESCs cannot be effectively obtained due to the limited proliferation capacity. In this study, we found that enhancement of liver transcription factor C/EBPβ during hepatic differentiation could not only significantly promote the expression of hepatic genes, such as albumin, alpha fetoprotein, and alpha-1 antitrypsin, but also dramatically reinforce proliferation-related phenotypes, including increasing the expression of proliferative genes, such as CDC25C, CDC45L, and PCNA, and the activation of cell cycle and DNA replication pathways. In addition, the analysis of CUT&Tag sequencing further revealed that C/EBPβ is directly bound to the promoter region of proliferating genes to promote cell proliferation; this interaction between C/EBPβ and DNA sequences of the promoters was verified by luciferase assay. On the contrary, the knockdown of C/EBPβ could significantly inhibit the expression of the aforementioned proliferative genes. RNA transcriptome analysis and GSEA enrichment indicated that the E2F family was enriched, and the expression of E2F2 was changed with the overexpression or knockdown of C/EBPβ. Moreover, the results of CUT&Tag sequencing showed that C/EBPβ also directly bound the promoter of E2F2, regulating E2F2 expression. Interestingly, Co-IP analysis exhibited a direct binding between C/EBPβ and E2F2 proteins, and this interaction between these two proteins was also verified in the LO2 cell line, a hepatic progenitor cell line. Thus, our results demonstrated that C/EBPβ first initiated E2F2 expression and then coupled with E2F2 to regulate the expression of proliferative genes in hepatocytes during the differentiation of hESCs. Therefore, our findings open a new avenue to provide an in vitro efficient approach to generate proliferative hepatocytes to potentially meet the demands for use in cell-based therapeutics as well as for pharmaceutical and toxicological studies.

Combined effect of polystyrene microplastics and bisphenol A on the human embryonic stem cells-derived liver organoids: The hepatotoxicity and lipid accumulation

Human are exposed to microplastics (MP) via inhalation or ingestion daily and inevitably. The liver is an important target organ of MP. Bisphenol A (BPA) is one of commonly used plasticizers. It is added in plastics, but also generally detected in the biological samples of human beings. However, the combined toxic effect of MP and BPA on human liver is unclear. In this study, a novel 3D in vitro model, the liver organoid (LO) derived from human-pluripotent stem cells, has been utilized to explore the 1 μm polystyrene (PS)-induced hepatotoxicity with BPA individually and jointly. Conclusively, all the changes in the cytotoxicity, cellular and molecular makers regarding the energy supplement, hepatic injury, oxidative stress, inflammatory response, disruption in the lipid accumulation, as well as epigenetics regulation induced by BPA or PS on the LOs individually were comparable to previous study. The BPA levels in the culture medium were declined by the added PS. The combined adverse effect of PS and BPA on the LOs was identified to be synergistic upon deteriorated hepatotoxicity and interfered the gene panels related to multiple processes of lipid metabolism, together with the proteins of HNF4A, CD36, ACC1, CPT1A, CYP2E1, ERα and ERβ. Specifically, PS didn't change the ERα or ERβ individually, but when the LOs were co-exposed to PS and BPA, the ERα further elevated significantly and synergistically. Our findings highlight the metabolic-related health risk due to co-exposure to MP and BPA, even at low-doses equivalent to human internal exposure level. Based on these findings, the potential adverse outcome pathway related to PS and BPA singly and jointly were proposed, predicting two possible outcomes to be hepatic steatosis. Moreover, the ERα and HNF4A were proposed to be potential candidate markers to investigate the "vector-like effect" of PS in the present of BPA.

ID1 and CEBPA coordinate epidermal progenitor cell differentiation

The regulatory circuits that coordinate epidermal differentiation during development are still not fully understood. Here, we report that the transcriptional regulator ID1 is enriched in mouse basal epidermal progenitor cells and find ID1 expression to be diminished upon differentiation. In utero silencing of Id1 impairs progenitor cell proliferation, leads to precocious delamination of targeted progenitor cells and enables differentiated keratinocytes to retain progenitor markers and characteristics. Transcriptional profiling suggests that ID1 acts by mediating adhesion to the basement membrane while inhibiting spinous layer differentiation. Co-immunoprecipitation reveals ID1 binding to transcriptional regulators of the class I bHLH family. We localize bHLH Tcf3, Tcf4 and Tcf12 to epidermal progenitor cells during epidermal stratification and establish TCF3 as a downstream effector of ID1-mediated epidermal proliferation. Finally, we identify crosstalk between CEBPA, a known mediator of epidermal differentiation, and Id1, and demonstrate that CEBPA antagonizes BMP-induced activation of Id1. Our work establishes ID1 as a key coordinator of epidermal development, acting to balance progenitor proliferation with differentiation and unveils how functional crosstalk between CEBPA and Id1 orchestrates epidermal lineage progression.

A single-nucleus and spatial transcriptomic atlas of the COVID-19 liver reveals topological, functional, and regenerative organ disruption in patients

The molecular underpinnings of organ dysfunction in acute COVID-19 and its potential long-term sequelae are under intense investigation. To shed light on these in the context of liver function, we performed single-nucleus RNA-seq and spatial transcriptomic profiling of livers from 17 COVID-19 decedents. We identified hepatocytes positive for SARS-CoV-2 RNA with an expression phenotype resembling infected lung epithelial cells. Integrated analysis and comparisons with healthy controls revealed extensive changes in the cellular composition and expression states in COVID-19 liver, reflecting hepatocellular injury, ductular reaction, pathologic vascular expansion, and fibrogenesis. We also observed Kupffer cell proliferation and erythrocyte progenitors for the first time in a human liver single-cell atlas, resembling similar responses in liver injury in mice and in sepsis, respectively. Despite the absence of a clinical acute liver injury phenotype, endothelial cell composition was dramatically impacted in COVID-19, concomitantly with extensive alterations and profibrogenic activation of reactive cholangiocytes and mesenchymal cells. Our atlas provides novel insights into liver physiology and pathology in COVID-19 and forms a foundational resource for its investigation and understanding.

A single-nucleus and spatial transcriptomic atlas of the COVID-19 liver reveals topological, functional, and regenerative organ disruption in patients

The molecular underpinnings of organ dysfunction in acute COVID-19 and its potential long-term sequelae are under intense investigation. To shed light on these in the context of liver function, we performed single-nucleus RNA-seq and spatial transcriptomic profiling of livers from 17 COVID-19 decedents. We identified hepatocytes positive for SARS-CoV-2 RNA with an expression phenotype resembling infected lung epithelial cells. Integrated analysis and comparisons with healthy controls revealed extensive changes in the cellular composition and expression states in COVID-19 liver, reflecting hepatocellular injury, ductular reaction, pathologic vascular expansion, and fibrogenesis. We also observed Kupffer cell proliferation and erythrocyte progenitors for the first time in a human liver single-cell atlas, resembling similar responses in liver injury in mice and in sepsis, respectively. Despite the absence of a clinical acute liver injury phenotype, endothelial cell composition was dramatically impacted in COVID-19, concomitantly with extensive alterations and profibrogenic activation of reactive cholangiocytes and mesenchymal cells. Our atlas provides novel insights into liver physiology and pathology in COVID-19 and forms a foundational resource for its investigation and understanding.

PIDDosome-SCAP crosstalk controls high-fructose-diet-dependent transition from simple steatosis to steatohepatitis

Sterol deficiency triggers SCAP-mediated SREBP activation, whereas hypernutrition together with ER stress activates SREBP1/2 via caspase-2. Whether these pathways interact and how they are selectively activated by different dietary cues are unknown. Here, we reveal regulatory crosstalk between the two pathways that controls the transition from hepatosteatosis to steatohepatitis. Hepatic ER stress elicited by NASH-inducing diets activates IRE1 and induces expression of the PIDDosome subunits caspase-2, RAIDD, and PIDD1, along with INSIG2, an inhibitor of SCAP-dependent SREBP activation. PIDDosome assembly activates caspase-2 and sustains IRE1 activation. PIDDosome ablation or IRE1 inhibition blunt steatohepatitis and diminish INSIG2 expression. Conversely, while inhibiting simple steatosis, SCAP ablation amplifies IRE1 and PIDDosome activation and liver damage in NASH-diet-fed animals, effects linked to ER disruption and preventable by IRE1 inhibition. Thus, the PIDDosome and SCAP pathways antagonistically modulate nutrient-induced hepatic ER stress to control non-linear transition from simple steatosis to hepatitis, a key step in NASH pathogenesis.

Dynamics of hepatocyte-cholangiocyte cell-fate decisions during liver development and regeneration

The immense regenerative potential of the liver is attributed to the ability of its two key cell types – hepatocytes and cholangiocytes – to trans-differentiate to one another either directly or through intermediate progenitor states. However, the dynamic features of decision-making between these cell-fates during liver development and regeneration remains elusive. Here, we identify a core gene regulatory network comprising c/EBPα, TGFBR2, and SOX9 which is multistable in nature, enabling three distinct cell states – hepatocytes, cholangiocytes, and liver progenitor cells (hepatoblasts/oval cells) – and stochastic switching among them. Predicted expression signature for these three states are validated through multiple bulk and single-cell transcriptomic datasets collected across developmental stages and injury-induced liver repair. This network can also explain the experimentally observed spatial organization of phenotypes in liver parenchyma and predict strategies for efficient cellular reprogramming. Our analysis elucidates how the emergent dynamics of underlying regulatory networks drive diverse cell-fate decisions in liver development and regeneration.

•

Identified minimal regulatory network to model liver development and regeneration

•

Changes in phenotypic landscapes by in-silico perturbations of regulatory networks

•

Ability to explain physiological spatial patterning of liver cell types

•

Decoded strategies for efficient reprogramming among liver cell phenotypes

The potential role of exosomal miRNAs and membrane proteins in acute HIV-infected people

Exosomes play an important role during human immunodeficiency virus (HIV) acute infection. Yet, information regarding its cargo and its association with HIV rapid progressors (RPs) and typical progressors (TPs) remain largely unknown. In this study, exosomal miRNAs sequencing and mass cytometry were used to identify differential exosomal miRNAs and membrane proteins that participate in the pathogenesis of TPs and RPs. We discovered that miR-144-5p, miR-1180-3p, miR-451a, miR-362-5p, and miR-625-5p are associated with the TPs and miR-362-5p with the RPs. Decreased autophagy, amino acid metabolism, immune response, and IL-6 are closely related to RPs. In addition, SP1 was selected as the most significant transcription factor (TF) associated with disease progression. CD49D, CD5, CCR5, CD40, CD14, and CD86 were selected as the differential exosomal membrane proteins between TPs and RPs. This study provides valuable information for clarifying the mechanism in people with acute HIV infection.

SimiC enables the inference of complex gene regulatory dynamics across cell phenotypes

Single-cell RNA-Sequencing has the potential to provide deep biological insights by revealing complex regulatory interactions across diverse cell phenotypes at single-cell resolution. However, current single-cell gene regulatory network inference methods produce a single regulatory network per input dataset, limiting their capability to uncover complex regulatory relationships across related cell phenotypes. We present SimiC, a single-cell gene regulatory inference framework that overcomes this limitation by jointly inferring distinct, but related, gene regulatory dynamics per phenotype. We show that SimiC uncovers key regulatory dynamics missed by previously proposed methods across a range of systems, both model and non-model alike. In particular, SimiC was able to uncover CAR T cell dynamics after tumor recognition and key regulatory patterns on a regenerating liver, and was able to implicate glial cells in the generation of distinct behavioral states in honeybees. SimiC hence establishes a new approach to quantitating regulatory architectures between distinct cellular phenotypes, with far-reaching implications for systems biology.

ConsensusPathDB 2022: molecular interactions update as a resource for network biology

Molecular interactions are key drivers of biological function. Providing interaction resources to the research community is important since they allow functional interpretation and network-based analysis of molecular data. ConsensusPathDB (http://consensuspathdb.org) is a meta-database combining interactions of diverse types from 31 public resources for humans, 16 for mice and 14 for yeasts. Using ConsensusPathDB, researchers commonly evaluate lists of genes, proteins and metabolites against sets of molecular interactions defined by pathways, Gene Ontology and network neighborhoods and retrieve complex molecular neighborhoods formed by heterogeneous interaction types. Furthermore, the integrated protein–protein interaction network is used as a basis for propagation methods. Here, we present the 2022 update of ConsensusPathDB, highlighting content growth, additional functionality and improved database stability. For example, the number of human molecular interactions increased to 859 848 connecting 200 499 unique physical entities such as genes/proteins, metabolites and drugs. Furthermore, we integrated regulatory datasets in the form of transcription factor–, microRNA– and enhancer–gene target interactions, thus providing novel functionality in the context of overrepresentation and enrichment analyses. We specifically emphasize the use of the integrated protein–protein interaction network as a scaffold for network inferences, present topological characteristics of the network and discuss strengths and shortcomings of such approaches.

Integration of chronological omics data reveals mitochondrial regulatory mechanisms during the development of hepatocellular carcinoma

Mitochondria participate in multiple functions in eukaryotic cells. Although disruption of mitochondrial function has been associated with energetic deregulation in cancer, the chronological changes in mitochondria during cancer development remain unclear. With the aim to assess the role of mitochondria throughout cancer development, we analyzed samples chronologically obtained from induced hepatocellular carcinoma (HCC) in rats. In our analyses, we integrated mitochondrial proteomic data, mitochondrial metabolomic data and nuclear genome transcriptomic data. We used pathway over-representation and weighted gene co-expression network analysis (WGCNA) to integrate expression profiles of genes, miRNAs, proteins and metabolite levels throughout HCC development. Our results show that mitochondria are dynamic organelles presenting specific modifications in different stages of HCC development. We also found that mitochondrial proteomic profiles from tissues adjacent to nodules or tumor are determined more by the stage of HCC development than by tissue type, and we evaluated two models to predict HCC stage of the samples using proteomic profiles. Finally, we propose an omics integration pipeline to massively identify molecular features that could be further evaluated as key regulators, biomarkers or therapeutic targets. As an example, we show a group of miRNAs and transcription factors as candidates, responsible for mitochondrial metabolic modification in HCC.

Genome-wide and sister chromatid-resolved profiling of protein occupancy in replicated chromatin with ChOR-seq and SCAR-seq

Elucidating the mechanisms underlying chromatin maintenance upon genome replication is critical for the understanding of how gene expression programs and cell identity are preserved across cell divisions. Here, we describe two recently developed techniques, chromatin occupancy after replication (ChOR)-seq and sister chromatids after replication (SCAR)-seq, that profile chromatin occupancy on newly replicated DNA in mammalian cells in 5 d of bench work. Both techniques share a common strategy that includes pulse labeling of newly synthesized DNA and chromatin immunoprecipitation (ChIP), followed by purification and high-throughput sequencing. Whereas ChOR-seq quantitatively profiles the post-replicative abundance of histone modifications and chromatin-associated proteins, SCAR-seq distinguishes chromatin occupancy between nascent sister chromatids. Together, these two complementary techniques have unraveled key mechanisms controlling the inheritance of modified histones during replication and revealed locus-specific dynamics of histone modifications across the cell cycle. Here, we provide the experimental protocols and bioinformatic pipelines for these methods.

PRISMA and BioID disclose a motifs-based interactome of the intrinsically disordered transcription factor C/EBPα

C/EBPα represents a paradigm intrinsically disordered transcription factor containing short linear motifs and post-translational modifications (PTM). Unraveling C/EBPα protein interaction networks is a prerequisite for understanding the multi-modal functions of C/EBPα in hematopoiesis and leukemia. Here, we combined arrayed peptide matrix screening (PRISMA) with BioID to generate an in vivo validated and isoform specific interaction map of C/EBPα. The myeloid C/EBPα interactome comprises promiscuous and PTM-regulated interactions with protein machineries involved in gene expression, epigenetics, genome organization, DNA replication, RNA processing, and nuclear transport. C/EBPα interaction hotspots coincide with homologous conserved regions of the C/EBP family that also score as molecular recognition features. PTMs alter the interaction spectrum of C/EBP-motifs to configure a multi-valent transcription factor hub that interacts with multiple co-regulatory components, including BAF/SWI-SNF or Mediator complexes. Combining PRISMA and BioID is a powerful strategy to systematically explore the PTM-regulated interactomes of intrinsically disordered transcription factors.

Cellular plasticity balances the metabolic and proliferation dynamics of a regenerating liver

The adult liver has an exceptional ability to regenerate, but how it maintains its specialized functions during regeneration is unclear. Here, we used partial hepatectomy (PHx) in tandem with single-cell transcriptomics to track cellular transitions and heterogeneities of ∼22,000 liver cells through the initiation, progression, and termination phases of mouse liver regeneration. Our results uncovered that, following PHx, a subset of hepatocytes transiently reactivates an early-postnatal-like gene expression program to proliferate, while a distinct population of metabolically hyperactive cells appears to compensate for any temporary deficits in liver function. Cumulative EdU labeling and immunostaining of metabolic, portal, and central vein-specific markers revealed that hepatocyte proliferation after PHx initiates in the midlobular region before proceeding toward the periportal and pericentral areas. We further demonstrate that portal and central vein proximal hepatocytes retain their metabolically active state to preserve essential liver functions while midlobular cells proliferate nearby. Through combined analysis of gene regulatory networks and cell-cell interaction maps, we found that regenerating hepatocytes redeploy key developmental regulons, which are guided by extensive ligand-receptor-mediated signaling events between hepatocytes and nonparenchymal cells. Altogether, our study offers a detailed blueprint of the intercellular crosstalk and cellular reprogramming that balances the metabolic and proliferative requirements of a regenerating liver.

Cellular plasticity balances the metabolic and proliferation dynamics of a regenerating liver

The adult liver has an exceptional ability to regenerate, but how it maintains its specialized functions during regeneration is unclear. Here, we used partial hepatectomy (PHx) in tandem with single-cell transcriptomics to track cellular transitions and heterogeneities of ∼22,000 liver cells through the initiation, progression, and termination phases of mouse liver regeneration. Our results uncovered that, following PHx, a subset of hepatocytes transiently reactivates an early-postnatal-like gene expression program to proliferate, while a distinct population of metabolically hyperactive cells appears to compensate for any temporary deficits in liver function. Cumulative EdU labeling and immunostaining of metabolic, portal, and central vein-specific markers revealed that hepatocyte proliferation after PHx initiates in the midlobular region before proceeding toward the periportal and pericentral areas. We further demonstrate that portal and central vein proximal hepatocytes retain their metabolically active state to preserve essential liver functions while midlobular cells proliferate nearby. Through combined analysis of gene regulatory networks and cell-cell interaction maps, we found that regenerating hepatocytes redeploy key developmental regulons, which are guided by extensive ligand-receptor-mediated signaling events between hepatocytes and nonparenchymal cells. Altogether, our study offers a detailed blueprint of the intercellular crosstalk and cellular reprogramming that balances the metabolic and proliferative requirements of a regenerating liver.

Short- and long-term effects of perinatal phthalate exposures on metabolic pathways in the mouse liver

Phthalates have been demonstrated to interfere with metabolism, presumably by interacting with peroxisome proliferator-activated receptors (PPARs). However, mechanisms linking developmental phthalate exposures to long-term metabolic effects have not yet been elucidated. We investigated the hypothesis that developmental phthalate exposure has long-lasting impacts on PPAR target gene expression and DNA methylation to influence hepatic metabolic profiles across the life course. We utilized an established longitudinal mouse model of perinatal exposures to diethylhexyl phthalate and diisononyl phthalate, and a mixture of diethylhexyl phthalate+diisononyl phthalate. Exposure was through the diet and spanned from 2 weeks before mating until weaning at postnatal day 21 (PND21). Liver tissue was analyzed from the offspring of exposed and control mice at PND21 and in another cohort of exposed and control mice at 10 months of age. RNA-seq and pathway enrichment analyses indicated that acetyl-CoA metabolic processes were altered in diisononyl phthalate-exposed female livers at both PND21 and 10 months (FDR = 0.0018). Within the pathway, all 13 significant genes were potential PPAR target genes. Promoter DNA methylation was altered at three candidate genes, but persistent effects were only observed for Fasn. Targeted metabolomics indicated that phthalate-exposed females had decreased acetyl-CoA at PND21 and increased acetyl-CoA and acylcarnitines at 10 months. Together, our data suggested that perinatal phthalate exposures were associated with short- and long-term activation of PPAR target genes, which manifested as increased fatty acid production in early postnatal life and increased fatty acid oxidation in adulthood. This presents a novel molecular pathway linking developmental phthalate exposures and metabolic health outcomes.

Senescence-activated enhancer landscape orchestrates the senescence-associated secretory phenotype in murine fibroblasts

The three-dimensional configuration of the chromatin architecture is known to be crucial for alterations in the transcriptional network; however, the underlying mechanisms of epigenetic control of senescence-related gene expression by modulating the chromatin architecture remain unknown. Here, we demonstrate frequent chromosomal compartment switching during mouse embryonic fibroblasts (MEFs) replicative senescence as characterized by senescence-inactivated (SIAEs) and -activated enhancers (SAEs) in topologically associated domains (TADs). Mechanistically, SAEs are closely correlated with senescence-associated secretory phenotype (SASP) genes, which are a key transcriptional feature of an aging microenvironment that contributes to tumor progression, aging acceleration, and immunoinflammatory responses. Moreover, SAEs can positively regulate robust changes in SASP expression. The transcription factor CCAAT/enhancer binding protein α (C/EBPα) is capable of enhancing SAE activity, which accelerates the emergence of SAEs flanking SASPs and the secretion of downstream factors, contributing to the progression of senescence. Our results provide novel insight into the TAD-related control of SASP gene expression by revealing hierarchical roles of the chromatin architecture, transcription factors, and enhancer activity in the regulation of cellular senescence.

Self-Organizing Human Induced Pluripotent Stem Cell Hepatocyte 3D Organoids Inform the Biology of the Pleiotropic TRIB1 Gene

Establishment of a physiologically relevant human hepatocyte‐like cell system for in vitro translational research has been hampered by the limited availability of cell models that accurately reflect human biology and the pathophysiology of human disease. Here we report a robust, reproducible, and scalable protocol for the generation of hepatic organoids from human induced pluripotent stem cells (hiPSCs) using short exposure to nonengineered matrices. These hepatic organoids follow defined stages of hepatic development and express higher levels of early (hepatocyte nuclear factor 4A [HNF4A], prospero‐related homeobox 1 [PROX1]) and mature hepatic and metabolic markers (albumin, asialoglycoprotein receptor 1 [ASGR1], CCAAT/enhancer binding protein α [C/EBPα]) than two‐dimensional (2D) hepatocyte‐like cells (HLCs) at day 20 of differentiation. We used this model to explore the biology of the pleiotropic TRIB1 (Tribbles‐1) gene associated with a number of metabolic traits, including nonalcoholic fatty liver disease and plasma lipids. We used genome editing to delete the TRIB1 gene in hiPSCs and compared TRIB1‐deleted iPSC‐HLCs to isogenic iPSC‐HLCs under both 2D culture and three‐dimensional (3D) organoid conditions. Under conventional 2D culture conditions, TRIB1‐deficient HLCs showed maturation defects, with decreased expression of late‐stage hepatic and lipogenesis markers. In contrast, when cultured as 3D hepatic organoids, the differentiation defects were rescued, and a clear lipid‐related phenotype was noted in the TRIB1‐deficient induced pluripotent stem cell HLCs. Conclusion: This work supports the potential of genome‐edited hiPSC‐derived hepatic 3D organoids in exploring human hepatocyte biology, including the functional interrogation of genes identified through human genetic investigation.

Identifying crucial genes for prognosis in septic patients: Gene integration study based on PRISMA guidelines

BACKGROUND

Sepsis is a serious clinical condition with a poor prognosis, despite improvements in diagnosis and treatment.Therefore, novel biomarkers are necessary that can help with estimating prognosis and improving clinical outcomes of patients with sepsis.

METHODS

The gene expression profiles GSE54514 and GSE63042 were downloaded from the GEO database. DEGs were screened by t test after logarithmization of raw data; then, the common DEGs between the 2 gene expression profiles were identified by up-regulation and down-regulation intersection. The DEGs were analyzed using bioinformatics, and a protein-protein interaction (PPI) survival network was constructed using STRING. Survival curves were constructed to explore the relationship between core genes and the prognosis of sepsis patients based on GSE54514 data.

RESULTS

A total of 688 common DEGs were identified between survivors and non-survivors of sepsis, and 96 genes were involved in survival networks. The crucial genes Signal transducer and activator of transcription 5A (STAT5A), CCAAT/enhancer-binding protein beta (CEBPB), Myc proto-oncogene protein (MYC), and REL-associated protein (RELA) were identified and showed increased expression in sepsis survivors. These crucial genes had a positive correlation with patients' survival time according to the survival analysis.

CONCLUSIONS

Our findings indicate that the genes STAT5A, CEBPB, MYC, and RELA may be important in predicting the prognosis of sepsis patients.

Molecular and Morphological Evidence of Hepatotoxicity after Silver Nanoparticle Exposure: A Systematic Review, In Silico, and Ultrastructure Investigation

Silver nanoparticles (AgNPs) have been widely used in a variety of applications in innovative development; consequently, people are more exposed to this particle. Growing concern about toxicity from AgNP exposure has attracted greater attention, while questions about nanosilver-responsive genes and consequences for human health remain unanswered. By considering early detection and prevention of nanotoxicology at the genetic level, this study aimed to identify 1) changes in gene expression levels that could be potential indicators for AgNP toxicity and 2) morphological phenotypes correlating to toxicity of HepG2 cells. To detect possible nanosilver-responsive genes in xenogenic targeted organs, a comprehensive systematic literature review of changes in gene expression in HepG2 cells after AgNP exposure and in silico method, connection up- and down-regulation expression analysis of microarrays (CU-DREAM), were performed. In addition, cells were extracted and processed for transmission electron microscopy to examine ultrastructural alterations. From the Gene Expression Omnibus (GEO) Series database, we selected genes that were up- and down-regulated in AgNPs, but not up- and down-regulated in silver ion exposed cells, as nanosilver-responsive genes. HepG2 cells in the AgNP-treated group showed distinct ultrastructural alterations. Our results suggested potential representative gene data after AgNPs exposure provide insight into assessment and prediction of toxicity from nanosilver exposure.

Structural Features of Transcription Factors Associating with Nucleosome Binding

Fate-changing transcription factors (TFs) scan chromatin to initiate new genetic programs during cell differentiation and reprogramming. Yet the protein structure domains that allow TFs to target nucleosomal DNA remain unexplored. We screened diverse TFs for binding to nucleosomes containing motif-enriched sequences targeted by pioneer factors in vivo. FOXA1, OCT4, ASCL1/E12α, PU1, CEBPα, and ZELDA display a range of nucleosome binding affinities that correlate with their cell reprogramming potential. We further screened 593 full-length human TFs on protein microarrays against different nucleosome sequences, followed by confirmation in solution, to distinguish among factors that bound nucleosomes, such as the neuronal AP-2α/β/γ, versus factors that only bound free DNA. Structural comparisons of DNA binding domains revealed that efficient nucleosome binders use short anchoring α helices to bind DNA, whereas weak nucleosome binders use unstructured regions and/or β sheets. Thus, specific modes of DNA interaction allow nucleosome scanning that confers pioneer activity to transcription factors.

Epigenetic Compensation Promotes Liver Regeneration

Two major functions of the epigenome are to regulate gene expression and to suppress transposons. It is unclear how these functions are balanced during physiological challenges requiring tissue regeneration, where exquisite coordination of gene expression is essential. Transcriptomic analysis of seven time points following partial hepatectomy identified the epigenetic regulator UHRF1, which is essential for DNA methylation, as dynamically expressed during liver regeneration in mice. UHRF1 deletion in hepatocytes (Uhrf1HepKO) caused genome-wide DNA hypomethylation but, surprisingly, had no measurable effect on gene or transposon expression or liver homeostasis. Partial hepatectomy of Uhrf1HepKO livers resulted in early and sustained activation of proregenerative genes and enhanced liver regeneration. This was attributed to redistribution of H3K27me3 from promoters to transposons, effectively silencing them and, consequently, alleviating repression of liver regeneration genes, priming them for expression in Uhrf1HepKO livers. Thus, epigenetic compensation safeguards the genome against transposon activation, indirectly affecting gene regulation.

Identification and Conservation Analysis of Cis-Regulatory Elements in Pig Liver

The liver plays a key role in metabolism and affects pig production. However, the functional annotation of noncoding regions of the pig liver remains poorly understood. We revealed the landscape of cis-regulatory elements and their functional characterization in pig liver. We identified 102,373 cis-regulatory elements in the pig liver, including enhancers, promoters, super-enhancers, and broad H3K4me3 domains, and highlighted 26 core transcription regulatory factors in the pig liver as well. We found similarity of cis-regulatory elements among those of pigs, humans, and cattle. Despite the low proportion of functionally conserved enhancers (~30%) between pig and human liver tissue, ~78% of the pig liver enhancer orthologues sequence could play an enhancer role in other human tissues. Additionally, we observed that the ratio of consistent super-enhancer-associated genes was significantly higher than the ratio of functionally conserved super-enhancers. Approximately 54% of the core regulation factors driven by super-enhancers were consistent across the liver from these three species. Our pig liver annotation and functional characterization studies provide a system and resource for noncoding annotation for future gene regulatory studies in pigs. Furthermore, our study also showed the high level functional conservation of cis-regulatory elements in mammals; it also improved our understanding of regulation function of mammal cis-regulatory elements.

Mechanisms involved in the activation of C/EBPα by small activating RNA in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is generally accompanied by high mortality and low cure rate. CCAAT enhancer-binding proteins (CEBPs) are transcriptional regulators that play a key role in maintaining liver function. Altered expression of C/EBPα and C/EBPβ occurs in many tumours including HCC. saRNAs are small double-stranded RNAs that enhance target gene expression at the transcriptional level. In this report, we activate CEPBA with saRNAs and suppress CEBPB with siRNAs in cells that represent three different degrees of HCC. We performed functional assays to investigate the effects of enhancing C/EBPα and its downstream targets, p21 and albumin across these lines. We also used Mass-spectrometry (MS) subsequent to a ChIP pull-down assay to characterise the components of the protein complex involved in regulating saRNA function. Putative saRNA interacting protein candidates that were identified by MS were knocked-down with siRNAs to investigate its impact on saRNA activity. We confirmed CEBPA-saRNA decreased proliferation and migration in the differentiated lines (HepG3/Hep3B). The undifferentiated line (PLCPRF5) showed saRNA-induced increase in CEBPA but with no loss in proliferation. This effect was reversed when CEBPB was suppressed with CEBPB-siRNA. When interrogating saRNA mode of action; three saRNA interacting proteins, CTR9, HnRNPA2/B1 and DDX5 were identified by MS. Targeted knock-down of these two proteins (by siRNA) abrogated saRNA activity. This study provides insight into how different HCC lines are affected by CEBPA-saRNAs and that endogenous abundance of CEBPB and saRNA accessory proteins may dictate efficacy of CEBPA-saRNA when used in a therapeutic context.

Aberrant enhancer hypomethylation contributes to hepatic carcinogenesis through global transcriptional reprogramming

Hepatocellular carcinomas (HCC) exhibit distinct promoter hypermethylation patterns, but the epigenetic regulation and function of transcriptional enhancers remain unclear. Here, our affinity- and bisulfite-based whole-genome sequencing analyses reveal global enhancer hypomethylation in human HCCs. Integrative epigenomic characterization further pinpoints a recurrent hypomethylated enhancer of CCAAT/enhancer-binding protein-beta (C/EBPβ) which correlates with C/EBPβ over-expression and poorer prognosis of patients. Demethylation of C/EBPβ enhancer reactivates a self-reinforcing enhancer-target loop via direct transcriptional up-regulation of enhancer RNA. Conversely, deletion of this enhancer via CRISPR/Cas9 reduces C/EBPβ expression and its genome-wide co-occupancy with BRD4 at H3K27ac-marked enhancers and super-enhancers, leading to drastic suppression of driver oncogenes and HCC tumorigenicity. Hepatitis B X protein transgenic mouse model of HCC recapitulates this paradigm, as C/ebpβ enhancer hypomethylation associates with oncogenic activation in early tumorigenesis. These results support a causal link between aberrant enhancer hypomethylation and C/EBPβ over-expression, thereby contributing to hepatocarcinogenesis through global transcriptional reprogramming.

There are distinct hypermethylation patterns in gene promoters in hepatocellular carcinomas (HCCs). Here, the authors show that the enhancer of C/EBPβ is recurrently hypomethylated in human HCCs, recapitulating this in a transgenic murine model and linking aberrant enhancer hypomethylation to hepatocarcinogenesis.

Accurate Recycling of Parental Histones Reproduces the Histone Modification Landscape during DNA Replication

Chromatin organization is disrupted genome-wide during DNA replication. On newly synthesized DNA, nucleosomes are assembled from new naive histones and old modified histones. It remains unknown whether the landscape of histone post-translational modifications (PTMs) is faithfully copied during DNA replication or the epigenome is perturbed. Here we develop chromatin occupancy after replication (ChOR-seq) to determine histone PTM occupancy immediately after DNA replication and across the cell cycle. We show that H3K4me3, H3K36me3, H3K79me3, and H3K27me3 positional information is reproduced with high accuracy on newly synthesized DNA through histone recycling. Quantitative ChOR-seq reveals that de novo methylation to restore H3K4me3 and H3K27me3 levels occurs across the cell cycle with mark- and locus-specific kinetics. Collectively, this demonstrates that accurate parental histone recycling preserves positional information and allows PTM transmission to daughter cells while modification of new histones gives rise to complex epigenome fluctuations across the cell cycle that could underlie cell-to-cell heterogeneity.

Chromosome 19 miRNA cluster and CEBPB expression specifically mark and potentially drive triple negative breast cancers

Triple negative breast cancers (TNBCs) are known to express low PGR, ESR1, and ERBB2, and high KRT5, KRT14, and KRT17. However, the reasons behind the increased expressions of KRT5, KRT14, KRT17 and decreased expressions of PGR, ESR1, and ERBB2 in TNBCs are not fully understood. Here we show that, expression of chromosome 19 miRNA cluster (C19MC) specifically marks human TNBCs. Low REST and high CEBPB correlate with expression of C19MC, KRT5, KRT14, and KRT17 and enhancers of these genes/cluster are regulated by CEBPB and REST binding sites. The C19MC miRNAs in turn can potentially target REST to offer a positive feedback loop, and might target PGR, ESR1, ERBB2, GATA3, SCUBE2, TFF3 mRNAs to contribute towards TNBC phenotype. Thus our study demonstrates that C19MC miRNA expression marks TNBCs and that C19MC miRNAs and CEBPB might together determine the TNBC marker expression pattern.

Cycles of gene expression and genome response during mammalian tissue regeneration

Background

Compensatory liver hyperplasia—or regeneration—induced by two-thirds partial hepatectomy (PH) permits the study of synchronized activation of mammalian gene expression, particularly in relation to cell proliferation. Here, we measured genomic transcriptional responses and mRNA accumulation changes after PH and sham surgeries.

Results

During the first 10–20 h, the PH- and sham-surgery responses were very similar, including parallel early activation of cell-division-cycle genes. After 20 h, however, whereas post-PH livers continued with a robust and coordinate cell-division-cycle gene-expression response before returning to the resting state by 1 week, sham-surgery livers returned directly to a resting gene-expression state. Localization of RNA polymerase II (Pol II), and trimethylated histone H3 lysine 4 (H3K4me3) and 36 (H3K36me3) on genes dormant in the resting liver and activated during the PH response revealed a general de novo promoter Pol II recruitment and H3K4me3 increase during the early 10–20 h phase followed by Pol II elongation and H3K36me3 accumulation in gene bodies during the later proliferation phase. H3K36me3, generally appearing at the first internal exon, was preceded 5′ by H3K36me2; 3′ of the first internal exon, in about half of genes H3K36me3 predominated and in the other half H3K36me2 and H3K36me3 co-existed. Further, we observed some unusual gene profiles with abundant Pol II but little evident H3K4me3 or H3K36me3 modification, indicating that these modifications are neither universal nor essential partners to Pol II transcription.

Conclusions

PH and sham surgical procedures on mice reveal striking early post-operatory gene expression similarities followed by synchronized mRNA accumulation and epigenetic histone mark changes specific to PH.

Electronic supplementary material

The online version of this article (10.1186/s13072-018-0222-0) contains supplementary material, which is available to authorized users.

MMARGE: Motif Mutation Analysis for Regulatory Genomic Elements

Cell-specific patterns of gene expression are determined by combinatorial actions of sequence-specific transcription factors at cis-regulatory elements. Studies indicate that relatively simple combinations of lineage-determining transcription factors (LDTFs) play dominant roles in the selection of enhancers that establish cell identities and functions. LDTFs require collaborative interactions with additional transcription factors to mediate enhancer function, but the identities of these factors are often unknown. We have shown that natural genetic variation between individuals has great utility for discovering collaborative transcription factors. Here, we introduce MMARGE (Motif Mutation Analysis of Regulatory Genomic Elements), the first publicly available suite of software tools that integrates genome-wide genetic variation with epigenetic data to identify collaborative transcription factor pairs. MMARGE is optimized to work with chromatin accessibility assays (such as ATAC-seq or DNase I hypersensitivity), as well as transcription factor binding data collected by ChIP-seq. Herein, we provide investigators with rationale for each step in the MMARGE pipeline and key differences for analysis of datasets with different experimental designs. We demonstrate the utility of MMARGE using mouse peritoneal macrophages, liver cells, and human lymphoblastoid cells. MMARGE provides a powerful tool to identify combinations of cell type-specific transcription factors while simultaneously interpreting functional effects of non-coding genetic variation.

Responses of Transgenic Melatonin-Enriched Goats on LPS Stimulation and the Proteogenomic Profiles of Their PBMCs

The anti-inflammatory activity of melatonin (MT) has been well documented; however, little is known regarding endogenously occurring MT in this respect, especially for large animals. In the current study, we created a MT-enriched animal model (goats) overexpressing the MT synthetase gene Aanat. The responses of these animals to lipopolysaccharide (LPS) stimulation were systematically studied. It was found that LPS treatment exacerbated the inflammatory response in wild-type (WT) goats and increased their temperature to 40 °C. In addition, their granulocyte counts were also significantly elevated. In contrast, these symptoms were not observed in transgenic goats with LPS treatment. The rescue study with MT injection into WT goats who were treated with LPS confirmed that the protective effects in transgenic goats against LPS were attributed to a high level of endogenously produced MT. The proteomic analysis in the peripheral blood mononuclear cells (PBMCs) isolated from the transgenic animals uncovered several potential mechanisms. MT suppressed the lysosome formation as well as its function by downregulation of the lysosome-associated genes Lysosome-associated membrane protein 2 (LAMP2), Insulin-like growth factor 2 receptor (IGF2R), and Arylsulfatase B (ARSB). A high level of MT enhanced the antioxidant capacity of these cells to reduce the cell apoptosis induced by the LPS. In addition, the results also uncovered previously unknown information that showed that MT may have protective effects on some human diseases, including tuberculosis, bladder cancer, and rheumatoid arthritis, by downregulation of these disease-associated genes. All these observations warranted further investigations.

Transcription factor profiling reveals molecular choreography and key regulators of human retrotransposon expression

Transposable elements (TEs) represent a substantial fraction of many eukaryotic genomes, and transcriptional regulation of these factors is important to determine TE activities in human cells. However, due to the repetitive nature of TEs, identifying transcription factor (TF)-binding sites from ChIP-sequencing (ChIP-seq) datasets is challenging. Current algorithms are focused on subtle differences between TE copies and thus bias the analysis to relatively old and inactive TEs. Here we describe an approach termed "MapRRCon" (mapping repeat reads to a consensus) which allows us to identify proteins binding to TE DNA sequences by mapping ChIP-seq reads to the TE consensus sequence after whole-genome alignment. Although this method does not assign binding sites to individual insertions in the genome, it provides a landscape of interacting TFs by capturing factors that bind to TEs under various conditions. We applied this method to screen TFs' interaction with L1 in human cells/tissues using ENCODE ChIP-seq datasets and identified 178 of the 512 TFs tested as bound to L1 in at least one biological condition with most of them (138) localized to the promoter. Among these L1-binding factors, we focused on Myc and CTCF, as they play important roles in cancer progression and 3D chromatin structure formation. Furthermore, we explored the transcriptomes of The Cancer Genome Atlas breast and ovarian tumor samples in which a consistent anti-/correlation between L1 and Myc/CTCF expression was observed, suggesting that these two factors may play roles in regulating L1 transcription during the development of such tumors.

Computational modeling of in vivo and in vitro protein-DNA interactions by multiple instance learning

Motivation

The study of transcriptional regulation is still difficult yet fundamental in molecular biology research. While the development of both in vivo and in vitro profiling techniques have significantly enhanced our knowledge of transcription factor (TF)-DNA interactions, computational models of TF-DNA interactions are relatively simple and may not reveal sufficient biological insight. In particular, supervised learning based models for TF-DNA interactions attempt to map sequence-level features ( k -mers) to binding event but usually ignore the location of k -mers, which can cause data fragmentation and consequently inferior model performance.

Results

Here, we propose a novel algorithm based on the so-called multiple-instance learning (MIL) paradigm. MIL breaks each DNA sequence into multiple overlapping subsequences and models each subsequence separately, therefore implicitly takes into consideration binding site locations, resulting in both higher accuracy and better interpretability of the models. The result from both in vivo and in vitro TF-DNA interaction data show that our approach significantly outperform conventional single-instance learning based algorithms. Importantly, the models learned from in vitro data using our approach can predict in vivo binding with very good accuracy. In addition, the location information obtained by our method provides additional insight for motif finding results from ChIP-Seq data. Finally, our approach can be easily combined with other state-of-the-art TF-DNA interaction modeling methods.

Availability and Implementation

http://www.cs.utsa.edu/∼jruan/MIL/.

Contact

jianhua.ruan@utsa.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

Genome-Wide Mapping of DNA Accessibility and Binding Sites for CREB and C/EBPβ in Vasopressin-Sensitive Collecting Duct Cells

Background Renal water excretion is controlled by vasopressin, in part through regulation of the transcription of the aquaporin-2 gene (Aqp2).Methods To identify enhancer regions likely to be involved in the regulation of Aqp2 and other principal cell-specific genes, we used several next generation DNA-sequencing techniques in a well characterized cultured cell model of collecting duct principal cells (mpkCCD). To locate enhancers, we performed the assay for transposase-accessible chromatin using sequencing (ATAC-Seq) to identify accessible regions of DNA and integrated the data with data generated by chromatin immunoprecipitation followed by next generation DNA-sequencing (ChIP-Seq) for CCCTC binding factor (CTCF) binding, histone H3 lysine-27 acetylation, and RNA polymerase II.Results We identified two high-probability enhancers centered 81 kb upstream and 5.8 kb downstream from the Aqp2 transcriptional start site. Motif analysis of these regions and the Aqp2 promoter identified several potential transcription factor binding sites, including sites for two b-ZIP transcription factors: CCAAT/enhancer binding protein-β (C/EBPβ) and cAMP-responsive element binding protein (CREB). To identify genomic binding sites for both, we conducted ChIP-Seq using well characterized antibodies. In the presence of vasopressin, C/EBPβ, a pioneer transcription factor critical to cell-specific gene expression, bound strongly at the identified enhancer downstream from Aqp2 However, over multiple replicates, we found no detectable CREB binding sites within 390 kb of Aqp2 Thus, any role for CREB in the regulation of Aqp2 gene transcription is likely to be indirect.Conclusions The analysis identified two enhancer regions pertinent to transcriptional regulation of the Aqp2 gene and showed C/EBPβ (but not CREB) binding.

A screening system to identify transcription factors that induce binding site-directed DNA demethylation

Background

DNA methylation is a fundamental epigenetic modification that is involved in many biological systems such as differentiation and disease. We and others recently showed that some transcription factors (TFs) are involved in the site-specific determination of DNA demethylation in a binding site-directed manner, although the reports of such TFs are limited.

Results

Here, we develop a screening system to identify TFs that induce binding site-directed DNA methylation changes. The system involves the ectopic expression of target TFs in model cells followed by DNA methylome analysis and overrepresentation analysis of the corresponding TF binding motif at differentially methylated regions. It successfully identified binding site-directed demethylation of SPI1, which is known to promote DNA demethylation in a binding site-directed manner. We extended our screening system to 15 master TFs involved in cellular differentiation and identified eight novel binding site-directed DNA demethylation-inducing TFs (RUNX3, GATA2, CEBPB, MAFB, NR4A2, MYOD1, CEBPA, and TBX5). Gene ontology and tissue enrichment analysis revealed that these TFs demethylate genomic regions associated with corresponding biological roles. We also describe the characteristics of binding site-directed DNA demethylation induced by these TFs, including the targeting of highly methylated CpGs, local DNA demethylation, and the overlap of demethylated regions between TFs of the same family.

Conclusions

Our results show the usefulness of the developed screening system for the identification of TFs that induce DNA demethylation in a site-directed manner.

Electronic supplementary material

The online version of this article (10.1186/s13072-017-0169-6) contains supplementary material, which is available to authorized users.