Post-GWAS multiomic functional investigation of the TNIP1 locus in Alzheimer's disease highlights a potential role for GPX3.

INTRODUCTION: Recent genome-wide association studies (GWAS) have reported a genetic association with Alzheimer's disease (AD) at the TNIP1/GPX3 locus, but the mechanism is unclear. METHODS: We used cerebrospinal fluid (CSF) proteomics data to test (n = 137) and replicate (n = 446) the association of glutathione peroxidase 3 (GPX3) with CSF biomarkers (including amyloid and tau) and the GWAS-implicated variants (rs34294852 and rs871269). RESULTS: CSF GPX3 levels decreased with amyloid and tau positivity (analysis of variance P = 1.5 × 10-5) and higher CSF phosphorylated tau (p-tau) levels (P = 9.28 × 10-7). The rs34294852 minor allele was associated with decreased GPX3 (P = 0.041). The replication cohort found associations of GPX3 with amyloid and tau positivity (P = 2.56 × 10-6) and CSF p-tau levels (P = 4.38 × 10-9). DISCUSSION: These results suggest variants in the TNIP1 locus may affect the oxidative stress response in AD via altered GPX3 levels. HIGHLIGHTS: Cerebrospinal fluid (CSF) glutathione peroxidase 3 (GPX3) levels decreased with amyloid and tau positivity and higher CSF phosphorylated tau. The minor allele of rs34294852 was associated with lower CSF GPX3. levels when also controlling for amyloid and tau category. GPX3 transcript levels in the prefrontal cortex were lower in Alzheimer's disease than controls. rs34294852 is an expression quantitative trait locus for GPX3 in blood, neutrophils, and microglia.

Evolution of diapause in the African turquoise killifish by remodeling the ancient gene regulatory landscape.

Suspended animation states allow organisms to survive extreme environments. The African turquoise killifish has evolved diapause as a form of suspended development to survive a complete drought. However, the mechanisms underlying the evolution of extreme survival states are unknown. To understand diapause evolution, we performed integrative multi-omics (gene expression, chromatin accessibility, and lipidomics) in the embryos of multiple killifish species. We find that diapause evolved by a recent remodeling of regulatory elements at very ancient gene duplicates (paralogs) present in all vertebrates. CRISPR-Cas9-based perturbations identify the transcription factors REST/NRSF and FOXOs as critical for the diapause gene expression program, including genes involved in lipid metabolism. Indeed, diapause shows a distinct lipid profile, with an increase in triglycerides with very-long-chain fatty acids. Our work suggests a mechanism for the evolution of complex adaptations and offers strategies to promote long-term survival by activating suspended animation programs in other species.

Deconvolution of polygenic risk score in single cells unravels cellular and molecular heterogeneity of complex human diseases.

Polygenic risk scores (PRSs) are commonly used for predicting an individual's genetic risk of complex diseases. Yet, their implication for disease pathogenesis remains largely limited. Here, we introduce scPRS, a geometric deep learning model that constructs single-cell-resolved PRS leveraging reference single-cell chromatin accessibility profiling data to enhance biological discovery as well as disease prediction. Real-world applications across multiple complex diseases, including type 2 diabetes (T2D), hypertrophic cardiomyopathy (HCM), and Alzheimer's disease (AD), showcase the superior prediction power of scPRS compared to traditional PRS methods. Importantly, scPRS not only predicts disease risk but also uncovers disease-relevant cells, such as hormone-high alpha and beta cells for T2D, cardiomyocytes and pericytes for HCM, and astrocytes, microglia and oligodendrocyte progenitor cells for AD. Facilitated by a layered multi-omic analysis, scPRS further identifies cell-type-specific genetic underpinnings, linking disease-associated genetic variants to gene regulation within corresponding cell types. We substantiate the disease relevance of scPRS-prioritized HCM genes and demonstrate that the suppression of these genes in HCM cardiomyocytes is rescued by Mavacamten treatment. Additionally, we establish a novel microglia-specific regulatory relationship between the AD risk variant rs7922621 and its target genes ANXA11 and TSPAN14. We further illustrate the detrimental effects of suppressing these two genes on microglia phagocytosis. Our work provides a multi-tasking, interpretable framework for precise disease prediction and systematic investigation of the genetic, cellular, and molecular basis of complex diseases, laying the methodological foundation for single-cell genetics.

