Dissection of pleiotropic effects in genome-wide association studies of phenotypes related to cardiometabolic health

In the past seven years, Genome-Wide Association Studies (GWAS) have identified hundreds of variants associated with cardiometabolic quantitative traits and diseases. Many genetic loci appear to harbour variants associated with multiple phenotypes (cross-phenotype associations, CP). CP associations highlight that phenotypes may share common underlying genetic mechanisms that might, or might not, be consistent with epidemiological expectations and, therefore, add complexity to the relationships between human phenotypes. Pleiotropy occurs when the same genetic causal element affects more than one phenotype “in parallel” and can explain the presence of CP associations. It can appear at a single variant level, where a single causal variant is related to multiple phenotypes, or at a locus level, that is when multiple variants in the same gene or locus are associated with different phenotypes by affecting the same functional element. However, other potential genetic mechanisms, that can explain CP associations, exist. Among them, mediation occurs when a genetic variant is directly associated with a phenotype and that phenotype is itself causal for a second phenotype or more phenotypes; multi-phenotype allelic heterogeneity is a phenomenon which involves independent uncorrelated variants within the same locus which cause changes in multiple phenotypes, by affecting them through independent pathways related to distinct functional elements. The identification and characterisation of CP associations across the genome may help uncovering the mechanistic basis of physiological processes that underlie variability of cardiometabolic quantitative traits, and of pathogenetic processes leading to metabolic disorders. The definition of specific patterns of effect combinations on cardiometabolic phenotypes will highlight novel biological pathways, targets for translational research, for therapeutic intervention, and for the understanding of the pathophysiology of human metabolism. Based on this hypothesis, and in collaboration with the Cross-Consortia pleiotropy group and with the European Network for Genetic and Genomic Epidemiology (ENGAGE) consortium, my PhD project focused on dissection of CP effects, pleiotropy in particular, at common variants across the genome in association with cardiometabolic phenotypes. The objective was to improve our understanding of the extent of shared genetics between cardiometabolic phenotypes and of the influences of DNA sequence variation on risk of metabolic diseases, considering phenotypes as a range of inter-related manifestations of biological mechanisms rather than as isolated events. My research has been divided into three sub-projects: Project 1: Clustering and pathway analysis of univariate GWAS results for the detection of pleiotropic effects. We explored multi-phenotype effects at hundreds of established cardiometabolic genetic variants from published univariate GWAS meta-analyses on more than 20 respective phenotypes, by defining clusters of loci with similar multiple effects, comparing them to known epidemiological expectations, and identifying enriched biological networks within the most interesting groups of loci. Our results highlighted that many variants at cardiometabolic loci have multiple associations that characterise different aspects of metabolism. Cardiometabolic loci can be grouped according to their shared multi-phenotype effects and metabolic syndrome represents just one possible combination; in fact, several other unexpected combinations might be observed, for example healthy obesity/unhealthy leanness. We also highlighted that genetic loci with similar cardiometabolic effects are involved in shared biological pathways. Some of these may be expected, for instance, regulation of lipids metabolism or cholesterol transport for groups of loci with strong effects on lipids, and circulatory system processes for genes near blood pressure-association signals. Sometimes groups of loci affected fundamental cell functions, such as regulation of cellular processes, for the loci with effects on obesity and anthropometric traits. The enriched connectivity within pathway networks revealed new potential candidate genes and tissues of action that are more likely to have causal effect on phenotypes. Project 2: Validating pleiotropy and analysis of locus architecture in potential pleiotropic regions. We aimed to dissect the architecture of established cardiometabolic loci showing multiple associations for a better definition of the underlying mechanisms of multi-phenotype effects and for the discernment of potential pleiotropy from allelic heterogeneity. To this aim, we applied an approximate conditional analysis, based on observed linkage disequilibrium patterns, which led us to the discovery of multiple associations at adjacent variants that underlie the same genetic cause for variability of different phenotypes. Our results also highlighted that a substantial proportion of metabolic loci incorporate complex patterns of multi-phenotype allelic heterogeneity, thus suggesting an important contribution of this mechanism into cross-phenotype effects. Project 3: Application of a multivariate statistical approach for the study of pleiotropy within cardiometabolic phenotypes. We developed and applied a statistical strategy for joint multivariate analysis of multiple correlated phenotypes using individual genetic data from the ENGAGE consortium to discover new uncovered multiple associations and to follow-up GWAS meta-analysis at two loci, FTO and FADS1. Using this approach we were able to take into account correlation between phenotypes, and we achieved a boost in power; moreover, we improved precision of parameter estimates and of the identification of novel candidate genes. Our results allowed us to identify several variants jointly associated with multiple lipid traits and body mass index. Our approach was useful for the identification of mediation: we, in fact, confirmed mediation underlying causal relationship between adiposity and other cardiometabolic phenotypes at the FTO locus. Additionally, we demonstrated that multiple effects on cardiometabolic phenotypes attributable to the FADS1 locus are mediated by its independent, thus pleiotropic, effect on total cholesterol and triglycerides. In conclusion, we applied several statistical approaches which allowed dissecting suggestive CP effects and their mechanisms, including pleiotropy, mediation and allelic heterogeneity. Our analyses have demonstrated the complexity of the relationships between cardiometabolic phenotypes related to the variability of both, underlying genetic mechanisms and genetic loci architecture

