Common, low-frequency, and rare genetic variants associated with lipoprotein subclasses and triglyceride measures in Finnish men from the METSIM study

PathGPS: discover shared genetic architecture using GWAS summary data

The increasing availability and scale of biobanks and "omic" datasets bring new horizons for understanding biological mechanisms. PathGPS is an exploratory data analysis tool to discover genetic architectures using Genome Wide Association Studies (GWAS) summary data. PathGPS is based on a linear structural equation model where traits are regulated by both genetic and environmental pathways. PathGPS decouples the genetic and environmental components by contrasting the GWAS associations of "signal" genes with those of "noise" genes. From the estimated genetic component, PathGPS then extracts genetic pathways via principal component and factor analysis, leveraging the low-rank and sparse properties. In addition, we provide a bootstrap aggregating ("bagging") algorithm to improve stability under data perturbation and hyperparameter tuning. When applied to a metabolomics dataset and the UK Biobank, PathGPS confirms several known gene-trait clusters and suggests multiple new hypotheses for future investigations.

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.

Familial LCAT Deficiency and Low HDL-C Levels: In silico Characterization of Two Rare LCAT Missense Mutations

Mutations in the lecithin-cholesterol acyltransferase (LCAT) gene, which catalyzes the esterification of cholesterol, result in two types of autosomal recessive disorders: Familial LCAT deficiency (FLD) and Fish Eye Disease (FED). While both phenotypes are characterized by corneal opacities and different forms of dyslipidemia, such as low levels of high-density lipoprotein-cholesterol (HDL-C), FLD exhibits more severe clinical manifestations like splenomegaly, anemia, and renal failure. We describe the first clinically and genetically confirmed case of FLD in Colombia which corresponds to a 46-year-old woman with corneal opacity, hypothyroidism, and dyslipidemia, who does not have any manifestations of renal failure, with two pathogenic heterozygous missense variants in the LCAT gene: LCAT (NM_000229.2):c.803G>A (p.Arg268His) and LCAT (NM_000229.2):c.368G>C (p.Arg123Pro). In silico analysis of the mutations predicted the physicochemical properties of the mutated protein, causing instability and potentially decreased LCAT function. These compound mutations highlight the clinical heterogeneity of the phenotypes associated with LCAT gene mutations.

ADAMTS7: a Novel Therapeutic Target in Atherosclerosis

PURPOSE OF REVIEW

Genome-wide association studies have repeatedly linked the metalloproteinase ADAMTS7 to coronary artery disease. Here we aim to highlight recent findings surrounding the human genetics of ADAMTS7, novel mouse models that investigate ADAMTS7 function, and potential substrates of ADAMTS7 cleavage.

RECENT FINDINGS

Recent genome-wide association studies in coronary artery disease have replicated the GWAS signal for ADAMTS7 and shown that the signal holds true even across different ethnic groups. However, the direction of effect in humans remains unclear. A recent novel mouse model revealed that the proatherogenicity of ADAMTS7 is derived from its catalytic functions, while at the translational level, vaccinating mice against ADAMTS7 reduced atherosclerosis. Finally, in vitro proteomics approaches have identified extracellular matrix proteins as candidate substrates that may be causal for the proatherogenicity of ADAMTS7. ADAMTS7 represents an enticing target for therapeutic intervention. The recent studies highlighted here have replicated prior findings, confirming the genetic link between ADAMTS7 and atherosclerosis, while providing further evidence in mice that ADAMTS7 is a targetable proatherogenic enzyme.

Genome-wide association study on coordination and agility in 461 Chinese Han males

There is growing evidence that genetic factors can influence human athletic performance. In many sports performances, excellent coordination and agility are the keys to mastery. However, few studies have been devoted to identifying genetic influences on athletic performance. Methods: We generated a derived measure of coordination and agility from the data of hexagonal jumps and T-runs and conducted genome-wide association and meta-analysis studies focused on coordination and agility. Results: The phenotypic correlation and genetic covariance analysis indicated that hexagonal jumps and T-runs were possibly influenced by the same set of genetic factors (R = 0.27, genetic covariance = 0.59). Meta-analysis identified rs117047321 genome-wide significant association (N = 143, P < 10E-5) with coordination and agility, and this association was replicated in the replication group (N = 318, P < 0.05). The CG genotype samples of this single nucleotide polymorphism (SNP) required a longer average movement time than the CC genotype samples, and the CG genotype only exists in Asia, which may belong to the East Asia-specific variation. This SNP is located on MYO5B, which is highly expressed in tissues such as the brain, heart, and muscle, suggesting that this locus might be a genetic factor related to human energy metabolism. Conclusion: Our study indicated that genetic factors can affect the athletic performance of coordination and agility. These findings may provide valuable insights for using genetic factors to evaluate sports characteristics.

Effects of protein-coding variants on blood metabolite measurements and clinical biomarkers in the UK Biobank

Genome-wide association studies (GWASs) have established the contribution of common and low-frequency variants to metabolic blood measurements in the UK Biobank (UKB). To complement existing GWAS findings, we assessed the contribution of rare protein-coding variants in relation to 355 metabolic blood measurements-including 325 predominantly lipid-related nuclear magnetic resonance (NMR)-derived blood metabolite measurements (Nightingale Health Plc) and 30 clinical blood biomarkers-using 412,393 exome sequences from four genetically diverse ancestries in the UKB. Gene-level collapsing analyses were conducted to evaluate a diverse range of rare-variant architectures for the metabolic blood measurements. Altogether, we identified significant associations (p < 1 × 10^-8) for 205 distinct genes that involved 1,968 significant relationships for the Nightingale blood metabolite measurements and 331 for the clinical blood biomarkers. These include associations for rare non-synonymous variants in PLIN1 and CREB3L3 with lipid metabolite measurements and SYT7 with creatinine, among others, which may not only provide insights into novel biology but also deepen our understanding of established disease mechanisms. Of the study-wide significant clinical biomarker associations, 40% were not previously detected on analyzing coding variants in a GWAS in the same cohort, reinforcing the importance of studying rare variation to fully understand the genetic architecture of metabolic blood measurements.