Personalized transcriptome signatures in a cardiomyopathy stem cell biobank.

BACKGROUND: There is growing evidence that pathogenic mutations do not fully explain hypertrophic (HCM) or dilated (DCM) cardiomyopathy phenotypes. We hypothesized that if a patient's genetic background was influencing cardiomyopathy this should be detectable as signatures in gene expression. We built a cardiomyopathy biobank resource for interrogating personalized genotype phenotype relationships in human cell lines. METHODS: We recruited 308 diseased and control patients for our cardiomyopathy stem cell biobank. We successfully reprogrammed PBMCs (peripheral blood mononuclear cells) into induced pluripotent stem cells (iPSCs) for 300 donors. These iPSCs underwent whole genome sequencing and were differentiated into cardiomyocytes for RNA-seq. In addition to annotating pathogenic variants, mutation burden in a panel of cardiomyopathy genes was assessed for correlation with echocardiogram measurements. Line-specific co-expression networks were inferred to evaluate transcriptomic subtypes. Drug treatment targeted the sarcomere, either by activation with omecamtiv mecarbil or inhibition with mavacamten, to alter contractility. RESULTS: We generated an iPSC biobank from 300 donors, which included 101 individuals with HCM and 88 with DCM. Whole genome sequencing of 299 iPSC lines identified 78 unique pathogenic or likely pathogenic mutations in the diseased lines. Notably, only DCM lines lacking a known pathogenic or likely pathogenic mutation replicated a finding in the literature for greater nonsynonymous SNV mutation burden in 102 cardiomyopathy genes to correlate with lower left ventricular ejection fraction in DCM. We analyzed RNA-sequencing data from iPSC-derived cardiomyocytes for 102 donors. Inferred personalized co-expression networks revealed two transcriptional subtypes of HCM. The first subtype exhibited concerted activation of the co-expression network, with the degree of activation reflective of the disease severity of the donor. In contrast, the second HCM subtype and the entire DCM cohort exhibited partial activation of the respective disease network, with the strength of specific gene by gene relationships dependent on the iPSC-derived cardiomyocyte line. ADCY5 was the largest hubnode in both the HCM and DCM networks and partially corrected in response to drug treatment. CONCLUSIONS: We have a established a stem cell biobank for studying cardiomyopathy. Our analysis supports the hypothesis the genetic background influences pathologic gene expression programs and support a role for ADCY5 in cardiomyopathy.

Molecular adaptations in response to exercise training are associated with tissue-specific transcriptomic and epigenomic signatures.

Regular exercise has many physical and brain health benefits, yet the molecular mechanisms mediating exercise effects across tissues remain poorly understood. Here we analyzed 400 high-quality DNA methylation, ATAC-seq, and RNA-seq datasets from eight tissues from control and endurance exercise-trained (EET) rats. Integration of baseline datasets mapped the gene location dependence of epigenetic control features and identified differing regulatory landscapes in each tissue. The transcriptional responses to 8 weeks of EET showed little overlap across tissues and predominantly comprised tissue-type enriched genes. We identified sex differences in the transcriptomic and epigenomic changes induced by EET. However, the sex-biased gene responses were linked to shared signaling pathways. We found that many G protein-coupled receptor-encoding genes are regulated by EET, suggesting a role for these receptors in mediating the molecular adaptations to training across tissues. Our findings provide new insights into the mechanisms underlying EET-induced health benefits across organs.

Genome-wide Cas9-mediated screening of essential non-coding regulatory elements via libraries of paired single-guide RNAs.

The functions of non-coding regulatory elements (NCREs), which constitute a major fraction of the human genome, have not been systematically studied. Here we report a method involving libraries of paired single-guide RNAs targeting both ends of an NCRE as a screening system for the Cas9-mediated deletion of thousands of NCREs genome-wide to study their functions in distinct biological contexts. By using K562 and 293T cell lines and human embryonic stem cells, we show that NCREs can have redundant functions, and that many ultra-conserved elements have silencer activity and play essential roles in cell growth and in cellular responses to drugs (notably, the ultra-conserved element PAX6_Tarzan may be critical for heart development, as removing it from human embryonic stem cells led to defects in cardiomyocyte differentiation). The high-throughput screen, which is compatible with single-cell sequencing, may allow for the identification of druggable NCREs.