Immunomediated and ischemia-independent inflammation of coronary microvessels in unstable angina.

This study investigated whether the myocardium is involved in the acute inflammatory reaction associated with bursts of unstable angina (UA). We looked for the presence of activated DR+ inflammatory cells and the expression patterns, localization, and immunostaining identification of genes for cytokines (IL-1beta, TNF-alpha, IL-6, and IFN-gamma), MCP-1, and iNOS in the left ventricle biopsies from 2-vessel disease anginal patients, 24 with UA and 12 with stable angina (SA), who underwent coronary bypass surgery. Biopsy specimens from 6 patients with mitral stenosis who underwent valve replacement were examined as control hearts (CHs). Plasma levels of IL-2 soluble receptor (sIL-2R) were measured as a marker of systemic immune reaction. In CHs, DR+ cells were undetectable, and cytokine and iNOS mRNA expression were negligible. UA patients had higher sIL-2R levels than SA patients (P<0.01), and their biopsy specimens showed both numerous DR+ cells identified as lymphocytes, macrophages, endothelial cells, and elevated expression levels of cytokine and iNOS genes (from 2.4- to 6.1-fold vs SA; P<0.01). Cytokine and iNOS genes and proteins were localized in endothelial cells without involvement of myocytes. IL-1beta and MCP-1 mRNAs were nearly undetectable. No significant differences were found in the number of DR+ cells, levels of cytokine, and iNOS genes between potentially ischemic and nonischemic left ventricle areas. In SA specimens, DR+ cells were very rare and only mRNAs for TNF-alpha and iNOS genes were overexpressed versus CHs. These results indicated that an acute immunomediated inflammatory reaction, essentially involving coronary microvessels, is demonstrable in UA patients

Genome-wide association study of offspring birth weight in 86 577 women identifies five novel loci and highlights maternal genetic effects that are independent of fetal genetics