Novel genetic loci associated with osteoarthritis in multi-ancestry analyses in the Million Veteran Program and UK Biobank

Osteoarthritis is a common progressive joint disease. As no effective medical interventions are available, osteoarthritis often progresses to the end stage, in which only surgical options such as total joint replacement are available. A more thorough understanding of genetic influences of osteoarthritis is essential to develop targeted personalized approaches to treatment, ideally long before the end stage is reached. To date, there have been no large multiancestry genetic studies of osteoarthritis. Here, we leveraged the unique resources of 484,374 participants in the Million Veteran Program and UK Biobank to address this gap. Analyses included participants of European, African, Asian and Hispanic descent. We discovered osteoarthritis-associated genetic variation at 10 loci and replicated findings from previous osteoarthritis studies. We also present evidence that some osteoarthritis-associated regions are robust to population ancestry. Drug repurposing analyses revealed enrichment of targets of several medication classes and provide potential insight into the etiology of beneficial effects of antiepileptics on osteoarthritis pain.

Spotlight on very-low-density lipoprotein as a driver of cardiometabolic disorders: Implications for disease progression and mechanistic insights

Very-low-density lipoprotein (VLDL) is the only lipoprotein containing apolipoprotein B that is secreted from the liver, where VLDL is assembled from apolipoproteins, cholesterol, and triglycerides. The primary function of VLDL is to transport cholesterol and other lipids to organs and cells for utilization. Apart from its role in normal biologic processes, VLDL is also known to contribute to the development of atherosclerotic cardiovascular disease. Large VLDL particles, which are subclassified according to their size by nuclear magnetic resonance spectrometry, are significantly correlated not only with atherosclerosis, but also with insulin resistance and diabetes incidence. VLDL can also be subclassified according to surface electrical charge by using anion-exchange chromatography. The most electronegative VLDL subclass is highly cytotoxic to endothelial cells and may contribute to coronary heart disease. In addition, electronegative VLDL contributes to the development of atrial remodeling, especially in patients with metabolic syndrome, which is an established risk factor for atrial fibrillation. In this review, we focus on the VLDL subclasses that are associated with apolipoprotein alterations and are involved in cardiometabolic disease. The postprandial enhancement of VLDL’s pathogenicity is a critical medical issue, especially in patients with metabolic syndrome. Therefore, the significance of the postprandial modification of VLDL’s chemical and functional properties is extensively discussed.

GWAS of longitudinal trajectories at biobank scale

Biobanks linked to massive, longitudinal electronic health record (EHR) data make numerous new genetic research questions feasible. One among these is the study of biomarker trajectories. For example, high blood pressure measurements over visits strongly predict stroke onset, and consistently high fasting glucose and Hb1Ac levels define diabetes. Recent research reveals that not only the mean level of biomarker trajectories but also their fluctuations, or within-subject (WS) variability, are risk factors for many diseases. Glycemic variation, for instance, is recently considered an important clinical metric in diabetes management. It is crucial to identify the genetic factors that shift the mean or alter the WS variability of a biomarker trajectory. Compared to traditional cross-sectional studies, trajectory analysis utilizes more data points and captures a complete picture of the impact of time-varying factors, including medication history and lifestyle. Currently, there are no efficient tools for genome-wide association studies (GWASs) of biomarker trajectories at the biobank scale, even for just mean effects. We propose TrajGWAS, a linear mixed effect model-based method for testing genetic effects that shift the mean or alter the WS variability of a biomarker trajectory. It is scalable to biobank data with 100,000 to 1,000,000 individuals and many longitudinal measurements and robust to distributional assumptions. Simulation studies corroborate that TrajGWAS controls the type I error rate and is powerful. Analysis of eleven biomarkers measured longitudinally and extracted from UK Biobank primary care data for more than 150,000 participants with 1,800,000 observations reveals loci that significantly alter the mean or WS variability.

Profile‐likelihood Bayesian model averaging for two‐sample summary data Mendelian randomization in the presence of horizontal pleiotropy

Two‐sample summary data Mendelian randomization is a popular method for assessing causality in epidemiology, by using genetic variants as instrumental variables. If genes exert pleiotropic effects on the outcome not entirely through the exposure of interest, this can lead to heterogeneous and (potentially) biased estimates of causal effect. We investigate the use of Bayesian model averaging to preferentially search the space of models with the highest posterior likelihood. We develop a Metropolis‐Hasting algorithm to perform the search using the recently developed MR‐RAPS as the basis for defining a posterior distribution that efficiently accounts for pleiotropic and weak instrument bias. We demonstrate how our general modeling approach can be extended from a standard one‐component causal model to a two‐component model, which allows a large proportion of SNPs to violate the InSIDE assumption. We use Monte Carlo simulations to illustrate our methods and compare it to several related approaches. We finish by applying our approach to investigate the causal role of cholesterol on the development age‐related macular degeneration.

New therapeutic horizons for plasma phospholipid transfer protein (PLTP): Targeting endotoxemia, infection and sepsis

Phospholipid Transfer Protein (PLTP) transfers amphiphilic lipids between circulating lipoproteins and between lipoproteins, cells and tissues. Indeed, PLTP is a major determinant of the plasma levels, turnover and functionality of the main lipoprotein classes: very low-density lipoproteins (VLDL), low-density lipoproteins (LDL) and high-density lipoproteins (HDL). To date, most attention has been focused on the role of PLTP in the context of cardiometabolic diseases, with additional insights in neurodegenerative diseases and immunity. Importantly, beyond its influence on plasma triglyceride and cholesterol transport, PLTP plays a key role in the modulation of the immune response, with immediate relevance to a wide range of inflammatory diseases including bacterial infection and sepsis. Indeed, emerging evidence supports the role of PLTP, in the context of its association with lipoproteins, in the neutralization and clearance of bacterial lipopolysaccharides (LPS) or endotoxins. LPS are amphipathic molecules originating from Gram-negative bacteria which harbor major pathogen-associated patterns, triggering an innate immune response in the host. Although the early inflammatory reaction constitutes a key step in the anti-microbial defense of the organism, it can lead to a dysregulated inflammatory response and to hemodynamic disorders, organ failure and eventually death. Moreover, and in addition to endotoxemia and acute inflammation, small amounts of LPS in the circulation can induce chronic, low-grade inflammation with long-term consequences in several metabolic disorders such as atherosclerosis, obesity and diabetes. After an updated overview of the role of PLTP in lipid transfer, lipoprotein metabolism and related diseases, current knowledge of its impact on inflammation, infection and sepsis is critically appraised. Finally, the relevance of PLTP as a new player and novel therapeutic target in the fight against inflammatory diseases is considered.

