Genomic and epigenomic active vitamin D responses in human colonic organoids

Silibinin, a commonly used therapeutic agent for non-alcohol fatty liver disease, functions through upregulating intestinal expression of fibroblast growth factor 15/19

BACKGROUND AND PURPOSE

Silibinin is used to treat non-alcohol fatty liver disease (NAFLD) despite having rapid liver metabolism. Therefore, we investigated the role of the intestine in silibinin mechanism of action.

EXPERIMENTAL APPROACH

NAFLD mice model was established by feeding them with a high-fat diet (HFD). Liver pathological were examined using H&E and oil red O staining. Tissue distribution of silibinin was detected by LC-MS/MS. SiRNA was employed for gene silencing and plasmid was used for gene overexpression. ChIP-qPCR assay was performed to detect the levels of histone acetylation. Recombinant adeno-associated virus 9-short hairpin-fibroblast growth factor (FGF)-15 and -farnesoid X receptor (FXR; NR1H4) were used to knockdown expression of FGF-15 and FXR.

KEY RESULTS

Oral silibinin significantly reversed NAFLD in mice, although liver concentration was insufficient for reduction of lipid accumulation in hepatocytes. Among endogenous factors capable of reversing NAFLD, the expression of Fgf-15 was selectively up-regulated by silibinin in ileum and colon of mice. When intestinal expression of Fgf-15 was knocked down, protection of silibinin against lipid accumulation and injury of livers nearly disappeared. Silibinin could reduce activity of histone deacetylase 2 (HDAC2), enhance histone acetylation in the promoter region of FXR and consequently increase intestinal expression of FGF-15/19.

CONCLUSION AND IMPLICATIONS

Oral silibinin selectively promotes expression of FGF-15/19 in ileum by enhancing transcription of FXR via reduction of HDAC2 activity, and FGF-15/19 enters into circulation to exert anti-NAFLD action. As the site of action is the intestine this would explain the discrepancy between pharmacodynamics and pharmacokinetics of silibinin.

Translational Utility of Organoid Models for Biomedical Research on Gastrointestinal Diseases

Organoid models have recently been utilized to study 3D human-derived tissue systems to uncover tissue architecture and adult stem cell biology. Patient-derived organoids unambiguously provide the most suitable in vitro system to study disease biology with the actual genetic background. With the advent of much improved and innovative approaches, patient-derived organoids can potentially be used in regenerative medicine. Various human tissues were explored to develop organoids due to their multifold advantage over the conventional in vitro cell line culture approach and in vivo models. Gastrointestinal (GI) tissues have been widely studied to establish organoids and organ-on-chip for screening drugs, nutraceuticals, and other small molecules having therapeutic potential. The function of channel proteins, transporters, and transmembrane proteins was also explained. The successful application of genome editing in organoids using the CRISPR-Cas approach has been reported recently. GI diseases such as Celiac disease (CeD), Inflammatory bowel disease (IBD), and common GI cancers have been investigated using several patient-derived organoid models. Recent advancements on organoid bio-banking and 3D bio-printing contributed significantly in personalized disease management and therapeutics. This article reviews the available literature on investigations and translational applications of patient-derived GI organoid models, notably on elucidating gut-microbial interaction and epigenetic modifications.

Genomic and non-genomic action of vitamin D on ion channels - Targeting mitochondria

Recent studies revealed that mitochondria are not only a place of vitamin D3 metabolism but also direct or indirect targets of its activities. This review summarizes current knowledge on the regulation of ion channels from plasma and mitochondrial membranes by the active form of vitamin D3 (1,25(OH)2D3). 1,25(OH)2D3, is a naturally occurring hormone with pleiotropic activities; implicated in the modulation of cell differentiation, and proliferation and in the prevention of various diseases, including cancer. Many experimental data indicate that 1,25(OH)2D3 deficiency induces ionic remodeling and 1,25(OH)2D3 regulates the activity of multiple ion channels. There are two main theories on how 1,25(OH)2D3 can modify the function of ion channels. First, describes the involvement of genomic pathways of response to 1,25(OH)2D3 in the regulation of the expression of the genes encoding channels, their auxiliary subunits, or additional regulators. Interestingly, intracellular ion channels, like mitochondrial, are encoded by the same genes as plasma membrane channels. Therefore, the comprehensive genomic regulation of the channels from these two different cellular compartments we analyzed using a bioinformatic approach. The second theory explores non-genomic pathways of vitamin D3 activities. It was shown, that 1,25(OH)2D3 indirectly regulates enzymes that impact ion channels, change membrane physical properties, or directly bind to channel proteins. In this article, the involvement of genomic and non-genomic pathways regulated by 1,25(OH)2D3 in the modulation of the levels and activity of plasma membrane and mitochondrial ion channels was investigated by an extensive review of the literature and analysis of the transcriptomic data using bioinformatics.

