The transcriptional landscape of age in human peripheral blood

Disease incidences increase with age, but the molecular characteristics of ageing that lead to increased disease susceptibility remain inadequately understood. Here we perform a whole-blood gene expression meta-analysis in 14,983 individuals of European ancestry (including replication) and identify 1,497 genes that are differentially expressed with chronological age. The age-associated genes do not harbor more age-associated CpG-methylation sites than other genes, but are instead enriched for the presence of potentially functional CpG-methylation sites in enhancer and insulator regions that associate with both chronological age and gene expression levels. We further used the gene expression profiles to calculate the ‘transcriptomic age' of an individual, and show that differences between transcriptomic age and chronological age are associated with biological features linked to ageing, such as blood pressure, cholesterol levels, fasting glucose, and body mass index. The transcriptomic prediction model adds biological relevance and complements existing epigenetic prediction models, and can be used by others to calculate transcriptomic age in external cohorts.

Ageing increases the risk of many diseases. Here the authors compare blood cell transcriptomes of over 14,000 individuals and identify a set of about 1,500 genes that are differently expressed with age, shedding light on transcriptional programs linked to the ageing process and age-associated diseases.

Characterization of whole-genome autosomal differences of DNA methylation between men and women

Background

Disease risk and incidence between males and females reveal differences, and sex is an important component of any investigation of the determinants of phenotypes or disease etiology. Further striking differences between men and women are known, for instance, at the metabolic level. The extent to which men and women vary at the level of the epigenome, however, is not well documented. DNA methylation is the best known epigenetic mechanism to date.

Results

In order to shed light on epigenetic differences, we compared autosomal DNA methylation levels between men and women in blood in a large prospective European cohort of 1799 subjects, and replicated our findings in three independent European cohorts. We identified and validated 1184 CpG sites to be differentially methylated between men and women and observed that these CpG sites were distributed across all autosomes. We showed that some of the differentially methylated loci also exhibit differential gene expression between men and women. Finally, we found that the differentially methylated loci are enriched among imprinted genes, and that their genomic location in the genome is concentrated in CpG island shores.

Conclusion

Our epigenome-wide association study indicates that differences between men and women are so substantial that they should be considered in design and analyses of future studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-015-0035-3) contains supplementary material, which is available to authorized users.

Extensive alterations of the whole-blood transcriptome are associated with body mass index: results of an mRNA profiling study involving two large population-based cohorts

Background

Obesity, defined as pathologically increased body mass index (BMI), is strongly related to an increased risk for numerous common cardiovascular and metabolic diseases. It is particularly associated with insulin resistance, hyperglycemia, and systemic oxidative stress and represents the most important risk factor for type 2 diabetes (T2D). However, the pathophysiological mechanisms underlying these associations are still not completely understood. Therefore, in order to identify potentially disease-relevant BMI-associated gene expression signatures, a transcriptome-wide association study (TWAS) on BMI was performed.

Methods

Whole-blood mRNA levels determined by array-based transcriptional profiling were correlated with BMI in two large independent population-based cohort studies (KORA F4 and SHIP-TREND) comprising a total of 1977 individuals.

Results

Extensive alterations of the whole-blood transcriptome were associated with BMI: More than 3500 transcripts exhibited significant positive or negative BMI-correlation. Three major whole-blood gene expression signatures associated with increased BMI were identified. The three signatures suggested: i) a ratio shift from mature erythrocytes towards reticulocytes, ii) decreased expression of several genes essentially involved in the transmission and amplification of the insulin signal, and iii) reduced expression of several key genes involved in the defence against reactive oxygen species (ROS).

Conclusions

Whereas the first signature confirms published results, the other two provide possible mechanistic explanations for well-known epidemiological findings under conditions of increased BMI, namely attenuated insulin signaling and increased oxidative stress. The putatively causative BMI-dependent down-regulation of the expression of numerous genes on the mRNA level represents a novel finding. BMI-associated negative transcriptional regulation of insulin signaling and oxidative stress management provide new insights into the pathogenesis of metabolic syndrome and T2D.

Electronic supplementary material

The online version of this article (doi:10.1186/s12920-015-0141-x) contains supplementary material, which is available to authorized users.

A Metabolome-Wide Association Study of Kidney Function and Disease in the General Population