Multi-Omic Approaches to Identify Genetic Factors in Metabolic Syndrome

Metabolic syndrome (MetS) is a highly heritable disease and a major public health burden worldwide. MetS diagnosis criteria are met by the simultaneous presence of any three of the following: high triglycerides, low HDL/high LDL cholesterol, insulin resistance, hypertension, and central obesity. These diseases act synergistically in people suffering from MetS and dramatically increase risk of morbidity and mortality due to stroke and cardiovascular disease, as well as certain cancers. Each of these component features is itself a complex disease, as is MetS. As a genetically complex disease, genetic risk factors for MetS are numerous, but not very powerful individually, often requiring specific environmental stressors for the disease to manifest. When taken together, all sequence variants that contribute to MetS disease risk explain only a fraction of the heritable variance, suggesting additional, novel loci have yet to be discovered. In this article, we will give a brief overview on the genetic concepts needed to interpret genome-wide association studies (GWAS) and quantitative trait locus (QTL) data, summarize the state of the field of MetS physiological genomics, and to introduce tools and resources that can be used by the physiologist to integrate genomics into their own research on MetS and any of its component features. There is a wealth of phenotypic and molecular data in animal models and humans that can be leveraged as outlined in this article. Integrating these multi-omic QTL data for complex diseases such as MetS provides a means to unravel the pathways and mechanisms leading to complex disease and promise for novel treatments. © 2022 American Physiological Society. Compr Physiol 12:1-40, 2022.

Pharmacogenomics of statins: lipid response and other outcomes in Brazilian cohorts

Statins are inhibitors of 3-hydroxy-3-methylglutaryl-CoA reductase, a key enzyme in cholesterol biosynthesis, that are highly effective in reducing plasma low-density lipoprotein (LDL) cholesterol and decreasing the risk of cardiovascular events. In recent years, a multitude of variants in genes involved in pharmacokinetics (PK) and pharmacodynamics (PD) have been suggested to influence the cholesterol-lowering response. However, the vast majority of studies have analyzed the pharmacogenetic associations in populations in Europe and the USA, whereas data in other populations, including Brazil, are mostly lacking. This narrative review provides an update of clinical studies on statin pharmacogenomics in Brazilian cohorts exploring lipid-lowering response, adverse events and pleiotropic effects. We find that variants in drug transporter genes (SLCO1B1 and ABCB1) positively impacted atorvastatin and simvastatin response, whereas variants in genes of drug metabolizing enzymes (CYP3A5) decreased response. Furthermore, multiple associations of variants in PD genes (HMGCR, LDLR and APOB) with statin response were identified. Few studies have explored statin-related adverse events, and only ABCB1 but not SLCO1B1 variants were robustly associated with increased risk in Brazil. Statin-related pleiotropic effects were shown to be influenced by variants in PD (LDLR, NR1H2) and antioxidant enzyme (NOS3, SOD2, MTHFR, SELENOP) genes. The findings of these studies indicate that statin pharmacogenomic associations are distinctly different in Brazil compared to other populations. This review also discusses the clinical implications of pharmacogenetic studies and the rising importance of investigating rare variants to explore their association with statin response.

The Genetic Basis of Hypertriglyceridemia

PURPOSE OF REVIEW

Hypertriglyceridemia is a common dyslipidemia associated with an increased risk of cardiovascular disease and pancreatitis. Severe hypertriglyceridemia may sometimes be a monogenic condition. However, in the vast majority of patients, hypertriglyceridemia is due to the cumulative effect of multiple genetic risk variants along with lifestyle factors, medications, and disease conditions that elevate triglyceride levels. In this review, we will summarize recent progress in the understanding of the genetic basis of hypertriglyceridemia.

RECENT FINDINGS

More than 300 genetic loci have been identified for association with triglyceride levels in large genome-wide association studies. Studies combining the loci into polygenic scores have demonstrated that some hypertriglyceridemia phenotypes previously attributed to monogenic inheritance have a polygenic basis. The new genetic discoveries have opened avenues for the development of more effective triglyceride-lowering treatments and raised interest towards genetic screening and tailored treatments against hypertriglyceridemia. The discovery of multiple genetic loci associated with elevated triglyceride levels has led to improved understanding of the genetic basis of hypertriglyceridemia and opened new translational opportunities.

Joo JW. MARS: leveraging allelic heterogeneity to increase power of association testing

In standard genome-wide association studies (GWAS), the standard association test is underpowered to detect associations between loci with multiple causal variants with small effect sizes. We propose a statistical method, Model-based Association test Reflecting causal Status (MARS), that finds associations between variants in risk loci and a phenotype, considering the causal status of variants, only requiring the existing summary statistics to detect associated risk loci. Utilizing extensive simulated data and real data, we show that MARS increases the power of detecting true associated risk loci compared to previous approaches that consider multiple variants, while controlling the type I error.

A Mendelian randomization study of the role of lipoprotein subfractions in coronary artery disease

Recent genetic data can offer important insights into the roles of lipoprotein subfractions and particle sizes in preventing coronary artery disease (CAD), as previous observational studies have often reported conflicting results. We used the LD score regression to estimate the genetic correlation of 77 subfraction traits with traditional lipid profile and identified 27 traits that may represent distinct genetic mechanisms. We then used Mendelian randomization (MR) to estimate the causal effect of these traits on the risk of CAD. In univariable MR, the concentration and content of medium high-density lipoprotein (HDL) particles showed a protective effect against CAD. The effect was not attenuated in multivariable analyses. Multivariable MR analyses also found that small HDL particles and smaller mean HDL particle diameter may have a protective effect. We identified four genetic markers for HDL particle size and CAD. Further investigations are needed to fully understand the role of HDL particle size.