Deep learning modeling of rare noncoding genetic variants in human motor neurons defines CCDC146 as a therapeutic target for ALS.

Amyotrophic lateral sclerosis (ALS) is a fatal and incurable neurodegenerative disease caused by the selective and progressive death of motor neurons (MNs). Understanding the genetic and molecular factors influencing ALS survival is crucial for disease management and therapeutics. In this study, we introduce a deep learning-powered genetic analysis framework to link rare noncoding genetic variants to ALS survival. Using data from human induced pluripotent stem cell (iPSC)-derived MNs, this method prioritizes functional noncoding variants using deep learning, links cis-regulatory elements (CREs) to target genes using epigenomics data, and integrates these data through gene-level burden tests to identify survival-modifying variants, CREs, and genes. We apply this approach to analyze 6,715 ALS genomes, and pinpoint four novel rare noncoding variants associated with survival, including chr7:76,009,472:C>T linked to CCDC146. CRISPR-Cas9 editing of this variant increases CCDC146 expression in iPSC-derived MNs and exacerbates ALS-specific phenotypes, including TDP-43 mislocalization. Suppressing CCDC146 with an antisense oligonucleotide (ASO), showing no toxicity, completely rescues ALS-associated survival defects in neurons derived from sporadic ALS patients and from carriers of the ALS-associated G4C2-repeat expansion within C9ORF72. ASO targeting of CCDC146 may be a broadly effective therapeutic approach for ALS. Our framework provides a generic and powerful approach for studying noncoding genetics of complex human diseases.

Longitudinal cytokine and multi-modal health data of an extremely severe ME/CFS patient with HSD reveals insights into immunopathology, and disease severity.

Introduction: Myalgic Encephalomyelitis/Chronic Fatigue Syndrome (ME/CFS) presents substantial challenges in patient care due to its intricate multisystem nature, comorbidities, and global prevalence. The heterogeneity among patient populations, coupled with the absence of FDA-approved diagnostics and therapeutics, further complicates research into disease etiology and patient managment. Integrating longitudinal multi-omics data with clinical, health,textual, pharmaceutical, and nutraceutical data offers a promising avenue to address these complexities, aiding in the identification of underlying causes and providing insights into effective therapeutics and diagnostic strategies. Methods: This study focused on an exceptionally severe ME/CFS patient with hypermobility spectrum disorder (HSD) during a period of marginal symptom improvements. Longitudinal cytokine profiling was conducted alongside the collection of extensive multi-modal health data to explore the dynamic nature of symptoms, severity, triggers, and modifying factors. Additionally, an updated severity assessment platform and two applications, ME-CFSTrackerApp and LexiTime, were introduced to facilitate real-time symptom tracking and enhance patient-physician/researcher communication, and evaluate response to medical intervention. Results: Longitudinal cytokine profiling revealed the significance of Th2-type cytokines and highlighted synergistic activities between mast cells and eosinophils, skewing Th1 toward Th2 immune responses in ME/CFS pathogenesis, particularly in cognitive impairment and sensorial intolerance. This suggests a potentially shared underlying mechanism with major ME/CFS comorbidities such as HSD, Mast cell activation syndrome, postural orthostatic tachycardia syndrome (POTS), and small fiber neuropathy. Additionally, the data identified potential roles of BCL6 and TP53 pathways in ME/CFS etiology and emphasized the importance of investigating adverse reactions to medication and supplements and drug interactions in ME/CFS severity and progression. Discussion: Our study advocates for the integration of longitudinal multi-omics with multi-modal health data and artificial intelligence (AI) techniques to better understand ME/CFS and its major comorbidities. These findings highlight the significance of dysregulated Th2-type cytokines in patient stratification and precision medicine strategies. Additionally, our results suggest exploring the use of low-dose drugs with partial agonist activity as a potential avenue for ME/CFS treatment. This comprehensive approach emphasizes the importance of adopting a patient-centered care approach to improve ME/CFS healthcare management, disease severity assessment, and personalized medicine. Overall, these findings contribute to our understanding of ME/CFS and offer avenues for future research and clinical practice.

Multi-omics in stress and health research: study designs that will drive the field forward.

