Obesity-associated variants within FTO form long-range functional connections with IRX3

Genome-Wide Association Study of Treatment-Resistant Depression: Shared Biology With Metabolic Traits

OBJECTIVE

Treatment-resistant depression (TRD) occurs in roughly one-third of all individuals with major depressive disorder (MDD). Although research has suggested a significant common variant genetic component of liability to TRD, with heritability estimated at 8% when compared with non-treatment-resistant MDD, no replicated genetic loci have been identified, and the genetic architecture of TRD remains unclear. A key barrier to this work has been the paucity of adequately powered cohorts for investigation, largely because of the challenge in prospectively investigating this phenotype. The objective of this study was to perform a well-powered genetic study of TRD.

METHODS

Using receipt of electroconvulsive therapy (ECT) as a surrogate for TRD, the authors applied standard machine learning methods to electronic health record data to derive predicted probabilities of receiving ECT. These probabilities were then applied as a quantitative trait in a genome-wide association study of 154,433 genotyped patients across four large biobanks.

RESULTS

Heritability estimates ranged from 2% to 4.2%, and significant genetic overlap was observed with cognition, attention deficit hyperactivity disorder, schizophrenia, alcohol and smoking traits, and body mass index. Two genome-wide significant loci were identified, both previously implicated in metabolic traits, suggesting shared biology and potential pharmacological implications.

CONCLUSIONS

This work provides support for the utility of estimation of disease probability for genomic investigation and provides insights into the genetic architecture and biology of TRD.

Single-nucleus RNA transcriptome profiling reveals murine adipose tissue endothelial cell proliferation gene networks involved in obesity development

Endothelial cells play an important role in the metabolism of adipose tissue (AT). This study aimed to analyze the changes that adipose tissue in AT endothelial cells undergo during the development of obesity, using single-nucleus RNA sequence (snRNA-seq). Mouse paraepididymal AT cells were subjected to snRNA-seq with the 10X Genomics platform. The cell types were then clustered using t-distributed stochastic neighbor embedding and unbiased computational informatics analyses. Protein-protein interactions network was established using the STRING database and visualized using Cytoscape. The dataset was subjected to differential gene enrichment analysis. In total, 21,333 cells acquired from 24 mouse paraepididymal AT samples were analyzed using snRNA-seq. This study identified 18 distinct clusters and annotated macrophages, fibroblasts, epithelial cells, T cells, endothelial cells, stem cells, neutrophil cells, and neutrophil cell types based on representative markers. Cluster 12 was defined as endothelial cells. The proportion of endothelial cells decreased with the development of obesity. Inflammatory factors, such as Vegfa and Prdm16 were upregulated in the medium obesity group but downregulated in the obesity group. Genes, such as Prox1, Erg, Flt4, Kdr, Flt1, and Pecam1 promoted the proliferation of AT endothelial cells and maintained the internal environment of AT. This study established a reference model and general framework for studying the mechanisms, biomarkers, and therapeutic targets of endothelial cell dysfunction-related diseases at the single-cell level.

Variant-to-function analysis of the childhood obesity chr12q13 locus implicates rs7132908 as a causal variant within the 3' UTR of FAIM2

The ch12q13 locus is among the most significant childhood obesity loci identified in genome-wide association studies. This locus resides in a non-coding region within FAIM2; thus, the underlying causal variant(s) presumably influence disease susceptibility via cis-regulation. We implicated rs7132908 as a putative causal variant by leveraging our in-house 3D genomic data and public domain datasets. Using a luciferase reporter assay, we observed allele-specific cis-regulatory activity of the immediate region harboring rs7132908. We generated isogenic human embryonic stem cell lines homozygous for either rs7132908 allele to assess changes in gene expression and chromatin accessibility throughout a differentiation to hypothalamic neurons, a key cell type known to regulate feeding behavior. The rs7132908 obesity risk allele influenced expression of FAIM2 and other genes and decreased the proportion of neurons produced by differentiation. We have functionally validated rs7132908 as a causal obesity variant that temporally regulates nearby effector genes and influences neurodevelopment and survival.

Efficacy of a Comprehensive Weight Reduction Intervention in Male Adolescents With Different FTO Genotypes

BACKGROUND

The FTO gene polymorphisms may influence the effects of lifestyle interventions on obesity. The present study aimed to assess the influence of the rs9930506 FTO gene polymorphism on the success of a comprehensive weight loss intervention in male adolescents with overweight and obesity.

METHODS

This study was carried out on 96 adolescent boys with overweight and obesity who were randomly assigned to the intervention (n = 53) and control (n = 43) groups. The blood samples of the participants were collected, and the FTO gene was genotyped for the rs9930506 polymorphism. A comprehensive lifestyle intervention including changes in diet and physical activity was performed for 8 weeks in the intervention group.

RESULTS

Following the lifestyle intervention, BMI and fat mass decreased significantly in the intervention group compared with the control group (both p < 0.05), while no change was found in weight, height or body muscle percentage between the groups. The participants in the intervention group with the AA/AG genotype and not in carriers of the GG genotype had a significantly higher reduction in BMI (-1.21 vs. 1.87 kg/m2, F = 4.07, p < 0.05) compared with the control group.

CONCLUSION

The intervention in individuals with the AA/AG genotype has been significantly effective in weight loss compared with the control group. The intervention had no association effect on anthropometric indices in adolescents with the GG genotype of the FTO rs9930506 polymorphism.

TRIAL REGISTRATION

Name of the registry: National Nutrition and Food Technology Research Institute; Trial registration number: IRCT2016020925699N2; Date of registration: 24/04/2016; URL of trial registry record: https://www.irct.ir/trial/21447.

Genetics, Nutrition, and Health: A New Frontier in Disease Prevention

The field of nutrition research has traditionally focused on the effects of macronutrients and micronutrients on the body. However, it has become evident that individuals have unique genetic makeups that influence their response to food. Nutritional genomics, which includes nutrigenetics and nutrigenomics, explores the interaction between an individual's genetic makeup, diet, and health outcomes. Nutrigenetics studies the impact of genetic variation on an individual's response to dietary nutrients, while nutrigenomics investigates how dietary components affect gene regulation and expression. These disciplines seek to understand the impact of diet on the genome, transcriptome, proteome, and metabolome. It provides insights into the mechanisms underlying the effect of diet on gene expression. Nutrients can cause the modification of genetic expression through epigenetic changes, such as DNA methylation and histone modifications. The aim of nutrigenomics is to create personalized diets based on the unique metabolic profile of an individual, gut microbiome, and genetic makeup to prevent diseases and promote health. Nutrigenomics has the potential to revolutionize the field of nutrition by combining the practicality of personalized nutrition with knowledge of genetic factors underlying health and disease. Thus, nutrigenomics offers a promising approach to improving health outcomes (in terms of disease prevention) through personalized nutrition strategies based on an individual's genetic and metabolic characteristics.

FTO genotype and body mass index reduction in childhood obesity interventions: A systematic review and meta-analysis

Numerous guidelines have called for personalized interventions to address childhood obesity. The role of fat mass and obesity-associated gene (FTO) in the risk of childhood obesity has been summarized. However, it remains unclear whether FTO could influence individual responses to obesity interventions, especially in children. To address this, we systematically reviewed 12,255 records across 10 databases/registers and included 13 lifestyle-based obesity interventions (3980 children with overweight/obesity) reporting changes in body mass index (BMI) Z-score, BMI, waist circumference, waist-to-hip ratio, and body fat percentage after interventions. These obesity-related outcomes were first compared between children carrying different FTO genotypes (rs9939609 or its proxy) and then synthesized by random-effect meta-analysis models. The results from single-group interventions showed no evidence of associations between FTO risk allele and changes in obesity-related outcomes after interventions (e.g., BMI Z-score: -0.01; 95% CI: -0.04, 0.01). The results from controlled trials showed that associations between the FTO risk allele and changes in obesity-related outcomes did not differ by intervention/control group. To conclude, the FTO risk allele might play a minor role in the response to obesity interventions among children. Future studies might pay more attention to the accumulation effect of multiple genes in the intervention process among children.

Investigation of inherited noncoding genetic variation impacting the pharmacogenomics of childhood acute lymphoblastic leukemia treatment

Defining genetic factors impacting chemotherapy failure can help to better predict response and identify drug resistance mechanisms. However, there is limited understanding of the contribution of inherited noncoding genetic variation on inter-individual differences in chemotherapy response in childhood acute lymphoblastic leukemia (ALL). Here we map inherited noncoding variants associated with treatment outcome and/or chemotherapeutic drug resistance to ALL cis-regulatory elements and investigate their gene regulatory potential and target gene connectivity using massively parallel reporter assays and three-dimensional chromatin looping assays, respectively. We identify 54 variants with transcriptional effects and high-confidence gene connectivity. Additionally, functional interrogation of the top variant, rs1247117, reveals changes in chromatin accessibility, PU.1 binding affinity and gene expression, and deletion of the genomic interval containing rs1247117 sensitizes cells to vincristine. Together, these data demonstrate that noncoding regulatory variants associated with diverse pharmacological traits harbor significant effects on allele-specific transcriptional activity and impact sensitivity to antileukemic agents.