Funding Information: Researchers were funded by investment from the European Regional Development Fund (ERDF) and the European Social Fund (ESF) Convergence Programme for Cornwall and the Isles of Scilly [J.T.]; European Research Council (ERC) [grant: SZ-245 50371-GLUCOSEGENES-FP7-IDEAS-ERC to T.M.F., A.R.W.], [ERC Consolidator Grant, ERC-2014-CoG-648916 to V.W.V.J.], [P.R.N.]; University of Bergen, KG Jebsen and Helse Vest [P.R.N.]; Wellcome Trust Senior Investigator Awards [A.T.H. (WT098395), M.I.M. (WT098381)]; National Institute for Health Research (NIHR) Senior Investigator Award (NF-SI-0611–10219); Sir Henry Dale Fellowship (Wellcome Trust and Royal Society grant: WT104150) [R.M.F., R.N.B.]; 4-year studentship (Grant Code: WT083431MF) [R.C.R]; the European Research Council under the European Union’s Seventh Framework Programme (FP/2007– 2013)/ERC Grant Agreement (grant number 669545; Develop Obese) [D.A.L.]; US National Institute of Health (grant: R01 DK10324) [D.A.L, C.L.R]; Wellcome Trust GWAS grant (WT088806) [D.A.L] and NIHR Senior Investigator Award (NF-SI-0611–10196) [D.A.L]; Wellcome Trust Institutional Strategic Support Award (WT097835MF) [M.A.T.]; The Diabetes Research and Wellness Foundation Non-Clinical Fellowship [J.T.]; Australian National Health and Medical Research Council Early Career Fellowship (APP1104818) [N.M.W.]; Daniel B. Burke Endowed Chair for Diabetes Research [S.F.A.G.]; UK Medical Research Council Unit grants MC_UU_12013_5 [R.C.R, L.P, S.R, C.L.R, D.M.E., D.A.L.] and MC_UU_12013_4 [D.M.E.]; Medical Research Council (grant: MR/M005070/1) [M.N.W., S.E.J.]; Australian Research Council Future Fellowship (FT130101709) [D.M.E] and (FT110100548) [S.E.M.]; NIHR Oxford Biomedical Research Centre (BRC); Oak Foundation Fellowship and Novo Nordisk Foundation (12955) [B.F.]; FRQS research scholar and Clinical Scientist Award by the Canadian Diabetes Association and the Maud Menten Award from the Institute of Genetics– Canadian Institute of Health Research (CIHR) [MFH]; CIHR— Frederick Banting and Charles Best Canada Graduate Scholarships [C.A.]; FRQS [L.B.]; Netherlands Organization for Health Research and Development (ZonMw–VIDI 016.136.361) [V.W.V.J.]; National Institute on Aging (R01AG29451) [J.M.M.]; 2010–2011 PRIN funds of the University of Ferrara—Holder: Prof. Guido Barbujani, Supervisor: Prof. Chiara Scapoli—and in part sponsored by the European Foundation for the Study of Diabetes (EFSD) Albert Renold Travel Fellowships for Young Scientists, ‘5 per mille’ contribution assigned to the University of Ferrara, income tax return year 2009 and the ENGAGE Exchange and Mobility Program for ENGAGE training funds, ENGAGE project, grant agreement HEALTH-F4–2007-201413 [L.M.]; ESRC (RES-060–23-0011) [C.L.R.]; National Institute of Health Research ([S.D., M.I.M.], Senior Investigator Award (NF-SI-0611–10196) [D.A.L]); Australian NHMRC Fellowships Scheme (619667) [G.W.M]. For study-specific funding, please see Supplementary Material. The views expressed are those of the authors and not necessarily those of the NHS, the NIHR or the Department of Health. Funding to pay the Open Access publication charges for this article was provided by the Charity Open Access Fund (COAF). Funding Information: We are extremely grateful to the participants and families who contributed to all of the studies and the teams of investigators involved in each one. These include interviewers, computer and laboratory technicians, clerical workers, research scientists, volunteers, managers, receptionists and nurses. This research has been conducted using the UK Biobank Resource (Application numbers 7036 and 12703). For additional study-specific acknowledgements, please see Supplementary Material. Conflict of Interest statement. D.A.L. has received support from Roche Diagnostics and Medtronic for biomarker research unrelated to the work presented here. Funding Researchers were funded by investment from the European Regional Development Fund (ERDF) and the European Social Fund (ESF) Convergence Programme for Cornwall and the Isles of Scilly [J.T.]; European Research Council (ERC) [grant: SZ-245 50371-GLUCOSEGENES-FP7-IDEAS-ERC to T.M.F., A.R.W.], [ERC Consolidator Grant, ERC-2014-CoG-648916 to V.W.V.J.], [P.R.N.]; University of Bergen, KG Jebsen and Helse Vest [P.R.N.]; Wellcome Trust Senior Investigator Awards [A.T.H. (WT098395), M.I.M. (WT098381)]; National Institute for Health Research (NIHR) Senior Investigator Award (NF-SI-0611-10219); Sir Henry Dale Fellowship (Wellcome Trust and Royal Society grant: WT104150) [R.M.F., R.N.B.]; 4-year studentship (Grant Code: WT083431MF) [R.C.R]; the European Research Council under the European Union's Seventh Framework Programme (FP/2007-2013)/ERC Grant Agreement (grant number 669545; Develop Obese) [D.A.L.]; US National Institute of Health (grant: R01 DK10324) [D.A.L, C.L.R]; Wellcome Trust GWAS grant (WT088806) [D.A.L] and NIHR Senior Investigator Award (NF-SI-0611-10196) [D.A.L]; Wellcome Trust Institutional Strategic Support Award (WT097835MF) [M.A.T.]; The Diabetes Research and Wellness Foundation Non-Clinical Fellowship [J.T.]; Australian National Health and Medical Research Council Early Career Fellowship (APP1104818) [N.M.W.]; Daniel B. Burke Endowed Chair for Diabetes Research [S.F.A.G.]; UK Medical Research Council Unit grants MC_UU_12013_5 [R.C.R, L.P, S.R, C.L.R, D.M.E., D.A.L.] and MC_UU_12013_4 [D.M.E.]; Medical Research Council (grant: MR/M005070/1) [M.N.W., S.E.J.]; Australian Research Council Future Fellowship (FT130101709) [D.M.E] and (FT110100548) [S.E.M.]; NIHR Oxford Biomedical Research Centre (BRC); Oak Foundation Fellowship and Novo Nordisk Foundation (12955) [B.F.]; FRQS research scholar and Clinical Scientist Award by the Canadian Diabetes Association and the Maud Menten Award from the Institute of Genetics-Canadian Institute of Health Research (CIHR) [MFH]; CIHR-Frederick Banting and Charles Best Canada Graduate Scholarships [C.A.]; FRQS [L.B.]; Netherlands Organization for Health Research and Development (ZonMw-VIDI 016.136.361) [V.W.V.J.]; National Institute on Aging (R01AG29451) [J.M.M.]; 2010-2011 PRIN funds of the University of Ferrara-Holder: Prof. Guido Barbujani, Supervisor: Prof. Chiara Scapoli-and in part sponsored by the European Foundation for the Study of Diabetes (EFSD) Albert Renold Travel Fellowships for Young Scientists, '5 per mille' contribution assigned to the University of Ferrara, income tax return year 2009 and the ENGAGE Exchange and Mobility Program for ENGAGE training funds, ENGAGE project, grant agreement HEALTH-F4-2007-201413 [L.M.]; ESRC (RES-060-23-0011) [C.L.R.]; National Institute of Health Research ([S.D., M.I.M.], Senior Investigator Award (NFSI-0611-10196) [D.A.L]); Australian NHMRC Fellowships Scheme (619667) [G.W.M]. For study-specific funding, please see Supplementary Material. The views expressed are those of the authors and not necessarily those of the NHS, the NIHR or the Department of Health. Funding to pay the Open Access publication charges for this article was provided by the Charity Open Access Fund (COAF). Publisher Copyright: © The Author(s) 2018.Genome-wide association studies of birth weight have focused on fetal genetics, whereas relatively little is known about the role of maternal genetic variation. We aimed to identify maternal genetic variants associated with birth weight that could highlight potentially relevant maternal determinants of fetal growth. We meta-analysed data on up to 8.7 million SNPs in up to 86 577 women of European descent from the Early Growth Genetics (EGG) Consortium and the UK Biobank. We used structural equation modelling (SEM) and analyses of mother-child pairs to quantify the separate maternal and fetal genetic effects. Maternal SNPs at 10 loci (MTNR1B, HMGA2, SH2B3, KCNAB1, L3MBTL3, GCK, EBF1, TCF7L2, ACTL9, CYP3A7) were associated with offspring birth weight at P<5 x 10(-8). In SEM analyses, at least 7 of the 10 associations were consistent with effects of the maternal genotype acting via the intrauterine environment, rather than via effects of shared alleles with the fetus. Variants, or correlated proxies, at many of the loci had been previously associated with adult traits, including fasting glucose (MTNR1B, GCK and TCF7L2) and sex hormone levels (CYP3A7), and one (EBF1) with gestational duration. The identified associations indicate that genetic effects on maternal glucose, cytochrome P450 activity and gestational duration, and potentially on maternal blood pressure and immune function, are relevant for fetal growth. Further characterization of these associations in mechanistic and causal analyses will enhance understanding of the potentially modifiable maternal determinants of fetal growth, with the goal of reducing the morbidity and mortality associated with low and high birth weights.Peer reviewe

