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

Fat mass and obesity associated gene and homeobox transcription factor iriquois-3 mRNA profiles in the metabolic tissues of zebrafish are modulated by feeding and food deprivation

Fat mass and obesity associated gene (FTO) has been strongly associated with obesity, and it is functionally linked to the homeobox transcription factor iriquois-3 (IRX3). In mammals, FTO and IRX3 are involved in the regulation of food intake and metabolism. This study aimed to determine whether FTO and IRX3are affected by feeding and food unavailability. FTO and IRX3 mRNA and protein were found widely distributed in all tissues examined, including the brain, muscle, gut, and liver. Postprandial increase in the abundance of FTO and IRX3 mRNAs was observed in metabolic tissues of both male and female zebrafish at 1 h post-feeding. Meanwhile, their expression in the brain and gut decreased at 3 h post-feeding, reaching preprandial levels. Additionally, FTO and IRX3 mRNA abundance in examined tissues increased after 7 days of food deprivation, but substantially decreased after refeeding for 24 h. In summary, we report that both FTO and IRX3 are meal-sensitive genes in zebrafish. The fasting-induced increase suggests a possible appetite regulatory role for FTO and IRX3 in zebrafish. These findings highlight the importance of FTO and IRX3 in appetite and metabolic regulation in zebrafish.

Depot-specific mRNA expression programs in human adipocytes suggest physiological specialization via distinct developmental programs

Adipose tissue is distributed in diverse locations throughout the human body. Not much is known about the extent to which anatomically distinct adipose depots are functionally distinct, specialized organs, nor whether depot-specific characteristics result from intrinsic developmental programs, as opposed to reversible physiological responses to differences in tissue microenvironment. We used DNA microarrays to compare mRNA expression patterns of isolated human adipocytes and cultured adipose stem cells, before and after ex vivo adipocyte differentiation, from seven anatomically diverse adipose tissue depots. Adipocytes from different depots display distinct gene expression programs, which are most closely shared with anatomically related depots. mRNAs whose expression differs between anatomically diverse groups of depots (e.g., subcutaneous vs. internal) suggest important functional specializations. These depot-specific differences in gene expression were recapitulated when adipocyte progenitor cells from each site were differentiated ex vivo, suggesting that progenitor cells from specific anatomic sites are deterministically programmed to differentiate into depot-specific adipocytes. Many developmental transcription factors show striking depot-specific patterns of expression, suggesting that adipocytes in each anatomic depot are programmed during early development in concert with anatomically related tissues and organs. Our results support the hypothesis that adipocytes from different depots are functionally distinct and that their depot-specific specialization reflects distinct developmental programs.

A novel genetic mouse model of osteoporosis with double heterozygosity of Irx3 and Irx5 characterizes sex-dependent phenotypes in bone homeostasis

Iroquois homeobox gene 3 (Irx3) and Irx5 encode transcription factors that play crucial roles in limb development and bone formation. Previous studies using knockout mice have revealed a role of Irx3 and Irx5 in osteogenesis in young adult mice. However, whether these genes are also essential for bone homeostasis in adulthood and contribute to bone diseases remain poorly understood. Osteoporosis is a disease characterized by lower bone mineral density and disrupted bone microarchitecture, typically occurs in postmenopausal women. Here, we demonstrate that Irx3/5dHet mice with a half-reduction of Irx3 and Irx5 dosage serve as a novel model of osteoporosis. By micro-computed tomography, we found that Irx3/5dHet mice exhibited sex-dependent bone loss patterns. While male Irx3/5dHet mice progressively lost trabecular microstructures with aging, female mutants exhibited lower bone mineral density (BMD) and bone volume fraction (BV/TV) at early adulthood (9-15 weeks old) but without further loss later at 1 year of age. Bone marrow adipocytes are known to be elevated at the expenses of lower osteogenesis in osteoporotic bone marrow. Surprisingly, we found sex-dependent changes in adipogenesis at the age of skeletal maturity that bone marrow adipocytes were reduced in female Irx3/5dHet mice along with deteriorated osteogenesis, while male mice exhibited elevated adipogenesis. In summary, we reported a novel genetic model for osteoporosis-like phenotypes, highlighting sex-dependent bone mineral density and bone marrow adipocyte characteristics.

Melanoma-specific mutation hotspots in distal, non-coding, promoter-interacting regions implicate novel candidate driver genes

BACKGROUND

To develop targeted treatments, it is crucial to identify the full spectrum of genetic drivers in melanoma, including those in non-coding regions. However, recent efforts to explore non-coding regions have primarily focused on gene-adjacent elements such as promoters and non-coding RNAs, leaving intergenic distal regulatory elements largely unexplored.

METHODS

We used Hi-C chromatin contact data from melanoma cells to map distal, non-coding, promoter-interacting regulatory elements genome-wide in melanoma. Using this "promoter-interaction network", alongside whole-genome sequence and gene expression data from the Pan Cancer Analysis of Whole Genomes, we developed multivariate linear regression models to identify distal somatic mutation hotspots that affect promoter activity.

RESULTS

We identified eight recurrently mutated hotspots that are novel, melanoma-specific, located in promoter-interacting distal regulatory elements, alter transcription factor binding motifs, and affect the expression of genes (e.g., HSPB7, CLDN1, ADCY9 and FDXR) previously implicated as tumour suppressors/oncogenes in various cancers.

CONCLUSIONS

Our study suggests additional non-coding drivers beyond the well-characterised TERT promoter in melanoma, offering new insights into the disruption of complex regulatory networks by non-coding mutations that may contribute to melanoma development. Furthermore, our study provides a framework for integrating multiple levels of biological data to uncover cancer-specific non-coding drivers.

Translation of genome-wide association study: from genomic signals to biological insights