Atp1a2 and Kcnj9 are candidate genes underlying oxycodone behavioral sensitivity and withdrawal in C57BL/6 substrains

Opioid use disorder is heritable, yet its genetic etiology is largely unknown. Analysis of addiction model traits in rodents (e.g., opioid behavioral sensitivity and withdrawal) can facilitate genetic and mechanistic discovery. C57BL/6J and C57BL/6NJ substrains have extremely limited genetic diversity, yet can show reliable phenotypic diversity which together, can facilitate gene discovery. The C57BL/6NJ substrain was less sensitive to oxycodone (OXY)-induced locomotor activity compared to the C57BL/6J substrain. Quantitative trait locus (QTL) mapping in an F2 cross identified a distal chromosome 1 QTL explaining 7-12% of the variance in OXY locomotor sensitivity and anxiety-like withdrawal in the elevated plus maze. We identified a second QTL for withdrawal on chromosome 5 near the candidate gene Gabra2 (alpha-2 subunit of GABA-A receptor) explaining 9% of the variance. Next, we generated recombinant lines from an F2 founder spanning the distal chromosome 1 locus (163-181 Mb), captured the QTL for OXY sensitivity and withdrawal, and fine-mapped a 2.45-Mb region (170.16-172.61 Mb). There were five striatal cis-eQTL transcripts in this region (Pcp4l1, Ncstn, Atp1a2, Kcnj9, Igsf9), two of which were confirmed at the protein level (KCNJ9, ATP1A2). Kcnj9, a.k.a., GIRK3, codes for a potassium channel that is a major effector of mu opioid receptor signaling. Atp1a2 codes for a subunit of a Na+/K+ ATPase enzyme that regulates neuronal excitability and shows adaptations following chronic opioid administration. To summarize, we identified genetic sources of opioid behavioral differences in C57BL/6 substrains, two of the most widely and often interchangeably used substrains in opioid addiction research.

A novel multi-omics mendelian randomization method for gene set enrichment and its application to psychiatric disorders

Genome-wide association studies (GWAS) of psychiatric disorders (PD) yield numerous loci with significant signals, but often do not implicate specific genes. Because GWAS risk loci are enriched in expression/protein/methylation quantitative loci (e/p/mQTL, hereafter xQTL), transcriptome/proteome/methylome-wide association studies (T/P/MWAS, hereafter XWAS) that integrate xQTL and GWAS information, can link GWAS signals to effects on specific genes. To further increase detection power, gene signals are aggregated within relevant gene sets (GS) by performing gene set enrichment (GSE) analyses. Often GSE methods test for enrichment of "signal" genes in curated GS while overlooking their linkage disequilibrium (LD) structure, allowing for the possibility of increased false positive rates. Moreover, no GSE tool uses xQTL information to perform mendelian randomization (MR) analysis. To make causal inference on association between PD and GS, we develop a novel MR GSE (MR-GSE) procedure. First, we generate a "synthetic" GWAS for each MSigDB GS by aggregating summary statistics for x-level (mRNA, protein or DNA methylation (DNAm) levels) from the largest xQTL studies available) of genes in a GS. Second, we use synthetic GS GWAS as exposure in a generalized summary-data-based-MR analysis of complex trait outcomes. We applied MR-GSE to GWAS of nine important PD. When applied to the underpowered opioid use disorder GWAS, none of the four analyses yielded any signals, which suggests a good control of false positive rates. For other PD, MR-GSE greatly increased the detection of GO terms signals (2,594) when compared to the commonly used (non-MR) GSE method (286). Some of the findings might be easier to adapt for treatment, e.g., our analyses suggest modest positive effects for supplementation with certain vitamins and/or omega-3 for schizophrenia, bipolar and major depression disorder patients. Similar to other MR methods, when applying MR-GSE researchers should be mindful of the confounding effects of horizontal pleiotropy on statistical inference.

Prioritizing susceptibility genes for the prognosis of male-pattern baldness with transcriptome-wide association study

BACKGROUND

Male-pattern baldness (MPB) is the most common cause of hair loss in men. It can be categorized into three types: type 2 (T2), type 3 (T3), and type 4 (T4), with type 1 (T1) being considered normal. Although various MPB-associated genetic variants have been suggested, a comprehensive study for linking these variants to gene expression regulation has not been performed to the best of our knowledge.

RESULTS

In this study, we prioritized MPB-related tissue panels using tissue-specific enrichment analysis and utilized single-tissue panels from genotype-tissue expression version 8, as well as cross-tissue panels from context-specific genetics. Through a transcriptome-wide association study and colocalization analysis, we identified 52, 75, and 144 MPB associations for T2, T3, and T4, respectively. To assess the causality of MPB genes, we performed a conditional and joint analysis, which revealed 10, 11, and 54 putative causality genes for T2, T3, and T4, respectively. Finally, we conducted drug repositioning and identified potential drug candidates that are connected to MPB-associated genes.

CONCLUSIONS

Overall, through an integrative analysis of gene expression and genotype data, we have identified robust MPB susceptibility genes that may help uncover the underlying molecular mechanisms and the novel drug candidates that may alleviate MPB.

Increased enhancer-promoter interactions during developmental enhancer activation in mammals

Remote enhancers are thought to interact with their target promoters via physical proximity, yet the importance of this proximity for enhancer function remains unclear. Here we investigate the three-dimensional (3D) conformation of enhancers during mammalian development by generating high-resolution tissue-resolved contact maps for nearly a thousand enhancers with characterized in vivo activities in ten murine embryonic tissues. Sixty-one percent of developmental enhancers bypass their neighboring genes, which are often marked by promoter CpG methylation. The majority of enhancers display tissue-specific 3D conformations, and both enhancer-promoter and enhancer-enhancer interactions are moderately but consistently increased upon enhancer activation in vivo. Less than 14% of enhancer-promoter interactions form stably across tissues; however, these invariant interactions form in the absence of the enhancer and are likely mediated by adjacent CTCF binding. Our results highlight the general importance of enhancer-promoter physical proximity for developmental gene activation in mammals.

Tissue-specific enhancer-gene maps from multimodal single-cell data identify causal disease alleles

Translating genome-wide association study (GWAS) loci into causal variants and genes requires accurate cell-type-specific enhancer-gene maps from disease-relevant tissues. Building enhancer-gene maps is essential but challenging with current experimental methods in primary human tissues. Here we developed a nonparametric statistical method, SCENT (single-cell enhancer target gene mapping), that models association between enhancer chromatin accessibility and gene expression in single-cell or nucleus multimodal RNA sequencing and ATAC sequencing data. We applied SCENT to 9 multimodal datasets including >120,000 single cells or nuclei and created 23 cell-type-specific enhancer-gene maps. These maps were highly enriched for causal variants in expression quantitative loci and GWAS for 1,143 diseases and traits. We identified likely causal genes for both common and rare diseases and linked somatic mutation hotspots to target genes. We demonstrate that application of SCENT to multimodal data from disease-relevant human tissue enables the scalable construction of accurate cell-type-specific enhancer-gene maps, essential for defining noncoding variant function.

p53 rapidly restructures 3D chromatin organization to trigger a transcriptional response

Activation of the p53 tumor suppressor triggers a transcriptional program to control cellular response to stress. However, the molecular mechanisms by which p53 controls gene transcription are not completely understood. Here, we uncover the critical role of spatio-temporal genome architecture in this process. We demonstrate that p53 drives direct and indirect changes in genome compartments, topologically associating domains, and DNA loops prior to one hour of its activation, which escort the p53 transcriptional program. Focusing on p53-bound enhancers, we report 340 genes directly regulated by p53 over a median distance of 116 kb, with 74% of these genes not previously identified. Finally, we showcase that p53 controls transcription of distal genes through newly formed and pre-existing enhancer-promoter loops in a cohesin dependent manner. Collectively, our findings demonstrate a previously unappreciated architectural role of p53 as regulator at distinct topological layers and provide a reliable set of new p53 direct target genes that may help designs of cancer therapies.

Comprehensive functional interrogation of susceptibility loci in GWASs identified KIAA0391 as a novel oncogenic driver via regulating pyroptosis in NSCLC

