The association between body mass index and metabolite response to a liquid mixed meal challenge: a Mendelian randomization study.

BACKGROUND: Metabolite abundance is a dynamic trait that varies in response to environmental stimuli and phenotypic traits, such as food consumption and body mass index (BMI, kg/m2). OBJECTIVES: In this study, we used the Netherlands Epidemiology of Obesity (NEO) study data to identify observational and causal associations between BMI and metabolite response to a liquid meal. METHODS: A liquid meal challenge was performed, and Nightingale Health metabolite profiles were collected in 5744 NEO participants. Observational and one-sample Mendelian randomization (MR) analysis were conducted to estimate the effect of BMI on metabolites (n = 229) in the fasting, postprandial, and response (or change in abundance) states. RESULTS: We observed 473 associations with BMI (175 fasting, 188 postprandial, and 110 response) in observational analyses. In MR analyses, we observed 20 metabolite traits (5 fasting, 12 postprandial, and 3 response) to be associated with BMI. MR associations included the glucogenic amino acid alanine, which was inversely associated with BMI in the response state (ß: -0.081; SE: 0.023; P = 5.91 × 10-4), suggesting that as alanine increased in postprandial abundance, that increase was attenuated with increasing BMI. CONCLUSIONS: Overall, this study showed that MR estimates were strongly correlated with observational effect estimates, suggesting that the broad associations seen between BMI and metabolite variation has a causal underpinning. Specific effects in previously unassessed postprandial and response states are detected, and these may likely mark novel life course risk exposures driven by regular nutrition.

Genome-wide characterization of circulating metabolic biomarkers.

Genome-wide association analyses using high-throughput metabolomics platforms have led to novel insights into the biology of human metabolism1-7. This detailed knowledge of the genetic determinants of systemic metabolism has been pivotal for uncovering how genetic pathways influence biological mechanisms and complex diseases8-11. Here we present a genome-wide association study for 233 circulating metabolic traits quantified by nuclear magnetic resonance spectroscopy in up to 136,016 participants from 33 cohorts. We identify more than 400 independent loci and assign probable causal genes at two-thirds of these using manual curation of plausible biological candidates. We highlight the importance of sample and participant characteristics that can have significant effects on genetic associations. We use detailed metabolic profiling of lipoprotein- and lipid-associated variants to better characterize how known lipid loci and novel loci affect lipoprotein metabolism at a granular level. We demonstrate the translational utility of comprehensively phenotyped molecular data, characterizing the metabolic associations of intrahepatic cholestasis of pregnancy. Finally, we observe substantial genetic pleiotropy for multiple metabolic pathways and illustrate the importance of careful instrument selection in Mendelian randomization analysis, revealing a putative causal relationship between acetone and hypertension. Our publicly available results provide a foundational resource for the community to examine the role of metabolism across diverse diseases.

Nutrition and Cardiovascular Disease: The Potential Role of Marine Bioactive Proteins and Peptides in Thrombosis Prevention.

Thrombus and cardiovascular diseases pose a significant health threat, and dietary interventions have shown promising potential in reducing the incidence of these diseases. Marine bioactive proteins and peptides have been extensively studied for their antithrombotic properties. They can inhibit platelet activation and aggregation by binding to key receptors on the platelet surface. Additionally, they can competitively anchor to critical enzyme sites, leading to the inhibition of coagulation factors. Marine microorganisms also offer alternative sources for the development of novel fibrinolytic proteins, which can help dissolve blood clots. The advancements in technologies, such as targeted hydrolysis, specific purification, and encapsulation, have provided a solid foundation for the industrialization of bioactive peptides. These techniques enable precise control over the production and delivery of bioactive peptides, enhancing their efficacy and safety. However, it is important to note that further research and clinical studies are needed to fully understand the mechanisms of action and therapeutic potential of marine bioactive proteins and peptides in mitigating thrombotic events. The challenges and future application perspectives of these bioactive peptides also need to be explored.

Association between venous thromboembolism-associated genetic variants, coagulation factor levels, and thrombin generation potential.