Since the turn of the 21st century, genome-wide association study (GWAS) have successfully identified genetic signals associated with a myriad of common complex traits and diseases. As we transition from establishing robust genetic associations with diverse phenotypes, the central challenge is now focused on characterizing the underlying functional mechanisms driving these signals. Previous GWAS efforts have revealed multiple variants, each conferring relatively subtle susceptibility, collectively contributing to the pathogenesis of various common diseases. Such variants can further exhibit associations with multiple other traits and differ across ancestries, plus disentangling causal variants from non-causal due to linkage disequilibrium complexities can lead to challenges in drawing direct biological conclusions. Combined with cellular context considerations, such challenges can reduce the capacity to definitively elucidate the biological significance of GWAS signals, limiting the potential to define mechanistic insights. This review will detail current and anticipated approaches for functional interpretation of GWAS signals, both in terms of characterizing the underlying causal variants and the corresponding effector genes.

Powerful mapping of cis-genetic effects on gene expression across diverse populations reveals novel disease-critical genes

While disease-associated variants identified by genome-wide association studies (GWAS) most likely regulate gene expression levels, linking variants to target genes is critical to determining the functional mechanisms of these variants. Genetic effects on gene expression have been extensively characterized by expression quantitative trait loci (eQTL) studies, yet data from non-European populations is limited. This restricts our understanding of disease to genes whose regulatory variants are common in European populations. While previous work has leveraged data from multiple populations to improve GWAS power and polygenic risk score (PRS) accuracy, multi-ancestry data has not yet been used to better estimate cis-genetic effects on gene expression. Here, we present a new method, Multi-Ancestry Gene Expression Prediction Regularized Optimization (MAGEPRO), which constructs robust genetic models of gene expression in understudied populations or cell types by fitting a regularized linear combination of eQTL summary data across diverse cohorts. In simulations, our tool generates more accurate models of gene expression than widely-used LASSO and the state-of-the-art multi-ancestry PRS method, PRS-CSx, adapted to gene expression prediction. We attribute this improvement to MAGEPRO's ability to more accurately estimate causal eQTL effect sizes (p < 3.98 × 10-4, two-sided paired t-test). With real data, we applied MAGEPRO to 8 eQTL cohorts representing 3 ancestries (average n = 355) and consistently outperformed each of 6 competing methods in gene expression prediction tasks. Integration with GWAS summary statistics across 66 complex traits (representing 22 phenotypes and 3 ancestries) resulted in 2,331 new gene-trait associations, many of which replicate across multiple ancestries, including PHTF1 linked to white blood cell count, a gene which is overexpressed in leukemia patients. MAGEPRO also identified biologically plausible novel findings, such as PIGB, an essential component of GPI biosynthesis, associated with heart failure, which has been previously evidenced by clinical outcome data. Overall, MAGEPRO is a powerful tool to enhance inference of gene regulatory effects in underpowered datasets and has improved our understanding of population-specific and shared genetic effects on complex traits.

Cross-species evidence for a developmental origin of adult hypersomnia with loss of synaptic adhesion molecules beat-Ia/CADM2

Idiopathic hypersomnia (IH) is a poorly-understood sleep disorder characterized by excessive daytime sleepiness despite normal nighttime sleep. Combining human genomics with behavioral and mechanistic studies in fish and flies, we uncover a role for beat-Ia/CADM2 , synaptic adhesion molecules of the immunoglobulin superfamily, in excessive sleepiness. Neuronal knockdown of Drosophila beat-Ia results in sleepy flies and loss of the vertebrate ortholog of beat-Ia , CADM2 , results in sleepy fish. We delineate a developmental function for beat-Ia in synaptic elaboration of neuropeptide F (NPF) neurites projecting to the suboesophageal zone (SEZ) of the fly brain. Brain connectome and experimental evidence demonstrate these NPF outputs synapse onto a subpopulation of SEZ GABAergic neurons to stabilize arousal. NPF is the Drosophila homolog of vertebrate neuropeptide Y (NPY), and an NPY receptor agonist restores sleep to normal levels in zebrafish lacking CADM2 . These findings point towards NPY modulation as a treatment target for human hypersomnia.

Dealkylation of Macromolecules by Eukaryotic α-Ketoglutarate-Dependent Dioxygenases from the AlkB-like Family

Alkylating modifications induced by either exogenous chemical agents or endogenous metabolites are some of the main types of damage to DNA, RNA, and proteins in the cell. Although research in recent decades has been almost entirely devoted to the repair of alkyl and in particular methyl DNA damage, more and more data lately suggest that the methylation of RNA bases plays an equally important role in normal functioning and in the development of diseases. Among the most prominent participants in the repair of methylation-induced DNA and RNA damage are human homologs of Escherichia coli AlkB, nonheme Fe(II)/α-ketoglutarate-dependent dioxygenases ABH1-8, and FTO. Moreover, some of these enzymes have been found to act on several protein targets. In this review, we present up-to-date data on specific features of protein structure, substrate specificity, known roles in the organism, and consequences of disfunction of each of the nine human homologs of AlkB. Special attention is given to reports about the effects of natural single-nucleotide polymorphisms on the activity of these enzymes and to potential consequences for carriers of such natural variants.

Atp1a2 and Kcnj9 are candidate genes underlying sensitivity to oxycodone-induced locomotor activation and withdrawal-induced anxiety-like behaviors in C57BL/6 substrains

