👤 J Kuusisto

Alternate splicing events can create isoforms that alter gene function, and genetic variants associated with alternate gene isoforms may reveal molecular mechanisms of disease. We used subcutaneous adipose tissue of 426 Finnish men from the METSIM study and identified splice junction quantitative trait loci (sQTLs) for 6,077 splice junctions (FDR < 1%). In the same individuals, we detected expression QTLs (eQTLs) for 59,443 exons and 15,397 genes (FDR < 1%). We identified 595 genes with an sQTL and exon eQTL but no gene eQTL, which could indicate potential isoform differences. Of the significant sQTL signals, 2,114 (39.8%) included at least one proxy variant (linkage disequilibrium r Show less

Lean body mass (LM) plays an important role in mobility and metabolic function. We previously identified five loci associated with LM adjusted for fat mass in kilograms. Such an adjustment may reduce the power to identify genetic signals having an association with both lean mass and fat mass. To determine the impact of different fat mass adjustments on genetic architecture of LM and identify additional LM loci. We performed genome-wide association analyses for whole-body LM (20 cohorts of European ancestry with n = 38,292) measured using dual-energy X-ray absorptiometry) or bioelectrical impedance analysis, adjusted for sex, age, age2, and height with or without fat mass adjustments (Model 1 no fat adjustment; Model 2 adjustment for fat mass as a percentage of body mass; Model 3 adjustment for fat mass in kilograms). Seven single-nucleotide polymorphisms (SNPs) in separate loci, including one novel LM locus (TNRC6B), were successfully replicated in an additional 47,227 individuals from 29 cohorts. Based on the strengths of the associations in Model 1 vs Model 3, we divided the LM loci into those with an effect on both lean mass and fat mass in the same direction and refer to those as "sumo wrestler" loci (FTO and MC4R). In contrast, loci with an impact specifically on LM were termed "body builder" loci (VCAN and ADAMTSL3). Using existing available genome-wide association study databases, LM increasing alleles of SNPs in sumo wrestler loci were associated with an adverse metabolic profile, whereas LM increasing alleles of SNPs in "body builder" loci were associated with metabolic protection. In conclusion, we identified one novel LM locus (TNRC6B). Our results suggest that a genetically determined increase in lean mass might exert either harmful or protective effects on metabolic traits, depending on its relation to fat mass. Show less

This manuscript has not been published before and is not currently being considered for publication elsewhere. Increased septal convexity of left ventricle has been described in subjects with hypertrophic cardiomyopathy (HCM) -causing mutations without left ventricular hypertrophy (LVH). Our objective was to study septal convexity by cardiac magnetic resonance (CMR) in subjects with the Finnish founder mutation Q1016X in the myosin-binding protein C gene (MYBPC3). Septal convexity was measured in end-diastolic 4-chamber CMR image in 67 study subjects (47 subjects with the MYBPC3-Q1061X mutation and 20 healthy relatives without the mutation). Septal convexity was significantly increased in subjects with the MYBPC3-Q1061X mutation and LVH (n = 32) compared to controls (11.4 ± 4.3 vs 2.7 ± 3.2 mm, P < 0.001). In mutation carriers without LVH, there was a trend for increased septal convexity compared to controls (4.9 ± 2.5 vs 2.7 ± 3.2 mm, P = 0.074). When indexed for BSA, septal convexity in mutation carriers without LVH was 2.8 ± 1.4 mm/m Show less

Nationwide large-scale genetic and outcome studies in cohorts with hypertrophic cardiomyopathy (HCM) have not been previously published. We sequenced 59 cardiomyopathy-associated genes in 382 unrelated Finnish patients with HCM and found 24 pathogenic or likely pathogenic mutations in six genes in 38.2% of patients. Most mutations were located in sarcomere genes (MYBPC3, MYH7, TPM1, and MYL2). Previously reported mutations by our study group (MYBPC3-Gln1061Ter, MYH7-Arg1053Gln, and TPM1-Asp175Asn) and a fourth major mutation MYH7-Val606Met accounted for 28.0% of cases. Mutations in GLA and PRKAG2 were found in three patients. Furthermore, we found 49 variants of unknown significance in 31 genes in 20.4% of cases. During a 6.7 ± 4.2 year follow-up, annual all-cause mortality in 482 index patients and their relatives with HCM was higher than that in the matched Finnish population (1.70 vs. 0.87%; P < 0.001). Sudden cardiac deaths were rare (n = 8). Systolic heart failure (hazard ratio 17.256, 95% confidence interval 3.266-91.170, P = 0.001) and maximal left ventricular wall thickness (hazard ratio 1.223, 95% confidence interval 1.098-1.363, P < 0.001) were independent predictors of HCM-related mortality and life-threatening cardiac events. The patients with a pathogenic or likely pathogenic mutation underwent an implantable cardioverter defibrillator implantation more often than patients without a pathogenic or likely pathogenic mutation (12.9 vs. 3.5%, P < 0.001), but there was no difference in all-cause or HCM-related mortality between the two groups. Mortality due to HCM during 10 year follow-up among the 5.2 million population of Finland was studied from death certificates of the National Registry, showing 269 HCM-related deaths, of which 32% were sudden. We identified pathogenic and likely pathogenic mutations in 38% of Finnish patients with HCM. Four major sarcomere mutations accounted for 28% of HCM cases, whereas HCM-related mutations in non-sarcomeric genes were rare. Mortality in patients with HCM exceeded that of the general population. Finally, among 5.2 million Finns, there were at least 27 HCM-related deaths annually. Show less