Association of structural variation with cardiometabolic traits in Finns

The contribution of genome structural variation (SV) to quantitative traits associated with cardiometabolic diseases remains largely unknown. Here, we present the results of a study examining genetic association between SVs and cardiometabolic traits in the Finnish population. We used sensitive methods to identify and genotype 129,166 high-confidence SVs from deep whole-genome sequencing (WGS) data of 4,848 individuals. We tested the 64,572 common and low-frequency SVs for association with 116 quantitative traits and tested candidate associations using exome sequencing and array genotype data from an additional 15,205 individuals. We discovered 31 genome-wide significant associations at 15 loci, including 2 loci at which SVs have strong phenotypic effects: (1) a deletion of the ALB promoter that is greatly enriched in the Finnish population and causes decreased serum albumin level in carriers (p = 1.47 × 10-54) and is also associated with increased levels of total cholesterol (p = 1.22 × 10-28) and 14 additional cholesterol-related traits, and (2) a multi-allelic copy number variant (CNV) at PDPR that is strongly associated with pyruvate (p = 4.81 × 10-21) and alanine (p = 6.14 × 10-12) levels and resides within a structurally complex genomic region that has accumulated many rearrangements over evolutionary time. We also confirmed six previously reported associations, including five led by stronger signals in single nucleotide variants (SNVs) and one linking recurrent HP gene deletion and cholesterol levels (p = 6.24 × 10-10), which was also found to be strongly associated with increased glycoprotein level (p = 3.53 × 10-35). Our study confirms that integrating SVs in trait-mapping studies will expand our knowledge of genetic factors underlying disease risk.

Genome-Wide Identification of Rare and Common Variants Driving Triglyceride Levels in a Nevada Population

Clinical conditions correlated with elevated triglyceride levels are well-known: coronary heart disease, hypertension, and diabetes. Underlying genetic and phenotypic mechanisms are not fully understood, partially due to lack of coordinated genotypic-phenotypic data. Here we use a subset of the Healthy Nevada Project, a population of 9,183 sequenced participants with longitudinal electronic health records to examine consequences of altered triglyceride levels. Specifically, Healthy Nevada Project participants sequenced by the Helix Exome+ platform were cross-referenced to their electronic medical records to identify: (1) rare and common single-variant genome-wide associations; (2) gene-based associations using a Sequence Kernel Association Test; (3) phenome-wide associations with triglyceride levels; and (4) pleiotropic variants linked to triglyceride levels. The study identified 549 significant single-variant associations (p < 8.75 × 10-9), many in chromosome 11's triglyceride hotspot: ZPR1, BUD13, APOC3, APOA5. A well-known protective loss-of-function variant in APOC3 (R19X) was associated with a 51% decrease in triglyceride levels in the cohort. Sixteen gene-based triglyceride associations were identified; six of these genes surprisingly did not include a single variant with significant associations. Results at the variant and gene level were validated with the UK Biobank. The combination of a single-variant genome-wide association, a gene-based association method, and phenome wide-association studies identified rare and common variants, genes, and phenotypes associated with elevated triglyceride levels, some of which may have been overlooked with standard approaches.

Contrasting effects of Western vs Mediterranean diets on monocyte inflammatory gene expression and social behavior in a primate model

Dietary changes associated with industrialization increase the prevalence of chronic diseases, such as obesity, type II diabetes, and cardiovascular disease. This relationship is often attributed to an 'evolutionary mismatch' between human physiology and modern nutritional environments. Western diets enriched with foods that were scarce throughout human evolutionary history (e.g. simple sugars and saturated fats) promote inflammation and disease relative to diets more akin to ancestral human hunter-gatherer diets, such as a Mediterranean diet. Peripheral blood monocytes, precursors to macrophages and important mediators of innate immunity and inflammation, are sensitive to the environment and may represent a critical intermediate in the pathway linking diet to disease. We evaluated the effects of 15 months of whole diet manipulations mimicking Western or Mediterranean diet patterns on monocyte polarization in a well-established model of human health, the cynomolgus macaque (Macaca fascicularis). Monocyte transcriptional profiles differed markedly between diets, with 40% of transcripts showing differential expression (FDR < 0.05). Monocytes from Western diet consumers were polarized toward a more proinflammatory phenotype. The Western diet shifted the co-expression of 445 gene pairs, including small RNAs and transcription factors associated with metabolism and adiposity in humans, and dramatically altered behavior. For example, Western-fed individuals were more anxious and less socially integrated. These behavioral changes were also associated with some of the effects of diet on gene expression, suggesting an interaction between diet, central nervous system activity, and monocyte gene expression. This study provides new molecular insights into an evolutionary mismatch and uncovers new pathways through which Western diets alter monocyte polarization toward a proinflammatory phenotype.

ILRUN, a Human Plasma Lipid GWAS Locus, Regulates Lipoprotein Metabolism in Mice

RATIONALE

Single-nucleotide polymorphisms near the ILRUN (inflammation and lipid regulator with ubiquitin-associated-like and NBR1 [next to BRCA1 gene 1 protein]-like domains) gene are genome-wide significantly associated with plasma lipid traits and coronary artery disease (CAD), but the biological basis of this association is unknown.

OBJECTIVE

To investigate the role of ILRUN in plasma lipid and lipoprotein metabolism.

METHODS AND RESULTS

ILRUN encodes a protein that contains a ubiquitin-associated-like domain, suggesting that it may interact with ubiquitinylated proteins. We generated mice globally deficient for Ilrun and found they had significantly lower plasma cholesterol levels resulting from reduced liver lipoprotein production. Liver transcriptome analysis uncovered altered transcription of genes downstream of lipid-related transcription factors, particularly PPARα (peroxisome proliferator-activated receptor alpha), and livers from Ilrun-deficient mice had increased PPARα protein. Human ILRUN was shown to bind to ubiquitinylated proteins including PPARα, and the ubiquitin-associated-like domain of ILRUN was found to be required for its interaction with PPARα.

CONCLUSIONS