Small molecules are extensively metabolized and cleared by the kidney. Changes in serum metabolite concentrations may result from impaired kidney function and can be used to estimate filtration (e.g., the established marker creatinine) or may precede and potentially contribute to CKD development. Here, we applied a nontargeted metabolomics approach using gas and liquid chromatography coupled to mass spectrometry to quantify 493 small molecules in human serum. The associations of these molecules with GFR estimated on the basis of creatinine (eGFRcr) and cystatin C levels were assessed in ≤1735 participants in the KORA F4 study, followed by replication in 1164 individuals in the TwinsUK registry. After correction for multiple testing, 54 replicated metabolites significantly associated with eGFRcr, and six of these showed pairwise correlation (r≥0.50) with established kidney function measures: C-mannosyltryptophan, pseudouridine, N-acetylalanine, erythronate, myo-inositol, and N-acetylcarnosine. Higher C-mannosyltryptophan, pseudouridine, and O-sulfo-L-tyrosine concentrations associated with incident CKD (eGFRcr <60 ml/min per 1.73 m(2)) in the KORA F4 study. In contrast with serum creatinine, C-mannosyltryptophan and pseudouridine concentrations showed little dependence on sex. Furthermore, correlation with measured GFR in 200 participants in the AASK study was 0.78 for both C-mannosyltryptophan and pseudouridine concentration, and highly significant associations of both metabolites with incident ESRD disappeared upon adjustment for measured GFR. Thus, these molecules may be alternative or complementary markers of kidney function. In conclusion, our study provides a comprehensive list of kidney function-associated metabolites and highlights potential novel filtration markers that may help to improve the estimation of GFR.

Childhood asthma is associated with mutations and gene expression differences of ORMDL genes that can interact

BACKGROUND

Genomewide association studies identified ORMDL3 as a plausible asthma candidate gene. ORMDL proteins regulate sphingolipid metabolism and ceramide homeostasis and participate in lymphocyte activation and eosinophil recruitment. Strong sequence homology between the three ORMDL genes and ORMDL protein conservation among different species suggest that they may have shared functions. We hypothesized that if single nucleotide polymorphisms (SNPs) in ORMDL3 alter its gene expression and play a role in asthma, variants in ORMDL1 and ORMDL2 might also be associated with asthma.

METHODS

Asthma associations of 44 genotyped SNPs were determined in at least 1303 subjects (651 asthmatics). ORMDL expression was evaluated in peripheral blood mononuclear cells (PBMC) from 55 subjects (eight asthmatics) before and after allergen stimulation, and in blood (n = 60, 5 asthmatics). Allele-specific cis-effects on ORMDL expression were assessed. Interactions between human ORMDL proteins were determined in living cells.

RESULTS

Sixteen SNPs in all three ORMDLs were associated with asthma (14 in ORMDL3). Baseline expression of ORMDL1 (P = 1.7 × 10(-6) ) and ORMDL2 (P = 4.9 × 10(-5) ) was significantly higher in PBMC from asthmatics, while induction of ORMDLs upon stimulation was stronger in nonasthmatics. Disease-associated alleles (rs8079416, rs4795405, rs3902920) alter ORMDL3 expression. ORMDL proteins formed homo- and heterooligomers and displayed similar patterns of interaction with SERCA2 and SPT1.

CONCLUSIONS

Polymorphisms in ORMDL genes are associated with asthma. Asthmatics exhibit increased ORMDL levels, suggesting that ORMDLs contribute to asthma. Formation of heterooligomers and similar interaction patterns with proteins involved in calcium homeostasis and sphingolipid metabolism could indicate shared biological roles of ORMDLs, influencing airway remodeling and hyperresponsiveness.

The Influence of Age and Sex on Genetic Associations with Adult Body Size and Shape: A Large-Scale Genome-Wide Interaction Study

Genome-wide association studies (GWAS) have identified more than 100 genetic variants contributing to BMI, a measure of body size, or waist-to-hip ratio (adjusted for BMI, WHR_adjBMI), a measure of body shape. Body size and shape change as people grow older and these changes differ substantially between men and women. To systematically screen for age- and/or sex-specific effects of genetic variants on BMI and WHR_adjBMI, we performed meta-analyses of 114 studies (up to 320,485 individuals of European descent) with genome-wide chip and/or Metabochip data by the Genetic Investigation of Anthropometric Traits (GIANT) Consortium. Each study tested the association of up to ~2.8M SNPs with BMI and WHR_adjBMI in four strata (men ≤50y, men >50y, women ≤50y, women >50y) and summary statistics were combined in stratum-specific meta-analyses. We then screened for variants that showed age-specific effects (G x AGE), sex-specific effects (G x SEX) or age-specific effects that differed between men and women (G x AGE x SEX). For BMI, we identified 15 loci (11 previously established for main effects, four novel) that showed significant (FDR<5%) age-specific effects, of which 11 had larger effects in younger (<50y) than in older adults (≥50y). No sex-dependent effects were identified for BMI. For WHR_adjBMI, we identified 44 loci (27 previously established for main effects, 17 novel) with sex-specific effects, of which 28 showed larger effects in women than in men, five showed larger effects in men than in women, and 11 showed opposite effects between sexes. No age-dependent effects were identified for WHR_adjBMI. This is the first genome-wide interaction meta-analysis to report convincing evidence of age-dependent genetic effects on BMI. In addition, we confirm the sex-specificity of genetic effects on WHR_adjBMI. These results may provide further insights into the biology that underlies weight change with age or the sexually dimorphism of body shape.