Genome-wide association studies (GWAS) have identified >250 loci for body mass index (BMI), implicating pathways related to neuronal biology. Most GWAS loci represent clusters of common, noncoding variants from which pinpointing causal genes remains challenging. Here we combined data from 718,734 individuals to discover rare and low-frequency (minor allele frequency (MAF) < 5%) coding variants associated with BMI. We identified 14 coding variants in 13 genes, of which 8 variants were in genes (ZBTB7B, ACHE, RAPGEF3, RAB21, ZFHX3, ENTPD6, ZFR2 and ZNF169) newly implicated in human obesity, 2 variants were in genes (MC4R and KSR2) previously observed to be mutated in extreme obesity and 2 variants were in GIPR. The effect sizes of rare variants are ~10 times larger than those of common variants, with the largest effect observed in carriers of an MC4R mutation introducing a stop codon (p.Tyr35Ter, MAF = 0.01%), who weighed ~7 kg more than non-carriers. Pathway analyses based on the variants associated with BMI confirm enrichment of neuronal genes and provide new evidence for adipocyte and energy expenditure biology, widening the potential of genetically supported therapeutic targets in obesity. Show less

The sensitivity and specificity of the conventional 12-lead ECG to identify carriers of hypertrophic cardiomyopathy (HCM) - causing mutations without left ventricular hypertrophy (LVH) has been limited. We assessed the ability of novel electrocardiographic parameters to improve the detection of HCM mutation carriers. We studied 140 carriers (G+) of the TPM1-Asp175Asn or MYBPC3-Gln1061X pathogenic variants for HCM: The G+/LVH+ group (n = 98) consisted of mutation carriers with LVH and the G+/LVH- group (n = 42) without LVH. The control group consisted of 30 subjects. The standard 12-lead ECG was comprehensively analyzed and two novel ECG variables were introduced: RV1RV3 and septal remodeling. A subset of 65 individuals underwent cardiac magnetic resonance imaging and 2D strain echocardiography. Conventional major ECG criteria were sensitive (90%) and specific (97%) in identifying G+/LVH+ subjects. RV1RV3 and septal remodeling were more prevalent in the G+/LVH- subjects compared to the control group (33% vs 3%, p = 0.005 and 45% vs 3%, p < 0.001, respectively). The combination of RV1RV3 and Q waves and repolarization abnormalities (QR) differentiated between the G+/LVH- subjects and the control group with a sensitivity of 52% and specificity of 97%. The combination of septal remodeling and QR differentiated between G+/LVH- subjects and the control group with a sensitivity of 64% and specificity of 97%. The novel ECG-parameters RV1RV3 and septal remodeling were effective in identifying G+/LVH- subjects and could be useful in the diagnostics of new suspected HCM patients and in the screening and follow-up of HCM families. Show less

Hypertrophic cardiomyopathy (HCM) is characterized by ventricular repolarization abnormalities and risk of ventricular arrhythmias. Our aim was to study the association between the phenotype and ventricular repolarization dynamics in HCM patients. HCM patients with either the MYBPC3-Q1061X or TPM1-D175N mutation (n = 46) and control subjects without mutation and hypertrophy (n = 35) were studied with 24-hr ambulatory ECG recordings by measuring time intervals of rate-adapted QT (QTe), maximal QT, and T-wave apex to wave end (TPE) intervals and the QTe/RR slope. Findings were correlated to specified echocardiographic and cardiac magnetic resonance imaging (CMRI) findings. Rate-adapted QTe interval was progressively longer in HCM patients with decreasing heart rates compared to control subjects (p = 0.020). The degree of hypertrophy correlated with measured QTe values. HCM patients with maximal wall thickness higher than the mean (20.6 mm) had longer maximum QTe and median TPE intervals compared to control subjects and HCM patients with milder hypertrophy (p < 0.001 and p = 0.014, respectively). HCM patients with late gadolinium enhancement (LGE) on CMRI had steeper QTe/RR slopes compared to HCM patients without LGE and control subjects (p = 0.044 and p = 0.001, respectively). LGE was an independent predictor of QTe/RR slope (p = 0.023, B = 0.043). Dynamics of ventricular repolarization in HCM are affected by hypertrophy and fibrosis. LGE may confer an independent effect on QT dynamics which may increase the arrhythmogenic potential in HCM. Show less

We screened variants on an exome-focused genotyping array in >300,000 participants (replication in >280,000 participants) and identified 444 independent variants in 250 loci significantly associated with total cholesterol (TC), high-density-lipoprotein cholesterol (HDL-C), low-density-lipoprotein cholesterol (LDL-C), and/or triglycerides (TG). At two loci (JAK2 and A1CF), experimental analysis in mice showed lipid changes consistent with the human data. We also found that: (i) beta-thalassemia trait carriers displayed lower TC and were protected from coronary artery disease (CAD); (ii) excluding the CETP locus, there was not a predictable relationship between plasma HDL-C and risk for age-related macular degeneration; (iii) only some mechanisms of lowering LDL-C appeared to increase risk for type 2 diabetes (T2D); and (iv) TG-lowering alleles involved in hepatic production of TG-rich lipoproteins (TM6SF2 and PNPLA3) tracked with higher liver fat, higher risk for T2D, and lower risk for CAD, whereas TG-lowering alleles involved in peripheral lipolysis (LPL and ANGPTL4) had no effect on liver fat but decreased risks for both T2D and CAD. Show less