These findings establish ILRUN as a novel regulator of lipid metabolism that promotes hepatic lipoprotein production. Our results also provide functional evidence that ILRUN may be the casual gene underlying the observed genetic associations with plasma lipids at 6p21 in human.

High-density lipoprotein characteristics and coronary artery disease: a Mendelian randomization study

BACKGROUND

To assess whether genetically determined quantitative and qualitative HDL characteristics were independently associated with coronary artery disease (CAD).

METHODS

We designed a two-sample multivariate Mendelian randomization study with available genome-wide association summary data. We identified genetic variants associated with HDL cholesterol and apolipoprotein A-I levels, HDL size, particle levels, and lipid content to define our genetic instrumental variables in one sample (Kettunen et al. study, n = 24,925) and analyzed their association with CAD risk in a different study (CARDIoGRAMplusC4D, n = 184,305). We validated these results by defining our genetic variables in another database (METSIM, n = 8372) and studied their relationship with CAD in the CARDIoGRAMplusC4D dataset. To estimate the effect size of the associations of interest adjusted for other lipoprotein traits and minimize potential pleiotropy, we used the Multi-trait-based Conditional & Joint analysis.

RESULTS

Genetically determined HDL cholesterol and apolipoprotein A-I levels were not associated with CAD. HDL mean diameter (β = 0.27 [95%CI = 0.19; 0.35]), cholesterol levels in very large HDLs (β = 0.29 [95%CI = 0.17; 0.40]), and triglyceride content in very large HDLs (β = 0.14 [95%CI = 0.040; 0.25]) were directly associated with CAD risk, whereas the cholesterol content in medium-sized HDLs (β = -0.076 [95%CI = -0.10; -0.052]) was inversely related to this risk. These results were validated in the METSIM-CARDIoGRAMplusC4D data.

CONCLUSIONS

Some qualitative HDL characteristics (related to size, particle distribution, and cholesterol and triglyceride content) are related to CAD risk while HDL cholesterol levels are not.

Adiponectin GWAS loci harboring extensive allelic heterogeneity exhibit distinct molecular consequences

Loci identified in genome-wide association studies (GWAS) can include multiple distinct association signals. We sought to identify the molecular basis of multiple association signals for adiponectin, a hormone involved in glucose regulation secreted almost exclusively from adipose tissue, identified in the Metabolic Syndrome in Men (METSIM) study. With GWAS data for 9,262 men, four loci were significantly associated with adiponectin: ADIPOQ, CDH13, IRS1, and PBRM1. We performed stepwise conditional analyses to identify distinct association signals, a subset of which are also nearly independent (lead variant pairwise r²<0.01). Two loci exhibited allelic heterogeneity, ADIPOQ and CDH13. Of seven association signals at the ADIPOQ locus, two signals colocalized with adipose tissue expression quantitative trait loci (eQTLs) for three transcripts: trait-increasing alleles at one signal were associated with increased ADIPOQ and LINC02043, while trait-increasing alleles at the other signal were associated with decreased ADIPOQ-AS1. In reporter assays, adiponectin-increasing alleles at two signals showed corresponding directions of effect on transcriptional activity. Putative mechanisms for the seven ADIPOQ signals include a missense variant (ADIPOQ G90S), a splice variant, a promoter variant, and four enhancer variants. Of two association signals at the CDH13 locus, the first signal consisted of promoter variants, including the lead adipose tissue eQTL variant for CDH13, while a second signal included a distal intron 1 enhancer variant that showed ~2-fold allelic differences in transcriptional reporter activity. Fine-mapping and experimental validation demonstrated that multiple, distinct association signals at these loci can influence multiple transcripts through multiple molecular mechanisms.

Many DNA variants affect common human traits, and distinct variants can have different effects on the function or expression level of the same gene. We identified variants associated with levels of adiponectin, a hormone involved in glucose regulation. Among these variants, we specifically studied the sets of variants located near two genes, ADIPOQ and CDH13, to determine how the variants affect gene expression or function. We focused on sets of variants that can be inherited together but are not always inherited together. Of the variants associated with adiponectin and located near ADIPOQ, one set were also associated with higher expression levels of the protein-coding ADIPOQ gene and a nearby non-coding gene, a second set of variants were associated with lower levels of the ADIPOQ-AS1 antisense gene, and additional variants changed the amino acid sequence or size of the adiponectin protein. Our examples show the benefits of identifying multiple sets of trait-associated variants in the same DNA region. These variants explain more trait variation, help identify genes that affect the trait, and guide studies of gene regulation and biological processes.

Systems Genetics for Mechanistic Discovery in Heart Diseases

Cardiovascular diseases are the leading cause of death worldwide. Complex diseases with highly heterogenous disease progression among patient populations, cardiovascular diseases feature multifactorial contributions from both genetic and environmental stressors. Despite significant effort utilizing multiple approaches from molecular biology to genome-wide association studies, the genetic landscape of cardiovascular diseases, particularly for the nonfamilial forms of heart failure, is still poorly understood. In the past decade, systems-level approaches based on omics technologies have become an important approach for the study of complex traits in large populations. These advances create opportunities to integrate genetic variation with other biological layers to identify and prioritize candidate genes, understand pathogenic pathways, and elucidate gene-gene and gene-environment interactions. In this review, we will highlight some of the recent progress made using systems genetics approaches to uncover novel mechanisms and molecular bases of cardiovascular pathophysiological manifestations. The key technology and data analysis platforms necessary to implement systems genetics will be described, and the current major challenges and future directions will also be discussed. For complex cardiovascular diseases, such as heart failure, systems genetics represents a powerful strategy to obtain mechanistic insights and to develop individualized diagnostic and therapeutic regiments, paving the way for precision cardiovascular medicine.

Genome-wide association study of adipocyte lipolysis in the GENetics of adipocyte lipolysis (GENiAL) cohort

OBJECTIVES

Lipolysis, hydrolysis of triglycerides to fatty acids in adipocytes, is tightly regulated, poorly understood, and, if perturbed, can lead to metabolic diseases including obesity and type 2 diabetes. The goal of this study was to identify the genetic regulators of lipolysis and elucidate their molecular mechanisms.

