Partitioning the Pleiotropy Between Coronary Artery Disease and Body Mass Index Reveals the Importance of Low Frequency Variants and Central Nervous System-Specific Functional Elements

Subgroup Identification with Gene Expression Profiles of Adipose Tissue in Patients with Coronary Artery Disease

Among many diseases, coronary artery disease (CAD) is the primary cause of mortality and morbidity worldwide. With the aim of revealing the underlying genetic characteristics of the CAD subtypes, we recruited patients with CAD and categorized them into subgroups according to the transcriptome expression profiles of the adipose tissue.With the removal of the batch effect, consensus clustering was employed to determine the subgroup numbers. Subgroup-specific genes were determined to conduct analysis of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). Weighted gene co-expression network analysis (WGCNA) revealed the subgroup-specific WGCNA modules. Moreover, gene set enrichment analysis (GSEA) was conducted. Overrepresentation enrichment analysis (OEA) of subgroup-specific signatures was also conducted to reveal the significant gene module associated with the corresponding clinical characteristics.After the removal of the batch effect, 77 CAD objects were divided into three subgroups. It was observed that the patients in subgroup III tended to be fat. After analyzing the dominant pathways of each subgroup, we discovered that the protein digestion and absorption pathway was specifically upregulated in subgroup I, which might result from the lowest proportion of the epicardial adipose tissue (EAT) sample. Moreover, subgroup II patients had genetic characteristics of high expression of complement and coagulation cascades and TNF signaling pathway. Furthermore, Th17 cell differentiation was significantly upregulated in subgroup III, indicating that Th17 cell differentiation is related to the clinical characteristics of body mass index (BMI).In conclusion, the genetic classification of CAD subjects indicated that subjects from different subgroups may exhibit specific gene expression patterns, suggesting that more personalized treatment should be applied to patients in each subgroup.

PDX1 and MC4R genetic polymorphisms are associated with type 2 diabetes mellitus risk in the Chinese Han population

BACKGROUND

Diabetes mellitus (DM) is a complex metabolic disease that is caused by a complex interplay between genetic and environmental factors. This research aimed to investigate the association of genetic polymorphisms in PDX1 and MC4R with T2DM risk.

METHODS

The genotypes of 10 selected SNPs in PDX1 and MC4R were identified using the Agena MassARRAY platform. We utilized odds ratio (OR) and 95% confidence intervals (CIs) to assess the correlation between genetic polymorphisms and T2DM risk.

RESULTS

We found that PDX1-rs9581943 decreased susceptibility to T2DM among in a Chinese Han population (OR = 0.76, p = 0.045). We also found that selected genetic polymorphisms in PDX1 and MC4R could modify the risk of T2DM, which might also be influenced by age, sex, BMI, smoking status, and drinking status (p < 0.05).

CONCLUSIONS

We concluded that PDX1 and MC4R genetic variants were significantly associated with T2DM risk in a Chinese Han population. These single polymorphic markers may be considered to be new targets in the assessment and prevention of T2DM among Chinese Han people.

Heterozygous Nme7 Mutation Affects Glucose Tolerance in Male Rats

Complex metabolic conditions such as type 2 diabetes and obesity result from the interaction of numerous genetic and environmental factors. While the family of Nme proteins has been connected so far mostly to development, proliferation, or ciliary functions, several lines of evidence from human and experimental studies point to the potential involvement of one of its members, NME7 (non-metastatic cells 7, nucleoside diphosphate kinase 7) in carbohydrate and lipid metabolism. As a complete lack of Nme7 is semilethal in rats, we compared morphometric, metabolic, and transcriptomic profiles of standard diet-fed heterozygous Nme7+/− on male rats vs. their wild-type Nme7+/+ controls. Nme7+/− animals showed increased body weight, adiposity, higher insulin levels together with decreased glucose tolerance. Moreover, they displayed pancreatic islet fibrosis and kidney tubular damage. Despite no signs of overt liver steatosis or dyslipidemia, we found significant changes in the hepatic transcriptome of Nme7+/− male rats with a concerted increase of expression of lipogenic enzymes including Scd1, Fads1, Dhcr7 and a decrease of Cyp7b1 and Nme7. Network analyses suggested possible links between Nme7 and the activation of Srebf1 and Srebf2 upstream regulators. These results further support the implication of NME7 in the pathogenesis of glucose intolerance and adiposity.

A unified framework identifies new links between plasma lipids and diseases from electronic medical records across large-scale cohorts