Approximately 51 non-small-cell lung cancer (NSCLC) risk loci have been identified by genome-wide association studies (GWASs). We conducted a high throughput RNA-interference (RNAi) screening to identify the candidate causal genes in NSCLC risk loci. KIAA0391 at 14q13.1 had the highest score and could promote proliferation and metastasis of NSCLC in vitro and in vivo. We next prioritized rs3783313 as a causal variant at 14q13.1, by integrating a large-scale population study consisting of 27,120 lung cancer cases and 27,355 controls, functional annotation, and expression quantitative trait locus (eQTL) analysis. Then we found that rs3783313 could facilitate a promoter-enhancer interaction to upregulate KIAA0391 expression by affecting the affinity of transcription factor NFYA. Mechanistically, KIAA0391 knockdown dramatically influenced pyroptosis-related pathways and increased the expression of CASP1. And KIAA0391 transcriptionally repressed CASP1 by binding to SMAD2 and induced an anti-pyroptosis phenotype, promoting tumorigenesis of NSCLC, which provides new insights and potential target for NSCLC.

Investigating the Influence of ANTXR2 Gene Mutations on Protective Antigen Binding for Heightened Anthrax Resistance

Bacillus anthracis is the bacterium responsible for causing the zoonotic disease called anthrax. The disease presents itself in different forms like gastrointestinal, inhalation, and cutaneous. Bacterial spores are tremendously adaptable, can persist for extended periods and occasionally endanger human health. The Anthrax Toxin Receptor-2 (ANTXR2) gene acts as membrane receptor and facilitates the entry of the anthrax toxin into host cells. Additionally, mutations in the ANTXR2 gene have been linked to various autoimmune diseases, including Hyaline Fibromatosis Syndrome (HFS), Ankylosing Spondylitis (AS), Juvenile Hyaline Fibromatosis (JHF), and Infantile Systemic Hyalinosis (ISH). This study delves into the genetic landscape of ANTXR2, aiming to comprehend its associations with diverse disorders, elucidate the impacts of its mutations, and pinpoint minimal non-pathogenic mutations capable of reducing the binding affinity of the ANTXR2 gene with the protective antigen. Recognizing the pivotal role of single-nucleotide polymorphisms (SNPs) in shaping genetic diversity, we conducted computational analyses to discern highly deleterious and tolerated non-synonymous SNPs (nsSNPs) in the ANTXR2 gene. The Mutpred2 server determined that the Arg465Trp alteration in the ANTXR2 gene leads to altered DNA binding (p = 0.22) with a probability of a deleterious mutation of 0.808; notably, among the identified deleterious SNPs, rs368288611 (Arg465Trp) stands out due to its significant impact on altering the DNA-binding ability of ANTXR2. We propose these SNPs as potential candidates for hypertension linked to the ANTXR2 gene, which is implicated in blood pressure regulation. Noteworthy among the tolerated substitutions is rs200536829 (Ala33Ser), recognized as less pathogenic; this highlights its potential as a valuable biomarker, potentially reducing side effects on the host while also reducing binding with the protective antigen protein. Investigating these SNPs holds the potential to correlate with several autoimmune disorders and mitigate the impact of anthrax disease in humans.

Are trait-associated genes clustered together in a gene network?

Genome-wide association studies (GWAS) have provided an abundance of information about the genetic variants and their loci that are associated to complex traits and diseases. However, due to linkage disequilibrium (LD) and noncoding regions of loci, it remains a challenge to pinpoint the causal genes. Gene network-based approaches, paired with network diffusion methods, have been proposed to prioritize causal genes and to boost statistical power in GWAS based on the assumption that trait-associated genes are clustered in a gene network. Due to the difficulty in mapping trait-associated variants to genes in GWAS, this assumption has never been directly or rigorously tested empirically. On the other hand, whole exome sequencing (WES) data focuses on the protein-coding regions, directly identifying trait-associated genes. In this study, we tested the assumption by leveraging the recently available exome-based association statistics from the UK Biobank WES data along with two types of networks. We found that almost all trait-associated genes were significantly more proximal to each other than randomly selected genes within both networks. These results support the assumption that trait-associated genes are clustered in gene networks, which can be further leveraged to boost the power of GWAS such as by introducing less stringent p value thresholds.

Variant to gene mapping for carpal tunnel syndrome risk loci implicates skeletal muscle regulatory elements

BACKGROUND

Carpal tunnel syndrome (CTS) is a common disorder caused by compression of the median nerve in the wrist, resulting in pain and numbness throughout the hand and forearm. While multiple behavioural and physiological factors influence CTS risk, a growing body of evidence supports a strong genetic contribution. Recent genome-wide association study (GWAS) efforts have reported 53 independent signals associated with CTS. While GWAS can identify genetic loci conferring risk, it does not determine which cell types drive the genetic aetiology of the trait, which variants are "causal" at a given signal, and which effector genes correspond to these non-coding variants. These obstacles limit interpretation of potential disease mechanisms.

METHODS

We analysed CTS GWAS findings in the context of chromatin conformation between gene promoters and accessible chromatin regions across cellular models of bone, skeletal muscle, adipocytes and neurons. We identified proxy variants in high LD with the lead CTS sentinel SNPs residing in promoter connected open chromatin in the skeletal muscle and bone contexts.

FINDINGS

We detected significant enrichment for heritability in skeletal muscle myotubes, as well as a weaker correlation in human mesenchymal stem cell-derived osteoblasts. In myotubes, our approach implicated 117 genes contacting 60 proxy variants corresponding to 20 of the 53 GWAS signals. In the osteoblast context we implicated 30 genes contacting 24 proxy variants coinciding with 12 signals, of which 19 genes shared. We subsequently prioritized BZW2 as a candidate effector gene in CTS and implicated it as novel gene that perturbs myocyte differentiation in vitro.

INTERPRETATION

Taken together our results suggest that the CTS genetic component influences the size, integrity, and organization of multiple tissues surrounding the carpal tunnel, in particular muscle and bone, to predispose the nerve to being compressed in this disease setting.

FUNDING

This work was supported by NIH Grant UM1 DK126194 (SFAG and WY), R01AG072705 (SFAG & KDH) and the Center for Spatial and Functional Genomics at CHOP (SFAG & ADW). SFAG is supported by the Daniel B. Burke Endowed Chair for Diabetes Research. WY is supported by the Perelman School of Medicine of the University of Pennsylvania.

Unraveling Adipose Tissue Dysfunction: Molecular Mechanisms, Novel Biomarkers, and Therapeutic Targets for Liver Fat Deposition

Adipose tissue (AT), once considered a mere fat storage organ, is now recognized as a dynamic and complex entity crucial for regulating human physiology, including metabolic processes, energy balance, and immune responses. It comprises mainly two types: white adipose tissue (WAT) for energy storage and brown adipose tissue (BAT) for thermogenesis, with beige adipocytes demonstrating the plasticity of these cells. WAT, beyond lipid storage, is involved in various metabolic activities, notably lipogenesis and lipolysis, critical for maintaining energy homeostasis. It also functions as an endocrine organ, secreting adipokines that influence metabolic, inflammatory, and immune processes. However, dysfunction in WAT, especially related to obesity, leads to metabolic disturbances, including the inability to properly store excess lipids, resulting in ectopic fat deposition in organs like the liver, contributing to non-alcoholic fatty liver disease (NAFLD). This narrative review delves into the multifaceted roles of WAT, its composition, metabolic functions, and the pathophysiology of WAT dysfunction. It also explores diagnostic approaches for adipose-related disorders, emphasizing the importance of accurately assessing AT distribution and understanding the complex relationships between fat compartments and metabolic health. Furthermore, it discusses various therapeutic strategies, including innovative therapeutics like adipose-derived mesenchymal stem cells (ADMSCs)-based treatments and gene therapy, highlighting the potential of precision medicine in targeting obesity and its associated complications.

Dietary macronutrient composition impacts gene regulation in adipose tissue

Diet is a key lifestyle component that influences metabolic health through several factors, including total energy intake and macronutrient composition. While the impact of caloric intake on gene expression and physiological phenomena in various tissues is well described, the influence of dietary macronutrient composition on these parameters is less well studied. Here, we use the Nutritional Geometry framework to investigate the role of macronutrient composition on metabolic function and gene regulation in adipose tissue. Using ten isocaloric diets that vary systematically in their proportion of energy from fat, protein, and carbohydrates, we find that gene expression and splicing are highly responsive to macronutrient composition, with distinct sets of genes regulated by different macronutrient interactions. Specifically, the expression of many genes associated with Bardet-Biedl syndrome is responsive to dietary fat content. Splicing and expression changes occur in largely separate gene sets, highlighting distinct mechanisms by which dietary composition influences the transcriptome and emphasizing the importance of considering splicing changes to more fully capture the gene regulation response to environmental changes such as diet. Our study provides insight into the gene regulation plasticity of adipose tissue in response to macronutrient composition, beyond the already well-characterized response to caloric intake.