METHODS

Adipocytes from abdominal subcutaneous adipose tissue biopsies were isolated and were incubated without (spontaneous lipolysis) or with a catecholamine (stimulated lipolysis) to analyze lipolysis. DNA was extracted and genome-wide genotyping and imputation conducted. After quality control, 939 samples with genetic and lipolysis data were available. Genome-wide association studies of spontaneous and stimulated lipolysis were conducted. Subsequent in vitro gene expression analyses were used to identify candidate genes and explore their regulation of adipose tissue biology.

RESULTS

One locus on chromosome 19 demonstrated genome-wide significance with spontaneous lipolysis. 60 loci showed suggestive associations with spontaneous or stimulated lipolysis, of which many influenced both traits. In the chromosome 19 locus, only HIF3A was expressed in the adipocytes and displayed genotype-dependent gene expression. HIF3A knockdown in vitro increased lipolysis and the expression of key lipolysis-regulating genes.

CONCLUSIONS

In conclusion, we identified a genetic regulator of spontaneous lipolysis and provided evidence of HIF3A as a novel key regulator of lipolysis in subcutaneous adipocytes as the mechanism through which the locus influences adipose tissue biology.

FAM13A affects body fat distribution and adipocyte function

Genetic variation in the FAM13A (Family with Sequence Similarity 13 Member A) locus has been associated with several glycemic and metabolic traits in genome-wide association studies (GWAS). Here, we demonstrate that in humans, FAM13A alleles are associated with increased FAM13A expression in subcutaneous adipose tissue (SAT) and an insulin resistance-related phenotype (e.g. higher waist-to-hip ratio and fasting insulin levels, but lower body fat). In human adipocyte models, knockdown of FAM13A in preadipocytes accelerates adipocyte differentiation. In mice, Fam13a knockout (KO) have a lower visceral to subcutaneous fat (VAT/SAT) ratio after high-fat diet challenge, in comparison to their wild-type counterparts. Subcutaneous adipocytes in KO mice show a size distribution shift toward an increased number of smaller adipocytes, along with an improved adipogenic potential. Our results indicate that GWAS-associated variants within the FAM13A locus alter adipose FAM13A expression, which in turn, regulates adipocyte differentiation and contribute to changes in body fat distribution.

The influence of rare variants in circulating metabolic biomarkers

Circulating metabolite levels are biomarkers for cardiovascular disease (CVD). Here we studied, association of rare variants and 226 serum lipoproteins, lipids and amino acids in 7,142 (discovery plus follow-up) healthy participants. We leveraged the information from multiple metabolite measurements on the same participants to improve discovery in rare variant association analyses for gene-based and gene-set tests by incorporating correlated metabolites as covariates in the validation stage. Gene-based analysis corrected for the effective number of tests performed, confirmed established associations at APOB, APOC3, PAH, HAL and PCSK (p<1.32x10-7) and identified novel gene-trait associations at a lower stringency threshold with ACSL1, MYCN, FBXO36 and B4GALNT3 (p<2.5x10-6). Regulation of the pyruvate dehydrogenase (PDH) complex was associated for the first time, in gene-set analyses also corrected for effective number of tests, with IDL and LDL parameters, as well as circulating cholesterol (pMETASKAT<2.41x10-6). In conclusion, using an approach that leverages metabolite measurements obtained in the same participants, we identified novel loci and pathways involved in the regulation of these important metabolic biomarkers. As large-scale biobanks continue to amass sequencing and phenotypic information, analytical approaches such as ours will be useful to fully exploit the copious amounts of biological data generated in these efforts.

Carotid Intima-Media Thickness: Novel Loci, Sex-Specific Effects, and Genetic Correlations With Obesity and Glucometabolic Traits in UK Biobank

OBJECTIVE

Atherosclerosis is the underlying cause of most cardiovascular disease, but mechanisms underlying atherosclerosis are incompletely understood. Ultrasound measurement of the carotid intima-media thickness (cIMT) can be used to measure vascular remodeling, which is indicative of atherosclerosis. Genome-wide association studies have identified many genetic loci associated with cIMT, but heterogeneity of measurements collected by many small cohorts have been a major limitation in these efforts. Here, we conducted genome-wide association analyses in UKB (UK Biobank; N=22 179), the largest single study with consistent cIMT measurements. Approach and Results: We used BOLT-LMM software to run linear regression of cIMT in UKB, adjusted for age, sex, and genotyping chip. In white British participants, we identified 5 novel loci associated with cIMT and replicated most previously reported loci. In the first sex-specific analyses of cIMT, we identified a locus on chromosome 5, associated with cIMT in women only and highlight VCAN as a good candidate gene at this locus. Genetic correlations with body mass index and glucometabolic traits were also observed. Two loci influenced risk of ischemic heart disease.

CONCLUSIONS

These findings replicate previously reported associations, highlight novel biology, and provide new directions for investigating the sex differences observed in cardiovascular disease presentation and progression.

A genome-wide search for gene-by-obesity interaction loci of dyslipidemia in Koreans shows diverse genetic risk alleles

Dyslipidemia is a well-established risk factor for CVD. Studies suggest that similar fat accumulation in a given population might result in different levels of dyslipidemia risk among individuals; for example, despite similar or leaner body composition compared with Caucasians, Asians of Korean descent experience a higher prevalence of dyslipidemia. These variations imply a possible role of gene-obesity interactions on lipid profiles. Genome-wide association studies have identified more than 500 loci regulating plasma lipids, but the interaction structure between genes and obesity traits remains unclear. We hypothesized that some loci modify the effects of obesity on dyslipidemia risk and analyzed extensive gene-environment interactions (G×Es) at genome-wide levels to search for replicated gene-obesity interactive SNPs. In four Korean cohorts (n = 18,025), we identified and replicated 20 gene-obesity interactions, including novel variants (SCN1A and SLC12A8) and known lipid-associated variants (APOA5, BUD13, ZNF259, and HMGCR). When we estimated the additional heritability of dyslipidemia by considering G×Es, the gain was substantial for triglycerides (TGs) but mild for LDL cholesterol (LDL-C) and total cholesterol (Total-C); the interaction explained up to 18.7% of TG, 2.4% of LDL-C, and 1.9% of Total-C heritability associated with waist-hip ratio. Our findings suggest that some individuals are prone to develop abnormal lipid profiles, particularly with regard to TGs, even with slight increases in obesity indices; ethnic diversities in the risk alleles might partly explain the differential dyslipidemia risk between populations. Research about these interacting variables may facilitate knowledge-based approaches to personalize health guidelines according to individual genetic profiles.