Plasma lipids are known heritable risk factors for cardiovascular disease, but increasing evidence also supports shared genetics with diseases of other organ systems. We devised a comprehensive three-phase framework to identify new lipid-associated genes and study the relationships among lipids, genotypes, gene expression and hundreds of complex human diseases from the Electronic Medical Records and Genomics (347 traits) and the UK Biobank (549 traits). Aside from 67 new lipid-associated genes with strong replication, we found evidence for pleiotropic SNPs/genes between lipids and diseases across the phenome. These include discordant pleiotropy in the HLA region between lipids and multiple sclerosis and putative causal paths between triglycerides and gout, among several others. Our findings give insights into the genetic basis of the relationship between plasma lipids and diseases on a phenome-wide scale and can provide context for future prevention and treatment strategies.

Convergence of biomarkers and risk factor trait loci of coronary artery disease at 3p21.31 and HLA region

Here we seek to identify molecular biomarkers that mediate the effect of risk factors on coronary artery disease (CAD). We perform a SNP-based multiomics data analysis to find biomarkers (probes) causally associated with the risk of CAD within known genomic loci for its risk factors. We identify 78 biomarkers, the majority (64%) of which are methylation probes. We detect the convergence of several CNS and lifestyle trait loci and their biomarkers at the 3p21.31 and human leukocyte antigen (HLA) regions. The 3p21.31 locus was the most populated region in the convergence of biomarkers and risk factors. In this region, we noted as the BSN gene becomes methylated the level of stomatin (STOM) in blood increases and this contributes to higher risk of CAD. In the HLA locus, we identify several methylation biomarkers associated with various CAD risk factors. SNPs in the CFB gene display a trans-regulatory impact on the GRIA4 protein level. A methylation site upstream of the APOE gene is associated with a higher protein level of S100A13 which in turn leads to higher LDL-C and greater CAD risk. We find UHRF1BP1 and ILRUN mediate the effect of obesity on CAD whereas methylation sites within NOS3 and CKM mediate the effect of their associated-risk factors on CAD. This study provides further insight into the biology of CAD and identifies a list of biomarkers that mediate the impact of risk factors on CAD. A SNP-based initiative can unite data from various fields of omics into a single network of knowledge.

RIPK1 gene variants associate with obesity in humans and can be therapeutically silenced to reduce obesity in mice

Obesity is a major public health burden worldwide and is characterized by chronic low-grade inflammation driven by the cooperation of the innate immune system and dysregulated metabolism in adipose tissue and other metabolic organs. Receptor-interacting serine/threonine-protein kinase 1 (RIPK1) is a central regulator of inflammatory cell function that coordinates inflammation, apoptosis and necroptosis in response to inflammatory stimuli. Here we show that genetic polymorphisms near the human RIPK1 locus associate with increased RIPK1 gene expression and obesity. We show that one of these single nucleotide polymorphisms is within a binding site for E4BP4 and increases RIPK1 promoter activity and RIPK1 gene expression in adipose tissue. Therapeutic silencing of RIPK1 in vivo in a mouse model of diet-induced obesity dramatically reduces fat mass, total body weight and improves insulin sensitivity, while simultaneously reducing macrophage and promoting invariant natural killer T cell accumulation in adipose tissue. These findings demonstrate that RIPK1 is genetically associated with obesity, and reducing RIPK1 expression is a potential therapeutic approach to target obesity and related diseases.

Multiomics Screening Identifies Molecular Biomarkers Causally Associated With the Risk of Coronary Artery Disease

BACKGROUND

In this study, we aimed to investigate functional mechanisms underlying coronary artery disease (CAD) loci and find molecular biomarkers for CAD.

METHODS

We devised a multiomics data analysis approach based on Mendelian randomization and utilized it to search for molecular biomarkers causally associated with the risk of CAD within genomic regions known to be associated with CAD.

RESULTS

Through our CAD-centered multiomics data analysis approach, we identified 33 molecular biomarkers (probes) that were causally associated with the risk of CAD. The majority of these (N=19) were methylation probes; moreover, methylation was often behind the causal effect of expression/protein probes. We identified a number of novel loci that have a causal impact on CAD including C5orf38, SF3A3, DHX36, and MRPL33. Furthermore, by integrating the risk factors of CAD in our analysis, we were able to investigate the clinical pathways whereby several of our probes exert their effect. We found that the SELE protein level in the blood is under the trans-regulatory impact of methylation sites within the ABO gene and that SELE exerts its effect on CAD through immune, glycemic, and lipid metabolism, making it a candidate of interest for therapeutic interventions. We found the methylation site, cg05126514 within the BSN gene exert its effect on CAD through central nervous system-lifestyle risk factors. Finally, genes with a transcriptional regulatory role (SF3A3, ILF3, and N4BP2L2) exert their effect on CAD through height.