Adult body size and body shape differ substantially between men and women and change over time. More than 100 genetic variants that influence body mass index (measure of body size) or waist-to-hip ratio (measure of body shape) have been identified. While there is evidence that some genetic loci affect body shape differently in men than in women, little is known about whether genetic effects differ in older compared to younger adults, and whether such changes differ between men and women. Therefore, we conducted a systematic genome-wide search, including 114 studies (>320,000 individuals), to specifically identify genetic loci with age- and or sex-dependent effects on body size and shape. We identified 15 loci of which the effect on BMI was different in older compared to younger adults, whereas we found no evidence for loci with different effects in men compared to women. The opposite was seen for body shape as we identified 44 loci of which the effect on waist-to-hip ratio differed between men and women, but no difference between younger and older adults were observed. Our observations may provide new insights into the biology that underlies weight change with age or the sexual dimorphism of body shape.

Effects of metformin on metabolite profiles and LDL cholesterol in patients with type 2 diabetes

OBJECTIVE

Metformin is used as a first-line oral treatment for type 2 diabetes (T2D). However, the underlying mechanism is not fully understood. Here, we aimed to comprehensively investigate the pleiotropic effects of metformin.

RESEARCH DESIGN AND METHODS

We analyzed both metabolomic and genomic data of the population-based KORA cohort. To evaluate the effect of metformin treatment on metabolite concentrations, we quantified 131 metabolites in fasting serum samples and used multivariable linear regression models in three independent cross-sectional studies (n = 151 patients with T2D treated with metformin [mt-T2D]). Additionally, we used linear mixed-effect models to study the longitudinal KORA samples (n = 912) and performed mediation analyses to investigate the effects of metformin intake on blood lipid profiles. We combined genotyping data with the identified metformin-associated metabolites in KORA individuals (n = 1,809) and explored the underlying pathways.

RESULTS

We found significantly lower (P < 5.0E-06) concentrations of three metabolites (acyl-alkyl phosphatidylcholines [PCs]) when comparing mt-T2D with four control groups who were not using glucose-lowering oral medication. These findings were controlled for conventional risk factors of T2D and replicated in two independent studies. Furthermore, we observed that the levels of these metabolites decreased significantly in patients after they started metformin treatment during 7 years' follow-up. The reduction of these metabolites was also associated with a lowered blood level of LDL cholesterol (LDL-C). Variations of these three metabolites were significantly associated with 17 genes (including FADS1 and FADS2) and controlled by AMPK, a metformin target.

CONCLUSIONS

Our results indicate that metformin intake activates AMPK and consequently suppresses FADS, which leads to reduced levels of the three acyl-alkyl PCs and LDL-C. Our findings suggest potential beneficial effects of metformin in the prevention of cardiovascular disease.

Spinocerebellar ataxia type 36 exists in diverse populations and can be caused by a short hexanucleotide GGCCTG repeat expansion

OBJECTIVE

Spinocerebellar ataxia 36 (SCA36) is an autosomal-dominant neurodegenerative disorder caused by a large (>650) hexanucleotide GGCCTG repeat expansion in the first intron of the NOP56 gene. The aim of this study is to clarify the prevalence, clinical and genetic features of SCA36.

METHODS

The expansion was tested in 676 unrelated SCA index cases and 727 controls from France, Germany and Japan. Clinical and neuropathological features were investigated in available family members.

RESULTS

Normal alleles ranged between 5 and 14 hexanucleotide repeats. Expansions were detected in 12 families in France (prevalence: 1.9% of all French SCAs) including one family each with Spanish, Portuguese or Chinese ancestry, in five families in Japan (1.5% of all Japanese SCAs), but were absent in German patients. All the 17 SCA36 families shared one common haplotype for a 7.5 kb pairs region flanking the expansion. While 27 individuals had typically long expansions, three affected individuals harboured small hexanucleotide expansions of 25, 30 and 31 hexanucleotide repeat-units, demonstrating that such a small expansion could cause the disease. All patients showed slowly progressive cerebellar ataxia frequently accompanied by hearing and cognitive impairments, tremor, ptosis and reduced vibration sense, with the age at onset ranging between 39 and 65 years, and clinical features were indistinguishable between individuals with short and typically long expansions. Neuropathology in a presymptomatic case disclosed that Purkinje cells and hypoglossal neurons are affected.

CONCLUSIONS

SCA36 is rare with a worldwide distribution. It can be caused by a short GGCCTG expansion and associates various extracerebellar symptoms.

Feasibility and quality development of biomaterials in the pretest studies of the German National Cohort

BACKGROUND

The German National Cohort (GNC) is designed to address research questions concerning a wide range of possible causes of major chronic diseases (e.g. cancer, diabetes, infectious, allergic, neurologic and cardiovascular diseases) as well as to identify risk factors and prognostic biomarkers for early diagnosis and prevention of these diseases. The collection of biomaterials in combination with extensive information from questionnaires and medical examinations represents one of the central study components.