3D genomic features across >50 diverse cell types reveal insights into the genomic architecture of childhood obesity

The prevalence of childhood obesity is increasing worldwide, along with the associated common comorbidities of type 2 diabetes and cardiovascular disease in later life. Motivated by evidence for a strong genetic component, our prior genome-wide association study (GWAS) efforts for childhood obesity revealed 19 independent signals for the trait; however, the mechanism of action of these loci remains to be elucidated. To molecularly characterize these childhood obesity loci we sought to determine the underlying causal variants and the corresponding effector genes within diverse cellular contexts. Integrating childhood obesity GWAS summary statistics with our existing 3D genomic datasets for 57 human cell types, consisting of high-resolution promoter-focused Capture-C/Hi-C, ATAC-seq, and RNA-seq, we applied stratified LD score regression and calculated the proportion of genome-wide SNP heritability attributable to cell type-specific features, revealing pancreatic alpha cell enrichment as the most statistically significant. Subsequent chromatin contact-based fine-mapping was carried out for genome-wide significant childhood obesity loci and their linkage disequilibrium proxies to implicate effector genes, yielded the most abundant number of candidate variants and target genes at the BDNF, ADCY3, TMEM18 and FTO loci in skeletal muscle myotubes and the pancreatic beta-cell line, EndoC-BH1. One novel implicated effector gene, ALKAL2 - an inflammation-responsive gene in nerve nociceptors - was observed at the key TMEM18 locus across multiple immune cell types. Interestingly, this observation was also supported through colocalization analysis using expression quantitative trait loci (eQTL) derived from the Genotype-Tissue Expression (GTEx) dataset, supporting an inflammatory and neurologic component to the pathogenesis of childhood obesity. Our comprehensive appraisal of 3D genomic datasets generated in a myriad of different cell types provides genomic insights into pediatric obesity pathogenesis.

1,2-bis(2,4,6-tribromophenoxy) ethane, a novel brominated flame retardant, disrupts intestinal barrier function via the IRX3/NOS2 axis in rat small intestine

Novel brominated flame retardants are widely used in electronics, textiles, furniture, and other products; they can enter the human body through ingestion and respiration and cause harm to the human body, and have been proven to have potential biological toxicity and accumulation effects. 1,2-bis(2,4,6-tribromophenoxy) ethane (BTBPE) is a widely used novel brominated flame retardant; however, there is a lack of research on its mechanism of toxicity, particularly that of intestinal toxicity. Currently, studies on the functionality of iroquois homeobox 3 (IRX3) are extremely limited. In our study, BTBPE was administered to Sprague-Dawley (SD) rats and rat small intestinal crypt epithelial cells (IEC6 cells) in vivo and in vitro, respectively, and hematoxylin and eosin (HE), immunohistochemical, Alcian blue-periodic acid-Schiff (AB-PAS), CCK8, acridine orange/ethidium bromide (AO/EB), fluorescent probes, qPCR, western blotting, and immunofluorescence analyses were performed. To explore the damage mechanism of BTBPE, we used siRNA to silence IRX3 and iNOs-IN-1 (yeast extract-peptone-wheat; YPW) to inhibit nitric oxide synthase 2 (NOS2). The results showed that BTBPE exposure caused inflammation and necroptosis in the jejunum and ileum, as well as destruction of the tight junctions and mucus layer. Moreover, BTBPE activated the IRX3/NOS2 axis both in vivo and in vitro. Silencing IRX3 or inhibiting NOS2 inhibits necroptosis and restores tight junctions in IEC6 cells. In conclusion, our study found that in the jejunum, ileum, and IEC6 cells, BTBPE exposure caused necroptosis and tight junction destruction by activating the IRX3/NOS2 axis. Blocking the IRX3/NOS2 axis can effectively inhibit necroptosis and restore tight junction. In addition, BTBPE exposure caused inflammation and loss of the mucous layer in the jejunum and ileum. Our study is the first to explore the mechanism of intestinal damage caused by BTBPE exposure and to discover new biological functions regulated by the IRX3/NOS2 axis, providing new research directions for necroptosis and tight junctions.

Hypothalamic FTO promotes high-fat diet-induced leptin resistance in mice through increasing CX3CL1 expression

Long-term consumption of a high-fat diet (HFD) disrupts energy homeostasis and leads to weight gain. The fat mass and obesity-associated (FTO) gene has been consistently identified to be associated with HFD-induced obesity. The hypothalamus is crucial for regulating energy balance, and HFD-induced hypothalamic leptin resistance contributes to obesity. FTO, an N6-methyladenosine (m6A) RNA methylation regulator, may be a key mediator of leptin resistance. However, the exact mechanisms remain unclear. Therefore, the present study aims to investigate the association between FTO and leptin resistance. After HFD or standard diet (SD) feeding in male mice for 22 weeks, m6A-sequencing and western blotting assays were used to identify target genes and assess protein level, and molecular interaction changes. CRISPR/Cas9 gene knockout system was employed to investigate the potential function of FTO in leptin resistance and obesity. Our data showed that chemokine (C-X3-C motif) ligand 1 (CX3CL1) was a direct downstream target of FTO-mediated m6A modification. Furthermore, upregulation of FTO/CX3CL1 and suppressor of cytokine signaling 3 (SOCS3) in the hypothalamus impaired leptin-signal transducer and activator of transcription 3 signaling, resulting in leptin resistance and obesity. Compared to wild-type (WT) mice, FTO deficiency in leptin receptor-expressing neurons of the hypothalamus significantly inhibited the upregulation of CX3CL1 and SOCS3, and partially ameliorating leptin resistance under HFD conditions. Our findings reveal that FTO involved in the hypothalamic leptin resistance and provides novel insight into the function of FTO in the contribution to hypothalamic leptin resistance and obesity.

Characterization of metabolic phenotypes and distinctive genes in mice with low-weight gain

Being overweight exacerbates various metabolic diseases, necessitating the identification of target molecules for obesity control. In the current study, we investigated common physiological features related to metabolism in mice with low weight gain: (1) G protein-coupled receptor, family C, group 5, member B-knockout; (2) gastric inhibitory polypeptide receptor-knockout; and (3) Iroquois-related homeobox 3-knockout. Moreover, we explored genes involved in metabolism by analyzing differentially expressed genes (DEGs) between low-weight gain mice and the respective wild-type control mice. The common characteristics of the low-weight gain mice were low inguinal white adipose tissue (iWAT) and liver weight despite similar food intake along with lower blood leptin levels and high energy expenditure. The DEGs of iWAT, epididymal (gonadal) WAT, brown adipose tissue, muscle, liver, hypothalamus, and hippocampus common to these low-weight gain mice were designated as candidate genes associated with metabolism. One such gene tetraspanin 7 (Tspan7) from the iWAT was validated using knockout and overexpressing mouse models. Mice with low Tspan7 expression gained more weight, while those with high Tspan7 expression gained less weight, confirming the involvement of the Tspan7 gene in weight regulation. Collectively, these findings suggest that the candidate gene list generated in this study contains potential target molecules for obesity regulation. Further validation and additional data from low-weight gain mice will aid in understanding the molecular mechanisms associated with obesity.

Prioritization of risk genes in colorectal cancer by integrative analysis of multi-omics data and gene networks

Genome-wide association studies (GWASs) have identified over 140 colorectal cancer (CRC)-associated loci; however, target genes at the majority of loci and underlying molecular mechanisms are poorly understood. Here, we utilized a Bayesian approach, integrative risk gene selector (iRIGS), to prioritize risk genes at CRC GWAS loci by integrating multi-omics data. As a result, a total of 105 high-confidence risk genes (HRGs) were identified, which exhibited strong gene dependencies for CRC and enrichment in the biological processes implicated in CRC. Among the 105 HRGs, CEBPB, located at the 20q13.13 locus, acted as a transcription factor playing critical roles in cancer. Our subsequent assays indicated the tumor promoter function of CEBPB that facilitated CRC cell proliferation by regulating multiple oncogenic pathways such as MAPK, PI3K-Akt, and Ras signaling. Next, by integrating a fine-mapping analysis and three independent case-control studies in Chinese populations consisting of 8,039 cases and 12,775 controls, we elucidated that rs1810503, a putative functional variant regulating CEBPB, was associated with CRC risk (OR=0.90, 95%CI=0.86-0.93, P=1.07×10-7). The association between rs1810503 and CRC risk was further validated in three additional multi-ancestry populations consisting of 24,254 cases and 58,741 controls. Mechanistically, the rs1810503 A to T allele change weakened the enhancer activity in an allele-specific manner to decrease CEBPB expression via long-range promoter-enhancer interactions, mediated by the transcription factor, REST, and thus decreased CRC risk. In summary, our study provides a genetic resource and a generalizable strategy for CRC etiology investigation, and highlights the biological implications of CEBPB in CRC tumorigenesis, shedding new light on the etiology of CRC.