Glucose-dependent insulinotropic polypeptide (GIP) is an incretin hormone with extrapancreatic effects beyond glycemic control. Here we demonstrate unexpected effects of GIP signaling in the vasculature. GIP induces the expression of the proatherogenic cytokine osteopontin (OPN) in mouse arteries via local release of endothelin-1 and activation of CREB. Infusion of GIP increases plasma OPN concentrations in healthy individuals. Plasma endothelin-1 and OPN concentrations are positively correlated in patients with critical limb ischemia. Fasting GIP concentrations are higher in individuals with a history of cardiovascular disease (myocardial infarction, stroke) when compared with control subjects. GIP receptor (GIPR) and OPN mRNA levels are higher in carotid endarterectomies from patients with symptoms (stroke, transient ischemic attacks, amaurosis fugax) than in asymptomatic patients, and expression associates with parameters that are characteristic of unstable and inflammatory plaques (increased lipid accumulation, macrophage infiltration, and reduced smooth muscle cell content). While GIPR expression is predominantly endothelial in healthy arteries from humans, mice, rats, and pigs, remarkable upregulation is observed in endothelial and smooth muscle cells upon culture conditions, yielding a "vascular disease-like" phenotype. Moreover, the common variant rs10423928 in the GIPR gene is associated with increased risk of stroke in patients with type 2 diabetes. Show less

Non-alcoholic steatohepatitis (NASH) is associated with changes in fatty acid (FA) metabolism. However, specific changes in metabolism and hepatic mRNA expression related to NASH independent of simple steatosis, obesity and diet are unknown. Liver histology, serum and liver FA composition and estimated enzyme activities based on the FA ratios in cholesteryl esters and triglycerides were assessed in 92 obese participants of the Kuopio Obesity Surgery Study (KOBS) divided to those with normal liver, steatosis or NASH (30 men and 62 women, age 46.8±9.5years (mean±SD), BMI 44.2±6.2kg/m(2)). Plasma FA composition was also investigated in the Metabolic Syndrome in Men (METSIM) Study (n=769), in which serum alanine aminotransferase (ALT) was used as a marker of liver disease. Obese individuals with NASH had higher activity of estimated activities of delta-6 desaturase (D6D, p<0.002) and stearoyl-CoA desaturase 1 (SCD1, p<0.002) and lower activity of delta-5 desaturase (D5D, p<0.002) when compared to individuals with normal liver. Estimated activities of D5D, D6D and SCD1 correlated positively between liver and serum indicating that serum estimates reflected liver metabolism. Accordingly, NASH was associated with higher hepatic mRNA expression of corresponding genes FADS1, FADS2 and SCD. Finally, differences in FA metabolism that associated with NASH in obese individuals were also associated with high ALT in the METSIM Study. We demonstrated alterations in FA metabolism and endogenous desaturase activities that associate with NASH, independent of obesity and diet. This suggests that changes in endogenous FA metabolism are related to NASH and that they may contribute to the progression of the disease. Show less

Limited information exists on how the relationship between dietary intake of fat and fatty acids in erythrocytes and plasma is modulated by polymorphisms in the FADS gene cluster. We examined gene-diet interaction of total marine PUFA intake with a known gene encoding Δ-5 desaturase enzyme (FADS1) variant (rs174550) for fatty acids in erythrocyte membranes and plasma phospholipids (PL), cholesteryl esters (CE), and triglycerides (TG). In this cross-sectional study, fatty acid compositions were measured using GC, and total intake of polyunsaturated fat from fish and fish oil was estimated using a food frequency questionnaire in a subsample (n = 962) of the Metabolic Syndrome in Men Study. We found nominally significant gene-diet interactions for eicosapentaenoic acid (EPA, 20:5n-3) in erythrocytes (pinteraction = 0.032) and for EPA in plasma PL (pinteraction = 0.062), CE (pinteraction = 0.035), and TG (pinteraction = 0.035), as well as for docosapentaenoic acid (22:5n-3) in PL (pinteraction = 0.007). After excluding omega-3 supplement users, we found a significant gene-diet interaction for EPA in erythrocytes (pinteraction < 0.003). In a separate cohort of the Kuopio Obesity Surgery Study, the same locus was strongly associated with hepatic mRNA expression of FADS1 (p = 1.5 × 10(-10) ). FADS1 variants may modulate the relationship between marine fatty acid intake and circulating levels of long-chain omega-3 fatty acids. Show less

Previous data suggest that mitral valve leaflets are elongated in hypertrophic cardiomyopathy (HCM), and mitral valve leaflet elongation may constitute a primary phenotypic expression of HCM. Our objective was to measure the length of mitral valve leaflets by cardiovascular magnetic resonance (CMR) in subjects with HCM caused by a Finnish founder mutation in the myosin-binding protein C gene (MYBPC3-Q1061X), carriers of the same mutation without left ventricular hypertrophy, as well as in unselected consecutive patients with HCM, and respective controls. Anterior mitral valve leaflet (AML) and posterior mitral valve leaflet (PML) lengths were measured by CMR in 47 subjects with the Q1061X mutation in the gene encoding MYBPC3 and in 20 healthy relatives without the mutation. In addition, mitral valve leaflet lengths were measured by CMR in 80 consecutive non-genotyped patients with HCM in CMR and 71 age- and gender-matched healthy subjects. Of the subjects with the MYBPC-Q1016X mutation, 32 had left ventricular hypertrophy (LVH, LV maximal wall thickness ≥ 13 mm in CMR) and 15 had no hypertrophy. PML was longer in patients with the MYBPC3-Q1061X mutation and LVH than in controls of the MYBPC group (12.8 ± 2.8 vs 10.6 ± 1.9 mm, P = 0.013), but the difference between the groups was not statistically significant when PML was indexed for BSA (P = 0.066), or when PML length was adjusted for BSA, age, gender, LV mass and ejection fraction (P = 0.195). There was no significant difference in the PML length in mutation carriers without LVH and controls (11.1 ± 3.4 vs 10.6 ± 1.9, P = 0.52). We found no difference in AML lengths between the MYBPC mutation carriers with or without hypertrophy and controls. In 80 consecutive non-genotyped patients with HCM, there was no difference either in AML or PML lengths in subjects with HCM compared to respective control subjects. In subjects with HCM caused by the Q1061X mutation in the MYBPC3 gene, the posterior mitral valve leaflets may be elongated, but mitral valve elongation does not constitute primary phenotypic expression of the disease. Instead, elongated mitral valve leaflets seem to be associated with body size and left ventricular remodeling. Show less