Opioid use disorder is heritable, yet its genetic etiology is largely unknown. C57BL/6J and C57BL/6NJ mouse substrains exhibit phenotypic diversity in the context of limited genetic diversity which together can facilitate genetic discovery. Here, we found C57BL/6NJ mice were less sensitive to oxycodone (OXY)-induced locomotor activation versus C57BL/6J mice in a conditioned place preference paradigm. Narrow-sense heritability was estimated at 0.22-0.31, implicating suitability for genetic analysis. Quantitative trait locus (QTL) mapping in an F2 cross identified a chromosome 1 QTL explaining 7-12% of the variance in OXY locomotion and anxiety-like withdrawal in the elevated plus maze. A second QTL for EPM withdrawal behavior on chromosome 5 near Gabra2 (alpha-2 subunit of GABA-A receptor) explained 9% of the variance. To narrow the chromosome 1 locus, we generated recombinant lines spanning 163-181 Mb, captured the QTL for OXY locomotor traits and withdrawal, and fine-mapped a 2.45-Mb region (170.16-172.61 Mb). Transcriptome analysis identified five, localized striatal cis-eQTL transcripts and two were confirmed at the protein level (KCNJ9, ATP1A2). Kcnj9 codes for a potassium channel (GIRK3) 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 functional adaptations following chronic opioid administration. To summarize, we identified two candidate genes underlying the physiological and behavioral properties of opioids, with direct preclinical relevance to investigators employing these widely used substrains and clinical relevance to human genetic studies of opioid use disorder.

Genotype inference from aggregated chromatin accessibility data reveals genetic regulatory mechanisms

Background

Understanding the genetic causes for variability in chromatin accessibility can shed light on the molecular mechanisms through which genetic variants may affect complex traits. Thousands of ATAC-seq samples have been collected that hold information about chromatin accessibility across diverse cell types and contexts, but most of these are not paired with genetic information and come from diverse distinct projects and laboratories.

Results

We report here joint genotyping, chromatin accessibility peak calling, and discovery of quantitative trait loci which influence chromatin accessibility (caQTLs), demonstrating the capability of performing caQTL analysis on a large scale in a diverse sample set without pre-existing genotype information. Using 10,293 profiling samples representing 1,454 unique donor individuals across 653 studies from public databases, we catalog 23,381 caQTLs in total. After joint discovery analysis, we cluster samples based on accessible chromatin profiles to identify context-specific caQTLs. We find that caQTLs are strongly enriched for annotations of gene regulatory elements across diverse cell types and tissues and are often strongly linked with genetic variation associated with changes in expression (eQTLs), indicating that caQTLs can mediate genetic effects on gene expression. We demonstrate sharing of causal variants for chromatin accessibility and diverse complex human traits, enabling a more complete picture of the genetic mechanisms underlying complex human phenotypes.

Conclusions

Our work provides a proof of principle for caQTL calling from previously ungenotyped samples, and represents one of the largest, most diverse caQTL resources currently available, informing mechanisms of genetic regulation of gene expression and contribution to disease.

Deciphering the impact of genomic variation on function

Our genomes influence nearly every aspect of human biology-from molecular and cellular functions to phenotypes in health and disease. Studying the differences in DNA sequence between individuals (genomic variation) could reveal previously unknown mechanisms of human biology, uncover the basis of genetic predispositions to diseases, and guide the development of new diagnostic tools and therapeutic agents. Yet, understanding how genomic variation alters genome function to influence phenotype has proved challenging. To unlock these insights, we need a systematic and comprehensive catalogue of genome function and the molecular and cellular effects of genomic variants. Towards this goal, the Impact of Genomic Variation on Function (IGVF) Consortium will combine approaches in single-cell mapping, genomic perturbations and predictive modelling to investigate the relationships among genomic variation, genome function and phenotypes. IGVF will create maps across hundreds of cell types and states describing how coding variants alter protein activity, how noncoding variants change the regulation of gene expression, and how such effects connect through gene-regulatory and protein-interaction networks. These experimental data, computational predictions and accompanying standards and pipelines will be integrated into an open resource that will catalyse community efforts to explore how our genomes influence biology and disease across populations.

The complex web of obesity: from genetics to precision medicine

INTRODUCTION

Obesity is a growing public health concern affecting both children and adults. Since it involves both genetic and environmental components, the management of obesity requires both, an understanding of the underlying genetics and changes in lifestyle. The knowledge of obesity genetics will enable the possibility of precision medicine in anti-obesity medications.

AREAS COVERED

Here, we explore health complications and the prevalence of obesity. We discuss disruptions in energy balance as a symptom of obesity, examining evolutionary theories, its multi-factorial origins, and heritability. Additionally, we discuss monogenic and polygenic obesity, the converging biological pathways, potential pharmacogenomics applications, and existing anti-obesity medications - specifically focussing on the leptin-melanocortin and incretin pathways. Comparisons between childhood and adult obesity genetics are made, along with insights into structural variants, epigenetic changes, and environmental influences on epigenetic signatures.

EXPERT OPINION

With recent advancements in anti-obesity drugs, genetic studies pinpoint new targets and allow for repurposing existing drugs. This creates opportunities for genotype-informed treatment options. Also, lifestyle interventions can help in the prevention and treatment of obesity by altering the epigenetic signatures. The comparison of genetic architecture in adults and children revealed a significant overlap. However, more robust studies with diverse ethnic representation is required in childhood obesity.

Variant-to-function mapping of late-onset Alzheimer's disease GWAS signals in human microglial cell models implicates RTFDC1 at the CASS4 locus