Recently three large meta-analyses of genome-wide association studies for venous thromboembolism (VTE) identified over 130 genetic variants. However, mechanisms by which newly identified and therefore underexplored VTE-associated genetic variants influence VTE remain unclear. To elucidate the mechanism, we investigated the association between 61 newly identified VTE-associated genetic variants and the levels of coagulation factor (F) VIII, FIX, FXI, and fibrinogen as well as thrombin generation parameters (lag time, peak, endogenous thrombin potential, time-to-peak, and velocity), which are well-known biological traits associated with VTE. This study was conducted on 5341 participants of the Netherlands Epidemiology of Obesity study. The associations between VTE-associated genetic variants and coagulation factor levels and thrombin generation parameters were examined using linear regression analyses, adjusted for age, sex, body mass index, oral contraceptive use, hormone replacement therapy, and menopausal status. Of 61 genetic variants, 33 were associated with one or more of the coagulation factor levels and thrombin generation parameters. Following multiple testing corrections, five genetic variants remained significant, of which MAP1A rs55707100 exhibited the most robust association with thrombin generation parameters and FXI levels (ß = -5.33%, 95% confidence interval: -8.44, -2.22). Our findings shed light on the underlying mechanisms by which these genetic variants influence the risk of VTE.

Partitioned polygenic risk scores identify distinct types of metabolic dysfunction-associated steatotic liver disease.

Metabolic dysfunction-associated steatotic liver disease (MASLD) encompasses an excess of triglycerides in the liver, which can lead to cirrhosis and liver cancer. While there is solid epidemiological evidence of MASLD coexisting with cardiometabolic disease, several leading genetic risk factors for MASLD do not increase the risk of cardiovascular disease, suggesting no causal relationship between MASLD and cardiometabolic derangement. In this work, we leveraged measurements of visceral adiposity and identified 27 novel genetic loci associated with MASLD. Among these loci, we replicated 6 in several independent cohorts. Next, we generated two partitioned polygenic risk scores (PRS) based on the mechanism of genetic association with MASLD encompassing intra-hepatic lipoprotein retention. The two PRS suggest the presence of at least two distinct types of MASLD, one confined to the liver resulting in a more aggressive liver disease and one that is systemic and results in a higher risk of cardiometabolic disease.

A genetic association study of circulating coagulation factor VIII and von Willebrand factor levels.

ABSTRACT: Coagulation factor VIII (FVIII) and its carrier protein von Willebrand factor (VWF) are critical to coagulation and platelet aggregation. We leveraged whole-genome sequence data from the Trans-Omics for Precision Medicine (TOPMed) program along with TOPMed-based imputation of genotypes in additional samples to identify genetic associations with circulating FVIII and VWF levels in a single-variant meta-analysis, including up to 45 289 participants. Gene-based aggregate tests were implemented in TOPMed. We identified 3 candidate causal genes and tested their functional effect on FVIII release from human liver endothelial cells (HLECs) and VWF release from human umbilical vein endothelial cells. Mendelian randomization was also performed to provide evidence for causal associations of FVIII and VWF with thrombotic outcomes. We identified associations (P < 5 × 10-9) at 7 new loci for FVIII (ST3GAL4, CLEC4M, B3GNT2, ASGR1, F12, KNG1, and TREM1/NCR2) and 1 for VWF (B3GNT2). VWF, ABO, and STAB2 were associated with FVIII and VWF in gene-based analyses. Multiphenotype analysis of FVIII and VWF identified another 3 new loci, including PDIA3. Silencing of B3GNT2 and the previously reported CD36 gene decreased release of FVIII by HLECs, whereas silencing of B3GNT2, CD36, and PDIA3 decreased release of VWF by HVECs. Mendelian randomization supports causal association of higher FVIII and VWF with increased risk of thrombotic outcomes. Seven new loci were identified for FVIII and 1 for VWF, with evidence supporting causal associations of FVIII and VWF with thrombotic outcomes. B3GNT2, CD36, and PDIA3 modulate the release of FVIII and/or VWF in vitro.