From molecular basis to clinical insights: a challenging future for the vitamin D endocrine system in colorectal cancer

Colorectal cancer (CRC) is one of the most life-threatening neoplasias in terms of incidence and mortality worldwide. Vitamin D deficiency has been associated with an increased risk of CRC. 1α,25-Dihydroxyvitamin D3 [1,25(OH)2D3], the most active vitamin D metabolite, is a pleiotropic hormone that, through its binding to a transcription factor of the nuclear receptor superfamily, is a major regulator of the human genome. 1,25(OH)2D3 acts on colon carcinoma and stromal cells and displays tumor protective actions. Here, we review the variety of molecular mechanisms underlying the effects of 1,25(OH)2D3 in CRC, which affect multiple processes that are dysregulated during tumor initiation and progression. Additionally, we discuss the epidemiological data that associate vitamin D deficiency and CRC, and the most relevant randomized controlled trials of vitamin D3 supplementation conducted in both healthy individuals and CRC patients.

Changes in chromatin accessibility and transcriptional landscape induced by HDAC inhibitors in TP53 mutated patient-derived colon cancer organoids

Here we present the generation and characterization of patient-derived organoids (PDOs) from colorectal cancer patients. PDOs derived from two patients with TP53 mutations were tested with two different HDAC inhibitors (SAHA and NKL54). Cell death induction, transcriptome, and chromatin accessibility changes were analyzed. HDACIs promote the upregulation of low expressed genes and the downregulation of highly expressed genes. A similar differential effect is observed at the level of chromatin accessibility. Only SAHA is a potent inducer of cell death, which is characterized by the upregulation of BH3-only genes BIK and BMF. Up-regulation of BIK is associated with increased accessibility in an intronic region that has enhancer properties. SAHA, but not NKL54, also causes downregulation of BCL2L1 and decreases chromatin accessibility in three distinct regions of the BCL2L1 locus. Both inhibitors upregulate the expression of innate immunity genes and members of the MHC family. In summary, our exploratory study indicates a mechanism of action for SAHA and demonstrate the low efficacy of NKL54 as a single agent for apoptosis induction, using two PDOs. These observations need to be validated in a larger cohort of PDOs.

Dissection of an impact of VDR and RXRA on the genomic activity of 1,25(OH)2D3 in A431 squamous cell carcinoma

BACKGROUND

Human skin is the natural source, place of metabolism, and target for vitamin D3. The classical active form of vitamin D3, 1,25(OH)2D3, expresses pluripotent properties and is intensively studied in cancer prevention and therapy. To define the specific role of vitamin D3 receptor (VDR) and its co-receptor retinoid X receptor alpha (RXRA) in genomic regulation, VDR or RXRA genes were silenced in the squamous cell carcinoma cell line A431 and treated with 1,25(OH)2D3 at long incubation time points 24 h/72 h. Extending the incubation time of A431 WT (wild-type) cells with 1,25(OH)2D3 resulted in a two-fold increase in DEGs (differentially expressed genes) and a change in the amount of downregulated from 37% to 53%. VDR knockout led to a complete loss of 1,25(OH)2D3-induced genome-wide gene regulation at 24 h time point, but after 72 h, 20 DEGs were found, of which 75% were downregulated, and most of them belonged to the gene ontology group "immune response". This may indicate the existence of an alternative, secondary response to 1,25(OH)2D3. In contrast, treatment of A431 ΔRXRA cells with 1,25(OH)2D3 for 24 h only partially affected DEGs, suggesting RXRA-independent regulation. Interestingly, overexpression of classic 1,25(OH)2D3 targets, like CYP24A1 (family 24 of subfamily A of cytochrome P450 member 1) or CAMP (cathelicidin antimicrobial peptide) was found to be RXRA-independent. Also, immunofluorescence staining of A431 WT cells revealed partial VDR/RXRA colocalization after 24 h and 72 h 1,25(OH)2D3 treatment. Comparison of transcriptome changes induced by 1,25(OH)2D3 in normal keratinocytes vs. cancer cells showed high cell type specific expression pattern with only a few genes commonly regulated by 1,25(OH)2D3. Activation of the genomic pathway at least partially reversed the expression of cancer-related genes, forming a basis for anti-cancer activates of 1,25(OH)2D3. In summary, VDR or RXRA independent genomic activities of 1,25(OH)2D3 suggest the involvement of alternative factors, opening new challenges in this field.