EnhancerTracker: Comparing cell-type-specific enhancer activity of DNA sequence triplets via an ensemble of deep convolutional neural networks

Motivation

Transcriptional enhancers - unlike promoters - are unrestrained by distance or strand orientation with respect to their target genes, making their computational identification a challenge. Further, there are insufficient numbers of confirmed enhancers for many cell types, preventing robust training of machine-learning-based models for enhancer prediction for such cell types.

Results

We present EnhancerTracker , a novel tool that leverages an ensemble of deep separable convolutional neural networks to identify cell-type-specific enhancers with the need of only two confirmed enhancers. EnhancerTracker is trained, validated, and tested on 52,789 putative enhancers obtained from the FANTOM5 Project and control sequences derived from the human genome. Unlike available tools, which accept one sequence at a time, the input to our tool is three sequences; the first two are enhancers active in the same cell type. EnhancerTracker outputs 1 if the third sequence is an enhancer active in the same cell type(s) where the first two enhancers are active. It outputs 0 otherwise. On a held-out set (15%), EnhancerTracker achieved an accuracy of 64%, a specificity of 93%, a recall of 35%, a precision of 84%, and an F1 score of 49%.

Availability and implementation

https://github.com/BioinformaticsToolsmith/EnhancerTracker.

Contact

hani.girgis@tamuk.edu.

Integrated omics analysis of coronary artery calcifications and myocardial infarction: the Framingham Heart Study

Gene function can be described using various measures. We integrated association studies of three types of omics data to provide insights into the pathophysiology of subclinical coronary disease and myocardial infarction (MI). Using multivariable regression models, we associated: (1) single nucleotide polymorphism, (2) DNA methylation, and (3) gene expression with coronary artery calcification (CAC) scores and MI. Among 3106 participants of the Framingham Heart Study, 65 (2.1%) had prevalent MI and 60 (1.9%) had incident MI, median CAC value was 67.8 [IQR 10.8, 274.9], and 1403 (45.2%) had CAC scores > 0 (prevalent CAC). Prevalent CAC was associated with AHRR (linked to smoking) and EXOC3 (affecting platelet function and promoting hemostasis). CAC score was associated with VWA1 (extracellular matrix protein associated with cartilage structure in endomysium). For prevalent MI we identified FYTTD1 (down-regulated in familial hypercholesterolemia) and PINK1 (linked to cardiac tissue homeostasis and ischemia-reperfusion injury). Incident MI was associated with IRX3 (enhancing browning of white adipose tissue) and STXBP3 (controlling trafficking of glucose transporter type 4 to plasma). Using an integrative trans-omics approach, we identified both putatively novel and known candidate genes associated with CAC and MI. Replication of findings is warranted.

Single-cell sequencing analysis related to sphingolipid metabolism guides immunotherapy and prognosis of skin cutaneous melanoma

Background

We explore sphingolipid-related genes (SRGs) in skin melanoma (SKCM) to develop a prognostic indicator for patient outcomes. Dysregulated lipid metabolism is linked to aggressive behavior in various cancers, including SKCM. However, the exact role and mechanism of sphingolipid metabolism in melanoma remain partially understood.

Methods

We integrated scRNA-seq data from melanoma patients sourced from the GEO database. Through the utilization of the Seurat R package, we successfully identified distinct gene clusters associated with patient survival in the scRNA-seq data. Key prognostic genes were identified through single-factor Cox analysis and used to develop a prognostic model using LASSO and stepwise regression algorithms. Additionally, we evaluated the predictive potential of these genes within the immune microenvironment and their relevance to immunotherapy. Finally, we validated the functional significance of the high-risk gene IRX3 through in vitro experiments.

Results

Analysis of scRNA-seq data identified distinct expression patterns of 4 specific genes (SRGs) in diverse cell subpopulations. Re-clustering cells based on increased SRG expression revealed 7 subgroups with significant prognostic implications. Using marker genes, lasso, and Cox regression, we selected 11 genes to construct a risk signature. This signature demonstrated a strong correlation with immune cell infiltration and stromal scores, highlighting its relevance in the tumor microenvironment. Functional studies involving IRX3 knockdown in A375 and WM-115 cells showed significant reductions in cell viability, proliferation, and invasiveness.

Conclusion

SRG-based risk signature holds promise for precise melanoma prognosis. An in-depth exploration of SRG characteristics offers insights into immunotherapy response. Therapeutic targeting of the IRX3 gene may benefit melanoma patients.

Trans -eQTL hotspots shape complex traits by modulating cellular states

Regulatory genetic variation shapes gene expression, providing an important mechanism connecting DNA variation and complex traits. The causal relationships between gene expression and complex traits remain poorly understood. Here, we integrated transcriptomes and 46 genetically complex growth traits in a large cross between two strains of the yeast Saccharomyces cerevisiae . We discovered thousands of genetic correlations between gene expression and growth, suggesting functional connections. Local regulatory variation was a minor source of these genetic correlations. Instead, genetic correlations tended to arise from multiple independent trans -acting regulatory loci. Trans -acting hotspots that affect the expression of numerous genes accounted for particularly large fractions of genetic growth variation and of genetic correlations between gene expression and growth. Genes with genetic correlations were enriched for similar biological processes across traits, but with heterogeneous direction of effect. Our results reveal how trans -acting regulatory hotspots shape complex traits by altering cellular states.

Single-nucleotide variants within heart enhancers increase binding affinity and disrupt heart development

Transcriptional enhancers direct precise gene expression patterns during development and harbor the majority of variants associated with phenotypic diversity, evolutionary adaptations, and disease. Pinpointing which enhancer variants contribute to changes in gene expression and phenotypes is a major challenge. Here, we find that suboptimal or low-affinity binding sites are necessary for precise gene expression during heart development. Single-nucleotide variants (SNVs) can optimize the affinity of ETS binding sites, causing gain-of-function (GOF) gene expression, cell migration defects, and phenotypes as severe as extra beating hearts in the marine chordate Ciona robusta. In human induced pluripotent stem cell (iPSC)-derived cardiomyocytes, a SNV within a human GATA4 enhancer increases ETS binding affinity and causes GOF enhancer activity. The prevalence of suboptimal-affinity sites within enhancers creates a vulnerability whereby affinity-optimizing SNVs can lead to GOF gene expression, changes in cellular identity, and organismal-level phenotypes that could contribute to the evolution of novel traits or diseases.

Studies on the fat mass and obesity-associated (FTO) gene and its impact on obesity-associated diseases

Obesity has become a major health crisis in the past ∼50 years. The fat mass and obesity-associated (FTO) gene, identified by genome-wide association studies (GWAS), was first reported to be positively associated with obesity in humans. Mice with more copies of the FTO gene were observed to be obese, while loss of the gene in mice was found to protect from obesity. Later, FTO was found to encode an m6A RNA demethylase and has a profound effect on many biological and metabolic processes. In this review, we first summarize recent studies that demonstrate the critical roles and regulatory mechanisms of FTO in obesity and metabolic disease. Second, we discuss the ongoing debates concerning the association between FTO polymorphisms and obesity. Third, since several small molecule drugs and micronutrients have been found to regulate metabolic homeostasis through controlling the expression or activity of FTO, we highlight the broad potential of targeting FTO for obesity treatment. Improving our understanding of FTO and the underlying mechanisms may provide new approaches for treating obesity and metabolic diseases.

Novel regulators of islet function identified from genetic variation in mouse islet Ca2+ oscillations

Insufficient insulin secretion to meet metabolic demand results in diabetes. The intracellular flux of Ca2+ into β-cells triggers insulin release. Since genetics strongly influences variation in islet secretory responses, we surveyed islet Ca2+ dynamics in eight genetically diverse mouse strains. We found high strain variation in response to four conditions: (1) 8 mM glucose; (2) 8 mM glucose plus amino acids; (3) 8 mM glucose, amino acids, plus 10 nM glucose-dependent insulinotropic polypeptide (GIP); and (4) 2 mM glucose. These stimuli interrogate β-cell function, α- to β-cell signaling, and incretin responses. We then correlated components of the Ca2+ waveforms to islet protein abundances in the same strains used for the Ca2+ measurements. To focus on proteins relevant to human islet function, we identified human orthologues of correlated mouse proteins that are proximal to glycemic-associated single-nucleotide polymorphisms in human genome-wide association studies. Several orthologues have previously been shown to regulate insulin secretion (e.g. ABCC8, PCSK1, and GCK), supporting our mouse-to-human integration as a discovery platform. By integrating these data, we nominate novel regulators of islet Ca2+ oscillations and insulin secretion with potential relevance for human islet function. We also provide a resource for identifying appropriate mouse strains in which to study these regulators.