Discovery and Fine-Mapping of Glycaemic and Obesity-Related Trait Loci Using High-Density Imputation

Reference panels from the 1000 Genomes (1000G) Project Consortium provide near complete coverage of common and low-frequency genetic variation with minor allele frequency ≥0.5% across European ancestry populations. Within the European Network for Genetic and Genomic Epidemiology (ENGAGE) Consortium, we have undertaken the fi

Pleiotropic genes for metabolic syndrome and inflammation

Metabolic syndrome (MetS) has become a health and financial burden worldwide. The MetS definition captures clustering of risk factors that predict higher risk for diabetes mellitus and cardiovascular disease. Our study hypothesis is that additional to genes influencing individual MetS risk factors, genetic variants exist that influence MetS and inflammatory markers forming a predisposing MetS genetic network. To test this hypothesis a staged approach was undertaken. (a) We analyzed 17 metabolic and inflammatory traits in more than 85,500 participants from 14 large epidemiological studies within the Cross Consortia Pleiotropy Group. Individuals classified with MetS (NCEP definition), versus those without, showed on average significantly different levels for most inflammatory markers studied. (b) Paired average correlations between 8 metabolic traits and 9 inflammatory markers from the same studies as above, estimated with two methods, and factor analyses on large simulated data, helped in identifying 8 combinations of traits for follow-up in meta-analyses, out of 130,305 possible combinations between metabolic traits and inflammatory markers studied. (c) We performed correlated meta-analyses for 8 metabolic traits and 6 inflammatory markers by using existing GWAS published genetic summary results, with about 2.5 million SNPs from twelve predominantly largest GWAS consortia. These analyses yielded 130 unique SNPs/genes with pleiotropic associations (a SNP/gene associating at least one metabolic trait and one inflammatory marker). Of them twenty-five variants (seven loci newly reported) are proposed as MetS candidates. They map to genes MACF1, KIAA0754, GCKR, GRB14, COBLL1, LOC646736-IRS1, SLC39A8, NELFE, SKIV2L, STK19, TFAP2B, BAZ1B, BCL7B, TBL2, MLXIPL, LPL, TRIB1, ATXN2, HECTD4, PTPN11, ZNF664, PDXDC1, FTO, MC4R and TOMM40. Based on large data evidence, we conclude that inflammation is a feature of MetS and several gene variants show pleiotropic genetic associations across phenotypes and might explain a part of MetS correlated genetic architecture. These findings warrant further functional investigation