We assessed the value of speckle tracking two-dimensional (2D) strain echocardiography (2DSE) measured mechanical dispersion (MD) with other imaging and electrocardiographic parameters in differentiating hypertrophic cardiomyopathy (HCM) patients with and without nonsustained ventricular tachycardia (NSVT) on 24-h ambulatory ECG monitoring. We studied 31 patients with HCM caused by the Finnish founder mutation MYBPC3-Q1061X and 20 control subjects with comprehensive 2DSE echocardiography and cardiac magnetic resonance imaging (CMRI). The presence of NSVT was assessed from ambulatory 24-h ECG monitoring. NSVT episodes were recorded in 11 (35%) patients with HCM. MD was significantly higher in HCM patients with NSVT (93 ± 41 ms) compared to HCM patients without NSVT (50 ± 18 ms, p = 0.012) and control subjects (41 ± 16 ms, p < 0.001). MD was the only variable independently associated with the presence of NSVT (OR: 1.60, 95% CI: 1.05-2.45, p = 0.030). Assessed by ROC curves, MD performed best in differentiating between HCM patients with and without NSVT (AUC = 0.81). Increased mechanical dispersion was associated with NSVT in HCM patients on 24-h ambulatory ECG monitoring. Key messages The prediction of sudden cardiac death in hypertrophic cardiomyopathy remains a challenge and novel imaging methods are required to identify individuals at risk of malignant ventricular arrhythmias. Mechanical dispersion by speckle tracking echocardiography is associated with NSVT on 24-h ambulatory ECG monitoring in patients with hypertrophic cardiomyopathy. Show less

Mutations in the cardiac myosin-binding protein C gene (MYBPC3) are the most common genetic cause of hypertrophic cardiomyopathy (HCM) worldwide. The molecular mechanisms leading to HCM are poorly understood. We investigated the metabolic profiles of mutation carriers with the HCM-causing MYBPC3-Q1061X mutation with and without left ventricular hypertrophy (LVH) and non-affected relatives, and the association of the metabolome to the echocardiographic parameters. 34 hypertrophic subjects carrying the MYBPC3-Q1061X mutation, 19 non-hypertrophic mutation carriers and 20 relatives with neither mutation nor hypertrophy were examined using comprehensive echocardiography. Plasma was analyzed for molecular lipids and polar metabolites using two metabolomics platforms. Concentrations of branched chain amino acids, triglycerides and ether phospholipids were increased in mutation carriers with hypertrophy as compared to controls and non-hypertrophic mutation carriers, and correlated with echocardiographic LVH and signs of diastolic and systolic dysfunction in subjects with the MYBPC3-Q1061X mutation. Our study implicates the potential role of branched chain amino acids, triglycerides and ether phospholipids in HCM, as well as suggests an association of these metabolites with remodeling and dysfunction of the left ventricle. Show less

In the nationwide FinHCM Study including 306 Finnish patients with hypertrophic cardiomyopathy (HCM), we have previously identified two founder mutations in the alpha-tropomyosin (TPM1-D175N) and myosin-binding protein C (MYBPC3-Q1061X) genes, accounting for 18% of all cases. Objective. To screen additional mutations, previously identified in eastern Finnish cohorts with HCM, in the FinHCM Study population. Ten mutations in the beta-myosin heavy chain gene (MYH7), TPM1, and MYBPC3 were screened. MYH7-R1053Q was found in 17 of 306 patients (5.6%). No carriers of MYH7-R719W or N696S were found. A novel TPM1-D175G mutation was found in a single patient. MYBPC3 mutations were found in 14 patients: IVS5-2A-C in two, IVS14-13G-A in two, K811del in six, and A851insT in four patients. Altogether, a HCM-causing mutation was identified in 32 patients, accounting for 10.5% of all cases. In addition, two MYBPC3 variants R326Q and V896M with uncertain pathogenicity were found in eight and in 10 patients, respectively. Combining the present findings with our previous results, a causative mutation was identified in 28% of the FinHCM cohort. MYH7-R1053Q was the third most common mutation, and should be screened in all new cases of HCM in Finland. Show less