Pathological Remodeling of the Gut Barrier as a Prodromal Event of High-Fat Diet-Induced Obesity

Intestinal barrier alterations represent a primum movens in obesity and related intestinal dysfunctions. However, whether gut barrier remodeling represents prodromal events in obesity before weight gain, metabolic alterations, and systemic inflammation remains unclear. Herein, we examined morphologic changes in the gut barrier in a mouse model of high-fat diet (HFD) since the earliest phases of diet assumption. C57BL/6J mice were fed with standard diet (SD) or HFD for 1, 2, 4, or 8 weeks. Remodeling of intestinal epithelial barrier, inflammatory infiltrate, and collagen deposition in the colonic wall was assessed by histochemistry and immunofluorescence analysis. Obese mice displayed increased body and epididymal fat weight along with increased plasma resistin, IL-1β, and IL-6 levels after 8 weeks of HFD. Starting from 1 week of HFD, mice displayed (1) a decreased claudin-1 expression in lining epithelial cells, (2) an altered mucus in goblet cells, (3) an increase in proliferating epithelial cells in colonic crypts, (4) eosinophil infiltration along with an increase in vascular P-selectin, and (5) deposition of collagen fibers. HFD intake is associated with morphologic changes in the large bowel at mucosal and submucosal levels. In particular, the main changes include alterations in the mucous layer and intestinal epithelial barrier integrity and activation of mucosal defense-enhanced fibrotic deposition. These changes represent early events occurring before the development of obesity condition that could contribute to compromising the intestinal mucosal barrier and functions, opening the way for systemic dissemination.

RNA N6-methyladenosine of DHAPAT and PAP involves in regulation of diapause of Bombyx mori via the lipid metabolism pathway

Environment-induced epigenetics are involved in diapause regulation, but the molecular mechanism that epigenetically couples nutrient metabolism to diapause regulation remains unclear. In this study, we paid special attention to the significant differences in the level of N6-adenosine methylation (m6A) of dihydroxyacetone phosphate acyltransferase (DHAPAT) and phosphatidate phosphatase (PAP) genes in the lipid metabolism pathway of the bivoltine silkworm (Bombyx mori) strain Qiufeng developed from eggs incubated at a normal temperature (QFHT, diapause egg producer) compared to those from eggs incubated at a low temperature (QFLT, non-diapause egg producer). We knocked down DHAPAT in the pupal stage of the QFLT group, resulting in the non-diapause destined eggs becoming diapausing eggs. In the PAP knockdown group, the colour of the non-diapause destined eggs changed from light yellow to pink 3 days after oviposition, but they hatched as normal. Moreover, we validated that YTHDF3 binds to m6A-modified DHAPAT and PAP mRNAs to promote their stability and translation. These results suggest that RNA m6A methylation participates in the diapause regulation of silkworm by changing the expression levels of DHAPAT and PAP and reveal that m6A epigenetic modification can be combined with a lipid metabolism signal pathway to participate in the regulation of insect diapause traits, which provides a clearer image for exploring the physiological basis of insect diapause.

Assessing efficiency of fine-mapping obesity-associated variants through leveraging ancestry architecture and functional annotation using PAGE and UKBB cohorts

Inadequate representation of non-European ancestry populations in genome-wide association studies (GWAS) has limited opportunities to isolate functional variants. Fine-mapping in multi-ancestry populations should improve the efficiency of prioritizing variants for functional interrogation. To evaluate this hypothesis, we leveraged ancestry architecture to perform comparative GWAS and fine-mapping of obesity-related phenotypes in European ancestry populations from the UK Biobank (UKBB) and multi-ancestry samples from the Population Architecture for Genetic Epidemiology (PAGE) consortium with comparable sample sizes. In the investigated regions with genome-wide significant associations for obesity-related traits, fine-mapping in our ancestrally diverse sample led to 95% and 99% credible sets (CS) with fewer variants than in the European ancestry sample. Lead fine-mapped variants in PAGE regions had higher average coding scores, and higher average posterior probabilities for causality compared to UKBB. Importantly, 99% CS in PAGE loci contained strong expression quantitative trait loci (eQTLs) in adipose tissues or harbored more variants in tighter linkage disequilibrium (LD) with eQTLs. Leveraging ancestrally diverse populations with heterogeneous ancestry architectures, coupled with functional annotation, increased fine-mapping efficiency and performance, and reduced the set of candidate variants for consideration for future functional studies. Significant overlap in genetic causal variants across populations suggests generalizability of genetic mechanisms underpinning obesity-related traits across populations.

Epromoters are new players in the regulatory landscape with potential pleiotropic roles

Precise spatiotemporal control of gene expression during normal development and cell differentiation is achieved by the combined action of proximal (promoters) and distal (enhancers) cis-regulatory elements. Recent studies have reported that a subset of promoters, termed Epromoters, works also as enhancers to regulate distal genes. This new paradigm opened novel questions regarding the complexity of our genome and raises the possibility that genetic variation within Epromoters has pleiotropic effects on various physiological and pathological traits by differentially impacting multiple proximal and distal genes. Here, we discuss the different observations pointing to an important role of Epromoters in the regulatory landscape and summarize the evidence supporting a pleiotropic impact of these elements in disease. We further hypothesize that Epromoter might represent a major contributor to phenotypic variation and disease.

The interaction between enhancer variants and environmental factors as an overlooked aetiological paradigm in human complex disease

The interactions between genetic and environmental risk factors contribute to the aetiology of complex human diseases. Genome-wide association studies (GWAS) have revealed that most of the genetic variants associated with complex diseases are located in the non-coding part of the genome, preferentially within enhancers. Enhancers are distal cis-regulatory elements composed of clusters of transcription factors binding sites that positively regulate the expression of their target genes. The generation of genome-wide maps for histone marks (e.g., H3K27ac), chromatin accessibility and transcription factor and coactivator (e.g., p300) binding profiles have enabled the identification of enhancers across many human cell types and tissues. Nonetheless, the functional and pathological consequences of the majority of disease-associated genetic variants located within enhancers seem to be rather minor under normal conditions, thus questioning their medical relevance. Here we propose that, due to the prevalence of enhancer redundancy, the pathological effects of many disease-associated non-coding genetic variants might be preferentially (or even only) manifested under environmental stress.

Low levels of circulating methylated IRX3 are related to worse outcome after transcatheter aortic valve implantation in patients with severe aortic stenosis

BACKGROUND

Aortic stenosis (AS) is one of the most common cardiac diseases and major cause of morbidity and mortality in the elderly. Transcatheter aortic valve implantation (TAVI) is performed in such patients with symptomatic severe AS and reduces mortality for the majority of these patients. However, a significant percentage dies within the first two years after TAVI, such that there is an interest to identify parameters, which predict outcome and could guide pre-TAVI patient selection. High levels of cardiac fibrosis have been identified as such independent predictor of cardiovascular mortality after TAVI. Promoter hypermethylation commonly leads to gene downregulation, and the Iroquois homeobox 3 (IRX3) gene was identified in a genome-wide transcriptome and methylome to be hypermethylated and downregulated in AS patients. In a well-described cohort of 100 TAVI patients in which cardiac fibrosis levels were quantified histologically in cardiac biopsies, and which had a follow-up of up to two years, we investigated if circulating methylated DNA of IRX3 in the peripheral blood is associated with cardiac fibrosis and/or mortality in AS patients undergoing TAVI and thus could serve as a biomarker to add information on outcome after TAVI.

RESULTS

Patients with high levels of methylation in circulating IRX3 show a significantly increased survival as compared to patients with low levels of IRX3 methylation indicating that high peripheral IRX3 methylation is associated with an improved outcome. In the multivariable setting, peripheral IRX3 methylation acts as an independent predictor of all-cause mortality. While there is no significant correlation of levels of IRX3 methylation with cardiac death, there is a significant but very weak inverse correlation between circulating IRX3 promoter methylation level and the amount of cardiac fibrosis. Higher levels of peripheral IRX3 methylation further correlated with decreased cardiac IRX3 expression and vice versa.

CONCLUSIONS

High levels of IRX3 methylation in the blood of AS patients at the time of TAVI are associated with better overall survival after TAVI and at least partially reflect myocardial IRX3 expression. Circulating methylated IRX3 might aid as a potential biomarker to help guide both pre-TAVI patient selection and post-TAVI monitoring.

Cell Type-Specific Whole-Genome Landscape of ΔFOSB Binding in the Nucleus Accumbens After Chronic Cocaine Exposure

BACKGROUND

The ability of neurons to respond to external stimuli involves adaptations of gene expression. Induction of the transcription factor ΔFOSB in the nucleus accumbens, a key brain reward region, is important for the development of drug addiction. However, a comprehensive map of ΔFOSB's gene targets has not yet been generated.

METHODS

