Zari Dastani

Causal Relationship Between Obesity and Vitamin D Status: Bi-Directional Mendelian Randomization Analysis of Multiple Cohorts

A mendelian randomization study based on data from multiple cohorts conducted by Karani Santhanakrishnan Vimaleswaran and colleagues re-examines the causal nature of the relationship between vitamin D levels and obesity.

Novel loci for adiponectin levels and their influence on type 2 diabetes and metabolic traits: a multi-ethnic meta-analysis of 45,891 individuals.

Circulating levels of adiponectin, a hormone produced predominantly by adipocytes, are highly heritable and are inversely associated with type 2 diabetes mellitus (T2D) and other metabolic traits. We conducted a meta-analysis of genome-wide association studies in 39,883 individuals of European ancestry to identify genes associated with metabolic disease. We identified 8 novel loci associated with adiponectin levels and confirmed 2 previously reported loci (P = 4.5×10−8–1.2×10−43). Using a novel method to combine data across ethnicities (N = 4,232 African Americans, N = 1,776 Asians, and N = 29,347 Europeans), we identified two additional novel loci. Expression analyses of 436 human adipocyte samples revealed that mRNA levels of 18 genes at candidate regions were associated with adiponectin concentrations after accounting for multiple testing (p<3×10−4). We next developed a multi-SNP genotypic risk score to test the association of adiponectin decreasing risk alleles on metabolic traits and diseases using consortia-level meta-analytic data. This risk score was associated with increased risk of T2D (p = 4.3×10−3, n = 22,044), increased triglycerides (p = 2.6×10−14, n = 93,440), increased waist-to-hip ratio (p = 1.8×10−5, n = 77,167), increased glucose two hours post oral glucose tolerance testing (p = 4.4×10−3, n = 15,234), increased fasting insulin (p = 0.015, n = 48,238), but with lower in HDL-cholesterol concentrations (p = 4.5×10−13, n = 96,748) and decreased BMI (p = 1.4×10−4, n = 121,335). These findings identify novel genetic determinants of adiponectin levels, which, taken together, influence risk of T2D and markers of insulin resistance.

The Empirical Power of Rare Variant Association Methods: Results from Sanger Sequencing in 1,998 Individuals

The role of rare genetic variation in the etiology of complex disease remains unclear. However, the development of next-generation sequencing technologies offers the experimental opportunity to address this question. Several novel statistical methodologies have been recently proposed to assess the contribution of rare variation to complex disease etiology. Nevertheless, no empirical estimates comparing their relative power are available. We therefore assessed the parameters that influence their statistical power in 1,998 individuals Sanger-sequenced at seven genes by modeling different distributions of effect, proportions of causal variants, and direction of the associations (deleterious, protective, or both) in simulated continuous trait and case/control phenotypes. Our results demonstrate that the power of recently proposed statistical methods depend strongly on the underlying hypotheses concerning the relationship of phenotypes with each of these three factors. No method demonstrates consistently acceptable power despite this large sample size, and the performance of each method depends upon the underlying assumption of the relationship between rare variants and complex traits. Sensitivity analyses are therefore recommended to compare the stability of the results arising from different methods, and promising results should be replicated using the same method in an independent sample. These findings provide guidance in the analysis and interpretation of the role of rare base-pair variation in the etiology of complex traits and diseases.

Genetic evidence for a normal-weight “metabolically obese” phenotype linking insulin resistance, hypertension, coronary artery disease and type 2 diabetes

The mechanisms that predispose to hypertension, coronary artery disease (CAD), and type 2 diabetes (T2D) in individuals of normal weight are poorly understood. In contrast, in monogenic primary lipodystrophy—a reduction in subcutaneous adipose tissue—it is clear that it is adipose dysfunction that causes severe insulin resistance (IR), hypertension, CAD, and T2D. We aimed to test the hypothesis that common alleles associated with IR also influence the wider clinical and biochemical profile of monogenic IR. We selected 19 common genetic variants associated with fasting insulin–based measures of IR. We used hierarchical clustering and results from genome-wide association studies of eight nondisease outcomes of monogenic IR to group these variants. We analyzed genetic risk scores against disease outcomes, including 12,171 T2D cases, 40,365 CAD cases, and 69,828 individuals with blood pressure measurements. Hierarchical clustering identified 11 variants associated with a metabolic profile consistent with a common, subtle form of lipodystrophy. A genetic risk score consisting of these 11 IR risk alleles was associated with higher triglycerides (β = 0.018; P = 4 × 10−29), lower HDL cholesterol (β = −0.020; P = 7 × 10−37), greater hepatic steatosis (β = 0.021; P = 3 × 10−4), higher alanine transaminase (β = 0.002; P = 3 × 10−5), lower sex-hormone-binding globulin (β = −0.010; P = 9 × 10−13), and lower adiponectin (β = −0.015; P = 2 × 10−26). The same risk alleles were associated with lower BMI (per-allele β = −0.008; P = 7 × 10−8) and increased visceral-to-subcutaneous adipose tissue ratio (β = −0.015; P = 6 × 10−7). Individuals carrying ≥17 fasting insulin–raising alleles (5.5% population) were slimmer (0.30 kg/m2) but at increased risk of T2D (odds ratio [OR] 1.46; per-allele P = 5 × 10−13), CAD (OR 1.12; per-allele P = 1 × 10−5), and increased blood pressure (systolic and diastolic blood pressure of 1.21 mmHg [per-allele P = 2 × 10−5] and 0.67 mmHg [per-allele P = 2 × 10−4], respectively) compared with individuals carrying ≤9 risk alleles (5.5% population). Our results provide genetic evidence for a link between the three diseases of the “metabolic syndrome” and point to reduced subcutaneous adiposity as a central mechanism.