Low-grade inflammation in obesity is associated with accumulation of the macrophage-derived cytokine osteopontin (OPN) in adipose tissue and induction of local as well as systemic insulin resistance. Since glucose-dependent insulinotropic polypeptide (GIP) is a strong stimulator of adipogenesis and may play a role in the development of obesity, we explored whether GIP directly would stimulate OPN expression in adipose tissue and thereby induce insulin resistance. GIP stimulated OPN protein expression in a dose-dependent fashion in rat primary adipocytes. The level of OPN mRNA was higher in adipose tissue of obese individuals (0.13 ± 0.04 vs. 0.04 ± 0.01, P < 0.05) and correlated inversely with measures of insulin sensitivity (r = -0.24, P = 0.001). A common variant of the GIP receptor (GIPR) (rs10423928) gene was associated with a lower amount of the exon 9-containing isoform required for transmembrane activity. Carriers of the A allele with a reduced receptor function showed lower adipose tissue OPN mRNA levels and better insulin sensitivity. Together, these data suggest a role for GIP not only as an incretin hormone but also as a trigger of inflammation and insulin resistance in adipose tissue. Carriers of the GIPR rs10423928 A allele showed protective properties via reduced GIP effects. Identification of this unprecedented link between GIP and OPN in adipose tissue might open new avenues for therapeutic interventions. Show less

Genome-wide association studies (GWAS) have identified ~100 loci associated with blood lipid levels, but much of the trait heritability remains unexplained, and at most loci the identities of the trait-influencing variants remain unknown. We conducted a trans-ethnic fine-mapping study at 18, 22, and 18 GWAS loci on the Metabochip for their association with triglycerides (TG), high-density lipoprotein cholesterol (HDL-C), and low-density lipoprotein cholesterol (LDL-C), respectively, in individuals of African American (n = 6,832), East Asian (n = 9,449), and European (n = 10,829) ancestry. We aimed to identify the variants with strongest association at each locus, identify additional and population-specific signals, refine association signals, and assess the relative significance of previously described functional variants. Among the 58 loci, 33 exhibited evidence of association at P<1 × 10(-4) in at least one ancestry group. Sequential conditional analyses revealed that ten, nine, and four loci in African Americans, Europeans, and East Asians, respectively, exhibited two or more signals. At these loci, accounting for all signals led to a 1.3- to 1.8-fold increase in the explained phenotypic variance compared to the strongest signals. Distinct signals across ancestry groups were identified at PCSK9 and APOA5. Trans-ethnic analyses narrowed the signals to smaller sets of variants at GCKR, PPP1R3B, ABO, LCAT, and ABCA1. Of 27 variants reported previously to have functional effects, 74% exhibited the strongest association at the respective signal. In conclusion, trans-ethnic high-density genotyping and analysis confirm the presence of allelic heterogeneity, allow the identification of population-specific variants, and limit the number of candidate SNPs for functional studies. Show less

Human complex metabolic traits are in part regulated by genetic determinants. Here we applied exome sequencing to identify novel associations of coding polymorphisms at minor allele frequencies (MAFs) >1% with common metabolic phenotypes. The study comprised three stages. We performed medium-depth (8×) whole exome sequencing in 1,000 cases with type 2 diabetes, BMI >27.5 kg/m(2) and hypertension and in 1,000 controls (stage 1). We selected 16,192 polymorphisms nominally associated (p < 0.05) with case-control status, from four selected annotation categories or from loci reported to associate with metabolic traits. These variants were genotyped in 15,989 Danes to search for association with 12 metabolic phenotypes (stage 2). In stage 3, polymorphisms showing potential associations were genotyped in a further 63,896 Europeans. Exome sequencing identified 70,182 polymorphisms with MAF >1%. In stage 2 we identified 51 potential associations with one or more of eight metabolic phenotypes covered by 45 unique polymorphisms. In meta-analyses of stage 2 and stage 3 results, we demonstrated robust associations for coding polymorphisms in CD300LG (fasting HDL-cholesterol: MAF 3.5%, p = 8.5 × 10(-14)), COBLL1 (type 2 diabetes: MAF 12.5%, OR 0.88, p = 1.2 × 10(-11)) and MACF1 (type 2 diabetes: MAF 23.4%, OR 1.10, p = 8.2 × 10(-10)). We applied exome sequencing as a basis for finding genetic determinants of metabolic traits and show the existence of low-frequency and common coding polymorphisms with impact on common metabolic traits. Based on our study, coding polymorphisms with MAF above 1% do not seem to have particularly high effect sizes on the measured metabolic traits. Show less

Hypertrophic cardiomyopathy (HCM) is predominantly caused by a large number of various mutations in the genes encoding sarcomeric proteins. However, two prevalent founder mutations for HCM in the alpha-tropomyosin (TPM1-D175N) and myosin-binding protein C (MYBPC3-Q1061X) genes have previously been identified in eastern Finland. To assess the prevalence of these founder mutations in a large population of patients with HCM from all over Finland. Patients and methods. We screened for two founder mutations (TPM1-D175N and MYBPC3-Q1061X) in 306 unrelated Finnish patients with HCM from the regions covering a population of ∼4,000,000. The TPM1-D175N mutation was found in 20 patients (6.5%) and the MYBPC3-Q1061X in 35 patients (11.4%). Altogether, the two mutations accounted for 17.9% of the HCM cases. In addition, 61 and 59 relatives of the probands were found to be carriers of TPM1-D175N and MYBPC3-Q1061X, respectively. The mutations showed regional clustering. TPM1-D175N was prevalent in central and western Finland, and MYBPC3-Q1061X in central and eastern Finland. The TPM1-D175N and MYBPC3-Q1061X mutations account for a substantial part of all HCM cases in the Finnish population, indicating that routine genetic screening of these mutations is warranted in Finnish patients with HCM. Show less