We used CUT&RUN (cleavage under targets and release using nuclease) to map the genome-wide changes in ΔFOSB binding in the 2 main types of nucleus accumbens neurons-D1 or D2 medium spiny neurons-after chronic cocaine exposure. To annotate genomic regions of ΔFOSB binding sites, we also examined the distributions of several histone modifications. Resulting datasets were leveraged for multiple bioinformatic analyses.

RESULTS

The majority of ΔFOSB peaks occur outside promoter regions, including intergenic regions, and are surrounded by epigenetic marks indicative of active enhancers. BRG1, the core subunit of the SWI/SNF chromatin remodeling complex, overlaps with ΔFOSB peaks, a finding consistent with earlier studies of ΔFOSB's interacting proteins. Chronic cocaine use induces broad changes in ΔFOSB binding in both D1 and D2 nucleus accumbens medium spiny neurons of male and female mice. In addition, in silico analyses predict that ΔFOSB cooperatively regulates gene expression with homeobox and T-box transcription factors.

CONCLUSIONS

These novel findings uncover key elements of ΔFOSB's molecular mechanisms in transcriptional regulation at baseline and in response to chronic cocaine exposure. Further characterization of ΔFOSB's collaborative transcriptional and chromatin partners specifically in D1 and D2 medium spiny neurons will reveal a broader picture of the function of ΔFOSB and the molecular basis of drug addiction.

GCPBayes pipeline: a tool for exploring pleiotropy at the gene level

Cross-phenotype association using gene-set analysis can help to detect pleiotropic genes and inform about common mechanisms between diseases. Although there are an increasing number of statistical methods for exploring pleiotropy, there is a lack of proper pipelines to apply gene-set analysis in this context and using genome-scale data in a reasonable running time. We designed a user-friendly pipeline to perform cross-phenotype gene-set analysis between two traits using GCPBayes, a method developed by our team. All analyses could be performed automatically by calling for different scripts in a simple way (using a Shiny app, Bash or R script). A Shiny application was also developed to create different plots to visualize outputs from GCPBayes. Finally, a comprehensive and step-by-step tutorial on how to use the pipeline is provided in our group's GitHub page. We illustrated the application on publicly available GWAS (genome-wide association studies) summary statistics data to identify breast cancer and ovarian cancer susceptibility genes. We have shown that the GCPBayes pipeline could extract pleiotropic genes previously mentioned in the literature, while it also provided new pleiotropic genes and regions that are worthwhile for further investigation. We have also provided some recommendations about parameter selection for decreasing computational time of GCPBayes on genome-scale data.

Cell Differentiation and Aging Lead To Up-Regulation of FTO, While the ALKBH5 Protein Level Was Stable During Aging but Up-Regulated During in vitro-Induced Cardiomyogenesis

FTO and ALKBH5 proteins are essential erasers of N6-adenosine methylation in RNA. We studied how levels of FTO and ALKBH5 proteins changed during mouse embryonic development, aging, cardiomyogenesis, and neuroectodermal differentiation. We observed that aging in male and female mice was associated with FTO up-regulation in mouse hearts, brains, lungs, and kidneys, while the ALKBH5 level remained stable. FTO and ALKBH5 proteins were up-regulated during experimentally induced cardiomyogenesis, but the level of ALKBH5 protein was not changed when neuroectodermal differentiation was induced. HDAC1 depletion in mouse ES cells caused FTO down-regulation. In these cells, mRNA, carrying information from genes that regulate histone signature, RNA processing, and cell differentiation, was characterized by a reduced level of N6-adenosine methylation in specific gene loci, primarily regulating cell differentiation into neuroectoderm. Together, when we compared both RNA demethylating proteins, the FTO protein level undergoes the most significant changes during cell differentiation and aging. Thus, we conclude that during aging and neuronal differentiation, m6A RNA demethylation is likely regulated by the FTO protein but not via the function of ALKBH5.

Identification and validation of key biomarkers based on RNA methylation genes in sepsis

Background

RNA methylation is closely involved in immune regulation, but its role in sepsis remains unknown. Here, we aim to investigate the role of RNA methylation-associated genes (RMGs) in classifying and diagnosing of sepsis.

Methods

Five types of RMGs (m1A, m5C, m6Am, m7G and Ψ) were used to identify sepsis subgroups based on gene expression profile data obtained from the GEO database (GSE57065, GSE65682, and GSE95233). Unsupervised clustering analysis was used to identify distinct RNA modification subtypes. The CIBERSORT, WGCNA, GO and KEGG analysis were performed to explore immune infiltration pattern and biological function of each cluster. RF, SVM, XGB, and GLM algorithm were applied to identify the diagnostic RMGs in sepsis. Finally, the expression levels of the five key RMGs were verified by collecting PBMCs from septic patients using qRT-PCR, and their diagnostic efficacy for sepsis was verified in combination with clinical data using ROC analysis.

Results

Sepsis was divided into three subtypes (cluster 1 to 3). Cluster 1 highly expressed NSUN7 and TRMT6, with the characteristic of neutrophil activation and upregulation of MAPK signaling pathways. Cluster 2 highly expressed NSUN3, and was featured by the regulation of mRNA stability and amino acid metabolism. NSUN5 and NSUN6 were upregulated in cluster 3 which was involved in ribonucleoprotein complex biogenesis and carbohydrate metabolism pathways. In addition, we identified that five RMGs (NSUN7, NOP2, PUS1, PUS3 and FTO) could function as biomarkers for clinic diagnose of sepsis. For validation, we determined that the relative expressions of NSUN7, NOP2, PUS1 and PUS3 were upregulated, while FTO was downregulated in septic patients. The area under the ROC curve (AUC) of NSUN7, NOP2, PUS1, PUS3 and FTO was 0.828, 0.707, 0.846, 0.834 and 0.976, respectively.

Conclusions

Our study uncovered that dysregulation of RNA methylation genes (m1A, m5C, m6Am, m7G and Ψ) was closely involved in the pathogenesis of sepsis, providing new insights into the classification of sepsis endotypes. We also revealed that five hub RMGs could function as novel diagnostic biomarkers and potential targets for treatment.

Variant-to-function analysis of the childhood obesity chr12q13 locus implicates rs7132908 as a causal variant within the 3' UTR of FAIM2

The ch12q13 obesity locus is among the most significant childhood obesity loci identified in genome-wide association studies. This locus resides in a non-coding region within FAIM2; thus, the underlying causal variant(s) presumably influence disease susceptibility via an influence on cis-regulation within the genomic region. We implicated rs7132908 as a putative causal variant at this locus leveraging a combination of our inhouse 3D genomic data, public domain datasets, and several computational approaches. Using a luciferase reporter assay in human primary astrocytes, we observed allele-specific cis-regulatory activity of the immediate region harboring rs7132908. Motivated by this finding, we went on to generate isogenic human embryonic stem cell lines homozygous for either rs7132908 allele with CRISPR-Cas9 homology-directed repair to assess changes in gene expression due to genotype and chromatin accessibility throughout a differentiation to hypothalamic neurons, a key cell type known to regulate feeding behavior. We observed that the rs7132908 obesity risk allele influenced the expression of FAIM2 along with other genes, decreased the proportion of neurons produced during differentiation, up-regulated cell death gene sets, and conversely down-regulated neuron differentiation gene sets. We have therefore functionally validated rs7132908 as a causal obesity variant which temporally regulates nearby effector genes at the ch12q13 locus and influences neurodevelopment and survival.

Perturbation of the insomnia WDR90 GWAS locus pinpoints rs3752495 as a causal variant influencing distal expression of neighboring gene, PIG-Q

Although genome wide association studies (GWAS) have been crucial for the identification of loci associated with sleep traits and disorders, the method itself does not directly uncover the underlying causal variants and corresponding effector genes. The overwhelming majority of such variants reside in non-coding regions and are therefore presumed to impact the activity of cis-regulatory elements, such as enhancers. Our previously reported 'variant-to-gene mapping' effort in human induced pluripotent stem cell (iPSC)-derived neural progenitor cells (NPCs), combined with validation in both Drosophila and zebrafish, implicated PIG-Q as a functionally relevant gene at the insomnia 'WDR90' locus. However, importantly that effort did not characterize the corresponding underlying causal variant at this GWAS signal. Specifically, our genome-wide ATAC-seq and high-resolution promoter-focused Capture C datasets generated in this cell setting brought our attention to a shortlist of three tightly neighboring single nucleotide polymorphisms (SNPs) in strong linkage disequilibrium in a candidate intronic enhancer region of WDR90 that contacted the open PIG-Q promoter. The objective of this study was to investigate the influence of the proxy SNPs collectively and then individually on PIG-Q modulation and to pinpoint the causal "regulatory" variant among the three SNPs. Starting at a gross level perturbation, deletion of the entire region harboring all three SNPs in human iPSC-derived neural progenitor cells via CRISPR-Cas9 editing and subsequent RNA sequencing revealed expression changes in specific PIG-Q transcripts. Results from more refined individual luciferase reporter assays for each of the three SNPs in iPSCs revealed that the intronic region with the rs3752495 risk allele induced a ~2.5-fold increase in luciferase expression (n=10). Importantly, rs3752495 also exhibited an allele specific effect, with the risk allele increasing the luciferase expression by ~2-fold compared to the non-risk allele. In conclusion, our variant-to-function approach and subsequent in vitro validation implicates rs3752495 as a causal insomnia risk variant embedded at the WDR90-PIG-Q locus.

