Induction of AML cell differentiation using HOXA9/DNA binding inhibitors as a potential therapeutic option for HOXA9-dependent AML.

The mainstay of acute myeloid leukemia (AML) treatment still relies on traditional chemotherapy, with a survival rate of approximately 30% for patients under 65 years of age and as low as 5% for those beyond. This unfavorable prognosis primarily stems from frequent relapses, resistance to chemotherapy, and limited approved targeted therapies for specific AML subtypes. Around 70% of all AML cases show overexpression of the transcription factor HOXA9, which is associated with a poor prognosis, increased chemoresistance, and higher relapse rates. However, direct targeting of HOXA9 in a clinical setting has not been achieved yet. The dysregulation caused by the leukemic HOXA9 transcription factor primarily results from its binding activity to DNA, leading to differentiation blockade. Our previous investigations have identified two HOXA9/DNA binding competitors, namely DB1055 and DB818. We assessed their antileukemic effects in comparison to HOXA9 knockdown or cytarabine treatment. Using human AML cell models, DB1055 and DB818 induced in vitro cell growth reduction, death, differentiation, and common transcriptomic deregulation but did not impact human CD34+ bone marrow cells. Furthermore, DB1055 and DB818 exhibited potent antileukemic activities in a human THP-1 AML in vivo model, leading to the differentiation of monocytes into macrophages. In vitro assays also demonstrated the efficacy of DB1055 and DB818 against AML blasts from patients, with DB1055 successfully reducing leukemia burden in patient-derived xenografts in NSG immunodeficient mice. Our findings indicate that inhibiting HOXA9/DNA interaction using DNA ligands may offer a novel differentiation therapy for the future treatment of AML patients dependent on HOXA9.

Performance of non-invasive tests for liver fibrosis resolution after bariatric surgery.

BACKGROUND & AIMS: The value of non-invasive tests for monitoring the resolution of significant liver fibrosis after treatment is poorly investigated. We compared the performances of six non-invasive tests to predict the resolution of significant fibrosis after bariatric surgery. METHODS: Participants were individuals with obesity submitted to needle liver biopsy at the time of bariatric surgery, and 12 and/or 60 months after surgery. We calculated the fibrosis-4 index (FIB-4), NAFLD fibrosis score (NFS), AST to platelet ratio index (APRI), Hepatic fibrosis score (HFS), Fibrotic NASH index (FNI), and Liver risk score (LRS) at each time point, and compared their performances for predicting significant fibrosis (F ≥ 2) and its resolution following surgery. RESULTS: At baseline, 2436 patients had liver biopsy, including 261 (10.7 %) with significant fibrosis. Overall, 672 patients had pre- and post-operative biopsies (564 at M12 and 328 at M60). The fibrosis stage decreased at M12 and M60 (p < 0.001 vs M0). Resolution of significant fibrosis occurred in 58/121 (47.9 %) at M12 and 32/50 (64 %) at M60. The mean value of all tests decreased after surgery, except for FIB-4. Performances for predicting fibrosis resolution was higher at M60 than at M12 for all tests, and maximal at M60 for FNI and LRS: area under the curve 0.843 (95%CI 0.71-0.95) and 0.92 (95%CI 0.84-1.00); positive likelihood ratio 3.75 (95 % CI 1.33-10.59) and 4.58 (95 % CI 1.65-12.70), respectively. CONCLUSIONS: Results showed the value and limits of non-invasive tests for monitoring the evolution of liver fibrosis after an intervention. Following bariatric surgery, the best performances to predict the resolution of significant fibrosis were observed at M60 with tests combining liver and metabolic traits, namely FNI and LRS.

Combining diabetes, sex, and menopause as meaningful clinical features associated with NASH and liver fibrosis in individuals with class II and III obesity: A retrospective cohort study.