Transcriptome-wide analysis of mRNA N6 -methyladenosine modification in the embryonic development of Spodoptera frugiperda

N6 -methyladenosine (m6 A) RNA is the most abundant modification of mRNA, and has been demonstrated in regulating various post-transcriptional processes. Many studies have shown that m6 A methylation plays key roles in sex determination, neuronal functions, and embryonic development in Drosophila and mammals. Here, we analyzed transcriptome-wide profile of m6 A modification in the embryonic development of the destructive agricultural pest Spodoptera frugiperda. We found that the 2 key mRNA m6 A methyltransferases SfrMETTL3 and SfrMETTL14 have high homologies with other insects and mammals, suggesting that SfrMETTL3 and SfrMETTL14 may have conserved function among different species. From methylated RNA immunoprecipitation sequencing analysis, we obtained 46 869 m6 A peaks representing 8 587 transcripts in the 2-h embryos after oviposition, and 41 389 m6 A peaks representing 9 230 transcripts in the 24-h embryos. In addition, 5 995 m6 A peaks were differentially expressed including 3 752 upregulated and 2243 downregulated peaks. Functional analysis with Gene Ontology and Kyoto Encyclopedia of Genes and Genomes suggested that differentially expressed m6 A peak-modified genes were enriched in cell and organ development between the 2- and 24-h embryos. By conjoint analysis of methylated RNA immunoprecipitation-seq and RNA-seq data, we found that RNA m6 A methylation may regulate the transcriptional levels of genes related to tissue and organ development from 2- to 24-h embryos. Our study reveals the role of RNA m6 A epigenetic regulation in the embryonic development of S. frugiperda, and provides new insights for the embryonic development of insects.

Genomically anchored vitamin D receptor mediates an abundance of bioprotective actions elicited by its 1,25-dihydroxyvitamin D hormonal ligand

The nuclear vitamin D receptor (VDR) mediates the actions of its physiologic 1,25-dihydroxyvitamin D3 (1,25D) ligand produced in kidney and at extrarenal sites during times of physiologic and cellular stress. The ligand-receptor complex transcriptionally controls genes encoding factors that regulate calcium and phosphate sensing/transport, bone remodeling, immune function, and nervous system maintenance. With the aid of parathyroid hormone (PTH) and fibroblast growth factor 23 (FGF23), 1,25D/VDR primarily participates in an intricate network of feedback controls that govern extracellular calcium and phosphate concentrations, mainly influencing bone formation and mineralization, ectopic calcification, and indirectly supporting many fundamental roles of calcium. Beyond endocrine and intracrine effects, 1,25D/VDR signaling impacts multiple biochemical phenomena that potentially affect human health and disease, including autophagy, carcinogenesis, cell growth/differentiation, detoxification, metabolic homeostasis, and oxidative stress mitigation. Several health advantages conferred by 1,25D/VDR appear to be promulgated by induction of klotho, an anti-aging renal peptide hormone which functions as a co-receptor for FGF23 and, like 1,25D, regulates nrf2, foxo, mTOR and other cellular protective pathways. Among hundreds of genes for which expression is modulated by 1,25D/VDR either primarily or secondarily in a cell-specific manner, the resulting gene products (in addition to those expressed in the classic skeletal mineral regulatory tissues kidney, intestine, and bone), fall into multiple biochemical categories including apoptosis, cholesterol homeostasis, glycolysis, hypoxia, inflammation, p53 signaling, unfolded protein response and xenobiotic metabolism. Thus, 1,25D/VDR is a bone mineral control instrument that also signals the maintenance of multiple cellular processes in the face of environmental and genetic challenges.

Gene expression in organoids: an expanding horizon

Recent development of human three-dimensional organoid cultures has opened new doors and opportunities ranging from modelling human development in vitro to personalised cancer therapies. These new in vitro systems are opening new horizons to the classic understanding of human development and disease. However, the complexity and heterogeneity of these models requires cutting-edge techniques to capture and trace global changes in gene expression to enable identification of key players and uncover the underlying molecular mechanisms. Rapid development of sequencing approaches made possible global transcriptome analyses and epigenetic profiling. Despite challenges in organoid culture and handling, these techniques are now being adapted to embrace organoids derived from a wide range of human tissues. Here, we review current state-of-the-art multi-omics technologies, such as single-cell transcriptomics and chromatin accessibility assays, employed to study organoids as a model for development and a platform for precision medicine.