A comprehensive study of metabolite genetics reveals strong pleiotropy and heterogeneity across time and context

Genetic studies of metabolites have identified thousands of variants, many of which are associated with downstream metabolic and obesogenic disorders. However, these studies have relied on univariate analyses, reducing power and limiting context-specific understanding. Here we aim to provide an integrated perspective of the genetic basis of metabolites by leveraging the Finnish Metabolic Syndrome In Men (METSIM) cohort, a unique genetic resource which contains metabolic measurements, mostly lipids, across distinct time points as well as information on statin usage. We increase effective sample size by an average of two-fold by applying the Covariates for Multi-phenotype Studies (CMS) approach, identifying 588 significant SNP-metabolite associations, including 228 new associations. Our analysis pinpoints a small number of master metabolic regulator genes, balancing the relative proportion of dozens of metabolite levels. We further identify associations to changes in metabolic levels across time as well as genetic interactions with statin at both the master metabolic regulator and genome-wide level.

Exome sequencing of Finnish isolates enhances rare-variant association power

Exome-sequencing studies have generally been underpowered to identify deleterious alleles with a large effect on complex traits as such alleles are mostly rare. Because the population of northern and eastern Finland has expanded considerably and in isolation following a series of bottlenecks, individuals of these populations have numerous deleterious alleles at a relatively high frequency. Here, using exome sequencing of nearly 20,000 individuals from these regions, we investigate the role of rare coding variants in clinically relevant quantitative cardiometabolic traits. Exome-wide association studies for 64 quantitative traits identified 26 newly associated deleterious alleles. Of these 26 alleles, 19 are either unique to or more than 20 times more frequent in Finnish individuals than in other Europeans and show geographical clustering comparable to Mendelian disease mutations that are characteristic of the Finnish population. We estimate that sequencing studies of populations without this unique history would require hundreds of thousands to millions of participants to achieve comparable association power.

Enhancer deletion and allelic effects define a regulatory molecular mechanism at the VLDLR cholesterol GWAS locus

Total cholesterol (TC) and low-density lipoprotein cholesterol (LDL-C) are heritable risk factors for cardiovascular disease, yet the molecular mechanisms underlying the majority of blood lipid-associated genome-wide association studies signals remain elusive. One association signal is located in intron 3 of VLDLR; rs3780181-A is a risk allele associated (P ≤ 2 × 10-9) with increased TC and LDL-C. We investigated variants, genes and mechanisms underlying this association signal. We used a functional genetic approach to show that the intronic region spanning rs3780181 exhibited 1.6-7.6-fold enhancer activity in human HepG2 hepatocyte, THP-1 monocyte and Simpson-Golabi-Behmel Syndrome (SGBS) preadipocyte cells and that the rs3780181-A risk allele showed significantly less enhancer activity compared with the G allele, consistent with the direction of an expression quantitative trait locus in liver. In addition, rs3780181 alleles showed differential binding to multiple nuclear proteins, including stronger IRF2 binding to the rs3780181 G allele. We used a CRISPR-cas9 approach to delete 475 and 663 bp of the putative enhancer element in HEK293T kidney cells; compared to expression of mock-edited cell lines, the homozygous enhancer deletion cell lines showed 1.2-fold significantly (P < 0.04) decreased expression of VLDLR, as well as 1.5-fold decreased expression of SMARCA2, located 388 kb away. Together, these results identify an enhancer of VLDLR expression and suggest that altered binding of one or more factors bound to rs3780181 alleles decreases enhancer activity and reduces at least VLDLR expression, leading to increased TC and LDL-C.

Identification of seven novel loci associated with amino acid levels using single-variant and gene-based tests in 8545 Finnish men from the METSIM study

Comprehensive metabolite profiling captures many highly heritable traits, including amino acid levels, which are potentially sensitive biomarkers for disease pathogenesis. To better understand the contribution of genetic variation to amino acid levels, we performed single variant and gene-based tests of association between nine serum amino acids (alanine, glutamine, glycine, histidine, isoleucine, leucine, phenylalanine, tyrosine, and valine) and 16.6 million genotyped and imputed variants in 8545 non-diabetic Finnish men from the METabolic Syndrome In Men (METSIM) study with replication in Northern Finland Birth Cohort (NFBC1966). We identified five novel loci associated with amino acid levels (P = < 5×10-8): LOC157273/PPP1R3B with glycine (rs9987289, P = 2.3×10-26); ZFHX3 (chr16:73326579, minor allele frequency (MAF) = 0.42%, P = 3.6×10-9), LIPC (rs10468017, P = 1.5×10-8), and WWOX (rs9937914, P = 3.8×10-8) with alanine; and TRIB1 with tyrosine (rs28601761, P = 8×10-9). Gene-based tests identified two novel genes harboring missense variants of MAF <1% that show aggregate association with amino acid levels: PYCR1 with glycine (Pgene = 1.5×10-6) and BCAT2 with valine (Pgene = 7.4×10-7); neither gene was implicated by single variant association tests. These findings are among the first applications of gene-based tests to identify new loci for amino acid levels. In addition to the seven novel gene associations, we identified five independent signals at established amino acid loci, including two rare variant signals at GLDC (rs138640017, MAF=0.95%, Pconditional = 5.8×10-40) with glycine levels and HAL (rs141635447, MAF = 0.46%, Pconditional = 9.4×10-11) with histidine levels. Examination of all single variant association results in our data revealed a strong inverse relationship between effect size and MAF (Ptrend<0.001). These novel signals provide further insight into the molecular mechanisms of amino acid metabolism and potentially, their perturbations in disease.

Genome-Wide Association Studies of Metabolite Concentrations (mGWAS): Relevance for Nephrology

