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.
Publications
2024
Paternal diet can influence the phenotype of the next generation, yet, the dietary components inducing specific responses in the offspring are not identified. Here, we use the Nutritional Geometry Framework to determine the effects of pre-conception paternal dietary macronutrient balance on offspring metabolic and behavioral traits in mice. Ten isocaloric diets varying in the relative proportion of protein, fats, and carbohydrates are fed to male mice prior to mating. Dams and offspring are fed standard chow and never exposed to treatment diets. Body fat in female offspring is positively associated with the paternal consumption of fat, while in male offspring, an anxiety-like phenotype is associated to paternal diets low in protein and high in carbohydrates. Our study uncovers that the nature and the magnitude of paternal effects are driven by interactions between macronutrient balance and energy intake and are not solely the result of over- or undernutrition.
BACKGROUND: Diabetes in pregnancy is associated with increased risk of long-term metabolic disease in the offspring, potentially mediated by in utero epigenetic variation. Previously, we identified multiple differentially methylated single CpG sites in offspring of women with gestational diabetes mellitus (GDM), but whether stretches of differentially methylated regions (DMRs) can also be identified in adolescent GDM offspring is unknown. Here, we investigate which DNA regions in adolescent offspring are differentially methylated in blood by exposure to diabetes in pregnancy. The secondary aim was to characterize the RNA expression of the identified DMR, which contained the nc886 non-coding RNA.
METHODS: To identify DMRs, we employed the bump hunter method in samples from young (9-16 yr, n = 92) offspring of women with GDM (O-GDM) and control offspring (n = 94). Validation by pyrosequencing was performed in an adult offspring cohort (age 28-33 years) consisting of O-GDM (n = 82), offspring exposed to maternal type 1 diabetes (O-T1D, n = 67) and control offspring (O-BP, n = 57). RNA-expression was measured using RT-qPCR in subcutaneous adipose tissue and skeletal muscle.
RESULTS: One significant DMR represented by 10 CpGs with a bimodal methylation pattern was identified, located in the nc886/VTRNA2-1 non-coding RNA gene. Low methylation status across all CpGs of the nc886 in the young offspring was associated with maternal GDM. While low methylation degree in adult offspring in blood, adipose tissue, and skeletal muscle was not associated with maternal GDM, adipose tissue nc886 expression was increased in O-GDM compared to O-BP, but not in O-T1D. In addition, adipose tissue nc886 expression levels were positively associated with maternal pre-pregnancy BMI (p = 0.006), but not with the offspring's own adiposity.
CONCLUSIONS: Our results highlight that nc886 is a metastable epiallele, whose methylation in young offspring is negatively correlated with maternal obesity and GDM status. The physiological effect of nc886 may be more important in adipose tissue than in skeletal muscle. Further research should aim to investigate how nc886 regulation in adipose tissue by exposure to GDM may contribute to development of metabolic disease.
Harsh environments in poorly perfused tumor regions may select for traits driving cancer aggressiveness. Here, we investigated whether tumor acidosis interacts with driver mutations to exacerbate cancer hallmarks. We adapted mouse organoids from normal pancreatic duct (mN10) and early pancreatic cancer (mP4, KRAS-G12D mutation, ± p53 knockout) from extracellular pH 7.4 to 6.7, representing acidic niches. Viability was increased by acid adaptation, a pattern most apparent in wild-type (WT) p53 organoids, and exacerbated upon return to pH 7.4. This led to increased survival of acid-adapted organoids treated with gemcitabine and/or erlotinib, and, in WT p53 organoids, acid-induced attenuation of drug effects. New genetic variants became dominant during adaptation, yet they were unlikely to be its main drivers. Transcriptional changes induced by acid and drug adaptation differed overall, but acid adaptation increased the expression of gemcitabine resistance genes. Thus, adaptation to acidosis increases cancer cell viability after chemotherapy.
The transcription factor MYC is overexpressed in most cancers, where it drives multiple hallmarks of cancer progression. MYC is known to promote oncogenic transcription by binding to active promoters. In addition, MYC has also been shown to invade distal enhancers when expressed at oncogenic levels, but this enhancer binding has been proposed to have low gene-regulatory potential. Here, we demonstrate that MYC directly regulates enhancer activity to promote cancer type-specific gene programs predictive of poor patient prognosis. MYC induces transcription of enhancer RNA through recruitment of RNA polymerase II (RNAPII), rather than regulating RNAPII pause-release, as is the case at promoters. This process is mediated by MYC-induced H3K9 demethylation and acetylation by GCN5, leading to enhancer-specific BRD4 recruitment through its bromodomains, which facilitates RNAPII recruitment. We propose that MYC drives prognostic cancer type-specific gene programs through induction of an enhancer-specific epigenetic switch, which can be targeted by BET and GCN5 inhibitors.
Genetic studies have identified numerous regions associated with plasma fibrinogen levels in Europeans, yet missing heritability and limited inclusion of non-Europeans necessitates further studies with improved power and sensitivity. Compared with array-based genotyping, whole-genome sequencing (WGS) data provide better coverage of the genome and better representation of non-European variants. To better understand the genetic landscape regulating plasma fibrinogen levels, we meta-analyzed WGS data from the National Heart, Lung, and Blood Institute's Trans-Omics for Precision Medicine (TOPMed) program (n = 32 572), with array-based genotype data from the Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium (n = 131 340) imputed to the TOPMed or Haplotype Reference Consortium panel. We identified 18 loci that have not been identified in prior genetic studies of fibrinogen. Of these, 4 are driven by common variants of small effect with reported minor allele frequency (MAF) at least 10 percentage points higher in African populations. Three signals (SERPINA1, ZFP36L2, and TLR10) contain predicted deleterious missense variants. Two loci, SOCS3 and HPN, each harbor 2 conditionally distinct, noncoding variants. The gene region encoding the fibrinogen protein chain subunits (FGG;FGB;FGA) contains 7 distinct signals, including 1 novel signal driven by rs28577061, a variant common in African ancestry populations but extremely rare in Europeans (MAFAFR = 0.180; MAFEUR = 0.008). Through phenome-wide association studies in the VA Million Veteran Program, we found associations between fibrinogen polygenic risk scores and thrombotic and inflammatory disease phenotypes, including an association with gout. Our findings demonstrate the utility of WGS to augment genetic discovery in diverse populations and offer new insights for putative mechanisms of fibrinogen regulation.
OBJECTIVES: Experimental and genetic studies suggest that fibroblast growth factor 21 (FGF21) modulates macronutrient and alcohol preferences, but evidence of such regulation in humans remains scarce. To address this gap in translation, we aimed to map the relationships between plasma FGF21 levels, FGF21 genetic variation and habitual macronutrient intake in a large human population.
METHODS: We fine-mapped and performed colocalization of the FGF21 genetic region in GWAS summary statistics of plasma FGF21 levels and macronutrient intake. UK Biobank data were used to investigate the associations between FGF21 genetic variants, plasma FGF21 protein levels, and macronutrient intake (including alcohol) assessed with repeated 24-hour recalls. One- and two-sample mendelian randomization were performed to estimate the effects of plasma FGF21 on macronutrient intake.
RESULTS: We show that the main macronutrient-associated variant rs838133 and the FGF21 cis-pQTL rs838131, both in the FGF21 gene, are distinct genetic signals. Effect directions also suggest that the influence of FGF21 variation on macronutrient intake appear more complex than by direct mediation through plasma FGF21. Only when considering this complexity at FGF21, is plasma FGF21 estimated to reduce alcohol and increase protein and fat intake using mendelian randomization. Importantly, plasma FGF21 levels also appear markedly elevated by primarily high alcohol and low protein intake.
CONCLUSIONS: These findings support the feedback diet-regulatory mechanism of FGF21 in humans, but highlights the need for mechanistic characterization of the complex FGF21 genetic region.
To date only a fraction of the genetic footprint of thyroid function has been clarified. We report a genome-wide association study meta-analysis of thyroid function in up to 271,040 individuals of European ancestry, including reference range thyrotropin (TSH), free thyroxine (FT4), free and total triiodothyronine (T3), proxies for metabolism (T3/FT4 ratio) as well as dichotomized high and low TSH levels. We revealed 259 independent significant associations for TSH (61% novel), 85 for FT4 (67% novel), and 62 novel signals for the T3 related traits. The loci explained 14.1%, 6.0%, 9.5% and 1.1% of the total variation in TSH, FT4, total T3 and free T3 concentrations, respectively. Genetic correlations indicate that TSH associated loci reflect the thyroid function determined by free T3, whereas the FT4 associations represent the thyroid hormone metabolism. Polygenic risk score and Mendelian randomization analyses showed the effects of genetically determined variation in thyroid function on various clinical outcomes, including cardiovascular risk factors and diseases, autoimmune diseases, and cancer. In conclusion, our results improve the understanding of thyroid hormone physiology and highlight the pleiotropic effects of thyroid function on various diseases.
CONTEXT: Insulin sensitivity (IS) is an important factor in type 2 diabetes (T2D) and can be estimated by many different indices.
OBJECTIVE: We aimed to compare the genetic components underlying IS indices obtained from fasting and oral glucose-stimulated plasma glucose and serum insulin levels.
METHODS: We computed 21 IS indices, classified as fasting, OGTT0,120, and OGTT0,30,120 indices, using fasting and oral glucose tolerance test (OGTT) data in 2 cohorts. We used data from a family cohort (n = 313) to estimate the heritability and the genetic and phenotypic correlations of IS indices. The population cohort, Inter99 (n = 5343), was used to test for associations between IS indices and 426 genetic variants known to be associated with T2D.
RESULTS: Heritability estimates of IS indices ranged between 19% and 38%. Fasting and OGTT0,30,120 indices had high genetic (ρG) and phenotypic (ρP) pairwise correlations (ρG and ρP: 0.88 to 1) The OGTT0,120 indices displayed a wide range of pairwise correlations (ρG: 0.17-1.00 and ρP: 0.13-0.97). We identified statistically significant associations between IS indices and established T2D-associated variants. The PPARG rs11709077 variant was associated only with fasting indices and PIK3R rs4976033 only with OGTT0,30,120 indices. The variants in FAM63A/MINDY1, GCK, C2CD4A/B, and FTO loci were associated only with OGTT0,120 indices.
CONCLUSION: Even though the IS indices mostly share a common genetic background, notable differences emerged between OGTT0,120 indices. The fasting and OGTT-based indices have distinct associations with T2D risk variants. This work provides a basis for future large-scale genetic investigations into the differences between IS indices.
In the Greenlandic Inuit population, 4% are homozygous carriers of a genetic nonsense TBC1D4 p.Arg684Ter variant leading to loss of the muscle-specific isoform of TBC1D4 and an approximately tenfold increased risk of type 2 diabetes1. Here we show the metabolic consequences of this variant in four female and four male homozygous carriers and matched controls. An extended glucose tolerance test reveals prolonged hyperglycaemia followed by reactive hypoglycaemia in the carriers. Whole-body glucose disposal is impaired during euglycaemic-hyperinsulinaemic clamp conditions and associates with severe insulin resistance in skeletal muscle only. Notably, a marked reduction in muscle glucose transporter GLUT4 and associated proteins is observed. While metabolic regulation during exercise remains normal, the insulin-sensitizing effect of a single exercise bout is compromised. Thus, loss of the muscle-specific isoform of TBC1D4 causes severe skeletal muscle insulin resistance without baseline hyperinsulinaemia. However, physical activity can ameliorate this condition. These observations offer avenues for personalized interventions and targeted preventive strategies.