OBJECTIVES

In two pretest studies of the German National Cohort conducted between 2011 and 2013, a range of biomaterials from a defined number of participants was collected. Ten study centres were involved in pretest 1 and 18 study centres were involved in pretest 2. Standard operation procedures (SOP) were developed and evaluated to minimize pre-analytical artefacts during biosample collection. Within the pretest studies different aspects concerning feasibility of sample collection/preparation [pretest 1 (a)] and quality control of biomarkers and proteome analyses were investigated [pretest 1 (b), (c)]. Additionally, recruitment of study participants for specific projects and examination procedures of all study centres in a defined time period according to common standards as well as transportation and decentralized storage of biological samples were tested (pretest 2). These analyses will serve as the basis for the biomaterial collection in the main study of the GNC starting in 2014.

MATERIALS AND METHODS

Participants, randomly chosen from the population (n = 1000 subjects recruited at ten study sites in pretest 1) were asked to donate blood, urine, saliva and stool samples. Additionally, nasal and oropharyngeal swabs were collected at the study sites and nasal swabs were collected by the participants at home. SOPs for sample collection, preparation, storage and transportation were developed and adopted for pretest 2. In pretest 2, 18 study sites (n = 599 subjects) collected biomaterials mostly identical to pretest 1. Biomarker analyses to test the quality of the biomaterials were performed.

RESULTS

In pretest 1 and 2, it was feasible to collect all biomaterials from nearly all invited participants without major problems. The mean response rate of the subjects was 95 %. As one important result we found for example that after blood draw the cellular fraction should be separated from the plasma and serum fractions during the first hour with no significant variation for up to 6 h at 4 ℃ for all analysed biomarkers. Moreover, quality control of samples using a proteomics approach showed no significant clustering of proteins according to different storage conditions. All developed SOPs were validated for use in the main study after some adaptation and modification. Additionally, electronic and paper documentation sheets were developed and tested to record time stamps, volumes, freezing times, and aliquot numbers of the collected biomaterials.

DISCUSSION

The collection of the biomaterials was feasible without major problems at all participating study sites. However, the processing times were in some cases too long. To avoid pre-analytical artefacts in sample collection, appropriate standardisation among the study sites is necessary. To achieve this, blood and urine collection will have to be adapted to specific conditions of usage of liquid handling robots, which will be available at all participating study centres in the main study of the GNC. Strict compliance with the SOPs, thorough training of the staff and accurate documentation are mandatory to obtain high sample quality for later analyses. The so obtained biomaterials represent a valuable resource for research on infectious and other common complex diseases in the GNC.

Kindlin 3 (FERMT3) is associated with unstable atherosclerotic plaques, anti-inflammatory type II macrophages and upregulation of beta-2 integrins in all major arterial beds

BACKGROUND

Kindlins (FERMT) are cytoplasmic proteins required for integrin (ITG) activation, leukocyte transmigration, platelet aggregation and thrombosis. Characterization of kindlins and their association with atherosclerotic plaques in human(s) is lacking.

METHODS AND RESULTS

Exploratory microarray (MA) was first performed followed by selective quantitative validation of robustly expressed genes with qRT-PCR low-density array (LDA). In LDA, ITGA1 (1.30-fold, p = 0.041) and ITGB3 (1.37-fold, p = 0.036) were upregulated in whole blood samples of patients with coronary artery disease (CAD) compared to healthy controls. In arterial plaques, both robustly expressed transcript variants of FERMT3 (MA: 5.90- and 3.4-fold; LDA: 3.99-fold, p < 0.0001 for all) and ITGB2 (MA: 4.81- and 4.92-fold; LDA: 5.29-fold, p < 0.0001 for all) were upregulated while FERMT2 was downregulated (MA: -1.61-fold; LDA: -2.88-fold, p < 0.0001 for both). The other integrins (ITGA1, ITGAV, ITGB3, ITGB5) were downregulated. All these results were replicated in at least one arterial bed. The latter FERMT3 transcript variant associated with unstable plaques (p = 0.0004). FERMT3 correlated with M2 macrophage markers and in hierarchical cluster analysis clustered with inflammatory and macrophage markers, while FERMT2 correlated with SMC-rich plaque markers and clustered with SMC markers. In confocal immunofluorescence analysis, FERMT3 protein colocalized with abundant CD68-positive cells of monocytic origin in the atherosclerotic plaques, while co-localization of FERMT3 with HHF35 indicative of smooth muscle cells was low.

CONCLUSIONS

Kindlin-3 (FERMT3) is upregulated in atherosclerotic, especially unstable plaques, mainly in cells of monocytic origin and of M2 type. Simultaneous upregulation of ITGB2 suggests a synergistic effect on leukocyte adherence and transmigration into the vessel wall.

Epigenome-wide association of DNA methylation markers in peripheral blood from Indian Asians and Europeans with incident type 2 diabetes: a nested case-control study

BACKGROUND