Late-onset Alzheimer's disease (LOAD) research has principally focused on neurons over the years due to their known role in the production of amyloid beta plaques and neurofibrillary tangles. In contrast, recent genomic studies of LOAD have implicated microglia as culprits of the prolonged inflammation exacerbating the neurodegeneration observed in patient brains. Indeed, recent LOAD genome-wide association studies (GWAS) have reported multiple loci near genes related to microglial function, including TREM2, ABI3, and CR1. However, GWAS alone cannot pinpoint underlying causal variants or effector genes at such loci, as most signals reside in non-coding regions of the genome and could presumably confer their influence frequently via long-range regulatory interactions. We elected to carry out a combination of ATAC-seq and high-resolution promoter-focused Capture-C in two human microglial cell models (iPSC-derived microglia and HMC3) in order to physically map interactions between LOAD GWAS-implicated candidate causal variants and their corresponding putative effector genes. Notably, we observed consistent evidence that rs6024870 at the GWAS CASS4 locus contacted the promoter of nearby gene, RTFDC1. We subsequently observed a directionallly consistent decrease in RTFDC1 expression with the the protective minor A allele of rs6024870 via both luciferase assays in HMC3 cells and expression studies in primary human microglia. Through CRISPR-Cas9-mediated deletion of the putative regulatory region harboring rs6024870 in HMC3 cells, we observed increased pro-inflammatory cytokine secretion and decreased DNA double strand break repair related, at least in part, to RTFDC1 expression levels. Our variant-to-function approach therefore reveals that the rs6024870-harboring regulatory element at the LOAD 'CASS4' GWAS locus influences both microglial inflammatory capacity and DNA damage resolution, along with cumulative evidence implicating RTFDC1 as a novel candidate effector gene.

Targeted inhibition of m6A demethylase FTO by FB23 attenuates allergic inflammation in the airway epithelium

Epithelial cells play a crucial role in asthma, contributing to chronic inflammation and airway hyperresponsiveness. m6A modification, which involves key proteins such as the demethylase fat mass and obesity-associated protein (FTO), is crucial in the regulation of various diseases, including asthma. However, the role of FTO in epithelial cells and the development of asthma remains unclear. In this study, we investigated the demethylase activity of FTO using a small-molecule inhibitor FB23 in epithelial cells and allergic inflammation in vivo and in vitro. We examined the FTO-regulated transcriptome-wide m6A profiling by methylated RNA immunoprecipitation sequencing (MeRIP-seq) and RNA-seq under FB23 treatment and allergic inflammation conditions. Immunofluorescence staining was performed to assess the tissue-specific expression of FTO in asthmatic bronchial mucosa. We demonstrated that FB23 alleviated allergic inflammation in IL-4/IL-13-treated epithelial cells and house dust mite (HDM)-induced allergic airway inflammation mouse model. The demethylase activity of FTO contributed to the regulation of TNF-α signaling via NF-κB and epithelial-mesenchymal transition-related pathways under allergic inflammation conditions in epithelial cells. FTO was expressed in epithelial, submucosal gland, and smooth muscle cells in human bronchial mucosa. In conclusion, FB23-induced inhibition of FTO alleviates allergic inflammation in epithelial cells and HDM-induced mice, potentially through diverse cellular processes and epithelial-mesenchymal transition signaling pathways, suggesting that FTO is a potential therapeutic target in asthma management.

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.

Insight into the regulatory mechanism of m6A modification: From MAFLD to hepatocellular carcinoma

In recent years, there has been a significant increase in the incidence of metabolic-associated fatty liver disease (MAFLD), which has been attributed to the increasing prevalence of type 2 diabetes mellitus (T2DM) and obesity. MAFLD affects more than one-third of adults worldwide, making it the most prevalent liver disease globally. Moreover, MAFLD is considered a significant risk factor for hepatocellular carcinoma (HCC), with MAFLD-related HCC cases increasing. Approximately 1 in 6 HCC patients are believed to have MAFLD, and nearly 40 % of these HCC patients do not progress to cirrhosis, indicating direct transformation from MAFLD to HCC. N6-methyladenosine (m6A) is commonly distributed in eukaryotic mRNA and plays a crucial role in normal development and disease progression, particularly in tumors. Numerous studies have highlighted the close association between abnormal m6A modification and cellular metabolic alterations, underscoring its importance in the onset and progression of MAFLD. However, the specific impact of m6A modification on the progression of MAFLD to HCC remains unclear. Can targeting m6A effectively halt the progression of MAFLD-related HCC? In this review, we investigated the pivotal role of abnormal m6A modification in the transition from MAFLD to HCC, explored the potential of m6A modification as a therapeutic target for MAFLD-related HCC, and proposed possible directions for future investigations.

Gapped-kmer sequence modeling robustly identifies regulatory vocabularies and distal enhancers conserved between evolutionarily distant mammals

Gene regulatory elements drive complex biological phenomena and their mutations are associated with common human diseases. The impacts of human regulatory variants are often tested using model organisms such as mice. However, mapping human enhancers to conserved elements in mice remains a challenge, due to both rapid enhancer evolution and limitations of current computational methods. We analyze distal enhancers across 45 matched human/mouse cell/tissue pairs from a comprehensive dataset of DNase-seq experiments, and show that while cell-specific regulatory vocabulary is conserved, enhancers evolve more rapidly than promoters and CTCF binding sites. Enhancer conservation rates vary across cell types, in part explainable by tissue specific transposable element activity. We present an improved genome alignment algorithm using gapped-kmer features, called gkm-align, and make genome wide predictions for 1,401,803 orthologous regulatory elements. We show that gkm-align discovers 23,660 novel human/mouse conserved enhancers missed by previous algorithms, with strong evidence of conserved functional activity.

Nucleic acid demethylase MpAlkB1 regulates the growth, development, and secondary metabolite biosynthesis in Monascus purpureus

Nucleic acid demethylases of α-ketoglutarate-dependent dioxygenase (AlkB) family can reversibly erase methyl adducts from nucleobases, thus dynamically regulating the methylation status of DNA/RNA and playing critical roles in multiple cellular processes. But little is known about AlkB demethylases in filamentous fungi so far. The present study reports that Monascus purpureus genomes contain a total of five MpAlkB genes. The MpAlkB1 gene was disrupted and complemented through homologous recombination strategy to analyze its biological functions in M. purpureus. MpAlkB1 knockout significantly accelerated the growth of strain, increased biomass, promoted sporulation and cleistothecia development, reduced the content of Monascus pigments (Mps), and strongly inhibited citrinin biosynthesis. The downregulated expression of the global regulator gene LaeA, and genes of Mps biosynthesis gene cluster (BGC) or citrinin BGC in MpAlkB1 disruption strain supported the pleiotropic trait changes caused by MpAlkB1 deletion. These results indicate that MpAlkB1-mediated demethylation of nucleic acid plays important roles in regulating the growth and development, and secondary metabolism in Monascus spp.