The incretin hormone GIP (glucose-dependent insulinotropic polypeptide) promotes pancreatic β-cell function by potentiating insulin secretion and β-cell proliferation. Recently, a combined analysis of several genome-wide association studies (Meta-analysis of Glucose and Insulin-Related Traits Consortium [MAGIC]) showed association to postprandial insulin at the GIP receptor (GIPR) locus. Here we explored mechanisms that could explain the protective effects of GIP on islet function. Associations of GIPR rs10423928 with metabolic and anthropometric phenotypes in both nondiabetic (N = 53,730) and type 2 diabetic individuals (N = 2,731) were explored by combining data from 11 studies. Insulin secretion was measured both in vivo in nondiabetic subjects and in vitro in islets from cadaver donors. Insulin secretion was also measured in response to exogenous GIP. The in vitro measurements included protein and gene expression as well as measurements of β-cell viability and proliferation. The A allele of GIPR rs10423928 was associated with impaired glucose- and GIP-stimulated insulin secretion and a decrease in BMI, lean body mass, and waist circumference. The decrease in BMI almost completely neutralized the effect of impaired insulin secretion on risk of type 2 diabetes. Expression of GIPR mRNA was decreased in human islets from carriers of the A allele or patients with type 2 diabetes. GIP stimulated osteopontin (OPN) mRNA and protein expression. OPN expression was lower in carriers of the A allele. Both GIP and OPN prevented cytokine-induced reduction in cell viability (apoptosis). In addition, OPN stimulated cell proliferation in insulin-secreting cells. These findings support β-cell proliferative and antiapoptotic roles for GIP in addition to its action as an incretin hormone. Identification of a link between GIP and OPN may shed new light on the role of GIP in preservation of functional β-cell mass in humans. Show less

We investigated the effects of 34 genetic risk variants for hyperglycemia/type 2 diabetes on lipoprotein subclasses and particle composition in a large population-based cohort. The study included 6,580 nondiabetic Finnish men from the population-based Metabolic Syndrome in Men (METSIM) study (aged 57 ± 7 years; BMI 26.8 ± 3.7 kg/m(2)). Genotyping of 34 single nucleotide polymorphism (SNPs) for hyperglycemia/type 2 diabetes was performed. Proton nuclear magnetic resonance spectroscopy was used to measure particle concentrations of 14 lipoprotein subclasses and their composition in native serum samples. The glucose-increasing allele of rs780094 in GCKR was significantly associated with low concentrations of VLDL particles (independently of their size) and small LDL and was nominally associated with low concentrations of intermediate-density lipoprotein, all LDL subclasses, and high concentrations of very large and large HDL particles. The glucose-increasing allele of rs174550 in FADS1 was significantly associated with high concentrations of very large and large HDL particles and nominally associated with low concentrations of all VLDL particles. SNPs rs10923931 in NOTCH2 and rs757210 in HNF1B genes showed nominal or significant associations with several lipoprotein traits. The genetic risk score of 34 SNPs was not associated with any of the lipoprotein subclasses. Four of the 34 risk loci for type 2 diabetes or hyperglycemia (GCKR, FADS1, NOTCH2, and HNF1B) were significantly associated with lipoprotein traits. A GCKR variant predominantly affected the concentration of VLDL, and the FADS1 variant affected very large and large HDL particles. Only a limited number of risk loci for hyperglycemia/type 2 diabetes significantly affect lipoprotein metabolism. Show less

Proinsulin is a precursor of mature insulin and C-peptide. Higher circulating proinsulin levels are associated with impaired β-cell function, raised glucose levels, insulin resistance, and type 2 diabetes (T2D). Studies of the insulin processing pathway could provide new insights about T2D pathophysiology. We have conducted a meta-analysis of genome-wide association tests of ∼2.5 million genotyped or imputed single nucleotide polymorphisms (SNPs) and fasting proinsulin levels in 10,701 nondiabetic adults of European ancestry, with follow-up of 23 loci in up to 16,378 individuals, using additive genetic models adjusted for age, sex, fasting insulin, and study-specific covariates. Nine SNPs at eight loci were associated with proinsulin levels (P < 5 × 10(-8)). Two loci (LARP6 and SGSM2) have not been previously related to metabolic traits, one (MADD) has been associated with fasting glucose, one (PCSK1) has been implicated in obesity, and four (TCF7L2, SLC30A8, VPS13C/C2CD4A/B, and ARAP1, formerly CENTD2) increase T2D risk. The proinsulin-raising allele of ARAP1 was associated with a lower fasting glucose (P = 1.7 × 10(-4)), improved β-cell function (P = 1.1 × 10(-5)), and lower risk of T2D (odds ratio 0.88; P = 7.8 × 10(-6)). Notably, PCSK1 encodes the protein prohormone convertase 1/3, the first enzyme in the insulin processing pathway. A genotype score composed of the nine proinsulin-raising alleles was not associated with coronary disease in two large case-control datasets. We have identified nine genetic variants associated with fasting proinsulin. Our findings illuminate the biology underlying glucose homeostasis and T2D development in humans and argue against a direct role of proinsulin in coronary artery disease pathogenesis. Show less