Maternal and fetal genetic effects on birth weight and their relevance to cardio-metabolic risk factors.

Birth weight variation is influenced by fetal and maternal genetic and non-genetic factors, and has been reproducibly associated with future cardio-metabolic health outcomes. In expanded genome-wide association analyses of own birth weight (n = 321,223) and offspring birth weight (n = 230,069 mothers), we identified 190 independent association signals (129 of which are novel). We used structural equation modeling to decompose the contributions of direct fetal and indirect maternal genetic effects, then applied Mendelian randomization to illuminate causal pathways. For example, both indirect maternal and direct fetal genetic effects drive the observational relationship between lower birth weight and higher later blood pressure: maternal blood pressure-raising alleles reduce offspring birth weight, but only direct fetal effects of these alleles, once inherited, increase later offspring blood pressure. Using maternal birth weight-lowering genotypes to proxy for an adverse intrauterine environment provided no evidence that it causally raises offspring blood pressure, indicating that the inverse birth weight-blood pressure association is attributable to genetic effects, and not to intrauterine programming.The Fenland Study is funded by the Medical Research Council (MC_U106179471) and Wellcome Trust

Genome-wide Association Study of Change in Fasting Glucose over time in 13,807 non-diabetic European Ancestry Individuals

Type 2 diabetes (T2D) affects the health of millions of people worldwide. The identification of genetic determinants associated with changes in glycemia over time might illuminate biological features that precede the development of T2D. Here we conducted a genome-wide association study of longitudinal fasting glucose changes in up to 13,807 non-diabetic individuals of European descent from nine cohorts. Fasting glucose change over time was defined as the slope of the line defined by multiple fasting glucose measurements obtained over up to 14 years of observation. We tested for associations of genetic variants with inverse-normal transformed fasting glucose change over time adjusting for age at baseline, sex, and principal components of genetic variation. We found no genome-wide significant association (P < 5 x 10(-8)) with fasting glucose change over time. Seven loci previously associated with T2D, fasting glucose or HbA1c were nominally (P < 0.05) associated with fasting glucose change over time. Limited power influences unambiguous interpretation, but these data suggest that genetic effects on fasting glucose change over time are likely to be small. A public version of the data provides a genomic resource to combine with future studies to evaluate shared genetic links with T2D and other metabolic risk traits.Peer reviewe