Metabolites are small molecules that are intermediates or products of metabolism, many of which are freely filtered by the kidneys. In addition, the kidneys have a central role in metabolite anabolism and catabolism, as well as in active metabolite reabsorption and/or secretion during tubular passage. This review article illustrates how the coupling of genomics and metabolomics in genome-wide association analyses of metabolites can be used to illuminate mechanisms underlying human metabolism, with a special focus on insights relevant to nephrology. First, genetic susceptibility loci for reduced kidney function and chronic kidney disease (CKD) were reviewed systematically for their associations with metabolite concentrations in metabolomics studies of blood and urine. Second, kidney function and CKD-associated metabolites reported from observational studies were interrogated for metabolite-associated genetic variants to generate and discuss complementary insights. Finally, insights originating from the simultaneous study of both blood and urine or by modeling intermetabolite relationships are summarized. We also discuss methodologic questions related to the study of metabolite concentrations in urine as well as among CKD patients. In summary, genome-wide association analyses of metabolites using metabolite concentrations quantified from blood and/or urine are a promising avenue of research to illuminate physiological and pathophysiological functions of the kidney.

Multi-ancestry study of blood lipid levels identifies four loci interacting with physical activity

Many genetic loci affect circulating lipid levels, but it remains unknown whether lifestyle factors, such as physical activity, modify these genetic effects. To identify lipid loci interacting with physical activity, we performed genome-wide analyses of circulating HDL cholesterol, LDL cholesterol, and triglyceride levels in up to 120,979 individuals of European, African, Asian, Hispanic, and Brazilian ancestry, with follow-up of suggestive associations in an additional 131,012 individuals. We find four loci, in/near CLASP1, LHX1, SNTA1, and CNTNAP2, that are associated with circulating lipid levels through interaction with physical activity; higher levels of physical activity enhance the HDL cholesterol-increasing effects of the CLASP1, LHX1, and SNTA1 loci and attenuate the LDL cholesterol-increasing effect of the CNTNAP2 locus. The CLASP1, LHX1, and SNTA1 regions harbor genes linked to muscle function and lipid metabolism. Our results elucidate the role of physical activity interactions in the genetic contribution to blood lipid levels.

Deciphering the Emerging Complexities of Molecular Mechanisms at GWAS Loci

Genome-wide association studies (GWASs) have identified thousands of loci associated with hundreds of complex diseases and traits, and progress is being made toward elucidating the causal variants and genes underlying these associations. Functional characterization of mechanisms at GWAS loci is a multi-faceted challenge. Challenges include linkage disequilibrium and allelic heterogeneity at each locus, the noncoding nature of most loci, and the time and cost needed for experimentally evaluating the potential mechanistic contributions of genes and variants. As GWAS sample sizes increase, more loci are identified, and the complexities of individual loci emerge. Loci can consist of multiple association signals, each of which can reflect the influence of multiple variants, inseparable by association analyses. Each signal within a locus can influence the same or different target genes. Experimental studies of genes and variants can differ on the basis of cell type, cellular environment, or other context-specific variables. In this review, we describe the complexity of mechanisms at GWAS loci-including multiple signals, multiple variants, and/or multiple genes-and the implications these complexities hold for experimental study design and interpretation of GWAS mechanisms.

Exome sequencing of an adolescent with nonalcoholic fatty liver disease identifies a clinically actionable case of Wilson disease

Diagnostic whole-exome sequencing has proven highly successful in a range of rare diseases, particularly early-onset genetic conditions. In more common conditions, however, exome sequencing for diagnostic purposes remains the exception. Here we describe a patient initially diagnosed with a common, complex liver disease, nonalcoholic fatty liver disease (NAFLD), who was determined to have Wilson disease (WD) upon research-related exome sequencing. The patient presented as a 14.5-yr-old adolescent with chronically elevated aminotransferases, normal ceruloplasmin, and histologic examination consistent with NAFLD with advanced fibrosis. He was enrolled in a large longitudinal study of patients with NAFLD and was found to have WD by exome sequencing performed 4 yr later. This new diagnosis, confirmed clinically by 24 h urine copper quantification, led to a change in the therapy from lifestyle counseling to directed treatment with d-penicillamine, a copper chelating agent. In this case, the likelihood of making the correct diagnosis and thereby choosing the appropriate treatment was increased by exome sequencing and careful interpretation. This example illustrates the utility of exome sequencing diagnostically in more common conditions not currently considered as targets for genome-wide evaluation and adds to a growing body of evidence that patients diagnosed with more common conditions often in fact have rarer genetically determined syndromes that have escaped clinical detection.

ADAMTS proteins in human disorders

ADAMTS proteins are a superfamily of 26 secreted molecules comprising two related, but distinct families. ADAMTS proteases are zinc metalloendopeptidases, most of whose substrates are extracellular matrix (ECM) components, whereas ADAMTS-like proteins lack a metalloprotease domain, reside in the ECM and have regulatory roles vis-à-vis ECM assembly and/or ADAMTS activity. Evolutionary conservation and expansion of ADAMTS proteins in mammals is suggestive of crucial embryologic or physiological roles in humans. Indeed, Mendelian disorders or birth defects resulting from naturally occurring ADAMTS2, ADAMTS3, ADAMTS10, ADAMTS13, ADAMTS17, ADAMTS20, ADAMTSL2 and ADAMTSL4 mutations as well as numerous phenotypes identified in genetically engineered mice have revealed ADAMTS participation in major biological pathways. Important roles have been identified in a few acquired conditions. ADAMTS5 is unequivocally implicated in pathogenesis of osteoarthritis via degradation of aggrecan, a major structural proteoglycan in cartilage. ADAMTS7 is strongly associated with coronary artery disease and promotes atherosclerosis. Autoantibodies to ADAMTS13 lead to a platelet coagulopathy, thrombotic thrombocytopenic purpura, which is similar to that resulting from ADAMTS13 mutations. ADAMTS proteins have numerous potential connections to other human disorders that were identified by genome-wide association studies. Here, we review inherited and acquired human disorders in which ADAMTS proteins participate, and discuss progress and prospects in therapeutics.