Obesity is globally prevalent and highly heritable, but its underlying genetic factors remain largely elusive. To identify genetic loci for obesity susceptibility, we examined associations between body mass index and ∼ 2.8 million SNPs in up to 123,865 individuals with targeted follow up of 42 SNPs in up to 125,931 additional individuals. We confirmed 14 known obesity susceptibility loci and identified 18 new loci associated with body mass index (P < 5 × 10⁻⁸), one of which includes a copy number variant near GPRC5B. Some loci (at MC4R, POMC, SH2B1 and BDNF) map near key hypothalamic regulators of energy balance, and one of these loci is near GIPR, an incretin receptor. Furthermore, genes in other newly associated loci may provide new insights into human body weight regulation. Show less

OBJECTIVE Recent genome-wide association studies have revealed loci associated with glucose and insulin-related traits. We aimed to characterize 19 such loci using detailed measures of insulin processing, secretion, and sensitivity to help elucidate their role in regulation of glucose control, insulin secretion and/or action. RESEARCH DESIGN AND METHODS We investigated associations of loci identified by the Meta-Analyses of Glucose and Insulin-related traits Consortium (MAGIC) with circulating proinsulin, measures of insulin secretion and sensitivity from oral glucose tolerance tests (OGTTs), euglycemic clamps, insulin suppression tests, or frequently sampled intravenous glucose tolerance tests in nondiabetic humans (n = 29,084). RESULTS The glucose-raising allele in MADD was associated with abnormal insulin processing (a dramatic effect on higher proinsulin levels, but no association with insulinogenic index) at extremely persuasive levels of statistical significance (P = 2.1 x 10(-71)). Defects in insulin processing and insulin secretion were seen in glucose-raising allele carriers at TCF7L2, SCL30A8, GIPR, and C2CD4B. Abnormalities in early insulin secretion were suggested in glucose-raising allele carriers at MTNR1B, GCK, FADS1, DGKB, and PROX1 (lower insulinogenic index; no association with proinsulin or insulin sensitivity). Two loci previously associated with fasting insulin (GCKR and IGF1) were associated with OGTT-derived insulin sensitivity indices in a consistent direction. CONCLUSIONS Genetic loci identified through their effect on hyperglycemia and/or hyperinsulinemia demonstrate considerable heterogeneity in associations with measures of insulin processing, secretion, and sensitivity. Our findings emphasize the importance of detailed physiological characterization of such loci for improved understanding of pathways associated with alterations in glucose homeostasis and eventually type 2 diabetes. Show less

Glucose levels 2 h after an oral glucose challenge are a clinical measure of glucose tolerance used in the diagnosis of type 2 diabetes. We report a meta-analysis of nine genome-wide association studies (n = 15,234 nondiabetic individuals) and a follow-up of 29 independent loci (n = 6,958-30,620). We identify variants at the GIPR locus associated with 2-h glucose level (rs10423928, beta (s.e.m.) = 0.09 (0.01) mmol/l per A allele, P = 2.0 x 10(-15)). The GIPR A-allele carriers also showed decreased insulin secretion (n = 22,492; insulinogenic index, P = 1.0 x 10(-17); ratio of insulin to glucose area under the curve, P = 1.3 x 10(-16)) and diminished incretin effect (n = 804; P = 4.3 x 10(-4)). We also identified variants at ADCY5 (rs2877716, P = 4.2 x 10(-16)), VPS13C (rs17271305, P = 4.1 x 10(-8)), GCKR (rs1260326, P = 7.1 x 10(-11)) and TCF7L2 (rs7903146, P = 4.2 x 10(-10)) associated with 2-h glucose. Of the three newly implicated loci (GIPR, ADCY5 and VPS13C), only ADCY5 was found to be associated with type 2 diabetes in collaborating studies (n = 35,869 cases, 89,798 controls, OR = 1.12, 95% CI 1.09-1.15, P = 4.8 x 10(-18)). Show less

Levels of circulating glucose are tightly regulated. To identify new loci influencing glycemic traits, we performed meta-analyses of 21 genome-wide association studies informative for fasting glucose, fasting insulin and indices of beta-cell function (HOMA-B) and insulin resistance (HOMA-IR) in up to 46,186 nondiabetic participants. Follow-up of 25 loci in up to 76,558 additional subjects identified 16 loci associated with fasting glucose and HOMA-B and two loci associated with fasting insulin and HOMA-IR. These include nine loci newly associated with fasting glucose (in or near ADCY5, MADD, ADRA2A, CRY2, FADS1, GLIS3, SLC2A2, PROX1 and C2CD4B) and one influencing fasting insulin and HOMA-IR (near IGF1). We also demonstrated association of ADCY5, PROX1, GCK, GCKR and DGKB-TMEM195 with type 2 diabetes. Within these loci, likely biological candidate genes influence signal transduction, cell proliferation, development, glucose-sensing and circadian regulation. Our results demonstrate that genetic studies of glycemic traits can identify type 2 diabetes risk loci, as well as loci containing gene variants that are associated with a modest elevation in glucose levels but are not associated with overt diabetes. Show less

Blood low-density lipoprotein (LDL) cholesterol, high-density lipoprotein (HDL) cholesterol and triglyceride levels are risk factors for cardiovascular disease. To dissect the polygenic basis of these traits, we conducted genome-wide association screens in 19,840 individuals and replication in up to 20,623 individuals. We identified 30 distinct loci associated with lipoprotein concentrations (each with P < 5 x 10(-8)), including 11 loci that reached genome-wide significance for the first time. The 11 newly defined loci include common variants associated with LDL cholesterol near ABCG8, MAFB, HNF1A and TIMD4; with HDL cholesterol near ANGPTL4, FADS1-FADS2-FADS3, HNF4A, LCAT, PLTP and TTC39B; and with triglycerides near AMAC1L2, FADS1-FADS2-FADS3 and PLTP. The proportion of individuals exceeding clinical cut points for high LDL cholesterol, low HDL cholesterol and high triglycerides varied according to an allelic dosage score (P < 10(-15) for each trend). These results suggest that the cumulative effect of multiple common variants contributes to polygenic dyslipidemia. Show less