Discovery and Fine-Mapping of Glycaemic and Obesity-Related Trait Loci Using High-Density Imputation

Reference panels from the 1000 Genomes (1000G) Project Consortium provide near complete coverage of common and low-frequency genetic variation with minor allele frequency >= 0.5% across European ancestry populations. Within the European Network for Genetic and Genomic Epidemiology (ENGAGE) Consortium, we have undertaken the first large-scale meta-analysis of genome-wide association studies (GWAS), supplemented by 1000G imputation, for four quantitative glycaemic and obesity-related traits, in up to 87,048 individuals of European ancestry. We identified two loci for body mass index (BMI) at genome-wide significance, and two for fasting glucose (FG), none of which has been previously reported in larger meta-analysis efforts to combine GWAS of European ancestry. Through conditional analysis, we also detected multiple distinct signals of association mapping to established loci for waist-hip ratio adjusted for BMI (RSPO3) and FG (GCK and G6PC2). The index variant for one association signal at the G6PC2 locus is a low-frequency coding allele, H177Y, which has recently been demonstrated to have a functional role in glucose regulation. Fine-mapping analyses revealed that the non-coding variants most likely to drive association signals at established and novel loci were enriched for overlap with enhancer elements, which for FG mapped to promoter and transcription factor binding sites in pancreatic islets, in particular. Our study demonstrates that 1000G imputation and genetic fine-mapping of common and low-frequency variant association signals at GWAS loci, integrated with genomic annotation in relevant tissues, can provide insight into the functional and regulatory mechanisms through which their effects on glycaemic and obesity-related traits are mediated.Peer reviewe

New genetic loci link adipose and insulin biology to body fat distribution.

Body fat distribution is a heritable trait and a well-established predictor of adverse metabolic outcomes, independent of overall adiposity. To increase our understanding of the genetic basis of body fat distribution and its molecular links to cardiometabolic traits, here we conduct genome-wide association meta-analyses of traits related to waist and hip circumferences in up to 224,459 individuals. We identify 49 loci (33 new) associated with waist-to-hip ratio adjusted for body mass index (BMI), and an additional 19 loci newly associated with related waist and hip circumference measures (P < 5 × 10(-8)). In total, 20 of the 49 waist-to-hip ratio adjusted for BMI loci show significant sexual dimorphism, 19 of which display a stronger effect in women. The identified loci were enriched for genes expressed in adipose tissue and for putative regulatory elements in adipocytes. Pathway analyses implicated adipogenesis, angiogenesis, transcriptional regulation and insulin resistance as processes affecting fat distribution, providing insight into potential pathophysiological mechanisms