OBJECTIVE: Steatotic liver disease (SLD) is frequent in individuals with obesity. In this study, type 2 diabetes (T2D), sex, and menopausal status were combined to refine the stratification of obesity regarding the risk of advanced SLD and gain further insight into disease physiopathology. METHODS: This study enrolled 1446 participants with obesity from the ABOS cohort (NCT01129297), who underwent extensive phenotyping, including liver histology and transcriptome profiling. Hierarchical clustering was applied to classify participants. The prevalence of metabolic disorders associated with steatohepatitis (NASH) and liver fibrosis (F ≥ 2) was determined within each identified subgroup and aligned to clinical and biological characteristics. RESULTS: The prevalence of NASH and F ≥ 2 was, respectively, 9.5% (N = 138/1446) and 11.7% (N = 159/1365) in the overall population, 20.3% (N = 107/726) and 21.1% (N = 106/502) in T2D patients, and 3.4% (N = 31/920) and 6.1% (N = 53/863) in non-T2D patients. NASH and F ≥ 2 prevalence was 15.4% (33/215) and 15.5% (32/206) among premenopausal women with T2D vs. 29.5% (33/112) and 30.3% (N = 36/119) in postmenopausal women with T2D (p < 0.01); and 21.0% (21/100) / 27.0% (24/89) in men with T2D ≥ age 50 years and 17.9% (17/95) / 18.5% (17/92) in men with T2D < age 50 years (NS). The distinct contribution of menopause was confirmed by the interaction between sex and age with respect to NASH among T2D patients (p = 0.048). Finally, several NASH-associated biological traits (lower platelet count; higher serum uric acid; gamma-glutamyl transferase; aspartate aminotransferase) and liver expressed genes AKR1B10 and CCL20 were significantly associated with menopause in women with T2D but not with age in men with T2D. CONCLUSIONS: This study unveiled a remarkably high prevalence of advanced SLD after menopause in women with T2D, associated with a dysfunctional biological liver profile.

Hepatitis C virus alters the morphology and function of peroxisomes.

Despite the introduction of effective treatments for hepatitis C in clinics, issues remain regarding the liver disease induced by chronic hepatitis C virus (HCV) infection. HCV is known to disturb the metabolism of infected cells, especially lipid metabolism and redox balance, but the mechanisms leading to HCV-induced pathogenesis are still poorly understood. In an APEX2-based proximity biotinylation screen, we identified ACBD5, a peroxisome membrane protein, as located in the vicinity of HCV replication complexes. Confocal microscopy confirmed the relocation of peroxisomes near HCV replication complexes and indicated that their morphology and number are altered in approximately 30% of infected Huh-7 cells. Peroxisomes are small versatile organelles involved among other functions in lipid metabolism and ROS regulation. To determine their importance in the HCV life cycle, we generated Huh-7 cells devoid of peroxisomes by inactivating the PEX5 and PEX3 genes using CRISPR/Cas9 and found that the absence of peroxisomes had no impact on replication kinetics or infectious titers of HCV strains JFH1 and DBN3a. The impact of HCV on peroxisomal functions was assessed using sub-genomic replicons. An increase of ROS was measured in peroxisomes of replicon-containing cells, correlated with a significant decrease of catalase activity with the DBN3a strain. In contrast, HCV replication had little to no impact on cytoplasmic and mitochondrial ROS, suggesting that the redox balance of peroxisomes is specifically impaired in cells replicating HCV. Our study provides evidence that peroxisome function and morphology are altered in HCV-infected cells.

An extended transcription factor regulatory network controls hepatocyte identity.

Cell identity is specified by a core transcriptional regulatory circuitry (CoRC), typically limited to a small set of interconnected cell-specific transcription factors (TFs). By mining global hepatic TF regulons, we reveal a more complex organization of the transcriptional regulatory network controlling hepatocyte identity. We show that tight functional interconnections controlling hepatocyte identity extend to non-cell-specific TFs beyond the CoRC, which we call hepatocyte identity (Hep-ID)CONNECT TFs. Besides controlling identity effector genes, Hep-IDCONNECT TFs also engage in reciprocal transcriptional regulation with TFs of the CoRC. In homeostatic basal conditions, this translates into Hep-IDCONNECT TFs being involved in fine tuning CoRC TF expression including their rhythmic expression patterns. Moreover, a role for Hep-IDCONNECT TFs in the control of hepatocyte identity is revealed in dedifferentiated hepatocytes where Hep-IDCONNECT TFs are able to reset CoRC TF expression. This is observed upon activation of NR1H3 or THRB in hepatocarcinoma or in hepatocytes subjected to inflammation-induced loss of identity. Our study establishes that hepatocyte identity is controlled by an extended array of TFs beyond the CoRC.

A Novel 8-Predictors Signature to Predict Complicated Disease Course in Pediatric-onset Crohn's Disease: A Population-based Study.