Whole genome analysis of plasma fibrinogen reveals population-differentiated genetic regulators with putative liver roles.

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 provides 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 NHLBI'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 (CHARGE) 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, four are driven by common variants of small effect with reported MAF at least 10% higher in African populations. Three ( SERPINA1, ZFP36L2 , and TLR10) signals contain predicted deleterious missense variants. Two loci, SOCS3 and HPN , each harbor two conditionally distinct, non-coding variants. The gene region encoding the protein chain subunits ( FGG;FGB;FGA ), contains 7 distinct signals, including one novel signal driven by rs28577061, a variant common (MAF=0.180) in African reference panels but extremely rare (MAF=0.008) in Europeans. 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. Key Points: Largest and most diverse genetic study of plasma fibrinogen identifies 54 regions (18 novel), housing 69 conditionally distinct variants (20 novel).Sufficient power achieved to identify signal driven by African population variant.Links to (1) liver enzyme, blood cell and lipid genetic signals, (2) liver regulatory elements, and (3) thrombotic and inflammatory disease.

Bone mineral density loci specific to the skull portray potential pleiotropic effects on craniosynostosis.

Skull bone mineral density (SK-BMD) provides a suitable trait for the discovery of key genes in bone biology, particularly to intramembranous ossification, not captured at other skeletal sites. We perform a genome-wide association meta-analysis (n ~ 43,800) of SK-BMD, identifying 59 loci, collectively explaining 12.5% of the trait variance. Association signals cluster within gene-sets involved in skeletal development and osteoporosis. Among the four novel loci (ZIC1, PRKAR1A, AZIN1/ATP6V1C1, GLRX3), there are factors implicated in intramembranous ossification and as we show, inherent to craniosynostosis processes. Functional follow-up in zebrafish confirms the importance of ZIC1 on cranial suture patterning. Likewise, we observe abnormal cranial bone initiation that culminates in ectopic sutures and reduced BMD in mosaic atp6v1c1 knockouts. Mosaic prkar1a knockouts present asymmetric bone growth and, conversely, elevated BMD. In light of this evidence linking SK-BMD loci to craniofacial abnormalities, our study provides new insight into the pathophysiology, diagnosis and treatment of skeletal diseases.

Sex-specific association between microvascular health and coagulation parameters: the Netherlands Epidemiology of Obesity study.

BACKGROUND: Microvascular dysfunction is a growing determinant of sex differences in coronary heart disease (CHD). Dysregulation of the coagulation system is involved in CHD pathogenesis and can be induced by endothelial glycocalyx (EG) perturbation. However, little is known about the link between EG function and coagulation parameters in population-based studies on sex specificity. OBJECTIVES: We sought to examine the sex differences in the relationship between EG function and coagulation parameters in a middle-aged Dutch population. METHODS: Using baseline measurements of 771 participants from the Netherlands Epidemiology of Obesity study (age, 56 years [IQR, 51-61 years]; 53% women; body mass index, 27.9 kg/m2 [IQR, 25.1-30.9 kg/m2]), associations between glycocalyx-related perfused boundary region (PBR) derived using sidestream dark-field imaging and coagulation parameters (factor [F]VIII/IX/XI; thrombin generation parameters; and fibrinogen) were investigated using linear regression analyses, adjusting for possible confounders (including C-reactive protein, leptin, and glycoprotein acetyls), followed by sex-stratified analyses. RESULTS: There was a sex difference in the associations between PBR and coagulation parameters. Particularly in women, 1-SD PBR (both total and feed vessel, indicating poorer glycocalyx status) was associated with higher FIX activity ([1.8%; 95% CI, 0.3%-3.3%] and [2.0%; 95% CI, 0.5%-3.4%], respectively) and plasma fibrinogen levels ([5.1 mg/dL; 95% CI, 0.4-9.9 mg/dL] and [5.8 mg/dL; 95% CI, 1.1-10.6 mg/dL], respectively). Furthermore, 1-SD PBRcapillary was associated with higher FVIII activity (3.5%; 95% CI, 0.4%-6.5%) and plasma fibrinogen levels (5.3 mg/dL; 95% CI, 0.6-10.0 mg/dL). CONCLUSION: We revealed a sex-specific association between microcirculatory health and procoagulant status, which suggests that microvascular health be considered during early development of CHD in women.