Mendelian Randomization Studies Do Not Support a Causal Role for Reduced Circulating Adiponectin Levels in Insulin Resistance and Type 2 Diabetes

Adiponectin is strongly inversely associated with insulin resistance and type 2 diabetes, but its causal role remains controversial. We used a Mendelian randomization approach to test the hypothesis that adiponectin causally influences insulin resistance and type 2 diabetes. We used genetic variants at the ADIPOQ gene as instruments to calculate a regression slope between adiponectin levels and metabolic traits (up to 31,000 individuals) and a combination of instrumental variables and summary statistics–based genetic risk scores to test the associations with gold-standard measures of insulin sensitivity (2,969 individuals) and type 2 diabetes (15,960 case subjects and 64,731 control subjects). In conventional regression analyses, a 1-SD decrease in adiponectin levels was correlated with a 0.31-SD (95% CI 0.26–0.35) increase in fasting insulin, a 0.34-SD (0.30–0.38) decrease in insulin sensitivity, and a type 2 diabetes odds ratio (OR) of 1.75 (1.47–2.13). The instrumental variable analysis revealed no evidence of a causal association between genetically lower circulating adiponectin and higher fasting insulin (0.02 SD; 95% CI −0.07 to 0.11; N = 29,771), nominal evidence of a causal relationship with lower insulin sensitivity (−0.20 SD; 95% CI −0.38 to −0.02; N = 1,860), and no evidence of a relationship with type 2 diabetes (OR 0.94; 95% CI 0.75–1.19; N = 2,777 case subjects and 13,011 control subjects). Using the ADIPOQ summary statistics genetic risk scores, we found no evidence of an association between adiponectin-lowering alleles and insulin sensitivity (effect per weighted adiponectin-lowering allele: −0.03 SD; 95% CI −0.07 to 0.01; N = 2,969) or type 2 diabetes (OR per weighted adiponectin-lowering allele: 0.99; 95% CI 0.95–1.04; 15,960 case subjects vs. 64,731 control subjects). These results do not provide any consistent evidence that interventions aimed at increasing adiponectin levels will improve insulin sensitivity or risk of type 2 diabetes.

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.

Adjusted Sequence Kernel Association Test for Rare Variants Controlling for Cryptic and Family Relatedness

Recent progress in sequencing technologies makes it possible to identify rare and unique variants that may be associated with complex traits. However, the results of such efforts depend crucially on the use of efficient statistical methods and study designs. Although family‐based designs might enrich a data set for familial rare disease variants, most existing rare variant association approaches assume independence of all individuals. We introduce here a framework for association testing of rare variants in family‐based designs. This framework is an adaptation of the sequence kernel association test (SKAT) which allows us to control for family structure. Our adjusted SKAT (ASKAT) combines the SKAT approach and the factored spectrally transformed linear mixed models (FaST‐LMMs) algorithm to capture family effects based on a LMM incorporating the realized proportion of the genome that is identical by descent between pairs of individuals, and using restricted maximum likelihood methods for estimation. In simulation studies, we evaluated type I error and power of this proposed method and we showed that regardless of the level of the trait heritability, our approach has good control of type I error and good power. Since our approach uses FaST‐LMM to calculate variance components for the proposed mixed model, ASKAT is reasonably fast and can analyze hundreds of thousands of markers. Data from the UK twins consortium are presented to illustrate the ASKAT methodology.

WW-domain-containing oxidoreductase is associated with low plasma HDL-C levels.