Differentiation block in acute myeloid leukemia regulated by intronic sequences of FTO

Iroquois transcription factor gene IRX3 is highly expressed in 20-30% of acute myeloid leukemia (AML) and contributes to the pathognomonic differentiation block. Intron 8 FTO sequences ∼220kB downstream of IRX3 exhibit histone acetylation, DNA methylation, and contacts with the IRX3 promoter, which correlate with IRX3 expression. Deletion of these intronic elements confirms a role in positively regulating IRX3. RNAseq revealed long non-coding (lnc) transcripts arising from this locus. FTO-lncAML knockdown (KD) induced differentiation of AML cells, loss of clonogenic activity, and reduced FTO intron 8:IRX3 promoter contacts. While both FTO-lncAML KD and IRX3 KD induced differentiation, FTO-lncAML but not IRX3 KD led to HOXA downregulation suggesting transcript activity in trans. FTO-lncAMLhigh AML samples expressed higher levels of HOXA and lower levels of differentiation genes. Thus, a regulatory module in FTO intron 8 consisting of clustered enhancer elements and a long non-coding RNA is active in human AML, impeding myeloid differentiation.

A lncRNA from the FTO locus acts as a suppressor of the m6A writer complex and p53 tumor suppression signaling

N6-methyladenosine (m6A) of mRNAs modulated by the METTL3-METTL14-WTAP-RBM15 methyltransferase complex and m6A demethylases such as FTO play important roles in regulating mRNA stability, splicing, and translation. Here, we demonstrate that FTO-IT1 long noncoding RNA (lncRNA) was upregulated and positively correlated with poor survival of patients with wild-type p53-expressing prostate cancer (PCa). m6A RIP-seq analysis revealed that FTO-IT1 knockout increased mRNA m6A methylation of a subset of p53 transcriptional target genes (e.g., FAS, TP53INP1, and SESN2) and induced PCa cell cycle arrest and apoptosis. We further showed that FTO-IT1 directly binds RBM15 and inhibits RBM15 binding, m6A methylation, and stability of p53 target mRNAs. Therapeutic depletion of FTO-IT1 restored mRNA m6A level and expression of p53 target genes and inhibited PCa growth in mice. Our study identifies FTO-IT1 lncRNA as a bona fide suppressor of the m6A methyltransferase complex and p53 tumor suppression signaling and nominates FTO-IT1 as a potential therapeutic target of cancer.

Genetic underpinning of the comorbidity between type 2 diabetes and osteoarthritis

Multimorbidity is a rising public health challenge with important implications for health management and policy. The most common multimorbidity pattern is the combination of cardiometabolic and osteoarticular diseases. Here, we study the genetic underpinning of the comorbidity between type 2 diabetes and osteoarthritis. We find genome-wide genetic correlation between the two diseases and robust evidence for association-signal colocalization at 18 genomic regions. We integrate multi-omics and functional information to resolve the colocalizing signals and identify high-confidence effector genes, including FTO and IRX3, which provide proof-of-concept insights into the epidemiologic link between obesity and both diseases. We find enrichment for lipid metabolism and skeletal formation pathways for signals underpinning the knee and hip osteoarthritis comorbidities with type 2 diabetes, respectively. Causal inference analysis identifies complex effects of tissue-specific gene expression on comorbidity outcomes. Our findings provide insights into the biological basis for the type 2 diabetes-osteoarthritis disease co-occurrence.

The rs1421085 variant within FTO promotes brown fat thermogenesis

One lead genetic risk signal of obesity-the rs1421085 T>C variant within the FTO gene-is reported to be functional in vitro but lacks evidence at an organism level. Here we recapitulate the homologous human variant in mice with global and brown adipocyte-specific variant knock-in and reveal that mice carrying the C-allele show increased brown fat thermogenic capacity and resistance to high-fat diet-induced adiposity, whereas the obesity-related phenotypic changes are blunted at thermoneutrality. Both in vivo and in vitro data reveal that the C-allele in brown adipocytes enhances the transcription of the Fto gene, which is associated with stronger chromatin looping linking the enhancer region and Fto promoter. Moreover, FTO knockdown or inhibition effectively eliminates the increased thermogenic ability of brown adipocytes carrying the C-allele. Taken together, these findings identify rs1421085 T>C as a functional variant promoting brown fat thermogenesis.

Breaking enhancers to gain insights into developmental defects

Despite ground-breaking genetic studies that have identified thousands of risk variants for developmental diseases, how these variants lead to molecular and cellular phenotypes remains a gap in knowledge. Many of these variants are non-coding and occur at enhancers, which orchestrate key regulatory programs during development. The prevailing paradigm is that non-coding variants alter the activity of enhancers, impacting gene expression programs, and ultimately contributing to disease risk. A key obstacle to progress is the systematic functional characterization of non-coding variants at scale, especially since enhancer activity is highly specific to cell type and developmental stage. Here, we review the foundational studies of enhancers in developmental disease and current genomic approaches to functionally characterize developmental enhancers and their variants at scale. In the coming decade, we anticipate systematic enhancer perturbation studies to link non-coding variants to molecular mechanisms, changes in cell state, and disease phenotypes.

Linkage disequilibrium block single-nucleotide polymorphisms in FTO alpha ketoglutarate dependent dioxygenase gene inference with breast cancer and Type II diabetes in Pakistani female population

Breast Cancer and Diabetes Mellitus are top ranked non-communicable life threatening diseases concerned in South Asia. The growing scientific clues witness the involvement of genetic variations which readily serve as risk factors for the disease onset. Comprehending the association of genetic predictors and risk factors in the conserved genome regions can help reveal underlying disease genetics and identify the sustained druggable targets. The present study aims to identify discrete inference of FTO alpha-ketoglutarate dependent dioxygenase gene linkage disequilibrium block SNPs rs9940128 and rs9939609 as prognostic genetic elements in defining the disease course either as BrC or NIDDM in Pakistani population. Clinically abreast female Breast Cancer and Type II (Non-Insulin Dependent) Diabetes Mellitus cases with the healthy controls participated in the study. The genomic study was established on the DNA of cases and controls through Tetra primers ARMS PCR and PCR-RFLP; data were analyzed statistically to reach comprehensive scientific annotation. Breast Cancer incidence was high in post menopause women. Fretful cholesterol, triglycerides, hypertension, sugar profiles and high incidence in females was evident in Type II (Non-Insulin Dependent) Diabetes Mellitus. BMI and family history were meager factor for either of the diseases. FTO gene alpha-ketoglutarate dependent dioxygenase linkage disequilibrium block Single-Nucleotide Polymorphism rs9939609 and rs9940128 threating inference was significant in the cancer and diabetes subjects correspondingly. The conclusion indicates serious clinical derailments in breast cancer and Type II (Non-Insulin Dependent) Diabetes Mellitus auxiliary to disease complication in genetically risk bearing FTO alpha-ketoglutarate dependent dioxygenase gene haplotype/linkage disequilibrium block SNPs prevailing in the affected Pakistani population. These clinical and genetic indicators could decisively be adopted in health care practice to intervene the sky rising disease incidence.

Systematic characterization of regulatory variants of blood pressure genes

High blood pressure (BP) is the major risk factor for cardiovascular disease. Genome-wide association studies have identified genetic variants for BP, but functional insights into causality and related molecular mechanisms lag behind. We functionally characterize 4,608 genetic variants in linkage with 135 BP loci in vascular smooth muscle cells and cardiomyocytes by massively parallel reporter assays. High densities of regulatory variants at BP loci (i.e., ULK4, MAP4, CFDP1, PDE5A) indicate that multiple variants drive genetic association. Regulatory variants are enriched in repeats, alter cardiovascular-related transcription factor motifs, and spatially converge with genes controlling specific cardiovascular pathways. Using heuristic scoring, we define likely causal variants, and CRISPR prime editing finally determines causal variants for KCNK9, SFXN2, and PCGF6, which are candidates for developing high BP. Our systems-level approach provides a catalog of functionally relevant variants and their genomic architecture in two trait-relevant cell lines for a better understanding of BP gene regulation.