BACKGROUND: The identification of patients at high risk of a disabling disease course would be invaluable in guiding initial therapy in Crohn's disease (CD). Our objective was to evaluate a combination of clinical, serological, and genetic factors to predict complicated disease course in pediatric-onset CD. METHODS: Data for pediatric-onset CD patients, diagnosed before 17 years of age between 1988 and 2004 and followed more than 5 years, were extracted from the population-based EPIMAD registry. The main outcome was defined by the occurrence of complicated behavior (stricturing or penetrating) and/or intestinal resection within the 5 years following diagnosis. Lasso logistic regression models were used to build a predictive model based on clinical data at diagnosis, serological data (ASCA, pANCA, anti-OmpC, anti-Cbir1, anti-Fla2, anti-Flax), and 369 candidate single nucleotide polymorphisms. RESULTS: In total, 156 children with an inflammatory (B1) disease at diagnosis were included. Among them, 35% (n = 54) progressed to a complicated behavior or an intestinal resection within the 5 years following diagnosis. The best predictive model (PREDICT-EPIMAD) included the location at diagnosis, pANCA, and 6 single nucleotide polymorphisms. This model showed good discrimination and good calibration, with an area under the curve of 0.80 after correction for optimism bias (sensitivity, 79%, specificity, 74%, positive predictive value, 61%, negative predictive value, 87%). Decision curve analysis confirmed the clinical utility of the model. CONCLUSIONS: A combination of clinical, serotypic, and genotypic variables can predict disease progression in this population-based pediatric-onset CD cohort. Independent validation is needed before it can be used in clinical practice.

Identification of patient subtypes based on protein expression for prediction of heart failure after myocardial infarction.

This study investigates the ability of high-throughput aptamer-based platform to identify circulating biomarkers able to predict occurrence of heart failure (HF), in blood samples collected during hospitalization of patients suffering from a first myocardial infarction (MI). REVE-1 (derivation) and REVE-2 (validation) cohorts included respectively 254 and 238 patients, followed up respectively 9 · 2 ± 4 · 8 and 7 · 6 ± 3 · 0 years. A blood sample collected during hospitalization was used for quantifying 4,668 proteins. Fifty proteins were significantly associated with long-term occurrence of HF with all-cause death as the competing event. k-means, an unsupervised clustering method, identified two groups of patients based on expression levels of the 50 proteins. Group 2 was significantly associated with a higher risk of HF in both cohorts. These results showed that a subset of 50 selected proteins quantified during hospitalization of MI patients is able to stratify and predict the long-term occurrence of HF.

Biomarkers of haemodynamic severity of systemic sclerosis-associated pulmonary arterial hypertension by serum proteome analysis.

OBJECTIVES: To mine the serum proteome of patients with systemic sclerosis-associated pulmonary arterial hypertension (SSc-PAH) and to detect biomarkers that may assist in earlier and more effective diagnosis and treatment. METHODS: Patients with limited cutaneous SSc, no extensive interstitial lung disease and no PAH-specific therapy were included. They were classified as cases if they had PAH confirmed by right heart catheterisation (RHC) and serum collected on the same day as RHC; and as controls if they had no clinical evidence of PAH. RESULTS: Patients were mostly middle-aged females with anticentromere-associated SSc. Among 1129 proteins assessed by a high-throughput proteomic assay (SOMAscan), only 2 were differentially expressed and correlated significantly with pulmonary vascular resistance (PVR) in SSc-PAH patients (n=15): chemerin (ρ=0.62, p=0.01) and SET (ρ=0.62, p=0.01). To validate these results, serum levels of chemerin were measured by ELISA in an independent cohort. Chemerin levels were confirmed to be significantly higher (p=0.01) and correlate with PVR (ρ=0.42, p=0.04) in SSc-PAH patients (n=24). Chemerin mRNA expression was detected in fibroblasts, pulmonary artery smooth muscle cells (PA-SMCs)/pericytes and mesothelial cells in SSc-PAH lungs by single-cell RNA-sequencing. Confocal immunofluorescence revealed increased expression of a chemerin receptor, CMKLR1, on SSc-PAH PA-SMCs. SSc-PAH serum seemed to induce higher PA-SMC proliferation than serum from SSc patients without PAH. This difference appeared neutralised when adding the CMKLR1 inhibitor α-NETA. CONCLUSION: Chemerin seems an interesting surrogate biomarker for PVR in SSc-PAH. Increased chemerin serum levels and CMKLR1 expression by PA-SMCs may contribute to SSc-PAH pathogenesis by inducing PA-SMC proliferation.

Effects of Immunoglobulins G From Systemic Sclerosis Patients in Normal Dermal Fibroblasts: A Multi-Omics Study.