Indian Asians, who make up a quarter of the world's population, are at high risk of developing type 2 diabetes. We investigated whether DNA methylation is associated with future type 2 diabetes incidence in Indian Asians and whether differences in methylation patterns between Indian Asians and Europeans are associated with, and could be used to predict, differences in the magnitude of risk of developing type 2 diabetes.

METHODS

We did a nested case-control study of DNA methylation in Indian Asians and Europeans with incident type 2 diabetes who were identified from the 8-year follow-up of 25 372 participants in the London Life Sciences Prospective Population (LOLIPOP) study. Patients were recruited between May 1, 2002, and Sept 12, 2008. We did epigenome-wide association analysis using samples from Indian Asians with incident type 2 diabetes and age-matched and sex-matched Indian Asian controls, followed by replication testing of top-ranking signals in Europeans. For both discovery and replication, DNA methylation was measured in the baseline blood sample, which was collected before the onset of type 2 diabetes. Epigenome-wide significance was set at p<1 × 10(-7). We compared methylation levels between Indian Asian and European controls without type 2 diabetes at baseline to estimate the potential contribution of DNA methylation to increased risk of future type 2 diabetes incidence among Indian Asians.

FINDINGS

1608 (11·9%) of 13 535 Indian Asians and 306 (4·3%) of 7066 Europeans developed type 2 diabetes over a mean of 8·5 years (SD 1·8) of follow-up. The age-adjusted and sex-adjusted incidence of type 2 diabetes was 3·1 times (95% CI 2·8-3·6; p<0·0001) higher among Indian Asians than among Europeans, and remained 2·5 times (2·1-2·9; p<0·0001) higher after adjustment for adiposity, physical activity, family history of type 2 diabetes, and baseline glycaemic measures. The mean absolute difference in methylation level between type 2 diabetes cases and controls ranged from 0·5% (SD 0·1) to 1·1% (0·2). Methylation markers at five loci were associated with future type 2 diabetes incidence; the relative risk per 1% increase in methylation was 1·09 (95% CI 1·07-1·11; p=1·3 × 10(-17)) for ABCG1, 0·94 (0·92-0·95; p=4·2 × 10(-11)) for PHOSPHO1, 0·94 (0·92-0·96; p=1·4 × 10(-9)) for SOCS3, 1·07 (1·04-1·09; p=2·1 × 10(-10)) for SREBF1, and 0·92 (0·90-0·94; p=1·2 × 10(-17)) for TXNIP. A methylation score combining results for the five loci was associated with future type 2 diabetes incidence (relative risk quartile 4 vs quartile 1 3·51, 95% CI 2·79-4·42; p=1·3 × 10(-26)), and was independent of established risk factors. Methylation score was higher among Indian Asians than Europeans (p=1 × 10(-34)).

INTERPRETATION

DNA methylation might provide new insights into the pathways underlying type 2 diabetes and offer new opportunities for risk stratification and prevention of type 2 diabetes among Indian Asians.

FUNDING

The European Union, the UK National Institute for Health Research, the Wellcome Trust, the UK Medical Research Council, Action on Hearing Loss, the UK Biotechnology and Biological Sciences Research Council, the Oak Foundation, the Economic and Social Research Council, Helmholtz Zentrum Munchen, the German Research Center for Environmental Health, the German Federal Ministry of Education and Research, the German Center for Diabetes Research, the Munich Center for Health Sciences, the Ministry of Science and Research of the State of North Rhine-Westphalia, and the German Federal Ministry of Health.

Polymorphisms related to ORMDL3 are associated with asthma susceptibility, alterations in transcriptional regulation of ORMDL3, and changes in TH2 cytokine levels

BACKGROUND

Chromosome 17q21, harboring the orosomucoid 1-like 3 (ORMDL3) gene, has been consistently associated with childhood asthma in genome-wide association studies.

OBJECTIVE

We investigated genetic variants in and around ORMDL3 that can change the function of ORMDL3 and thus contribute to asthma susceptibility.

METHODS

We performed haplotype analyses and fine mapping of the ORMDL3 locus in a cross-sectional (International Study of Asthma and Allergies in Childhood Phase II, n = 3557 total subjects, n = 281 asthmatic patients) and case-control (Multicenter Asthma Genetics in Childhood Study/International Study of Asthma and Allergies in Childhood Phase II, n = 1446 total subjects, n = 763 asthmatic patients) data set to identify putative causal single nucleotide polymorphisms (SNPs) in the locus. Top asthma-associated polymorphisms were analyzed for allele-specific effects on transcription factor binding and promoter activity in vitro and gene expression in PBMCs after stimulation ex vivo.

RESULTS