Novel loci linked to serum lipid traits are identified in a genome-wide association study of a highly admixed Brazilian population - the 2015 ISA Nutrition

BACKGROUND

Cardiovascular diseases (CVDs) comprise major causes of death worldwide, leading to extensive burden on populations and societies. Alterations in normal lipid profiles, i.e., dyslipidemia, comprise important risk factors for CVDs. However, there is lack of comprehensive evidence on the genetic contribution to dyslipidemia in highly admixed populations. The identification of single nucleotide polymorphisms (SNPs) linked to blood lipid traits in the Brazilian population was based on genome-wide associations using data from the São Paulo Health Survey with Focus on Nutrition (ISA-Nutrition).

METHODS

A total of 667 unrelated individuals had genetic information on 330,656 SNPs available, and were genotyped with Axiom™ 2.0 Precision Medicine Research Array. Genetic associations were tested at the 10- 5 significance level for the following phenotypes: low-density lipoprotein cholesterol (LDL-c), very low-density lipoprotein cholesterol (VLDL-c), high-density lipoprotein cholesterol (HDL-c), HDL-c/LDL-c ratio, triglycerides (TGL), total cholesterol, and non-HDL-c.

RESULTS

There were 19 significantly different SNPs associated with lipid traits, the majority of which corresponding to intron variants, especially in the genes FAM81A, ZFHX3, PTPRD, and POMC. Three variants (rs1562012, rs16972039, and rs73401081) and two variants (rs8025871 and rs2161683) were associated with two and three phenotypes, respectively. Among the subtypes, non-HDL-c had the highest proportion of associated variants.

CONCLUSIONS

The results of the present genome-wide association study offer new insights into the genetic structure underlying lipid traits in underrepresented populations with high ancestry admixture. The associations were robust across multiple lipid phenotypes, and some of the phenotypes were associated with two or three variants. In addition, some variants were present in genes that encode ncRNAs, raising important questions regarding their role in lipid metabolism.

CREdb: A comprehensive database of Cis-Regulatory Elements and their activity in human cells and tissues

BACKGROUND

Cis-regulatory elements (CREs) play a pivotal role in gene expression regulation, allowing cells to serve diverse functions and respond to external stimuli. Understanding CREs is essential for personalized medicine and disease research, as an increasing number of genetic variants associated with phenotypes and diseases overlap with CREs. However, existing databases often focus on subsets of regulatory elements and present each identified instance of element individually, confounding the effort to obtain a comprehensive view. To address this gap, we have created CREdb, a comprehensive database with over 10 million human regulatory elements across 1,058 cell types and 315 tissues harmonized from different data sources. We curated and aligned the cell types and tissues to standard ontologies for efficient data query.

RESULTS

Data from 11 sources were curated and mapped to standard ontological terms. 11,223,434 combined elements are present in the final database, and these were merged into 5,666,240 consensus elements representing the combined ranges of the individual elements informed by their overlap. Each consensus element contains curated metadata including the number of elements supporting it and a hash linking to the source databases. The inferred activity of each consensus element in various cell-type and tissue context is also provided. Examples presented here show the potential utility of CREdb in annotating non-coding genetic variants and informing chromatin accessibility profiling analysis.

CONCLUSIONS

We developed CREdb, a comprehensive database of CREs, to simplify the analysis of CREs by providing a unified framework for researchers. CREdb compiles consensus ranges for each element by integrating the information from all instances identified across various source databases. This unified database facilitates the functional annotation of non-coding genetic variants and complements chromatin accessibility profiling analysis. CREdb will serve as an important resource in expanding our knowledge of the epigenome and its role in human diseases.

CoPheScan: phenome-wide association studies accounting for linkage disequilibrium

Phenome-wide association studies (PheWAS) facilitate the discovery of associations between a single genetic variant with multiple phenotypes. For variants which impact a specific protein, this can help identify additional therapeutic indications or on-target side effects of intervening on that protein. However, PheWAS is restricted by an inability to distinguish confounding due to linkage disequilibrium (LD) from true pleiotropy. Here we describe CoPheScan (Coloc adapted Phenome-wide Scan), a Bayesian approach that enables an intuitive and systematic exploration of causal associations while simultaneously addressing LD confounding. We demonstrate its performance through simulation, showing considerably better control of false positive rates than a conventional approach not accounting for LD. We used CoPheScan to perform PheWAS of protein-truncating variants and fine-mapped variants from disease and pQTL studies, in 2275 disease phenotypes from the UK Biobank. Our results identify the complexity of known pleiotropic genes such as APOE, and suggest a new causal role for TGM3 in skin cancer.

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