Autoantibodies (Aabs) are frequent in systemic sclerosis (SSc). Although recognized as potent biomarkers, their pathogenic role is debated. This study explored the effect of purified immunoglobulin G (IgG) from SSc patients on protein and mRNA expression of dermal fibroblasts (FBs) using an innovative multi-omics approach. Dermal FBs were cultured in the presence of sera or purified IgG from patients with diffuse cutaneous SSc (dcSSc), limited cutaneous SSc or healthy controls (HCs). The FB proteome and transcriptome were explored using liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS) and microarray assays, respectively. Proteomic analysis identified 3,310 proteins. SSc sera and purified IgG induced singular protein profile patterns. These FB proteome changes depended on the Aab serotype, with a singular effect observed with purified IgG from anti-topoisomerase-I autoantibody (ATA) positive patients compared to HC or other SSc serotypes. IgG from ATA positive SSc patients induced enrichment in proteins involved in focal adhesion, cadherin binding, cytosolic part, or lytic vacuole. Multi-omics analysis was performed in two ways: first by restricting the analysis of the transcriptomic data to differentially expressed proteins; and secondly, by performing a global statistical analysis integrating proteomics and transcriptomics. Transcriptomic analysis distinguished 764 differentially expressed genes and revealed that IgG from dcSSc can induce extracellular matrix (ECM) remodeling changes in gene expression profiles in FB. Global statistical analysis integrating proteomics and transcriptomics confirmed that IgG from SSc can induce ECM remodeling and activate FB profiles. This effect depended on the serotype of the patient, suggesting that SSc Aab might play a pathogenic role in some SSc subsets.

TgAP2IX-5 is a key transcriptional regulator of the asexual cell cycle division in Toxoplasma gondii.

Apicomplexan parasites have evolved efficient and distinctive strategies for intracellular replication where the timing of emergence of the daughter cells (budding) is a decisive element. However, the molecular mechanisms that provide the proper timing of parasite budding remain unknown. Using Toxoplasma gondii as a model Apicomplexan, we identified a master regulator that controls the timing of the budding process. We show that an ApiAP2 transcription factor, TgAP2IX-5, controls cell cycle events downstream of centrosome duplication. TgAP2IX-5 binds to the promoter of hundreds of genes and controls the activation of the budding-specific cell cycle expression program. TgAP2IX-5 regulates the expression of specific transcription factors that are necessary for the completion of the budding cycle. Moreover, TgAP2IX-5 acts as a limiting factor that ensures that asexual proliferation continues by promoting the inhibition of the differentiation pathway. Therefore, TgAP2IX-5 is a master regulator that controls both cell cycle and developmental pathways.

Circulating proteomic signature of early death in heart failure patients with reduced ejection fraction.

Heart failure (HF) remains a main cause of mortality worldwide. Risk stratification of patients with systolic chronic HF is critical to identify those who may benefit from advanced HF therapies. The aim of this study is to identify plasmatic proteins that could predict the early death (within 3 years) of HF patients with reduced ejection fraction hospitalized in CHRU de Lille. The subproteome targeted by an aptamer-based technology, the Slow Off-rate Modified Aptamer (SOMA) scan assay of 1310 proteins, was profiled in blood samples from 168 HF patients, and 203 proteins were significantly modulated between patients who died of cardiovascular death and patients who were alive after 3 years of HF evaluation (Wilcoxon test, FDR 5%). A molecular network was built using these 203 proteins, and the resulting network contained 2281 molecules assigned to 34 clusters annotated to biological pathways by Gene Ontology. This network model highlighted extracellular matrix organization as the main mechanism involved in early death in HF patients. In parallel, an adaptive Least Absolute Shrinkage and Selection Operator (LASSO) was performed on these 203 proteins, and six proteins were selected as candidates to predict early death in HF patients: complement C3, cathepsin S and F107B were decreased and MAPK5, MMP1 and MMP7 increased in patients who died of cardiovascular causes compared with patients living 3 years after HF evaluation. This proteomic signature of 6 circulating plasma proteins allows the identification of systolic HF patients with a risk of early death.

Transient Receptor Potential Channel Expression Signatures in Tumor-Derived Endothelial Cells: Functional Roles in Prostate Cancer Angiogenesis.