Genomic and epigenomic responses to aspirin in human colonic organoids

Aspirin (ASA) is a proven chemoprotective agent for colorectal cancer, though mechanisms underlying these effects are incompletely understood. Human organoids are an ideal system to study genomic and epigenomic host-environment interactions. We use human colonic organoids to profile ASA responses on genome-wide gene expression and chromatin accessibility. Human colonic organoids from one individual were cultured and treated in triplicate with 3 mM ASA or vehicle control (DMSO) for 24 h. Gene expression and chromatin accessibility were measured using RNA- and ATAC-sequencing, respectively. Differentially expressed genes were analyzed using DESeq2. Top genes were validated by qPCR. Gene set enrichment was performed by SetRank. Differentially accessible peaks were analyzed using DiffBind and edgeR. Peak annotation and differential transcription factor motifs were determined by HOMER and diffTF. The results showed robust transcriptional responses to ASA with significant enrichment for fatty acid oxidation and peroxisome proliferator-activated receptor (PPAR) signaling that were validated in independent organoid lines. A large number of differentially accessible chromatin regions were found in response to ASA with significant enrichment for Fos, Jun, and Hnf transcription factor motifs. Integrated analysis of epigenomic and genomic treatment responses highlighted gene regions that could mediate ASA's specific effects in the colon including those involved in chemoprotection and/or toxicity. Assessment of chromatin accessibility and transcriptional responses to ASA yielded new observations about genome-wide effects in the colon facilitated by application of human colonic organoids. This framework can be applied to study colonic ASA responses between individuals and populations in future studies.

Patient-Derived Organoid Model in the Prediction of Chemotherapeutic Drug Response in Colorectal Cancer

As an emerging technology in precision medicine, the patient-derived organoid (PDO) technology has been indicated to provide novel modalities to judge the sensitivity of individual tumors to cancer drugs. In this work, an in vitro model of colorectal cancer (CRC) was established using the PDO culture, and it is demonstrated that the PDO samples preserved, to a great extent, the histologic features and marker expression of the original tumor tissues. Subsequently, cancer drugs 5-FU, oxaliplatin, and irinotecan were selected and screened on five CRC PDO samples, while the patient-derived organoid xenograft (PDOX) model was applied for comparison. The receiver operating characteristic (ROC) curve was drawn according to the IC₅₀ data from the PDO model and the relative tumor proliferation rate (T/C%) from PDOX. Interestingly, the area under the ROC curve was 0.84 (95% CI, 0.64-1.04, P value = 0.028), which suggested that the IC₅₀ of cancer drugs from the PDO model was strongly correlated with PDOX responses. In addition, the optimal sensitivity cutoff value for drug screening in CRC PDOs was identified at 10.35 μM, which could act as a reference value for efficacy evaluation of 5-FU, oxaliplatin, and irinotecan in the colorectal cancer drug screening. Since there are no unified criteria to judge the sensitivity of drugs in vitro, our work provides a method for establishing in vitro evaluation criteria via PDO and PDOX model using the patient tissues received from local hospitals, exhibiting potential in clinical cancer therapy and precision medicine.

Oct4 dependent chromatin activation is required for chicken primordial germ cell migration

Primordial germ cells (PGCs) are the undifferentiated progenitors of the gametes. Unlike the poor maintenance of cultured mammalian PGCs, the avian PGCs can be expanded in vitro indefinitely while preserving pluripotency and germline competence. In mammals, the Oct4 is the master transcription factor that ensures the stemness of pluripotent cells such as PGCs, but the specific function of Oct4 in chicken PGCs remains unclear. As expected, the loss of Oct4 in chicken PGCs reduced the expression of key pluripotency factors and promoted the genes involved in endoderm and ectoderm differentiation. Furthermore, the global active chromatin was reduced as shown by the depletion of the H3K27ac upon Oct4 suppression. Interestingly, the de-activated chromatin caused the down-regulation of adjacent genes which are mostly known regulators of cell junction, chemotaxis and cell migration. Consequently, the Oct4-deficient PGCs show impaired cell migration and could not colonize the gonads when re-introduced into the bloodstream of the embryo. We propose that, in addition to maintaining pluripotency, the Oct4 mediated chromatin activation is dictating chicken PGC migration.