Altered high-density lipoprotein composition is associated with risk for complications in type 2 diabetes mellitus in South Asian descendants: A cross-sectional, case-control study on lipoprotein subclass profiling.

BACKGROUND: Composition of high-density lipoproteins (HDL) is emerging as an important determinant in the development of microvascular complications in type 2 diabetes mellitus (T2DM). Dutch South Asian (DSA) individuals with T2DM display an increased risk of microvascular complications compared with Dutch white Caucasian (DwC) individuals with T2DM. In this study, we aimed to investigate whether changes in HDL composition associate with increased microvascular risk in this ethnic group and lead to new lipoprotein biomarkers. MATERIALS AND METHODS: Using 1 H nuclear magnetic resonance spectroscopy and Bruker IVDr Lipoprotein Subclass Analysis (B.I.LISA) software, plasma lipoprotein changes were determined in 51 healthy individuals (30 DwC, 21 DSA) and 92 individuals with T2DM (45 DwC, 47 DSA) in a cross-sectional, case-control study. Differential HDL subfractions were investigated using multinomial logistic regression analyses, adjusting for possible confounders including BMI and diabetes duration. RESULTS: We identified HDL compositional differences between healthy and diabetic individuals in both ethnic groups. Specifically, levels of apolipoprotein A2 and HDL-4 subfractions were lower in DSA compared with DwC with T2DM. Apolipoprotein A2 and HDL-4 subfractions also negatively correlated with waist circumference, waist-to-hip ratio, haemoglobin A1c, glucose levels and disease duration in DSA with T2DM, and associated with increased incidence of microvascular complications. CONCLUSION: While HDL composition differed between controls and T2DM in both ethnic groups, the lower levels of lipid content in the smallest HDL subclass (HDL-4) in DSA with T2DM appeared to be more clinically relevant, with higher odds of having diabetes-related pan-microvascular complications such as retinopathy and neuropathy. These typical differences in HDL could be used as ethnicity-specific T2DM biomarkers.

Hepatic triglyceride content is intricately associated with numerous metabolites and biochemical pathways.

BACKGROUND AND AIMS: Non-alcoholic fatty liver disease (NAFLD) is characterized by the pathological accumulation of triglycerides in hepatocytes and is associated with insulin resistance, atherogenic dyslipidaemia and cardiometabolic diseases. Thus far, the extent of metabolic dysregulation associated with hepatic triglyceride accumulation has not been fully addressed. In this study, we aimed to identify metabolites associated with hepatic triglyceride content (HTGC) and map these associations using network analysis. METHODS: To gain insight in the spectrum of metabolites associated with hepatic triglyceride accumulation, we performed a comprehensive plasma metabolomics screening of 1363 metabolites in apparently healthy middle aged (age 45-65) individuals (N = 496) in whom HTGC was measured by proton magnetic resonance spectroscopy. An atlas of metabolite-HTGC associations, based on univariate results, was created using correlation-based Gaussian graphical model (GGM) and genome scale metabolic model network analyses. Pathways associated with the clinical prognosis marker fibrosis 4 (FIB-4) index were tested using a closed global test. RESULTS: Our analyses revealed that 118 metabolites were univariately associated with HTGC (p-value <6.59 × 10-5 ), including 106 endogenous, 1 xenobiotic and 11 partially characterized/uncharacterized metabolites. These associations were mapped to several biological pathways including branched amino acids (BCAA), diglycerols, sphingomyelin, glucosyl-ceramide and lactosyl-ceramide. We also identified a novel possible HTGC-related pathway connecting glutamate, metabolonic lactone sulphate and X-15245 using the GGM network. These pathways were confirmed to be associated with the FIB-4 index as well. The full interactive metabolite-HTGC atlas is provided online: https://tofaquih.github.io/AtlasLiver/. CONCLUSIONS: The combined network and pathway analyses indicated extensive associations between BCAA and the lipids pathways with HTGC and the FIB-4 index. Moreover, we report a novel pathway glutamate-metabolonic lactone sulphate-X-15245 with a potential strong association with HTGC. These findings can aid elucidating HTGC metabolomic profiles and provide insight into novel drug targets for fibrosis-related outcomes.