Background: Transient receptor potential (TRP) channels control multiple processes involved in cancer progression by modulating cell proliferation, survival, invasion and intravasation, as well as, endothelial cell (EC) biology and tumor angiogenesis. Nonetheless, a complete TRP expression signature in tumor vessels, including in prostate cancer (PCa), is still lacking. Methods: In the present study, we profiled by qPCR the expression of all TRP channels in human prostate tumor-derived ECs (TECs) in comparison with TECs from breast and renal tumors. We further functionally characterized the role of the 'prostate-associated' channels in proliferation, sprout formation and elongation, directed motility guiding, as well as in vitro and in vivo morphogenesis and angiogenesis. Results: We identified three 'prostate-associated' genes whose expression is upregulated in prostate TECs: TRPV2 as a positive modulator of TEC proliferation, TRPC3 as an endothelial PCa cell attraction factor and TRPA1 as a critical TEC angiogenic factor in vitro and in vivo. Conclusions: We provide here the full TRP signature of PCa vascularization among which three play a profound effect on EC biology. These results contribute to explain the aggressive phenotype previously observed in PTEC and provide new putative therapeutic targets.

Metabolic and Innate Immune Cues Merge into a Specific Inflammatory Response via the UPR.

Innate immune responses are intricately linked with intracellular metabolism of myeloid cells. Toll-like receptor (TLR) stimulation shifts intracellular metabolism toward glycolysis, while anti-inflammatory signals depend on enhanced mitochondrial respiration. How exogenous metabolic signals affect the immune response is unknown. We demonstrate that TLR-dependent responses of dendritic cells (DCs) are exacerbated by a high-fatty-acid (FA) metabolic environment. FAs suppress the TLR-induced hexokinase activity and perturb tricarboxylic acid cycle metabolism. These metabolic changes enhance mitochondrial reactive oxygen species (mtROS) production and, in turn, the unfolded protein response (UPR), leading to a distinct transcriptomic signature with IL-23 as hallmark. Interestingly, chemical or genetic suppression of glycolysis was sufficient to induce this specific immune response. Conversely, reducing mtROS production or DC-specific deficiency in XBP1 attenuated IL-23 expression and skin inflammation in an IL-23-dependent model of psoriasis. Thus, fine-tuning of innate immunity depends on optimization of metabolic demands and minimization of mtROS-induced UPR.

The miRNA-targeted transcriptome of porcine alveolar macrophages upon infection with Porcine Reproductive and Respiratory Syndrome Virus.

Host miRNAs are known to modulate the cell response to virus infections. We characterized the miRNA-targeted transcriptome of porcine alveolar macrophages (PAMs) at early times after infection with a subtype 1.1 strain of PRRSV (Porcine Reproductive and Respiratory Syndrome Virus). We performed the immunoprecipitation of RISC (RNA-induced Silencing Complex) followed by microarray analysis of the RISC-bound miRNA targets (RIP-Chip) to evaluate the relative enrichment or depletion of expressed genes in RISC. The miRNA-mediated regulation occurred early after PRRSV infection and decreased fast (1,241 and 141 RISC-bound genes at 7 h and 10 h post-infection, respectively); it affected several cell functions with evidence of miRNA buffering of upregulated interferon-related genes. Eight miRNAs were highly enriched in RISC of both control and infected cells with no evidence of differential expression. Although miR-335-5p was the miRNA with most predicted targets among enriched RISC-bound genes, no effects on surface markers, cytokine expression and PRRSV replication were detected upon miR-335-5p mimics of primary PAMs. Our results do not point to specific miRNA-driven mechanisms regulating the early response to infection with this PRRSV 1.1 strain and indicate that the miRNome expressed by steady-state PAMs reacts promptly to counterbalance PRRSV infection by a pervasive modulation of host functions.

SMAGEXP: a galaxy tool suite for transcriptomics data meta-analysis.

Background: With the proliferation of available microarray and high-throughput sequencing experiments in the public domain, the use of meta-analysis methods increases. In these experiments, where the sample size is often limited, meta-analysis offers the possibility to considerably enhance the statistical power and give more accurate results. For those purposes, it combines either effect sizes or results of single studies in an appropriate manner. R packages metaMA and metaRNASeq perform meta-analysis on microarray and next generation sequencing (NGS) data, respectively. They are not interchangeable as they rely on statistical modeling specific to each technology. Results: SMAGEXP (Statistical Meta-Analysis for Gene EXPression) integrates metaMA and metaRNAseq packages into Galaxy. We aim to propose a unified way to carry out meta-analysis of gene expression data, while taking care of their specificities. We have developed this tool suite to analyze microarray data from the Gene Expression Omnibus database or custom data from Affymetrix© microarrays. These data are then combined to carry out meta-analysis using metaMA package. SMAGEXP also offers to combine raw read counts from NGS experiments using DESeq2 and metaRNASeq package. In both cases, key values, independent from the technology type, are reported to judge the quality of the meta-analysis. These tools are available on the Galaxy main tool shed. A dockerized instance of galaxy containing SMAGEXP and its dependencies is available on Docker hub. Source code, help, and installation instructions are available on GitHub. Conclusion: The use of Galaxy offers an easy-to-use gene expression meta-analysis tool suite based on the metaMA and metaRNASeq packages.