Despite decades of stress research, there still exist substantial gaps in our understanding of how social, environmental, and biological factors interact and combine with developmental stressor exposures, cognitive appraisals of stressors, and psychosocial coping processes to shape individuals' stress reactivity, health, and disease risk. Relatively new biological profiling approaches, called multi-omics, are helping address these issues by enabling researchers to quantify thousands of molecules from a single blood or tissue sample, thus providing a panoramic snapshot of the molecular processes occurring in an organism from a systems perspective. In this review, we summarize two types of research designs for which multi-omics approaches are best suited, and describe how these approaches can help advance our understanding of stress processes and the development, prevention, and treatment of stress-related pathologies. We first discuss incorporating multi-omics approaches into theory-rich, intensive longitudinal study designs to characterize, in high-resolution, the transition to stress-related multisystem dysfunction and disease throughout development. Next, we discuss how multi-omics approaches should be incorporated into intervention research to better understand the transition from stress-related dysfunction back to health, which can help inform novel precision medicine approaches to managing stress and fostering biopsychosocial resilience. Throughout, we provide concrete recommendations for types of studies that will help advance stress research, and translate multi-omics data into better health and health care.

Longitudinal profiling of the microbiome at four body sites reveals core stability and individualized dynamics during health and disease.

To understand the dynamic interplay between the human microbiome and host during health and disease, we analyzed the microbial composition, temporal dynamics, and associations with host multi-omics, immune, and clinical markers of microbiomes from four body sites in 86 participants over 6 years. We found that microbiome stability and individuality are body-site specific and heavily influenced by the host. The stool and oral microbiome are more stable than the skin and nasal microbiomes, possibly due to their interaction with the host and environment. We identify individual-specific and commonly shared bacterial taxa, with individualized taxa showing greater stability. Interestingly, microbiome dynamics correlate across body sites, suggesting systemic dynamics influenced by host-microbial-environment interactions. Notably, insulin-resistant individuals show altered microbial stability and associations among microbiome, molecular markers, and clinical features, suggesting their disrupted interaction in metabolic disease. Our study offers comprehensive views of multi-site microbial dynamics and their relationship with host health and disease.

Immunotherapeutic IL-6R and targeting the MCT-1/IL-6/CXCL7/PD-L1 circuit prevent relapse and metastasis of triple-negative breast cancer.

Rationale: Multiple copies in T-cell malignancy 1 (MCT-1) is a prognostic biomarker for aggressive breast cancers. Overexpressed MCT-1 stimulates the IL-6/IL-6R/gp130/STAT3 axis, which promotes epithelial-to-mesenchymal transition and cancer stemness. Because cancer stemness largely contributes to the tumor metastasis and recurrence, we aimed to identify whether the blockade of MCT-1 and IL-6R can render these effects and to understand the underlying mechanisms that govern the process. Methods: We assessed primary tumor invasion, postsurgical local recurrence and distant metastasis in orthotopic syngeneic mice given the indicated immunotherapy and MCT-1 silencing (shMCT-1). Results: We found that shMCT-1 suppresses the transcriptomes of the inflammatory response and metastatic signaling in TNBC cells and inhibits tumor recurrence, metastasis and mortality in xenograft mice. IL-6R immunotherapy and shMCT-1 combined further decreased intratumoral M2 macrophages and T regulatory cells (Tregs) and avoided postsurgical TNBC expansion. shMCT-1 also enhances IL-6R-based immunotherapy effectively in preventing postsurgical TNBC metastasis, recurrence and mortality. Anti-IL-6R improved helper T, cytotoxic T and natural killer (NK) cells in the lymphatic system and decreased Tregs in the recurrent and metastatic tumors. Combined IL-6R and PD-L1 immunotherapies abridged TNBC cell stemness and M2 macrophage activity to a greater extent than monotherapy. Sequential immunotherapy of PD-L1 and IL-6R demonstrated the best survival outcome and lowest postoperative recurrence and metastasis compared with synchronized therapy, particularly in the shMCT-1 context. Multiple positive feedforward loops of the MCT-1/IL-6/IL-6R/CXCL7/PD-L1 axis were identified in TNBC cells, which boosted metastatic niches and immunosuppressive microenvironments. Clinically, MCT-1high/PD-L1high/CXCL7high and CXCL7high/IL-6high/IL-6Rhigh expression patterns predict worse prognosis and poorer survival of breast cancer patients. Conclusion: Systemic targeting the MCT-1/IL-6/IL-6R/CXCL7/PD-L1 interconnections enhances immune surveillance that inhibits the aggressiveness of TNBC.