Normal range CAG repeat size variations in the HTT gene are associated with an adverse lipoprotein profile partially mediated by body mass index.

Tandem cytosine-adenine-guanine (CAG) repeat sizes of 36 or more in the huntingtin gene (HTT) cause Huntington's disease (HD). Apart from neuropsychiatric complications, the disease is also accompanied by metabolic dysregulation and weight loss, which contribute to a progressive functional decline. Recent studies also reported an association between repeats below the pathogenic threshold (<36) for HD and body mass index (BMI), suggesting that HTT repeat sizes in the non-pathogenic range are associated with metabolic dysregulation. In this study, we hypothesized that HTT repeat sizes < 36 are associated with metabolite levels, possibly mediated through reduced BMI. We pooled data from three European cohorts (n = 10 228) with genotyped HTT CAG repeat size and metabolomic measurements. All 145 metabolites were measured on the same targeted platform in all studies. Multilevel mixed-effects analysis using the CAG repeat size in HTT identified 67 repeat size metabolite associations. Overall, the metabolomic profile associated with larger CAG repeat sizes in HTT were unfavorable-similar to those of higher risk of coronary artery disease and type 2 diabetes-and included elevated levels of amino acids, fatty acids, low-density lipoprotein (LDL)-, very low-density lipoprotein- and intermediate density lipoprotein (IDL)-related metabolites while with decreased levels of very large high-density lipoprotein (HDL)-related metabolites. Furthermore, the associations of 50 metabolites, in particular, specific very large HDL-related metabolites, were mediated by lower BMI. However, no mediation effect was found for 17 metabolites related to LDL and IDL. In conclusion, our findings indicate that large non-pathogenic CAG repeat sizes in HTT are associated with an unfavorable metabolomic profile despite their association with a lower BMI.

The Roles of Gut Microbiome and Plasma Metabolites in the Associations between ABO Blood Groups and Insulin Homeostasis: The Microbiome and Insulin Longitudinal Evaluation Study (MILES).

Non-O blood groups are associated with decreased insulin sensitivity and risk of type 2 diabetes. A recent study pinpointed the associations between ABO blood groups and gut microbiome, which may serve as potential mediators for the observed increased disease risks. We aimed to characterize associations between ABO haplotypes and insulin-related traits as well as potential mediating pathways. We assessed insulin homeostasis in African Americans (AAs; n = 109) and non-Hispanic whites (n = 210) from the Microbiome and Insulin Longitudinal Evaluation Study. The ABO haplotype was determined by six SNPs located in the ABO gene. Based on prior knowledge, we included 21 gut bacteria and 13 plasma metabolites for mediation analysis. In the white study cohort (60 ± 9 years, 42% male), compared to the O1 haplotype, A1 was associated with a higher Matsuda insulin sensitivity index, while a lower relative abundance of Bacteroides massiliensis and lactate levels. Lactate was a likely mediator of this association but not Bacteroides massiliensis. In the AAs group (57 ± 8 years, 33% male), we found no association between any haplotype and insulin-related traits. In conclusion, the A1 haplotype may promote healthy insulin sensitivity in non-Hispanic whites and lactate likely play a role in this process but not selected gut bacteria.

Cross-Ancestry Investigation of Venous Thromboembolism Genomic Predictors.