Pre-malignant transformation by senescence evasion is prevented by the PERK and ATF6alpha branches of the Unfolded Protein Response.

The incidence of carcinomas highly increases with age. However, the initial steps of the age-related molecular carcinogenic processes remain poorly characterized. We previously showed that normal human epidermal keratinocytes spontaneously and systematically escape from senescence to give rise to preneoplastic emerging cells through a process called post-senescence neoplastic emergence (PSNE). To identify molecular pathways involved in the switch from senescence to pre-transformation, we performed Connectivity Map analyses and DAVID functional annotations followed by hierarchical clustering and multidimensional scaling of the gene expression signature of PSNE cells. We identified endoplasmic reticulum stress related pathways as key regulators of PSNE. Invalidation by RNA interference of the UPR sensors PERK, ATF6α, but not IRE1α, delayed the occurrence of senescence when performed in pre-senescent cells, and increased the PSNE frequency when performed in already senescent cells. Conversely, endoplasmic reticulum stress inducers applied to already senescent cells decreased the frequency of PSNE. In conclusion, these results indicate that the activation of the UPR could protect from the early carcinogenic steps by senescence evasion. This opens new avenues to explore therapeutics that could be useful in decreasing the age-associated tumor incidence.

BACH2 promotes indolent clinical presentation in Waldenström macroglobulinemia.

Approximately 30% of the patients who fulfil the criteria of Waldenström's macroglobulinemia (WM) are diagnosed while asymptomatic (indolent), and will not require immediate therapy. Conversely, patients with a disease-related event will be considered for therapy. The physiopathology of these 2 groups remains unclear, and the mechanisms of progression from indolent to symptomatic WM have yet to be fully understood. Seventeen patients diagnosed with WM were included in this study, 8 asymptomatic WM (A-WM) and 9 symptomatic WM (S-WM). A differential analysis was performed on a first series of 11 patients and identified 48 genes whose expression separated samples from A- to S-WM. This gene signature was then confirmed on a second independent validation set of 6 WM. Within this expression profile, BACH2, a B-cell transcription factor known to be a tumor suppressor gene, was found to be over-expressed in A-MW relatively to S-MW. We specifically over-expressed BACH2 in a WM-related cell line and observed a significant reduction of the clonogenic activity. To the best of our knowledge, we report for the first time a specific gene expression signature that differentiates A-WM and S-WM. Within this expression profile, BACH2 was identified as a candidate gene that may help to understand better the behavior of tumor cells in indolent WM.

The logic of transcriptional regulator recruitment architecture at cis-regulatory modules controlling liver functions.

Control of gene transcription relies on concomitant regulation by multiple transcriptional regulators (TRs). However, how recruitment of a myriad of TRs is orchestrated at cis-regulatory modules (CRMs) to account for coregulation of specific biological pathways is only partially understood. Here, we have used mouse liver CRMs involved in regulatory activities of the hepatic TR, NR1H4 (FXR; farnesoid X receptor), as our model system to tackle this question. Using integrative cistromic, epigenomic, transcriptomic, and interactomic analyses, we reveal a logical organization where trans-regulatory modules (TRMs), which consist of subsets of preferentially and coordinately corecruited TRs, assemble into hierarchical combinations at hepatic CRMs. Different combinations of TRMs add to a core TRM, broadly found across the whole landscape of CRMs, to discriminate promoters from enhancers. These combinations also specify distinct sets of CRM differentially organized along the genome and involved in regulation of either housekeeping/cellular maintenance genes or liver-specific functions. In addition to these TRMs which we define as obligatory, we show that facultative TRMs, such as one comprising core circadian TRs, are further recruited to selective subsets of CRMs to modulate their activities. TRMs transcend TR classification into ubiquitous versus liver-identity factors, as well as TR grouping into functional families. Hence, hierarchical superimpositions of obligatory and facultative TRMs bring about independent transcriptional regulatory inputs defining different sets of CRMs with logical connection to regulation of specific gene sets and biological pathways. Altogether, our study reveals novel principles of concerted transcriptional regulation by multiple TRs at CRMs.