Two haplotypes (H1 and H2) were significantly associated with asthma in the cross-sectional (P = 9.9 × 10(-5) and P = .0035, respectively) and case-control (P = 3.15 × 10(-8) and P = .0021, respectively) populations. Polymorphisms rs8076131 and rs4065275 were identified to drive these effects. For rs4065275, a quantitative difference in transcription factor binding was found, whereas for rs8076131, changes in upstream stimulatory factor 1 and 2 transcription factor binding were observed in vitro by using different cell lines and PBMCs. This might contribute to detected alterations in luciferase activity paralleled with changes in ORMDL3 gene expression and IL-4 and IL-13 cytokine levels ex vivo in response to innate and adaptive stimuli in an allele-specific manner. Both SNPs were in strong linkage disequilibrium with asthma-associated 17q21 SNPs previously related to altered ORMDL3 gene expression.

CONCLUSION

Polymorphisms in a putative promoter region of ORMDL3, which are associated with childhood asthma, alter transcriptional regulation of ORMDL3, correlate with changes in TH2 cytokines levels, and therefore might contribute to the childhood asthma susceptibility signal from 17q21.

A genome-wide expression quantitative trait loci analysis of proprotein convertase subtilisin/kexin enzymes identifies a novel regulatory gene variant for FURIN expression and blood pressure

Proprotein convertase subtilisin/kexin (PCSK) enzymes cleave and convert their immature substrates into biologically active forms. Polymorphisms in the PCSK genes have been reported to associate with human diseases and phenotypes, including hypercholesterolemia and blood pressure (BP), and targeting PCSKs is considered a promising future form of drug therapy. PCSK processing is readily induced upon upregulation of the enzyme, but the genetic factors contributing to PCSK expression have not been thoroughly characterized. To gain a comprehensive understanding of the genetic regulation of PCSK expression, we performed, for the first time, a genome-wide expression quantitative trait loci (eQTL) analysis using mRNA expression in >1400 human peripheral blood samples from the Cardiovascular Risk in Young Finns Study and ca. ten million single-nucleotide polymorphisms (SNPs). The expression data showed clear expression for FURIN, PCSK5, PCSK7 and MBTPS1 (membrane-bound transcription factor peptidase, site 1) mRNAs in virtually all tested samples. A discovery analysis demonstrated a genome-wide significant (p < 5 × 10(-8)) association with the selected PCSK probes for 1024 variants, which were located at ten independent loci. Of these loci, 5/10 could be confirmed to regulate PCSK expression in two additional and independent sample sets. Finally, a phenotypic analysis demonstrated that a novel cis-eQTL SNP rs4702 for FURIN is strongly associated with both diastolic (p = 0.012) and systolic (p = 0.035) BP levels, as well as peripheral vascular resistance (p = 0.003). These findings indicate that the expression of the PCSK enzymes is regulated by genetic factors, which have biological roles in health and disease.

Multi-omic signature of body weight change: results from a population-based cohort study

BACKGROUND

Excess body weight is a major risk factor for cardiometabolic diseases. The complex molecular mechanisms of body weight change-induced metabolic perturbations are not fully understood. Specifically, in-depth molecular characterization of long-term body weight change in the general population is lacking. Here, we pursued a multi-omic approach to comprehensively study metabolic consequences of body weight change during a seven-year follow-up in a large prospective study.

METHODS

We used data from the population-based Cooperative Health Research in the Region of Augsburg (KORA) S4/F4 cohort. At follow-up (F4), two-platform serum metabolomics and whole blood gene expression measurements were obtained for 1,631 and 689 participants, respectively. Using weighted correlation network analysis, omics data were clustered into modules of closely connected molecules, followed by the formation of a partial correlation network from the modules. Association of the omics modules with previous annual percentage weight change was then determined using linear models. In addition, we performed pathway enrichment analyses, stability analyses, and assessed the relation of the omics modules with clinical traits.

RESULTS

Four metabolite and two gene expression modules were significantly and stably associated with body weight change (P-values ranging from 1.9 × 10(-4) to 1.2 × 10(-24)). The four metabolite modules covered major branches of metabolism, with VLDL, LDL and large HDL subclasses, triglycerides, branched-chain amino acids and markers of energy metabolism among the main representative molecules. One gene expression module suggests a role of weight change in red blood cell development. The other gene expression module largely overlaps with the lipid-leukocyte (LL) module previously reported to interact with serum metabolites, for which we identify additional co-expressed genes. The omics modules were interrelated and showed cross-sectional associations with clinical traits. Moreover, weight gain and weight loss showed largely opposing associations with the omics modules.

CONCLUSIONS

Long-term weight change in the general population globally associates with serum metabolite concentrations. An integrated metabolomics and transcriptomics approach improved the understanding of molecular mechanisms underlying the association of weight gain with changes in lipid and amino acid metabolism, insulin sensitivity, mitochondrial function as well as blood cell development and function.

Metabolite profiling and cardiovascular event risk: a prospective study of 3 population-based cohorts

BACKGROUND

High-throughput profiling of circulating metabolites may improve cardiovascular risk prediction over established risk factors.

METHODS AND RESULTS