BACKGROUND: Venous thromboembolism (VTE) is a life-threatening vascular event with environmental and genetic determinants. Recent VTE genome-wide association studies (GWAS) meta-analyses involved nearly 30 000 VTE cases and identified up to 40 genetic loci associated with VTE risk, including loci not previously suspected to play a role in hemostasis. The aim of our research was to expand discovery of new genetic loci associated with VTE by using cross-ancestry genomic resources. METHODS: We present new cross-ancestry meta-analyzed GWAS results involving up to 81 669 VTE cases from 30 studies, with replication of novel loci in independent populations and loci characterization through in silico genomic interrogations. RESULTS: In our genetic discovery effort that included 55 330 participants with VTE (47 822 European, 6320 African, and 1188 Hispanic ancestry), we identified 48 novel associations, of which 34 were replicated after correction for multiple testing. In our combined discovery-replication analysis (81 669 VTE participants) and ancestry-stratified meta-analyses (European, African, and Hispanic), we identified another 44 novel associations, which are new candidate VTE-associated loci requiring replication. In total, across all GWAS meta-analyses, we identified 135 independent genomic loci significantly associated with VTE risk. A genetic risk score of the significantly associated loci in Europeans identified a 6-fold increase in risk for those in the top 1% of scores compared with those with average scores. We also identified 31 novel transcript associations in transcriptome-wide association studies and 8 novel candidate genes with protein quantitative-trait locus Mendelian randomization analyses. In silico interrogations of hemostasis and hematology traits and a large phenome-wide association analysis of the 135 GWAS loci provided insights to biological pathways contributing to VTE, with some loci contributing to VTE through well-characterized coagulation pathways and others providing new data on the role of hematology traits, particularly platelet function. Many of the replicated loci are outside of known or currently hypothesized pathways to thrombosis. CONCLUSIONS: Our cross-ancestry GWAS meta-analyses identified new loci associated with VTE. These findings highlight new pathways to thrombosis and provide novel molecules that may be useful in the development of improved antithrombosis treatments.

Genetic and environmental contributions to stability and change in social inhibition across the adolescent and adult life span.

Feeling inhibited and socially not at ease is reflected in the trait social inhibition (SI). SI is associated with psychopathology that arises in young adulthood, such as anxiety. We aim for a better insight into the genetic and environmental contributions to SI across the life span, and as such examine their contributions to SI stability and change across adolescent and adult life span. We analyzed cohort-sequential longitudinal data from the Netherlands Twin Register (NTR), spanning a period of 25 years (Men (N, %): 17855, 37.4%; Age (Median, IQR): 19 years, 16-26 years; 7474 complete MZ twins and 8799 complete DZ twins). The data were organized into 7 age groups: < 14 (preadolescence), 15-16 (early adolescence), 17-18 (mid adolescence), 19-20 (late adolescence), 21-30 (young adulthood), 31-40 (adulthood), 41 + (middle-age-older adulthood). SI was assessed with the ASEBA-based proxy questionnaire. Phenotypic stability was established across the entire age range. Next, a longitudinal genetic simplex model was fitted to estimate the genetic and environmental contributions to the observed phenotypic stability. Results showed SI correlated well across follow-up of a single decade (.44 ≤ r ≤ .59) and moderately across the 25 years (.23 - .32) from adolescence to middle-age and older. Broad-sense heritability (h²) was between 40 and 48% across the 7 age groups. Additive and nonadditive genetic effects together explained most of the stability of SI across the life span (about 60-90% of the phenotypic correlation between ages), whereas environmental effects played a lesser role (about 10-40%). Concluding, SI, known to increase the risk of internalizing psychopathology, is phenotypically stable across the life span, which is largely attributable to genetic contributions to individual differences in SI. (PsycInfo Database Record (c) 2022 APA, all rights reserved).

Depression and Inflammatory Bowel Disease: A Bidirectional Two-sample Mendelian Randomization Study.