CONCLUSIONS

We demonstrate that multiomics data analysis is a powerful approach to unravel the functional mechanisms underlying CAD loci and to identify novel molecular biomarkers. Our results indicate epigenetic modifications are important in the pathogenesis of CAD and identifying and targeting these sites is of potential therapeutic interest to address the detrimental effects of both environmental and genetic factors.

SGCG rs679482 Associates With Weight Loss Success in Response to an Intensively Supervised Outpatient Program

Weight loss in response to energy restriction is highly variable, and identification of genetic contributors can provide insights into underlying biology. Leveraging 1000 Genomes imputed genotypes, we carried out genome-wide association study (GWAS) analysis in 551 unrelated obese subjects of European ancestry who participated in an intensively supervised weight loss program with replication of promising signals in an independent sample of 1,331 obese subjects who completed the program at a later date. By single nucleotide polymorphism-based and sib-pair analysis, we show that that weight loss is a heritable trait, with estimated heritability (h ² = 0.49) within the range reported for obesity. We find rs679482, intronic to SGCG (sarcoglycan γ), highly expressed in skeletal muscle, to concordantly associate with weight loss in discovery and replication samples reaching GWAS significance in the combined meta-analysis (β = -0.35, P = 1.7 × 10^-12). Located in a region of open chromatin, rs679482 is predicted to bind DMRT2, and allele-specific transcription factor binding analysis indicates preferential binding of DMRT2 to rs679482-A. Concordantly, rs679482-A impairs native repressor activity and increases basal and DMRT2-mediated enhancer activity. These findings confirm that weight loss is a heritable trait and provide evidence by which a novel variant in SGCG, rs679482, leads to impaired diet response.

Association of MC4R (rs17782313) with diabetes and cardiovascular disease in Korean men and women

Background

Diabetes is mostly assessed by the fasting glucose level. Several studies reported that serum fasting glucose levels and cardiovascular disease are associated with MC4R.

Methods

A total of 4294 subjects participated in this study. There were 1810 subjects with cardiovascular disease among the 4294 subjects. We used multivariate linear regression models and multiple logistic regression analysis.

Results

Individuals with the TC/CC genotype had a 1.29-fold higher risk of diabetes than did those with the TT genotype when adjusting for age, sex, and BMI (OR, 1.29; 95% CI, 1.04–1.60). For healthy subjects, the association was significant in women (OR, 1.99; 95% CI, 1.01–3.93). Men with the TC/CC genotype had a 1.21-fold higher risk of cardiovascular disease than did those with the TT genotype when adjusting for age, sex, and BMI (OR, 1.21; 95% CI, 1.04–1.41). The relationship between MC4R and cardiovascular disease was stronger in lean men (OR, 1.40; 95% CI, 1.12–1.74, p = 0.0028) than in overweight men.

Conclusions

This study suggests that the rs17782313 SNP in MC4R is related to diabetes and the SNP is also associated with cardiovascular disease in lean men.

A genome-wide scan for pleiotropy between bone mineral density and nonbone phenotypes

Osteoporosis is the most common metabolic bone disorder globally and is characterized by skeletal fragility and microarchitectural deterioration. Genetic pleiotropy occurs when a single genetic element is associated with more than one phenotype. We aimed to identify pleiotropic loci associated with bone mineral density (BMD) and nonbone phenotypes in genome-wide association studies. In the discovery stage, the NHGRI-EBI Catalog was searched for genome-wide significant associations (P value < 5 × 10-8), excluding bone-related phenotypes. SNiPA was used to identify proxies of the significantly associated single nucleotide polymorphisms (SNPs) (r 2 = 1). We then assessed putative genetic associations of this set of SNPs with femoral neck (FN) and lumbar spine (LS) BMD data from the GEFOS Consortium. Pleiotropic variants were claimed at a false discovery rate < 1.4 × 10-3 for FN-BMD and < 1.5 × 10-3 for LS-BMD. Replication of these genetic markers was performed among more than 400 000 UK Biobank participants of European ancestry with available genetic and heel bone ultrasound data. In the discovery stage, 72 BMD-related pleiotropic SNPs were identified, and 12 SNPs located in 11 loci on 8 chromosomes were replicated in the UK Biobank. These SNPs were associated, in addition to BMD, with 14 different phenotypes. Most pleiotropic associations were exhibited by rs479844 (AP5B1, OVOL1 genes), which was associated with dermatological and allergic diseases, and rs4072037 (MUC1 gene), which was associated with magnesium levels and gastroenterological cancer. In conclusion, 12 BMD-related genome-wide significant SNPs showed pleiotropy with nonbone phenotypes. Pleiotropic associations can deepen the genetic understanding of bone-related diseases by identifying shared biological mechanisms with other diseases or traits.