Miscarriage risk assessment: a bioinformatic approach to identifying candidate lethal genes and variants.

PURPOSE: Miscarriage, often resulting from a variety of genetic factors, is a common pregnancy outcome. Preconception genetic carrier screening (PGCS) identifies at-risk partners for newborn genetic disorders; however, PGCS panels currently lack miscarriage-related genes. In this study, we evaluated the potential impact of both known and candidate genes on prenatal lethality and the effectiveness of PGCS in diverse populations. METHODS: We analyzed 125,748 human exome sequences and mouse and human gene function databases. Our goals were to identify genes crucial for human fetal survival (lethal genes), to find variants not present in a homozygous state in healthy humans, and to estimate carrier rates of known and candidate lethal genes in various populations and ethnic groups. RESULTS: This study identified 138 genes in which heterozygous lethal variants are present in the general population with a frequency of 0.5% or greater. Screening for these 138 genes could identify 4.6% (in the Finnish population) to 39.8% (in the East Asian population) of couples at risk of miscarriage. This explains the cause of pregnancy loss in approximately 1.1-10% of cases affected by biallelic lethal variants. CONCLUSION: This study has identified a set of genes and variants potentially associated with lethality across different ethnic backgrounds. The variation of these genes across ethnic groups underscores the need for a comprehensive, pan-ethnic PGCS panel that includes genes related to miscarriage.

Harnessing human genetics and stem cells for precision cardiovascular medicine.

Human induced pluripotent stem cell (iPSC) platforms are valuable for biomedical and pharmaceutical research by providing tissue-specific human cells that retain patients' genetic integrity and display disease phenotypes in a dish. Looking forward, combining iPSC phenotyping platforms with genomic and screening technologies will continue to pave new directions for precision medicine, including genetic prediction, visualization, and treatment of heart disease. This review summarizes the recent use of iPSC technology to unpack the influence of genetic variants in cardiovascular pathology. We focus on various state-of-the-art genomic tools for cardiovascular therapies-including the expansion of genetic toolkits for molecular interrogation, in vitro population studies, and function-based drug screening-and their current applications in patient- and genome-edited iPSC platforms that are heralding new avenues for cardiovascular research.

Short-chain fatty acids propionate and butyrate control growth and differentiation linked to cellular metabolism.

The short-chain fatty acids (SCFA) propionate and butyrate are produced in large amounts by microbial metabolism and have been identified as unique acyl lysine histone marks. In order to better understand the function of these modifications we used ChIP-seq to map the genome-wide location of four short-chain acyl histone marks H3K18pr/bu and H4K12pr/bu in treated and untreated colorectal cancer (CRC) and normal cells, as well as in mouse intestines in vivo. We correlate these marks with open chromatin regions along with gene expression to access the function of the target regions. Our data demonstrate that propionate and butyrate act as promoters of growth, differentiation as well as ion transport. We propose a mechanism involving direct modification of specific genomic regions, resulting in increased chromatin accessibility, and in case of butyrate, opposing effects on the proliferation of normal versus CRC cells.

Validation of biomarkers of aging.

The search for biomarkers that quantify biological aging (particularly 'omic'-based biomarkers) has intensified in recent years. Such biomarkers could predict aging-related outcomes and could serve as surrogate endpoints for the evaluation of interventions promoting healthy aging and longevity. However, no consensus exists on how biomarkers of aging should be validated before their translation to the clinic. Here, we review current efforts to evaluate the predictive validity of omic biomarkers of aging in population studies, discuss challenges in comparability and generalizability and provide recommendations to facilitate future validation of biomarkers of aging. Finally, we discuss how systematic validation can accelerate clinical translation of biomarkers of aging and their use in gerotherapeutic clinical trials.

The Importance, Challenges, and Possible Solutions for Sharing Proteomics Data While Safeguarding Individuals' Privacy.