We applied quantitative nuclear magnetic resonance metabolomics to identify the biomarkers for incident cardiovascular disease during long-term follow-up. Biomarker discovery was conducted in the National Finnish FINRISK study (n=7256; 800 events). Replication and incremental risk prediction was assessed in the Southall and Brent Revisited (SABRE) study (n=2622; 573 events) and British Women's Health and Heart Study (n=3563; 368 events). In targeted analyses of 68 lipids and metabolites, 33 measures were associated with incident cardiovascular events at P<0.0007 after adjusting for age, sex, blood pressure, smoking, diabetes mellitus, and medication. When further adjusting for routine lipids, 4 metabolites were associated with future cardiovascular events in meta-analyses: higher serum phenylalanine (hazard ratio per standard deviation, 1.18; 95% confidence interval, 1.12-1.24; P=4×10(-10)) and monounsaturated fatty acid levels (1.17; 1.11-1.24; P=1×10(-8)) were associated with increased cardiovascular risk, while higher omega-6 fatty acids (0.89; 0.84-0.94; P=6×10(-5)) and docosahexaenoic acid levels (0.90; 0.86-0.95; P=5×10(-5)) were associated with lower risk. A risk score incorporating these 4 biomarkers was derived in FINRISK. Risk prediction estimates were more accurate in the 2 validation cohorts (relative integrated discrimination improvement, 8.8% and 4.3%), albeit discrimination was not enhanced. Risk classification was particularly improved for persons in the 5% to 10% risk range (net reclassification, 27.1% and 15.5%). Biomarker associations were further corroborated with mass spectrometry in FINRISK (n=671) and the Framingham Offspring Study (n=2289).

CONCLUSIONS

Metabolite profiling in large prospective cohorts identified phenylalanine, monounsaturated fatty acids, and polyunsaturated fatty acids as biomarkers for cardiovascular risk. This study substantiates the value of high-throughput metabolomics for biomarker discovery and improved risk assessment.

Rationale, design and objectives of ARegPKD, a European ARPKD registry study

BACKGROUND

Autosomal recessive polycystic kidney disease (ARPKD) is a rare but frequently severe disorder that is typically characterized by cystic kidneys and congenital hepatic fibrosis but displays pronounced phenotypic heterogeneity. ARPKD is among the most important causes for pediatric end stage renal disease and a leading reason for liver-, kidney- or combined liver kidney transplantation in childhood. The underlying pathophysiology, the mechanisms resulting in the observed clinical heterogeneity and the long-term clinical evolution of patients remain poorly understood. Current treatment approaches continue to be largely symptomatic and opinion-based even in most-advanced medical centers. While large clinical trials for the frequent and mostly adult onset autosomal dominant polycystic kidney diseases have recently been conducted, therapeutic initiatives for ARPKD are facing the challenge of small and clinically variable cohorts for which reliable end points are hard to establish.

METHODS/DESIGN

ARegPKD is an international, mostly European, observational study to deeply phenotype ARPKD patients in a pro- and retrospective fashion. This registry study is conducted with the support of the German Society for Pediatric Nephrology (GPN) and the European Study Consortium for Chronic Kidney Disorders Affecting Pediatric Patients (ESCAPE Network). ARegPKD clinically characterizes long-term ARPKD courses by a web-based approach that uses detailed basic data questionnaires in combination with yearly follow-up visits. Clinical data collection is accompanied by associated biobanking and reference histology, thus setting roots for future translational research.

DISCUSSION

The novel registry study ARegPKD aims to characterize miscellaneous subcohorts and to compare the applied treatment options in a large cohort of deeply characterized patients. ARegPKD will thus provide evidence base for clinical treatment decisions and contribute to the pathophysiological understanding of this severe inherited disorder.

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

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

DNA methylation of lipid-related genes affects blood lipid levels

BACKGROUND

Epigenetic mechanisms might be involved in the regulation of interindividual lipid level variability and thus may contribute to the cardiovascular risk profile. The aim of this study was to investigate the association between genome-wide DNA methylation and blood lipid levels high-density lipoprotein cholesterol, low-density lipoprotein cholesterol, triglycerides, and total cholesterol. Observed DNA methylation changes were also further analyzed to examine their relationship with previous hospitalized myocardial infarction.

METHODS AND RESULTS

Genome-wide DNA methylation patterns were determined in whole blood samples of 1776 subjects of the Cooperative Health Research in the Region of Augsburg F4 cohort using the Infinium HumanMethylation450 BeadChip (Illumina). Ten novel lipid-related CpG sites annotated to various genes including ABCG1, MIR33B/SREBF1, and TNIP1 were identified. CpG cg06500161, located in ABCG1, was associated in opposite directions with both high-density lipoprotein cholesterol (β coefficient=-0.049; P=8.26E-17) and triglyceride levels (β=0.070; P=1.21E-27). Eight associations were confirmed by replication in the Cooperative Health Research in the Region of Augsburg F3 study (n=499) and in the Invecchiare in Chianti, Aging in the Chianti Area study (n=472). Associations between triglyceride levels and SREBF1 and ABCG1 were also found in adipose tissue of the Multiple Tissue Human Expression Resource cohort (n=634). Expression analysis revealed an association between ABCG1 methylation and lipid levels that might be partly mediated by ABCG1 expression. DNA methylation of ABCG1 might also play a role in previous hospitalized myocardial infarction (odds ratio, 1.15; 95% confidence interval=1.06-1.25).