Although genome-wide association studies (GWAS) have identified loci for sleep-related traits, they do not directly uncover the underlying causal variants and corresponding effector genes. The majority of such variants reside in non-coding regions and are therefore presumed to impact cis-regulatory elements. 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 phosphatidyl inositol glycan (PIG)-Q as a functionally relevant gene at the insomnia "WDR90" GWAS locus. However, importantly that effort did not characterize the corresponding underlying causal variant. Specifically, our previous 3D genomic datasets nominated a shortlist of three neighboring single nucleotide polymorphisms (SNPs) in strong linkage disequilibrium within an intronic enhancer region of WDR90 that contacted the open PIG-Q promoter. We sought to investigate the influence of these SNPs collectively and then individually on PIG-Q modulation to pinpoint the causal "regulatory" variant. Starting with gross level perturbation, deletion of the entire region in NPCs via CRISPR-Cas9 editing and subsequent RNA sequencing revealed expression changes in specific PIG-Q transcripts. Results from individual luciferase reporter assays for each SNP in iPSCs revealed that the region with the rs3752495 risk allele (RA) induced a ~2.5-fold increase in luciferase expression. Importantly, rs3752495 also exhibited an allele-specific effect, with the RA increasing the luciferase expression by ~2-fold versus the non-RA. In conclusion, our variant-to-function approach and in vitro validation implicate rs3752495 as a causal insomnia variant embedded within WDR90 while modulating the expression of the distally located PIG-Q.

Partitioning and aggregating cross-tissue and tissue-specific genetic effects to identify gene-trait associations

TWAS have shown great promise in extending GWAS loci to a functional understanding of disease mechanisms. In an effort to fully unleash the TWAS and GWAS information, we propose MTWAS, a statistical framework that partitions and aggregates cross-tissue and tissue-specific genetic effects in identifying gene-trait associations. We introduce a non-parametric imputation strategy to augment the inaccessible tissues, accommodating complex interactions and non-linear expression data structures across various tissues. We further classify eQTLs into cross-tissue eQTLs and tissue-specific eQTLs via a stepwise procedure based on the extended Bayesian information criterion, which is consistent under high-dimensional settings. We show that MTWAS significantly improves the prediction accuracy across all 47 tissues of the GTEx dataset, compared with other single-tissue and multi-tissue methods, such as PrediXcan, TIGAR, and UTMOST. Applying MTWAS to the DICE and OneK1K datasets with bulk and single-cell RNA sequencing data on immune cell types showcases consistent improvements in prediction accuracy. MTWAS also identifies more predictable genes, and the improvement can be replicated with independent studies. We apply MTWAS to 84 UK Biobank GWAS studies, which provides insights into disease etiology.

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 potential 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-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.

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.

Genetic Variants in the Fat Mass and Obesity-Associated Gene and Risk of Obesity/Overweight in Children and Adolescents: A Systematic Review and Meta-Analysis

OBJECTIVE

The variations in the single-nucleotide polymorphisms (SNPs) of the fat mass and obesity (FTO)-associated gene have been linked to being overweight or obese in children. In this research a thorough examination was performed to elucidate the connection between various FTO gene SNPs and overweight or obesity in children and adolescents.

METHOD

We searched PubMed, Google scholar, Web of Science and Scopus until January 2024 to find studies that investigate the association between different SNPs of FTO gene and the risk of overweight/obesity in children and adolescents. After filtering the relevant studies, meta-analysis was used to quantify the association of FTO gene SNPs within different genetic inheritance models.

RESULTS

We have identified 32 eligible studies with 14,930 obese/overweight cases and 24,765 healthy controls. Our recessive model showed a significant association with rs9939609 (OR: 1.56, 95% CI: 1.20; 2.02, p < 0.01) and rs1421085 (OR: 1.77, 95% CI: 1.14; 2.75, p < 0.01). Besides, in the homozygote model, rs1421085 showed the highest association (OR: 2.32, 95% CI: 1.38; 3.89, p < 0.01) with the risk of obesity in a population of children and adolescents. Moreover, there are other SNPs of FTO genes, such as rs9921255, rs9928094 and rs9930333, which showed a positive association with obesity and overweight. However, their effects were evaluated in very few numbers of studies.

CONCLUSION

In this study, we have found that the FTO rs9939609 and rs1421085 are associated to an increased risk of obesity among children and adolescents. Besides, the findings of this study further reaffirmed the established link between rs9939609 and obesity in children and adolescents.

Conserved Noncoding Cis-Elements Associated with Hibernation Modulate Metabolic and Behavioral Adaptations in Mice

Our study elucidates functional roles for conserved cis-elements associated with the evolution of mammalian hibernation. Genomic analyses found topologically associated domains (TADs) that disproportionately accumulated convergent genomic changes in hibernators, including the TAD for the Fat Mass & Obesity (Fto) locus. Some hibernation-linked cis-elements in this TAD form regulatory contacts with multiple neighboring genes. Knockout mice for these cis-elements exhibit Fto, Irx3, and Irx5 gene expression changes, impacting hundreds of genes downstream. Profiles of pre-torpor, torpor, and post-torpor phenotypes found distinct roles for each cis-element in metabolic control, while a high caloric diet uncovered different obesogenic effects. One cis-element promoting a lean phenotype influences foraging behaviors throughout life, affecting specific behavioral sequences. Thus, convergent evolution in hibernators pinpoints functional genetic mechanisms of mammalian metabolic control.

Genetic association of FTO gene polymorphisms with obesity and its related phenotypes: A case-control study

Introduction

FTO gene belongs to the non-heme Fe (II) and 2 oxoglutarate-dependent dioxygenase superfamily. Polymorphisms within the first intron of the FTO gene have been examined across various populations, yielding disparate findings.The present study aimed to determine the impact of two intronic polymorphisms FTO 30685T/G (rs17817449) and -23525T/A (rs9939609) on the risk of obesity in Punjab, India.

Methods

Genotypic and biochemical analysis were done for 671 unrelated participants (obese=333 and non-obese=338) (age≥18 years). Genotyping of the polymorphisms was done by PCR-RFLP method. However, 50% of the samples were sequenced by Sanger sequencing.

Results