Low serum HDL-cholesterol (HDL-C) is a major risk factor for coronary artery disease. We performed targeted genotyping of a 12.4 Mb linked region on 16q to test for association with low HDL-C by using a regional-tag SNP strategy. We identified one SNP, rs2548861, in the WW-domain-containing oxidoreductase (WWOX) gene with region-wide significance for low HDL-C in dyslipidemic families of Mexican and European descent and in low-HDL-C cases and controls of European descent (p = 6.9 x 10(-7)). We extended our investigation to the population level by using two independent unascertained population-based Finnish cohorts, the cross-sectional METSIM cohort of 4,463 males and the prospective Young Finns cohort of 2,265 subjects. The combined analysis provided p = 4 x 10(-4) to 2 x 10(-5). Importantly, in the prospective cohort, we observed a significant longitudinal association of rs2548861 with HDL-C levels obtained at four different time points over 21 years (p = 0.003), and the T risk allele explained 1.5% of the variance in HDL-C levels. The rs2548861 resides in a highly conserved region in intron 8 of WWOX. Results from our in vitro reporter assay and electrophoretic mobility-shift assay demonstrate that this region functions as a cis-regulatory element whose associated rs2548861 SNP has a specific allelic effect and that the region forms an allele-specific DNA-nuclear-factor complex. In conclusion, analyses of 9,798 subjects show significant association between HDL-C and a WWOX variant with an allele-specific cis-regulatory function.

Identification of a novel C5L2 variant (S323I) in a french canadian family with familial combined hyperlipemia

Objective—A functional acylation stimulating protein (ASP) receptor, C5L2, has been recently identified in ASP-responsive cells. Impaired ASP-mediated triglyceride synthesis has previously been described in a subset of hyperapolipoprotein B/familial combined hyperlipidemia subjects. Methods and Results—DNA sequencing of C5L2 coding region in 61 unrelated probands identified a heterozygous variant (G968→T) in 1 subject, resulting in Ser323→Ile substitution in the carboxyl terminal region. This variant was not detected in 2176 additional chromosomes by restriction fragment length polymorphism or fluorescence polarization genotyping. Eight family members of the proband were identified with one altered (+/−)C5L2 allele. Nine other family members had the wild-type (+/+)C5L2 sequence. The abnormal allele was associated with increased plasma triglyceride, plasma cholesterol, low-density lipoprotein (LDL) cholesterol, apolipoprotein B and ASP. Of 23 subjects tested in cell-based ASP bioactivity assays, those with C5L2(+/−) variant (n=2) had a 50% reduction in ASP-stimulated triglyceride synthesis, glucose transport and marked reduction in maximal binding (Bmax). By contrast, a C5L2(+/+) family member responded normally, as did hyperapolipoprotein B normal ASP subjects compared with C5L2(+/+) controls (n=6). Conclusion—The S323I variant may alter C5L2 function and might be one molecular basis contributing to familial combined hyperlipidemia.

In Healthy Adults, Biological Activity of Vitamin D, as Assessed by Serum PTH, Is Largely Independent of DBP Concentrations

Vitamin D insufficiency, as measured by 25‐hydroxyvitamin D (25[OH]D) levels, has been associated with important health outcomes. The majority of vitamin D in circulation is bound to vitamin D–binding protein (DBP) and albumin, and recent genetic studies have demonstrated that serum DBP is a major determinant of 25(OH)D concentrations in adults. The impact of circulating DBP levels on vitamin Ds biologic action, is unclear, but is of particular relevance to vitamin D epidemiology, because a lack of control for DBP levels could strongly influence the association of vitamin D with disease. Serum parathyroid hormone (PTH) levels can act as a biological readout of 25(OH)D activity. We therefore assessed the relationship between serum total and free 25(OH)D and PTH with and without adjusting for DBP, in 2073 subjects of European descent. Total 25(OH)D levels correlated positively (r = 0.19, p = 1.8 × 10−17) with DBP, whereas the free 25(OH)D correlated negatively (r = −0.14, p = 5.0 × 10−12). Total and free 25(OH)D levels correlated negatively with PTH (r = −0.29, p = 1.3 × 10−39; r = −0.26, p = 1.9 × 10−33, respectively). Including age, body mass index (BMI), sex, estimated glomerular filtration rate, calcium, and season of blood draw as covariates, total 25(OH)D levels were significantly associated with log‐transformed PTH (lnPTH) levels (linear term: β = −0.010, p < 0.0001, squared term: β = 0.00004, p < 0.0001) and this association was not changed by adjusting for DBP. These findings provide evidence that in a largely vitamin D–sufficient cohort, the biological effect of vitamin D on PTH levels is mainly independent of DBP concentration. Accordingly, this study may provide useful information for studies investigating the relationship between vitamin D, DBP, and disease.