BACKGROUND AND AIMS: Observational studies have suggested a bidirectional association between depression and inflammatory bowel disease [IBD], including Crohn's disease [CD] and ulcerative colitis [UC]. However, it remains unclear whether the observed associations are causal due to the difficulties of determining sequential temporality. We investigated the association between depression and IBD by using bidirectional two-sample Mendelian randomization [MR]. METHODS: Independent genetic variants for depression and IBD were selected as instruments from published genome-wide association studies [GWAS] among individuals of predominantly European ancestry. Summary statistics for instrument-outcome associations were retrieved from three separate databases for both depression [Psychiatric Genomics Consortium, FinnGen and UK Biobank] and IBD [the largest GWAS meta-analysis, FinnGen and UK Biobank], respectively. MR analyses included the inverse-variance-weighted method, weighted-median estimator, MR-Egger regression, and sensitivity analyses of Steiger filtering and MR PRESSO. From either direction, analyses were performed per outcome database and were subsequently meta-analysed using a fixed-effect model. RESULTS: Genetically predicted depression [per log-odds ratio increase] was associated with a higher risk of IBD; odds ratios [95% confidence interval] for IBD, CD and UC were 1.20 [1.05, 1.36], 1.29 [1.07, 1.56] and 1.22 [1.01, 1.47] in a combined sample size of 693 183 [36 507 IBD cases], 212 172 [13 714 CD cases] and 219 686 [15 691 UC cases] individuals, respectively. In contrast, no association was observed between genetically influenced IBD and depression in 534 635 individuals [71 466 depression cases]. CONCLUSIONS: Our findings corroborated a causal association of depression on IBD, which may impact the clinical decision on the management of depression in patients with IBD. Though our results did not support a causal effect of IBD on depression, further investigations are needed to clarify the effect of IBD activity on depression [with different symptomology].

Higher thyrotropin leads to unfavorable lipid profile and somewhat higher cardiovascular disease risk: evidence from multi-cohort Mendelian randomization and metabolomic profiling.

BACKGROUND: Observational studies suggest interconnections between thyroid status, metabolism, and risk of coronary artery disease (CAD), but causality remains to be proven. The present study aimed to investigate the potential causal relationship between thyroid status and cardiovascular disease and to characterize the metabolomic profile associated with thyroid status. METHODS: Multi-cohort two-sample Mendelian randomization (MR) was performed utilizing genome-wide significant variants as instruments for standardized thyrotropin (TSH) and free thyroxine (fT4) within the reference range. Associations between TSH and fT4 and metabolic profile were investigated in a two-stage manner: associations between TSH and fT4 and the full panel of 161 metabolomic markers were first assessed hypothesis-free, then directional consistency was assessed through Mendelian randomization, another metabolic profile platform, and in individuals with biochemically defined thyroid dysfunction. RESULTS: Circulating TSH was associated with 52/161 metabolomic markers, and fT4 levels were associated with 21/161 metabolomic markers among 9432 euthyroid individuals (median age varied from 23.0 to 75.4 years, 54.5% women). Positive associations between circulating TSH levels and concentrations of very low-density lipoprotein subclasses and components, triglycerides, and triglyceride content of lipoproteins were directionally consistent across the multivariable regression, MR, metabolomic platforms, and for individuals with hypo- and hyperthyroidism. Associations with fT4 levels inversely reflected those observed with TSH. Among 91,810 CAD cases and 656,091 controls of European ancestry, per 1-SD increase of genetically determined TSH concentration risk of CAD increased slightly, but not significantly, with an OR of 1.03 (95% CI 0.99-1.07; p value 0.16), whereas higher genetically determined fT4 levels were not associated with CAD risk (OR 1.00 per SD increase of fT4; 95% CI 0.96-1.04; p value 0.59). CONCLUSIONS: Lower thyroid status leads to an unfavorable lipid profile and a somewhat increased cardiovascular disease risk.

Genetic Studies of Metabolomics Change After a Liquid Meal Illuminate Novel Pathways for Glucose and Lipid Metabolism.