Application of Genome-Wide Association Studies in Coronary Artery Disease

Coronary artery disease (CAD) is a complex disease caused by the combination of environmental and genetic factors. It is one of the leading causes of death and disability in the world. Much research has been focussed on CAD genetic mechanism. In recent years, genome-wide association study (GWAS) has developed rapidly around the world. Medical researchers around the world have successfully discovered a series of CAD genetic susceptibility genes or susceptible loci using medical research strategies, leading CAD research toward a new stage. This paper briefly summarizes the important progress made by GWAS for CAD in the world in recent years, and then analyzes the challenges faced by GWAS at this stage and the development trend of future research, to promote the transformation of genetic research results into clinical practice and provide guidance for further exploration of the genetic mechanism of CAD.

Phenome-wide screening for traits causally associated with the risk of coronary artery disease

Using two independent approaches, Mendelian randomization and Polygenic risk score in a sample of 6194 CAD cases and 4287 controls of European ancestry, we did a comprehensive phenome-wide search (PheWAS) for traits that causally associated with the risk of CAD. We found 46 risk factors that represented diverse categories including cardiovascular, CNS (central nervous system), diabetes, lipids, immune, anthropometry, and life style features; moreover, we noted numerous evidences of genetic correlations and causal associations between risk factors from different categories. Among the identified risk factors, CAD showed highest genetic relatedness with thrombotic conditions. The most represented category was life style features (29%) with evidence of strong genetic overlap with CNS traits. Genetic variants associated with higher cognition were associated with life style characteristics and cardiometabolic features that lower the risk of CAD. Conditional analysis indicated this trend is in part attributed to higher age of first sexual intercourse (AFS) in those with higher cognition. Lower AFS was concordantly associated with higher risk of CAD in males, females, and the combined sample; furthermore, lower AFS was causally associated with several CAD-risk factors including, higher fasting insulin, fasting glucose, LDL, immature reticulocyte fraction, HbA1c levels, as well as, higher risk of T2D and pulmonary embolism but lower levels of HDL. These results indicate CAD is the outcome of several phenotypically distinct but genetically interrelated sources; moreover, we identified lower AFS as an independent causal risk factor of CAD and revealed its role in mediating the effect of other risk factors.

Genome-wide identification of circulating-miRNA expression quantitative trait loci reveals the role of several miRNAs in the regulation of cardiometabolic phenotypes

Aims

To identify genetic variants that have a regulatory impact on circulating microRNAs (miRNAs) and to connect genetic risk to blood traits/biomarkers through the circulating miRNAs.

Methods and results

Leveraging miRNA-Seq data and the 1000 Genomes imputed genotypes, we carried out genome-wide association analysis for SNPs that regulate the expression of circulating miRNAs in a sample of 710 unrelated subjects of European ancestry. Wherever possible, we used data from the Framingham and the Geuvadis studies to replicate our findings. We found at least one genome-wide significant (P < 5e−8) miRNA-eQTL (mirQTL) for 143 circulating miRNAs. Overall each mirQTL explained a small portion (<1%) of variation in miRNA levels; however, we identified a few mirQTLs that explained 4% to 20% of variation in miRNA levels in plasma. Unlike trans-mirQTLs (P = 0.7), cis-mirQTLs tend to be also associated with their counterpart mature miRNAs (P < 0.0001), this suggests trans-mirQTLs exert their effect through processes that affect the stability of mature miRNAs; whereas, cis-mirQTLs mainly regulate the expression of primary-miRNAs. Next, we used the identified mirQTLs to investigate the links between circulating miRNAs with blood traits/biomarkers through Mendelian randomization analysis. We found miR-1908-5p plays an important role in regulating low-density lipoprotein (LDL), total cholesterol (TC), fasting glucose, HbA1c, and several lipid-metabolites in blood, whereas, miR-10b-5p mediates the trans-regulatory effect of the ABO locus on several blood proteins, coronary artery disease, and TC. Moreover, we demonstrated that a higher plasma level of miR-199a is causally associated with lower levels of LDL and TC. Finally, we found miR-143-3p and miR-145-5p are functionally related and mediate the effect of ZFPM2 on a number of its protein targets in blood including VEGFA, SERPINE1, and PDGFs.

Conclusions

This study identifies SNPs that have a regulatory impact on circulating miRNAs, and underlines the role of several circulating miRNAs in mediating the effect of a number of GWAS loci on cardiometabolic phenotypes.