Proteomics data sharing has profound benefits at the individual level as well as at the community level. While data sharing has increased over the years, mostly due to journal and funding agency requirements, the reluctance of researchers with regard to data sharing is evident as many shares only the bare minimum dataset required to publish an article. In many cases, proper metadata is missing, essentially making the dataset useless. This behavior can be explained by a lack of incentives, insufficient awareness, or a lack of clarity surrounding ethical issues. Through adequate training at research institutes, researchers can realize the benefits associated with data sharing and can accelerate the norm of data sharing for the field of proteomics, as has been the standard in genomics for decades. In this article, we have put together various repository options available for proteomics data. We have also added pros and cons of those repositories to facilitate researchers in selecting the repository most suitable for their data submission. It is also important to note that a few types of proteomics data have the potential to re-identify an individual in certain scenarios. In such cases, extra caution should be taken to remove any personal identifiers before sharing on public repositories. Data sets that will be useless without personal identifiers need to be shared in a controlled access repository so that only authorized researchers can access the data and personal identifiers are kept safe.

Longitudinal profiling of the microbiome at four body sites reveals core stability and individualized dynamics during health and disease.

To understand dynamic interplay between the human microbiome and host during health and disease, we analyzed the microbial composition, temporal dynamics, and associations with host multi-omics, immune and clinical markers of microbiomes from four body sites in 86 participants over six years. We found that microbiome stability and individuality are body-site-specific and heavily influenced by the host. The stool and oral microbiome were more stable than the skin and nasal microbiomes, possibly due to their interaction with the host and environment. Also, we identified individual-specific and commonly shared bacterial taxa, with individualized taxa showing greater stability. Interestingly, microbiome dynamics correlated across body sites, suggesting systemic coordination influenced by host-microbial-environment interactions. Notably, insulin-resistant individuals showed altered microbial stability and associations between microbiome, molecular markers, and clinical features, suggesting their disrupted interaction in metabolic disease. Our study offers comprehensive views of multi-site microbial dynamics and their relationship with host health and disease. Study Highlights: The stability of the human microbiome varies among individuals and body sites.Highly individualized microbial genera are more stable over time.At each of the four body sites, systematic interactions between the environment, the host and bacteria can be detected.Individuals with insulin resistance have lower microbiome stability, a more diversified skin microbiome, and significantly altered host-microbiome interactions.

Rare and common genetic determinants of mitochondrial function determine severity but not risk of amyotrophic lateral sclerosis.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease involving selective vulnerability of energy-intensive motor neurons (MNs). It has been unclear whether mitochondrial function is an upstream driver or a downstream modifier of neurotoxicity. We separated upstream genetic determinants of mitochondrial function, including genetic variation within the mitochondrial genome or autosomes; from downstream changeable factors including mitochondrial DNA copy number (mtCN). Across three cohorts including 6,437 ALS patients, we discovered that a set of mitochondrial haplotypes, chosen because they are linked to measurements of mitochondrial function, are a determinant of ALS survival following disease onset, but do not modify ALS risk. One particular haplotype appeared to be neuroprotective and was significantly over-represented in two cohorts of long-surviving ALS patients. Causal inference for mitochondrial function was achievable using mitochondrial haplotypes, but not autosomal SNPs in traditional Mendelian randomization (MR). Furthermore, rare loss-of-function genetic variants within, and reduced MN expression of, ACADM and DNA2 lead to ∼50 % shorter ALS survival; both proteins are implicated in mitochondrial function. Both mtCN and cellular vulnerability are linked to DNA2 function in ALS patient-derived neurons. Finally, MtCN responds dynamically to the onset of ALS independently of mitochondrial haplotype, and is correlated with disease severity. We conclude that, based on the genetic measures we have employed, mitochondrial function is a therapeutic target for amelioration of disease severity but not prevention of ALS.

Short-chain fatty acids propionate and butyrate control growth and differentiation linked to cellular metabolism.

The short-chain fatty acids (SCFA) propionate and butyrate are produced in large amounts by microbial metabolism and have been identified as unique acyl lysine histone marks. In order to better understand the function of these modifications we used ChIP-seq to map the genome-wide location of four short-chain acyl histone marks H3K18pr/bu and H4K12pr/bu in treated and untreated colorectal cancer (CRC) and normal cells, as well as in mouse intestines in vivo . We correlate these marks with open chromatin regions along with gene expression to access the function of the target regions. Our data demonstrate that propionate and butyrate act as promoters of growth, differentiation as well as ion transport. We propose a mechanism involving direct modification of specific genomic regions, resulting in increased chromatin accessibility, and in case of butyrate, opposing effects on the proliferation of normal versus CRC cells.