Humans spend the greater part of the day in a postprandial state. However, the genetic basis of postprandial blood measures is relatively uncharted territory. We examined the genetics of variation in concentrations of postprandial metabolites (t = 150 min) in response to a liquid mixed meal through genome-wide association studies (GWAS) performed in the Netherlands Epidemiology of Obesity (NEO) study (n = 5,705). The metabolite response GWAS identified an association between glucose change and rs10830963:G in the melatonin receptor 1B (ß [SE] -0.23 [0.03], P = 2.15 × 10-19). In addition, the ANKRD55 locus led by rs458741:C showed strong associations with extremely large VLDL (XXLVLDL) particle response (XXLVLDL total cholesterol: ß [SE] 0.17 [0.03], P = 5.76 × 10-10; XXLVLDL cholesterol ester: ß [SE] 0.17 [0.03], P = 9.74 × 10-10), which also revealed strong associations with body composition and diabetes in the UK Biobank (P < 5 × 10-8). Furthermore, the associations between XXLVLDL response and insulinogenic index, HOMA-ß, Matsuda insulin sensitivity index, and HbA1c in the NEO study implied the role of chylomicron synthesis in diabetes (with false discovery rate-corrected q <0.05). To conclude, genetic studies of metabolomics change after a liquid meal illuminate novel pathways for glucose and lipid metabolism. Further studies are warranted to corroborate biological pathways of the ANKRD55 locus underlying diabetes.

Rare and low-frequency exonic variants and gene-by-smoking interactions in pulmonary function.

Genome-wide association studies have identified numerous common genetic variants associated with spirometric measures of pulmonary function, including forced expiratory volume in one second (FEV1), forced vital capacity, and their ratio. However, variants with lower minor allele frequencies are less explored. We conducted a large-scale gene-smoking interaction meta-analysis on exonic rare and low-frequency variants involving 44,429 individuals of European ancestry in the discovery stage and sought replication in the UK BiLEVE study with 45,133 European ancestry samples and UK Biobank study with 59,478 samples. We leveraged data on cigarette smoking, the major environmental risk factor for reduced lung function, by testing gene-by-smoking interaction effects only and simultaneously testing the genetic main effects and interaction effects. The most statistically significant signal that replicated was a previously reported low-frequency signal in GPR126, distinct from common variant associations in this gene. Although only nominal replication was obtained for a top rare variant signal rs142935352 in one of the two studies, interaction and joint tests for current smoking and PDE3B were significantly associated with FEV1. This study investigates the utility of assessing gene-by-smoking interactions and underscores their effects on potential pulmonary function.

Sugar-Sweetened Beverage Consumption May Modify Associations Between Genetic Variants in the CHREBP (Carbohydrate Responsive Element Binding Protein) Locus and HDL-C (High-Density Lipoprotein Cholesterol) and Triglyceride Concentrations.

BACKGROUND: ChREBP (carbohydrate responsive element binding protein) is a transcription factor that responds to sugar consumption. Sugar-sweetened beverage (SSB) consumption and genetic variants in the CHREBP locus have separately been linked to HDL-C (high-density lipoprotein cholesterol) and triglyceride concentrations. We hypothesized that SSB consumption would modify the association between genetic variants in the CHREBP locus and dyslipidemia. METHODS: Data from 11 cohorts from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium (N=63 599) and the UK Biobank (N=59 220) were used to quantify associations of SSB consumption, genetic variants, and their interaction on HDL-C and triglyceride concentrations using linear regression models. A total of 1606 single nucleotide polymorphisms within or near CHREBP were considered. SSB consumption was estimated from validated questionnaires, and participants were grouped by their estimated intake. RESULTS: In a meta-analysis, rs71556729 was significantly associated with higher HDL-C concentrations only among the highest SSB consumers (ß, 2.12 [95% CI, 1.16-3.07] mg/dL per allele; P<0.0001), but not significantly among the lowest SSB consumers (P=0.81; PDiff <0.0001). Similar results were observed for 2 additional variants (rs35709627 and rs71556736). For triglyceride, rs55673514 was positively associated with triglyceride concentrations only among the highest SSB consumers (ß, 0.06 [95% CI, 0.02-0.09] ln-mg/dL per allele, P=0.001) but not the lowest SSB consumers (P=0.84; PDiff=0.0005). CONCLUSIONS: Our results identified genetic variants in the CHREBP locus that may protect against SSB-associated reductions in HDL-C and other variants that may exacerbate SSB-associated increases in triglyceride concentrations. Registration: URL: https://www.clinicaltrials.gov; Unique identifier: NCT00005133, NCT00005121, NCT00005487, and NCT00000479.