Both the FTO variants 30685 (TT vs GG: odds ratio (OR), 2.30; 95% confidence interval (CI), 1.39-3.79) and -23525 (TT vs AA: odds ratio (OR), 2.78; 95% confidence interval (CI), 1.37-5.64) showed substantial risk towards obesity by conferring it 2 times and 3 times, respectively. The analysis by logistic regression showed a significant association for both the variants 30685T/G (rs17817449) and -23525T/A (rs9939609) (OR=2.29; 95%CI: 1.47-3.57) and (OR=5.25; 95% CI: 2.68-10.28) under the recessive genetic model, respectively. The haplotype combination TA (30685; -23525) develops a 4 times risk for obesity (P=0.0001). Among obese, the G allele of 30685T/G and A- allele of -23525T/A showed variance in Body mass index (BMI), waist circumference (WC), waist-to-height ratio(WHtR), systolic blood pressure (SBP), diastolic blood pressure (DBP) and triglyceride(TG).

Conclusion

The present investigation indicated that both the FTO 30685T/G (rs17817449) and -23525T/A (rs9939609) polymorphisms have a key impact on an individual's vulnerability to obesity in this population.

Integrating Genetic and Transcriptomic Data to Identify Genes Underlying Obesity Risk Loci

Genome-wide association studies (GWAS) have identified numerous body mass index (BMI) loci. However, most underlying mechanisms from risk locus to BMI remain unknown. Leveraging omics data through integrative analyses could provide more comprehensive views of biological pathways on BMI. We analyzed genotype and blood gene expression data in up to 5,619 samples from the Framingham Heart Study (FHS). Using 3,992 single nucleotide polymorphisms (SNPs) at 97 BMI loci and 20,692 transcripts within 1 Mb, we performed separate association analyses of transcript with BMI and SNP with transcript (PBMI and PSNP, respectively) and then a correlated meta-analysis between the full summary data sets (PMETA). We identified transcripts that met Bonferroni-corrected significance for each omic, were more significant in the correlated meta-analysis than each omic, and were at least nominally associated with BMI in FHS data. Among 308 significant SNP-transcript-BMI associations, we identified seven genes (NT5C2, GSTM3, SNAPC3, SPNS1, TMEM245, YPEL3, and ZNF646) in five association regions. Using an independent sample of blood gene expression data, we validated results for SNAPC3 and YPEL3. We tested for generalization of these associations in hypothalamus, nucleus accumbens, and liver and observed significant (PMETA<0.05 & PMETA Collapse

Key Words Collapse

MESH Headings Collapse

Grants

• R01 HL143885 NHLBI NIH HHS
• R01 HL142302 NHLBI NIH HHS
• R01 HL151152 NHLBI NIH HHS
• HHSN268201500001I NHLBI NIH HHS
• N01HC25195 NHLBI NIH HHS
• R01 HG010297 NHGRI NIH HHS
• HHSN268201500001C NHLBI NIH HHS
• R01 HL163262 NHLBI NIH HHS
• R01 HL092577 NHLBI NIH HHS
• P20 MD000170 NIMHD NIH HHS
• 75N92019D00031 NHLBI NIH HHS
• R01 HD057194 NICHD NIH HHS
• UL1 TR000371 NCATS NIH HHS
• R01 DK122503 NIDDK NIH HHS

Collapse

Affiliation(s)

Hanfei Xu Department of Biostatistics, School of Public Health, Boston University, 801 Massachusettes Ave, Boston, MA, 02118, USA
Shreyash Gupta Department of Population Health Sciences, Geisinger, 100 N. Academy Ave., Danville, PA, 17822, USA
Ian Dinsmore Department of Genomic Health, Geisinger, 100 N. Academy Ave., Danville, PA, 17822, USA
Abbey Kollu Department of Psychology and Neuroscience, University of North Carolina, 235 E. Cameron Avenue, Chapel Hill, NC, 27599, USA
Anne Marie Cawley Marsico Lung Institute, University of North Carolina, 125 Mason Farm Rd, Chapel Hill, NC, 27599, USA
Mohammad Y. Anwar Department of Epidemiology, Gillings School of Global Public Health, University of North Carolina, 135 Dauer Drive, Chapel Hill, NC, 27599, USA
Hung-Hsin Chen Institute of Biomedical Sciences, Academia Sinica, No. 128, Section 2, Academia Rd., Taipei, Nangang District, 115201, Taiwan Vanderbilt Genetics Institute, Division of Genetic Medicine, Vanderbilt University Medical Center, 1161 21st Ave S, Nashville, TN, 37232, USA
Lauren E. Petty Vanderbilt Genetics Institute, Division of Genetic Medicine, Vanderbilt University Medical Center, 1161 21st Ave S, Nashville, TN, 37232, USA
Sudha Seshadri Department of Biostatistics, School of Public Health, Boston University, 801 Massachusettes Ave, Boston, MA, 02118, USA Department of Neurology, School of Medicine, Boston University, 85 East Concord Street, Boston, MA, 02118, USA Glenn Biggs Institute for Alzheimer’s & Neurodegenerative Diseases, UT Health San Antonio, 8300 Floyd Curl Drive, San Antonio, TX, 78229, USA
Misa Graff Department of Epidemiology, Gillings School of Global Public Health, University of North Carolina, 135 Dauer Drive, Chapel Hill, NC, 27599, USA
Piper Below Vanderbilt Genetics Institute, Division of Genetic Medicine, Vanderbilt University Medical Center, 1161 21st Ave S, Nashville, TN, 37232, USA
Jennifer A. Brody Department of Medicine, Cardiovascular Health Research Unit, University of Washington, 1730 Minor Ave, Seattle, WA, 98101, USA
Geetha Chittoor Department of Population Health Sciences, Geisinger, 100 N. Academy Ave., Danville, PA, 17822, USA
Susan P. Fisher-Hoch Department of Epidemiology, School of Public Health, UT Health Houston, Regional Academic Health Center, One West University Blvd, Brownsville, TX, 78520, USA
Nancy L. Heard-Costa Framingham Heart Study, 73 Mt Wayte Ave, Framingham, MA, 01702, USA Department of Neurology, Chobanian & Avedisian School of Medicine, Boston University, 72 E Concord St, Boston, MA, 02118, USA
Daniel Levy Population Sciences Branch, National Heart, Lung, and Blood Institute of the National Institutes of Health, 6701 Rockledge Drive, Bethesda, MD, 20892, USA
Honghuang Lin Department of Medicine, University of Massachusetts Chan Medical School, 55 N Lake Ave, Worcester, MA, 01655, USA
Ruth JF. Loos Charles Bronfman Institute for Personalized Medicine at Mount Sinai, Icahn School of Medicine at Mount Sinai, 1 Gustave L. Levy Pl, New York, NY, 10029, USA Novo Nordisk Foundation Center for Basic Metabolic Research, Faculty of Health and Medical Sciences, University of Copenhagen, Blegdamsvej 3A, 2200, Copenhagen, Denmark
Joseph B. Mccormick Department of Epidemiology, School of Public Health, UT Health Houston, Regional Academic Health Center, One West University Blvd, Brownsville, TX, 78520, USA
Jerome I. Rotter Department of Pediatrics, The Institute for Translational Genomics and Population Sciences, The Lundquist Institute for Biomedical Innovation at Harbor-UCLA Medical Center, 1124 West Carson Street, Torrance, CA, 90502, USA
Tooraj Mirshahi Department of Genomic Health, Geisinger, 100 N. Academy Ave., Danville, PA, 17822, USA
Christopher D. Still Center for Obesity and Metabolic Health, Geisinger, 100 N. Academy Ave., Danville, PA, 17822, USA
Anita Destefano Department of Biostatistics, School of Public Health, Boston University, 801 Massachusettes Ave, Boston, MA, 02118, USA Department of Neurology, School of Medicine, Boston University, 85 East Concord Street, Boston, MA, 02118, USA
L. Adrienne Cupples Department of Biostatistics, School of Public Health, Boston University, 801 Massachusettes Ave, Boston, MA, 02118, USA
Karen L Mohlke Department of Genetics, School of Medicine, University of North Carolina, Chapel Hill, NC, 27599, USA
Kari E. North Department of Epidemiology, Gillings School of Global Public Health, University of North Carolina, 135 Dauer Drive, Chapel Hill, NC, 27599, USA
Anne E. Justice Department of Population Health Sciences, Geisinger, 100 N. Academy Ave., Danville, PA, 17822, USA
Ching-Ti Liu Department of Biostatistics, School of Public Health, Boston University, 801 Massachusettes Ave, Boston, MA, 02118, USA

Collapse

Genetics, pharmacotherapy, and dietary interventions in childhood obesity

Childhood obesity has emerged as a major global health issue, contributing to the increased prevalence of chronic conditions and adversely affecting the quality of life and future prospects of affected individuals, thereby presenting a substantial societal challenge. This complex condition, influenced by the interplay of genetic predispositions and environmental factors, is characterized by excessive energy intake due to uncontrolled appetite regulation and a Westernized diet. Managing obesity in childhood requires specific considerations compared with adulthood, given the vulnerability of the critical juvenile-adolescent period to toxicity and developmental defects. Consequently, common treatment options for adult obesity may not directly apply to younger populations. Therefore, research on childhood obesity has focused on genetic defects in regulating energy intake, alongside pharmacotherapy and dietary interventions as management approaches, with an emphasis on safety concerns. This review aims to summarize canonical knowledge and recent findings on genetic factors contributing to childhood obesity. Additionally, it assesses the efficacy and safety of existing pharmacotherapies and dietary interventions and suggests future research directions. By providing a comprehensive understanding of the complex dynamics of childhood obesity, this review aims to offer insights into more targeted and effective strategies for addressing this condition, including personalized healthcare solutions.

The Impact of Haplotypes of the FTO Gene, Lifestyle, and Dietary Patterns on BMI and Metabolic Syndrome in Polish Young Adult Men

BACKGROUND

Variants in fat mass and the obesity-associated protein (FTO) gene have long been recognized as the most significant genetic predictors of body fat mass and obesity. Nevertheless, despite the overall evidence, there are conflicting reports regarding the correlation between different polymorphisms of the FTO gene and body mass index (BMI). Additionally, it is unclear whether FTO influences metabolic syndrome (MetS) through mechanisms other than BMI's impact. In this work, we aimed to analyze the impact of the following FTO polymorphisms on the BMI as well as MetS components in a population of young adult men.

METHODS

The patient group consisted of 279 Polish young adult men aged 28.92 (4.28) recruited for the MAGNETIC trial. The single-nucleotide polymorphisms (SNPs), located in the first intron of the FTO gene, were genotyped, and the results were used to identify "protective" and "risk" haplotypes and diplotypes based on the literature data. Laboratory, as well as anthropometric measurements regarding MetS, were performed. Measured MetS components included those used in the definition in accordance with the current guidelines. Data regarding dietary patterns were also collected, and principal components of the dietary patterns were identified.

RESULTS

No statistically significant correlations were identified between the analyzed FTO diplotypes and BMI (p = 0.53) or other MetS components (waist circumference p = 0.55; triglycerides p = 0.72; HDL cholesterol p = 0.33; blood glucose p = 0.20; systolic blood pressure p = 0.06; diastolic blood pressure p = 0.21). Stratification by the level of physical activity or adherence to the dietary patterns also did not result in any statistically significant result.

CONCLUSIONS

Some studies have shown that FTO SNPs such as rs1421085, rs1121980, rs8050136, rs9939609, and rs9930506 have an impact on the BMI or other MetS components; nevertheless, this was not replicated in this study of Polish young adult males.

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.

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.

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.