CONCLUSIONS

Epigenetic modifications of the newly identified loci might regulate disturbed blood lipid levels and thus contribute to the development of complex lipid-related diseases.

Genome-wide comparative analysis of atopic dermatitis and psoriasis gives insight into opposing genetic mechanisms

Atopic dermatitis and psoriasis are the two most common immune-mediated inflammatory disorders affecting the skin. Genome-wide studies demonstrate a high degree of genetic overlap, but these diseases have mutually exclusive clinical phenotypes and opposing immune mechanisms. Despite their prevalence, atopic dermatitis and psoriasis very rarely co-occur within one individual. By utilizing genome-wide association study and ImmunoChip data from >19,000 individuals and methodologies developed from meta-analysis, we have identified opposing risk alleles at shared loci as well as independent disease-specific loci within the epidermal differentiation complex (chromosome 1q21.3), the Th2 locus control region (chromosome 5q31.1), and the major histocompatibility complex (chromosome 6p21-22). We further identified previously unreported pleiotropic alleles with opposing effects on atopic dermatitis and psoriasis risk in PRKRA and ANXA6/TNIP1. In contrast, there was no evidence for shared loci with effects operating in the same direction on both diseases. Our results show that atopic dermatitis and psoriasis have distinct genetic mechanisms with opposing effects in shared pathways influencing epidermal differentiation and immune response. The statistical analysis methods developed in the conduct of this study have produced additional insight from previously published data sets. The approach is likely to be applicable to the investigation of the genetic basis of other complex traits with overlapping and distinct clinical features.

Gender-specific pathway differences in the human serum metabolome

The susceptibility for various diseases as well as the response to treatments differ considerably between men and women. As a basis for a gender-specific personalized healthcare, an extensive characterization of the molecular differences between the two genders is required. In the present study, we conducted a large-scale metabolomics analysis of 507 metabolic markers measured in serum of 1756 participants from the German KORA F4 study (903 females and 853 males). One-third of the metabolites show significant differences between males and females. A pathway analysis revealed strong differences in steroid metabolism, fatty acids and further lipids, a large fraction of amino acids, oxidative phosphorylation, purine metabolism and gamma-glutamyl dipeptides. We then extended this analysis by a network-based clustering approach. Metabolite interactions were estimated using Gaussian graphical models to get an unbiased, fully data-driven metabolic network representation. This approach is not limited to possibly arbitrary pathway boundaries and can even include poorly or uncharacterized metabolites. The network analysis revealed several strongly gender-regulated submodules across different pathways. Finally, a gender-stratified genome-wide association study was performed to determine whether the observed gender differences are caused by dimorphisms in the effects of genetic polymorphisms on the metabolome. With only a single genome-wide significant hit, our results suggest that this scenario is not the case. In summary, we report an extensive characterization and interpretation of gender-specific differences of the human serum metabolome, providing a broad basis for future analyses.

No evidence for genome-wide interactions on plasma fibrinogen by smoking, alcohol consumption and body mass index: results from meta-analyses of 80,607 subjects

Plasma fibrinogen is an acute phase protein playing an important role in the blood coagulation cascade having strong associations with smoking, alcohol consumption and body mass index (BMI). Genome-wide association studies (GWAS) have identified a variety of gene regions associated with elevated plasma fibrinogen concentrations. However, little is yet known about how associations between environmental factors and fibrinogen might be modified by genetic variation. Therefore, we conducted large-scale meta-analyses of genome-wide interaction studies to identify possible interactions of genetic variants and smoking status, alcohol consumption or BMI on fibrinogen concentration. The present study included 80,607 subjects of European ancestry from 22 studies. Genome-wide interaction analyses were performed separately in each study for about 2.6 million single nucleotide polymorphisms (SNPs) across the 22 autosomal chromosomes. For each SNP and risk factor, we performed a linear regression under an additive genetic model including an interaction term between SNP and risk factor. Interaction estimates were meta-analysed using a fixed-effects model. No genome-wide significant interaction with smoking status, alcohol consumption or BMI was observed in the meta-analyses. The most suggestive interaction was found for smoking and rs10519203, located in the LOC123688 region on chromosome 15, with a p value of 6.2×10⁻⁸. This large genome-wide interaction study including 80,607 participants found no strong evidence of interaction between genetic variants and smoking status, alcohol consumption or BMI on fibrinogen concentrations. Further studies are needed to yield deeper insight in the interplay between environmental factors and gene variants on the regulation of fibrinogen concentrations.