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... Show moreTo 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. Show less
We evaluate the shared genetic regulation of mRNA molecules, proteins and metabolites derived from whole blood from 3029 human donors. We find abundant allelic heterogeneity, where multiple... Show moreWe evaluate the shared genetic regulation of mRNA molecules, proteins and metabolites derived from whole blood from 3029 human donors. We find abundant allelic heterogeneity, where multiple variants regulate a particular molecular phenotype, and pleiotropy, where a single variant associates with multiple molecular phenotypes over multiple genomic regions. The highest proportion of share genetic regulation is detected between gene expression and proteins (66.6%), with a further median shared genetic associations across 49 different tissues of 78.3% and 62.4% between plasma proteins and gene expression. We represent the genetic and molecular associations in networks including 2828 known GWAS variants, showing that GWAS variants are more often connected to gene expression in trans than other molecular phenotypes in the network. Our work provides a roadmap to understanding molecular networks and deriving the underlying mechanism of action of GWAS variants using different molecular phenotypes in an accessible tissue. Show less
Vegte, Y. van de; Eppinga, R.N.; Ende, M.Y. van der; Hagemeijer, Y.; Mahendran, Y.V.; Salfati, E.Y.; ... ; DCCT EDIC Res Grp 2023
The genetics and clinical consequences of resting heart rate (RHR) remain incompletely understood. Here, the authors discover new genetic variants associated with RHR and find that higher... Show moreThe genetics and clinical consequences of resting heart rate (RHR) remain incompletely understood. Here, the authors discover new genetic variants associated with RHR and find that higher genetically predicted RHR decreases risk of atrial fibrillation and ischemic stroke.Resting heart rate is associated with cardiovascular diseases and mortality in observational and Mendelian randomization studies. The aims of this study are to extend the number of resting heart rate associated genetic variants and to obtain further insights in resting heart rate biology and its clinical consequences. A genome-wide meta-analysis of 100 studies in up to 835,465 individuals reveals 493 independent genetic variants in 352 loci, including 68 genetic variants outside previously identified resting heart rate associated loci. We prioritize 670 genes and in silico annotations point to their enrichment in cardiomyocytes and provide insights in their ECG signature. Two-sample Mendelian randomization analyses indicate that higher genetically predicted resting heart rate increases risk of dilated cardiomyopathy, but decreases risk of developing atrial fibrillation, ischemic stroke, and cardio-embolic stroke. We do not find evidence for a linear or non-linear genetic association between resting heart rate and all-cause mortality in contrast to our previous Mendelian randomization study. Systematic alteration of key differences between the current and previous Mendelian randomization study indicates that the most likely cause of the discrepancy between these studies arises from false positive findings in previous one-sample MR analyses caused by weak-instrument bias at lower P-value thresholds. The results extend our understanding of resting heart rate biology and give additional insights in its role in cardiovascular disease development. Show less
This case study explores the applicability of transcriptome data to characterize a common mechanism of action within groups of short-chain aliphatic alpha-, beta-, and gamma-diketones. Human... Show moreThis case study explores the applicability of transcriptome data to characterize a common mechanism of action within groups of short-chain aliphatic alpha-, beta-, and gamma-diketones. Human reference in vivo data indicate that the alpha-diketone diacetyl induces bronchiolitis obliterans in workers involved in the preparation of microwave popcorn. The other three alpha-diketones induced inflammatory responses in preclinical in vivo animal studies, whereas beta and gamma diketones in addition caused neuronal effects. We investigated early transcriptional responses in primary human bronchiolar (PBEC) cell cultures after 24 h and 72 h of air-liquid exposure. Differentially expressed genes (DEGs) were assessed based on transcriptome data generated with the EUToxRisk gene panel of Temp-O-Seq (R). For each individual substance, genes were identified displaying a consistent differential expression across dose and exposure duration. The log fold change values of the DEG profiles indicate that alpha- and beta-diketones are more active compared to gamma-diketones. alpha-diketones in particular showed a highly concordant expression pattern, which may serve as a first indication of the shared mode of action. In order to gain a better mechanistic understanding, the resultant DEGs were submitted to a pathway analysis using ConsensusPathDB. The four alpha-diketones showed very similar results with regard to the number of activated and shared pathways. Overall, the number of signaling pathways decreased from alpha-to beta-to gamma-diketones. Additionally, we reconstructed networks of genes that interact with one another and are associated with different adverse outcomes such as fibrosis, inflammation or apoptosis using the TRANSPATH-database. Transcription factor enrichment and upstream analyses with the geneXplain platform revealed highly interacting gene products (called master regulators, MRs) per case study compound. The mapping of the resultant MRs on the reconstructed networks, visualized similar gene regulation with regard to fibrosis, inflammation and apoptosis. This analysis showed that transcriptome data can strengthen the similarity assessment of compounds, which is of particular importance, e.g., in read-across approaches. It is one important step towards grouping of compounds based on biological profiles. Show less
The 3-dimensional spatial and 2-dimensional frontal QRS-T angles are measures derived from the vectorcardiogram. They are independent risk predictors for arrhythmia, but the underlying biology is... Show moreThe 3-dimensional spatial and 2-dimensional frontal QRS-T angles are measures derived from the vectorcardiogram. They are independent risk predictors for arrhythmia, but the underlying biology is unknown. Using multi-ancestry genome-wide association studies we identify 61 (58 previously unreported) loci for the spatial QRS-T angle (N=118,780) and 11 for the frontal QRS-T angle (N=159,715). Seven out of the 61 spatial QRS-T angle loci have not been reported for other electrocardiographic measures. Enrichments are observed in pathways related to cardiac and vascular development, muscle contraction, and hypertrophy. Pairwise genome-wide association studies with classical ECG traits identify shared genetic influences with PR interval and QRS duration. Phenome-wide scanning indicate associations with atrial fibrillation, atrioventricular block and arterial embolism and genetically determined QRS-T angle measures are associated with fascicular and bundle branch block (and also atrioventricular block for the frontal QRS-T angle). We identify potential biology involved in the QRS-T angle and their genetic relationships with cardiovascular traits and diseases, may inform future research and risk prediction. The spatial and frontal QRS-T angles are electrocardiographic (ECG) predictors for arrhythmia. This work used genetic analyses to identify associated loci and pathways, and explore their relationships with other ECG traits and cardiovascular disease. Show less
Hengst, J.M.A. van; Hellemons, R.J.C.; Remmerswaal, W.A.; Vrande, K.N.A. van de; Hansen, T.; Vorm, S. van der; ... ; Codee, J.D.C. 2023
The reactivity of the acceptor alcohol can have a tremendous influence on the outcome of a glycosylation reaction, both in terms of yield and stereoselectivity. Through a systematic survey of 67... Show moreThe reactivity of the acceptor alcohol can have a tremendous influence on the outcome of a glycosylation reaction, both in terms of yield and stereoselectivity. Through a systematic survey of 67 acceptor alcohols in glycosylation reactions with two glucosyl donors we here reveal how the reactivity of a carbohydrate acceptor depends on its configuration and substitution pattern. The study shows how the functional groups flanking the acceptor alcohol influence the reactivity of the alcohol and show that both the nature and relative orientation play an essential role. The empiric acceptor reactivity guidelines revealed here will aid in the rational optimization of glycosylation reactions and be an important tool in the assembly of oligosaccharides. Show less
The application of multiple omics technologies in biomedical cohorts has the potential to reveal patient-level disease characteristics and individualized response to treatment. However, the scale... Show moreThe application of multiple omics technologies in biomedical cohorts has the potential to reveal patient-level disease characteristics and individualized response to treatment. However, the scale and heterogeneous nature of multi-modal data makes integration and inference a non-trivial task. We developed a deep-learning-based framework, multi-omics variational autoencoders (MOVE), to integrate such data and applied it to a cohort of 789 people with newly diagnosed type 2 diabetes with deep multi-omics phenotyping from the DIRECT consortium. Using in silico perturbations, we identified drug-omics associations across the multi-modal datasets for the 20 most prevalent drugs given to people with type 2 diabetes with substantially higher sensitivity than univariate statistical tests. From these, we among others, identified novel associations between metformin and the gut microbiota as well as opposite molecular responses for the two statins, simvastatin and atorvastatin. We used the associations to quantify drug-drug similarities, assess the degree of polypharmacy and conclude that drug effects are distributed across the multi-omics modalities. Show less
Background: Genetic variants within nearly 1000 loci are known to contribute to modulation of blood lipid levels. However, the biological pathways underlying these associations are frequently... Show moreBackground: Genetic variants within nearly 1000 loci are known to contribute to modulation of blood lipid levels. However, the biological pathways underlying these associations are frequently unknown, limiting understanding of these findings and hindering downstream translational efforts such as drug target discovery. Results: To expand our understanding of the underlying biological pathways and mechanisms controlling blood lipid levels, we leverage a large multi-ancestry meta-analysis (N=1,654,960) of blood lipids to prioritize putative causal genes for 2286 lipid associations using six gene prediction approaches. Using phenome-wide association (PheWAS) scans, we identify relationships of genetically predicted lipid levels to other diseases and conditions. We confirm known pleiotropic associations with cardiovascular phenotypes and determine novel associations, notably with cholelithiasis risk. We perform sex-stratified GWAS meta-analysis of lipid levels and show that 3-5% of autosomal lipid-associated loci demonstrate sex-biased effects. Finally, we report 21 novel lipid loci identified on the X chromosome. Many of the sex-biased autosomal and X chromosome lipid loci show pleiotropic associations with sex hormones, emphasizing the role of hormone regulation in lipid metabolism. Conclusions: Taken together, our findings provide insights into the biological mechanisms through which associated variants lead to altered lipid levels and potentially cardiovascular disease risk. Show less
Oxocarbenium ions are key reactive intermediates in organic chemistry. To generate a series of structure-reactivity-stereoselectivity principles for these species, we herein investigated the... Show moreOxocarbenium ions are key reactive intermediates in organic chemistry. To generate a series of structure-reactivity-stereoselectivity principles for these species, we herein investigated the bimolecular electrophilic substitution reactions (SE2') between allyltrimethylsilane and a series of archetypal six-membered ring oxocarbenium ions using a combined density functional theory (DFT) and coupled-cluster theory approach. These reactions preferentially proceed following a reaction path where the oxocarbenium ion transforms from a half chair (3H4 or 4H3) to a chair conformation. The introduction of alkoxy substituents on six-membered ring oxocarbenium ions, dramatically influences the conformational preference of the canonical 3H4 and 4H3 conformers, and thereby the stereochemical outcome of the SE2' reaction. In general, we find that the stereoselectivity in the reactions correlates to the "intrinsic preference" of the cations, as dictated by their shape. However, for the C5-CH2OMe substituent, steric factors override the "intrinsic preference", showing a more selective reaction than expected based on the shape of the ion. Our SE2' energetics correlate well with experimentally observed stereoselectivity, and the use of the activation strain model has enabled us to quantify important interactions and structural features that occur in the transition state of the reactions to precisely understand the relative energy barriers of the diastereotopic addition reactions. The fundamental mechanistic insight provided in this study will aid in understanding the reactivity of more complex glycosyl cations featuring multiple substituents and will facilitate our general understanding of glycosylation reactions. Show less
Uronic acids are carbohydrates carrying a terminal carboxylic acid and have a unique reactivity in stereoselective glycosylation reactions. Herein, the competing intramolecular stabilization of... Show moreUronic acids are carbohydrates carrying a terminal carboxylic acid and have a unique reactivity in stereoselective glycosylation reactions. Herein, the competing intramolecular stabilization of uronic acid cations by the C-5 carboxylic acid or the C-4 acetyl group was studied with infrared ion spectroscopy (IRIS). IRIS reveals that a mixture of bridged ions is formed, in which the mixture is driven towards the C-1,C-5 dioxolanium ion when the C-5,C-2-relationship is cis, and towards the formation of the C-1,C-4 dioxepanium ion when this relation is trans. Isomer-population analysis and interconversion barrier computations show that the two bridged structures are not in dynamic equilibrium and that their ratio parallels the density functional theory computed stability of the structures. These studies reveal how the intrinsic interplay of the different functional groups influences the formation of the different regioisomeric products. Show less
Demkiw, K.M.; Remmerswaal, W.A.; Hansen, T.; Marel, G.A. van der; Codee, J.D.C.; Woerpel, K.A. 2022
Uronic acids are carbohydrates carrying a terminal carboxylic acid and have a unique reactivity in stereoselective glycosylation reactions. Herein, the competing intramolecular stabilization of... Show moreUronic acids are carbohydrates carrying a terminal carboxylic acid and have a unique reactivity in stereoselective glycosylation reactions. Herein, the competing intramolecular stabilization of uronic acid cations by the C-5 carboxylic acid or the C-4 acetyl group was studied with infrared ion spectroscopy (IRIS). IRIS reveals that a mixture of bridged ions is formed, in which the mixture is driven towards the C-1,C-5 dioxolanium ion when the C-5,C-2-relationship is cis, and towards the formation of the C-1,C-4 dioxepanium ion when this relation is trans. Isomer-population analysis and interconversion barrier computations show that the two bridged structures are not in dynamic equilibrium and that their ratio parallels the density functional theory computed stability of the structures. These studies reveal how the intrinsic interplay of the different functional groups influences the formation of the different regioisomeric products. Show less
A major challenge of genome-wide association studies (GWASs) is to translate phenotypic associations into biological insights. Here, we integrate a large GWAS on blood lipids involving 1.6 million... Show moreA major challenge of genome-wide association studies (GWASs) is to translate phenotypic associations into biological insights. Here, we integrate a large GWAS on blood lipids involving 1.6 million individuals from five ancestries with a wide array of functional genomic datasets to discover regulatory mechanisms underlying lipid associations. We first prioritize lipid-associated genes with expression quantitative trait locus (eQTL) colocalizations and then add chromatin interaction data to narrow the search for functional genes. Polygenic enrichment analysis across 697 annotations from a host of tissues and cell types confirms the central role of the liver in lipid levels and highlights the selective enrichment of adipose-specific chromatin marks in high-density lipoprotein cholesterol and triglycerides. Overlapping transcription factor (TF) binding sites with lipid-associated loci identifies TFs relevant in lipid biology. In addition, we present an integrative framework to prioritize causal variants at GWAS loci, producing a comprehensive list of candidate causal genes and variants with multiple layers of functional evidence. We highlight two of the prioritized genes, CREBRF and RRBP1, which show convergent evidence across functional datasets supporting their roles in lipid biology. Show less
The development of small molecule covalent inhibitors and probes continuously pushes the rapidly evolving field of chemical biology forward. A key element in these molecular tool compounds is the ... Show moreThe development of small molecule covalent inhibitors and probes continuously pushes the rapidly evolving field of chemical biology forward. A key element in these molecular tool compounds is the "electrophilic trap" that allows for a covalent linkage with the target enzyme. The reactivity of this entity needs to be well balanced to effectively trap the desired enzyme, while not being attacked by off-target nucleophiles. We here investigate the intrinsic reactivity of substrates containing a class of widely used electrophilic traps, the three-membered heterocycles with an N- (aziridine), P- (phosphirane), O- (epoxide) and S-atom (thiirane) as heteroatom. Using quantum chemical approaches, we studied the conformational flexibility and nucleophilic ring-opening reaction of a series of model substrates, in which these electrophilic traps are mounted on a cyclohexene scaffold (C6H10Y with Y = NH, PH, O, S). It is revealed that the activation energy of the ring-opening reaction does not necessarily follow the trend that is expected from C-Y leaving-group bond strength, but steeply decreases from NH, to PH, to O, to S. We illustrate that the HOMONu-LUMOSubstrate interaction is an all-important factor for the observed reactivity. In addition, we show that the activation energy of aziridines and phosphiranes can be tuned far below that of the corresponding epoxides and thiiranes by the addition of proper electron-withdrawing ring substituents. Our results provide mechanistic insights to rationally tune the reactivity of this class of popular electrophilic traps and can guide the experimental design of covalent inhibitors and probes for enzymatic activity. Show less
Remmerswaal, W.A.; Houthuijs, K.J.; Ven, R. van de; Elferink, H.; Hansen, T.; Berden, G.; ... ; Codee, J.D.C. 2022
The stereoselective introduction of glycosidic bonds is of paramount importance to oligosaccharide synthesis. Among the various chemical strategies to steer stereoselectivity, participation by... Show moreThe stereoselective introduction of glycosidic bonds is of paramount importance to oligosaccharide synthesis. Among the various chemical strategies to steer stereoselectivity, participation by either neighboring or distal acyl groups is used particularly often. Recently, the use of the 2,2-dimethyl-2-(ortho-nitrophenyl)acetyl (DMNPA) protection group was shown to offer enhanced stereoselective steering compared to other acyl groups. Here, we investigate the origin of the stereoselectivity induced by the DMNPA group through systematic glycosylation reactions and infrared ion spectroscopy (IRIS) combined with techniques such as isotopic labeling of the anomeric center and isomer population analysis. Our study indicates that the origin of the DMNPA stereoselectivity does not lie in the direct participation of the nitro moiety but in the formation of a dioxolenium ion that is strongly stabilized by the nitro group. Show less
We have studied the palladium-mediated activation of C(sp(n))-X bonds (n = 1-3 and X = H, CH3, Cl) in archetypal model substrates H3C-CH2-X, H2C=CH-X and HC equivalent to C-X by catalysts PdLn with... Show moreWe have studied the palladium-mediated activation of C(sp(n))-X bonds (n = 1-3 and X = H, CH3, Cl) in archetypal model substrates H3C-CH2-X, H2C=CH-X and HC equivalent to C-X by catalysts PdLn with L-n = no ligand, Cl-, and (PH3)(2), using relativistic density functional theory at ZORA-BLYP/TZ2P. The oxidative addition barrier decreases along this series, even though the strength of the bonds increases going from C(sp(3))-X, to C(sp(2))-X, to C(sp)-X. Activation strain and matching energy decomposition analyses reveal that the decreased oxidative addition barrier going from sp(3), to sp(2), to sp, originates from a reduction in the destabilizing steric (Pauli) repulsion between catalyst and substrate. This is the direct consequence of the decreasing coordination number of the carbon atom in C(sp(n))-X, which goes from four, to three, to two along this series. The associated net stabilization of the catalyst-substrate interaction dominates the trend in strain energy which indeed becomes more destabilizing along this same series as the bond becomes stronger from C(sp(3))-X to C(sp)-X. Show less
The presentation and underlying pathophysiology of type 2 diabetes (T2D) is complex and heterogeneous. Recent studies attempted to stratify T2D into distinct subgroups using data-driven approaches,... Show moreThe presentation and underlying pathophysiology of type 2 diabetes (T2D) is complex and heterogeneous. Recent studies attempted to stratify T2D into distinct subgroups using data-driven approaches, but their clinical utility may be limited if categorical representations of complex phenotypes are suboptimal. We apply a soft-clustering (archetype) method to characterize newly diagnosed T2D based on 32 clinical variables. We assign quantitative clustering scores for individuals and investigate the associations with glycemic deterioration, genetic risk scores, circulating omics biomarkers, and phenotypic stability over 36 months. Four archetype profiles represent dysfunction patterns across combinations of T2D etiological processes and correlate with multiple circulating biomarkers. One archetype associated with obesity, insulin resistance, dyslipidemia, and impaired 1 beta cell glucose sensitivity corresponds with the fastest disease progression and highest demand for anti-diabetic treatment. We demonstrate that clinical heterogeneity in T2D can be mapped to heterogeneity in individual etiological processes, providing a potential route to personalized treatments. Show less