We investigated the presence of left ventricular hypertrophy (LVH) and features of diastolic dysfunction in genotype-confirmed children from families with hypertrophic cardiomyopathy (HCM) and healthy control children. In subjects with HCM-causing mutations, LVH usually does not evolve until adolescence. Diastolic dysfunction has not been systematically evaluated in children carrying HCM-causing mutations. All children (aged 1.5-16.7 years) from 14 HCM families with identified disease-causing mutations (the Arg719Trp mutation in the beta-myosin heavy chain gene [MYH7], the Asp175Asn mutation in the alpha-tropomyosin gene [TPM1], the Gln1061X mutation in the myosin-binding protein C gene [MYBPC3], and the IVS5-2A-->C mutation in the MYBPC3 gene) and 53 matched control children were examined with electrocardiography and 2- and 3-dimensional echocardiography (2DE and 3DE). Natriuretic peptides were measured in children from HCM families and 67 control children. Of 53 children from HCM families, 27 (51%) had a disease-causing mutation (G+). G+ children had slightly thicker septum on 2DE compared with the control children (P = .004), but only 3 (11%) of 27 G+ children exceeded the 95th percentile values of the body surface area-adjusted maximal LV thickness of healthy children (the major echocardiographic criterion for HCM). However, prolonged isovolumetric relaxation time, increased left atrial volume on 3DE, or increased levels of NT-proANP, all features suggestive of diastolic dysfunction, were found in 14 (52%) of 27 G+ children. In children with HCM-causing mutations, signs of diastolic dysfunction are found in about half of the cases, as LVH is present only in small percentage of these children. Show less

Hypertrophic cardiomyopathy (HCM) is a genetically and clinically heterogeneous myocardial disease caused by mutations in genes encoding sarcomeric proteins. To assess the genetic background and phenotypic expression of HCM in eastern Finland, we screened 35 unrelated patients with HCM from the Kuopio University Hospital area for variants in 9 genes encoding sarcomeric proteins with the PCR-SSCP method. We herewith describe our previous findings in five sarcomeric genes and also report hitherto unpublished data on four additional sarcomeric genes. Mutations in the cardiac myosin-binding protein C gene (MYBPC3) were most frequent, accounting for 26% of cases. A novel mutation (Gln1061X) in this gene was the most common mutation, found in 6 of 35 families and accounting for 17% of all cases. Other novel mutations in MYBPC3 (IVS5-2A --> C, IVS14-13G --> A, and Ex25deltaLys) were found in one family each. A previously described alpha-tropomyosin (TPM1) mutation (Asp175Asn) was found in 11% of cases. Haplotype analysis suggested that the two most common variants (MYBPC3-Gln1061X and TPM1-Asp175Asn) were founder mutations. Only one mutation (Arg719Trp) in the beta-myosin heavy chain gene (MYH7) was found in one family, and no disease-causing mutations were found in the genes encoding alpha-actin, cardiac troponin I, T, C, or myosin essential and regulatory light chains. Altogether, the aforementioned 6 mutations found in MYBPC3, TPM1, and MYH7 accounted for 61% of familial and 40% of all HCM cases. The mutations were associated mostly with benign or intermediary phenotypes with only few HCM-related deaths. We conclude that the genetic profile of HCM in eastern Finland is unique, characterized by few founder mutations with benign or intermediary phenotypes. Show less

Hypertrophic cardiomyopathy (HCM) is a genetic disorder characterized by cardiac hypertrophy caused by mutations in genes encoding sarcomere proteins. This study screened all patients with HCM from the Kuopio University Hospital region in eastern Finland for variants in the cardiac myosin-binding protein C gene ( MYBPC3). All 35 exons of MYBPC3 were screened by the single-strand conformation polymorphism method in 37 unrelated patients with HCM. In MYBPC3 we identified seven novel (Gln1061X, IVS5-2A-->C, IVS14-13G-->A, Ex25DeltaLys, Pro147Leu, Ser236Gly, and Arg1138His) and two previously reported (Arg326Gln, Val896Met) variants, all of which are predicted to affect the structure of the encoded protein. Four of the nine variants, a nonsense mutation Gln1061X, a splice acceptor mutation (IVS5-2A-->C), a novel substitution in intron 14 (IVS14-13G-->A), and a novel 3-bp deletion in exon 25 (Ex25DeltaLys) were concluded to be disease-causing mutations because they cosegregated with the HCM phenotype or were absent in more than 200 normal chromosomes, or both. The mutation Gln1061X was found most frequently, being present in 6 families (23 subjects) while the other three mutations were found in single families each. Haplotype analysis indicated a likely founder effect among the families carrying the Gln1061X mutation. We found four novel mutations in MYBPC3, accounting for approx. 38% of familial and 24% of all cases of HCM. In our previous and unpublished studies no more frequent cause of HCM has been found in genetic analyses of other eight sarcomeric proteins. Consequently MYBPC3 is the predominant gene for HCM in eastern Finland. In addition, several amino acid substitutions in MYBPC3 suspected to be not associated with HCM were identified, indicating that some of the missense variants found in MYBPC3 are possibly not disease-causing mutations. Show less