Novel associations for coronary artery disease derived from genome wide association studies are not associated with increased carotid intima-media thickness, suggesting they do not act via early atherosclerosis or vessel remodeling

BACKGROUND

Recent genome-wide association studies (GWAS) have identified associations with myocardial infarction and coronary artery disease (CAD), but the mechanisms underlying these associations remain largely unclear. Carotid intima-media thickness (IMT) is a measure of early arterial remodeling and arteriosclerosis. Therefore, if CAD associated SNPs are also associated with carotid IMT; it suggests that they are acting via the early stages of the atherosclerotic process.

METHODS

In three large community based independent populations (CAPS, KORA and Young Finns) of European ancestry in which common carotid IMT had been measured (total 4961 individuals), we determined whether SNPs that have been associated with CAD in GWAS studies are also associated with carotid IMT. Associations with plaque were not examined.

RESULTS

We identified 11 SNPs and one haplotype previously associated with CAD. None of these were associated with common carotid IMT.

CONCLUSIONS

We found no evidence that SNPs associated with CAD on GWAS are also associated with carotid IMT. This suggests these genetic associations are not acting via early vessel remodeling or early arteriosclerosis.

Discovery of sexual dimorphisms in metabolic and genetic biomarkers

Metabolomic profiling and the integration of whole-genome genetic association data has proven to be a powerful tool to comprehensively explore gene regulatory networks and to investigate the effects of genetic variation at the molecular level. Serum metabolite concentrations allow a direct readout of biological processes, and association of specific metabolomic signatures with complex diseases such as Alzheimer's disease and cardiovascular and metabolic disorders has been shown. There are well-known correlations between sex and the incidence, prevalence, age of onset, symptoms, and severity of a disease, as well as the reaction to drugs. However, most of the studies published so far did not consider the role of sexual dimorphism and did not analyse their data stratified by gender. This study investigated sex-specific differences of serum metabolite concentrations and their underlying genetic determination. For discovery and replication we used more than 3,300 independent individuals from KORA F3 and F4 with metabolite measurements of 131 metabolites, including amino acids, phosphatidylcholines, sphingomyelins, acylcarnitines, and C6-sugars. A linear regression approach revealed significant concentration differences between males and females for 102 out of 131 metabolites (p-values<3.8×10(-4); Bonferroni-corrected threshold). Sex-specific genome-wide association studies (GWAS) showed genome-wide significant differences in beta-estimates for SNPs in the CPS1 locus (carbamoyl-phosphate synthase 1, significance level: p<3.8×10(-10); Bonferroni-corrected threshold) for glycine. We showed that the metabolite profiles of males and females are significantly different and, furthermore, that specific genetic variants in metabolism-related genes depict sexual dimorphism. Our study provides new important insights into sex-specific differences of cell regulatory processes and underscores that studies should consider sex-specific effects in design and interpretation.

Dopamine-related genes and spontaneous smoking cessation in ever-heavy smokers

Several studies have provided evidence for associations of polymorphisms located in and near dopamine-related genes and nicotine dependence and other smoking-related phenotypes, including pharmacogenetic interactions.

AIM

The purpose of the present work was to examine the association of SNPs in the DOPA decarboxylase (DDC), dopamine receptor D2 (DRD2) and dopamine transporter (SLC6A3) genes with smoking cessation in a large retrospective study featuring approximately 900 cessation events.

MATERIALS & METHODS

Data originated from the enrollment questionnaire of the epidemiological ESTHER study of community-dwelling adults aged 50-74 years, conducted in the German state of Saarland between July 2000 and December 2002. Restricting the analyses to subjects who reported to have regularly smoked > 20 cigarettes per day at some point in their life, we used survival analysis methods to model the time from initiation of regular smoking to cessation (defined as quitting with abstinence lasting until enrollment) and its relation with eight polymorphisms in the aforementioned genes (five in DDC, two in DRD2 and one in SLC6A3) in 1446 participants.

RESULTS

Neither individual variants nor DDC haplotypes were associated with the probability of overcoming nicotine dependence in this cohort.

CONCLUSION

The repeated suggestion of associations between the variants examined and nicotine dependence in previous reports seems to contrast the negative results in the present study. This would appear consistent with the hypothesis that the establishment of regular heavy smoking might abolish associations between genetic determinants of nicotine dependence and nicotine dependence-related phenotypes, in particular the probability of successful smoking cessation.

miR-21, miR-210, miR-34a, and miR-146a/b are up-regulated in human atherosclerotic plaques in the Tampere Vascular Study

OBJECTIVE

MicroRNAs are small non-coding RNAs that inversely regulate their target gene expression. The whole miRNA profile of human atherosclerotic plaques has not been studied previously. The aim of this study was to investigate the miRNA expression profile in human atherosclerotic plaques as compared to non-atherosclerotic left internal thoracic arteries (LITA), and to connect this expression to the processes in atherosclerosis.

METHODS

The miRNA expression profiles of six LITAs and 12 atherosclerotic plaques obtained from aortic, carotid, and femoral atherosclerotic arteries from Tampere Vascular Study were analyzed. The analyses were performed with Agilent's miRNA Microarray. The expression levels of over 4-fold up-regulated miRNAs were verified with qRT-PCR from a larger population (n=50). Messenger RNA levels were analyzed with Illumina's Expression BeadChip to study miRNA target expression.

RESULTS

Ten miRNAs were found to be differently expressed in atherosclerotic plaques when compared to controls (p<0.05). The expression of miR-21, -34a, -146a, -146b-5p, and -210 was verified and found to be significantly up-regulated in atherosclerotic arteries versus LITAs (p<0.001, fold changes 4.61, 2.55, 2.87, 2.82, and 3.92, respectively). Several predicted targets of these miRNAs were down-regulated, and gene set enrichment analysis showed several pathways which could be differently expressed due to this miRNA profile.

CONCLUSIONS

The microRNA expression profile differs significantly between atherosclerotic plaques and control arteries. The most up-regulated miRNAs are involved in processes known to be connected to atherosclerosis. Interfering with the miRNA expression in the artery wall is a potential way to affect atherosclerotic plaque and cardiovascular disease development.

Differences between human plasma and serum metabolite profiles

BACKGROUND

Human plasma and serum are widely used matrices in clinical and biological studies. However, different collecting procedures and the coagulation cascade influence concentrations of both proteins and metabolites in these matrices. The effects on metabolite concentration profiles have not been fully characterized.

METHODOLOGY/PRINCIPAL FINDINGS

We analyzed the concentrations of 163 metabolites in plasma and serum samples collected simultaneously from 377 fasting individuals. To ensure data quality, 41 metabolites with low measurement stability were excluded from further analysis. In addition, plasma and corresponding serum samples from 83 individuals were re-measured in the same plates and mean correlation coefficients (r) of all metabolites between the duplicates were 0.83 and 0.80 in plasma and serum, respectively, indicating significantly better stability of plasma compared to serum (p = 0.01). Metabolite profiles from plasma and serum were clearly distinct with 104 metabolites showing significantly higher concentrations in serum. In particular, 9 metabolites showed relative concentration differences larger than 20%. Despite differences in absolute concentration between the two matrices, for most metabolites the overall correlation was high (mean r = 0.81±0.10), which reflects a proportional change in concentration. Furthermore, when two groups of individuals with different phenotypes were compared with each other using both matrices, more metabolites with significantly different concentrations could be identified in serum than in plasma. For example, when 51 type 2 diabetes (T2D) patients were compared with 326 non-T2D individuals, 15 more significantly different metabolites were found in serum, in addition to the 25 common to both matrices.

CONCLUSIONS/SIGNIFICANCE

Our study shows that reproducibility was good in both plasma and serum, and better in plasma. Furthermore, as long as the same blood preparation procedure is used, either matrix should generate similar results in clinical and biological studies. The higher metabolite concentrations in serum, however, make it possible to provide more sensitive results in biomarker detection.

Association of single nucleotide polymorphisms in the diamine oxidase gene with diamine oxidase serum activities

BACKGROUND

Histamine intolerance (HIT) is associated with an excess of histamine because of an impaired function of the histamine-degrading enzyme diamine oxidase (DAO). The genetic background of HIT is unknown yet.

METHODS

Case-control association study of all haplotype tagging and four previously reported DAO SNPs and one HNMT Single nucleotide polymorphism with symptoms of HIT and DAO serum activity in 484 German individuals including 285 patients with clinical symptoms of HIT and 199 controls.

RESULTS

Diamine oxidase serum activity was significantly associated with seven SNPs within the DAO gene. The minor allele at rs2052129, rs2268999, rs10156191 and rs1049742 increased the risk for a reduced DAO activity whereas showing a moderate protective effect at rs2071514, rs1049748 and rs2071517 in the genotypic (P = 2.1 × 10(-8) , 7.6 × 10(-10) , 8.3 × 10(-10) , 0.009, 0.005, 0.00001, 0.006, respectively) and allelic genetic model (P = 2.5 × 10(-11) , 5.4 × 10(-13) , 8.9 × 10(-13) , 0.00002, 0.006, 0.0003, 0.005, respectively). Reporter gene assays at rs2052129 revealed a lower promoter activity (P = 0.016) of the minor allele. DAO mRNA expression in peripheral blood mononuclear cells of homozygous carriers of the minor allele at rs2052129, rs2268999, rs10156191 was lower (P = 0.002) than homozygous carriers of the major allele. Diamine oxidase variants were not associated with the HIT phenotype per se, only with DAO activity alone and the subgroup of HIT patients displaying a reduced DAO activity.

CONCLUSIONS

DAO gene variants strongly influence DAO expression and activity but alone are not sufficient to fully effectuate the potentially associated disease state of HIT, suggesting an interplay of genetic and environmental factors.

Multiple loci are associated with white blood cell phenotypes

White blood cell (WBC) count is a common clinical measure from complete blood count assays, and it varies widely among healthy individuals. Total WBC count and its constituent subtypes have been shown to be moderately heritable, with the heritability estimates varying across cell types. We studied 19,509 subjects from seven cohorts in a discovery analysis, and 11,823 subjects from ten cohorts for replication analyses, to determine genetic factors influencing variability within the normal hematological range for total WBC count and five WBC subtype measures. Cohort specific data was supplied by the CHARGE, HeamGen, and INGI consortia, as well as independent collaborative studies. We identified and replicated ten associations with total WBC count and five WBC subtypes at seven different genomic loci (total WBC count—6p21 in the HLA region, 17q21 near ORMDL3, and CSF3; neutrophil count—17q21; basophil count- 3p21 near RPN1 and C3orf27; lymphocyte count—6p21, 19p13 at EPS15L1; monocyte count—2q31 at ITGA4, 3q21, 8q24 an intergenic region, 9q31 near EDG2), including three previously reported associations and seven novel associations. To investigate functional relationships among variants contributing to variability in the six WBC traits, we utilized gene expression- and pathways-based analyses. We implemented gene-clustering algorithms to evaluate functional connectivity among implicated loci and showed functional relationships across cell types. Gene expression data from whole blood was utilized to show that significant biological consequences can be extracted from our genome-wide analyses, with effect estimates for significant loci from the meta-analyses being highly corellated with the proximal gene expression. In addition, collaborative efforts between the groups contributing to this study and related studies conducted by the COGENT and RIKEN groups allowed for the examination of effect homogeneity for genome-wide significant associations across populations of diverse ancestral backgrounds.

WBC traits are highly variable, moderately heritable, and commonly assayed as part of clinical complete blood count (CBC) examinations. The counts of constituent cell subtypes comprising the WBC count measure are assayed as part of a standard clinical WBC differential test. In this study we employed meta-analytic techniques and identified ten associations with WBC measures at seven genomic loci in a large sample set of over 31,000 participants. Cohort specific data was supplied by the CHARGE, HeamGen, and INGI consortia, as well as independent collaborative studies. We confirm previous associations of WBC traits with three loci and identified seven novel loci. We also utilize a number of additional analytic methods to infer the functional relatedness of independently implicated loci across WBC phenotypes, as well as investigate direct functional consequences of these loci through analyses of genomic variation affecting the expression of proximal genes in samples of whole blood. In addition, subsequent collaborative efforts with studies of WBC traits in African-American and Japanese cohorts allowed for the investigation of the effects of these genomic variants across populations of diverse continental ancestries.

CUBN is a gene locus for albuminuria

Identification of genetic risk factors for albuminuria may alter strategies for early prevention of CKD progression, particularly among patients with diabetes. Little is known about the influence of common genetic variants on albuminuria in both general and diabetic populations. We performed a meta-analysis of data from 63,153 individuals of European ancestry with genotype information from genome-wide association studies (CKDGen Consortium) and from a large candidate gene study (CARe Consortium) to identify susceptibility loci for the quantitative trait urinary albumin-to-creatinine ratio (UACR) and the clinical diagnosis microalbuminuria. We identified an association between a missense variant (I2984V) in the CUBN gene, which encodes cubilin, and both UACR (P = 1.1 × 10(-11)) and microalbuminuria (P = 0.001). We observed similar associations among 6981 African Americans in the CARe Consortium. The associations between this variant and both UACR and microalbuminuria were significant in individuals of European ancestry regardless of diabetes status. Finally, this variant associated with a 41% increased risk for the development of persistent microalbuminuria during 20 years of follow-up among 1304 participants with type 1 diabetes in the prospective DCCT/EDIC Study. In summary, we identified a missense CUBN variant that associates with levels of albuminuria in both the general population and in individuals with diabetes.

Reliability of serum metabolite concentrations over a 4-month period using a targeted metabolomic approach

Metabolomics is a promising tool for discovery of novel biomarkers of chronic disease risk in prospective epidemiologic studies. We investigated the between- and within-person variation of the concentrations of 163 serum metabolites over a period of 4 months to evaluate the metabolite reliability expressed by the intraclass-correlation coefficient (ICC: the ratio of between-person variance and total variance). The analyses were performed with the BIOCRATES AbsoluteIDQ™ targeted metabolomics technology, including acylcarnitines, amino acids, glycerophospholipids, sphingolipids and hexose in 100 healthy individuals from the European Prospective Investigation into Cancer and Nutrition (EPIC)-Potsdam study who had provided two fasting blood samples 4 months apart. Overall, serum reliability of metabolites over a 4-month period was good. The median ICC of the 163 metabolites was 0.57. The highest ICC was observed for hydroxysphingomyelin C14:1 (ICC = 0.85) and the lowest was found for acylcarnitine C3:1 (ICC = 0). Reliability was high for hexose (ICC = 0.76), sphingolipids (median ICC = 0.66; range: 0.24-0.85), amino acids (median ICC = 0.58; range: 0.41-0.72) and glycerophospholipids (median ICC = 0.58; range: 0.03-0.81). Among acylcarnitines, reliability of short and medium chain saturated compounds was good to excellent (ICC range: 0.50-0.81). Serum reliability was lower for most hydroxyacylcarnitines and monounsaturated acylcarnitines (ICC range: 0.11-0.45 and 0.00-0.63, respectively). For most of the metabolites a single measurement may be sufficient for risk assessment in epidemiologic studies with healthy subjects.

A genome-wide association study identifies LIPA as a susceptibility gene for coronary artery disease

BACKGROUND

eQTL analyses are important to improve the understanding of genetic association results. We performed a genome-wide association and global gene expression study to identify functionally relevant variants affecting the risk of coronary artery disease (CAD).

METHODS AND RESULTS

In a genome-wide association analysis of 2078 CAD cases and 2953 control subjects, we identified 950 single-nucleotide polymorphisms (SNPs) that were associated with CAD at P<10(-3). Subsequent in silico and wet-laboratory replication stages and a final meta-analysis of 21 428 CAD cases and 38 361 control subjects revealed a novel association signal at chromosome 10q23.31 within the LIPA (lysosomal acid lipase A) gene (P=3.7×10(-8); odds ratio, 1.1; 95% confidence interval, 1.07 to 1.14). The association of this locus with global gene expression was assessed by genome-wide expression analyses in the monocyte transcriptome of 1494 individuals. The results showed a strong association of this locus with expression of the LIPA transcript (P=1.3×10(-96)). An assessment of LIPA SNPs and transcript with cardiovascular phenotypes revealed an association of LIPA transcript levels with impaired endothelial function (P=4.4×10(-3)).

CONCLUSIONS

The use of data on genetic variants and the addition of data on global monocytic gene expression led to the identification of the novel functional CAD susceptibility locus LIPA, located on chromosome 10q23.31. The respective eSNPs associated with CAD strongly affect LIPA gene expression level, which was related to endothelial dysfunction, a precursor of CAD.

Polymorphisms in MUC1, MUC2, MUC5B and MUC6 genes are not associated with the risk of chronic atrophic gastritis

Mucins represent major components of the mucous layer in the stomach, protecting the underlying epithelium from acid, mechanical trauma, proteases and pathogenic bacteria. Previous studies have shown an association of neoplastic transformation in the stomach with aberrant mucin levels, suggesting a potential role of genetic variation in mucin genes in the development of gastric cancer (GC). We assessed the association of genetic variation in candidate single nucleotide polymorphisms (SNPs) in mucin genes with the risk of chronic atrophic gastritis (CAG), a well-established precursor of GC in the German population-based ESTHER study. We genotyped MUC1 T31T, MUC2 L58P, MUC2 V116M, MUC5B E34G, MUC5B R51W, MUC5B rs2014486 (intronic) and MUC6 V619M for 533 serologically defined CAG cases and 1054 age- and sex-matched controls. None of the analysed SNPs was associated with CAG. However, large studies are needed to disclose or exclude potential weak associations of these SNPs with CAG risk.

Effect of genome-wide simultaneous hypotheses tests on the discovery rate

An increasing number of genome-wide association studies are being performed in hundreds of thousands of single nucleotide polymorphisms (SNPs). Many of such studies carry on a second stage in which a selected number of SNPs are genotyped in new individuals in order to validate genome-wide findings. Unfortunately, a large proportion of such studies have been unable to validate the genome-wide findings. In this study we aim to better understand how to distinguish the truly associated features from the false positives in genome-wide scans. In order to achieve this goal we use empirical data to look at three aspects that may play a key role in determining which features are called to be associated with the phenotype. First, we examine the usual assumption of a uniform distribution on null p-values and assess whether or not it affects which features are called significant and the number of significant features. Second, we compare the global behavior of the p-value distribution genome-wide with the local behavior at regions such as chromosomes. Third, we look at the effect of minor allele frequency in the p-value distribution. We show empirically that the uniform distribution is not a generally valid assumption and we find that as a consequence strikingly different conclusions can be drawn regarding what we call significant associations and the number of significant findings. We propose that in order to better assign significance to potential associations one needs to estimate the true distribution of null and non-null p-values.

Genetic predictors of fibrin D-dimer levels in healthy adults

BACKGROUND

Fibrin fragment D-dimer, one of several peptides produced when crosslinked fibrin is degraded by plasmin, is the most widely used clinical marker of activated blood coagulation. To identity genetic loci influencing D-dimer levels, we performed the first large-scale, genome-wide association search.

METHODS AND RESULTS

A genome-wide investigation of the genomic correlates of plasma D-dimer levels was conducted among 21 052 European-ancestry adults. Plasma levels of D-dimer were measured independently in each of 13 cohorts. Each study analyzed the association between ≈2.6 million genotyped and imputed variants across the 22 autosomal chromosomes and natural-log–transformed D-dimer levels using linear regression in additive genetic models adjusted for age and sex. Among all variants, 74 exceeded the genome-wide significance threshold and marked 3 regions. At 1p22, rs12029080 (P=6.4×10(-52)) was 46.0 kb upstream from F3, coagulation factor III (tissue factor). At 1q24, rs6687813 (P=2.4×10(-14)) was 79.7 kb downstream of F5, coagulation factor V. At 4q32, rs13109457 (P=2.9×10(-18)) was located between 2 fibrinogen genes: 10.4 kb downstream from FGG and 3.0 kb upstream from FGA. Variants were associated with a 0.099-, 0.096-, and 0.061-unit difference, respectively, in natural-log–transformed D-dimer and together accounted for 1.8% of the total variance. When adjusted for nonsynonymous substitutions in F5 and FGA loci known to be associated with D-dimer levels, there was no evidence of an additional association at either locus.

CONCLUSIONS

Three genes were associated with fibrin D-dimer levels. Of these 3, the F3 association was the strongest, and has not been previously reported.

Large-scale association analysis identifies 13 new susceptibility loci for coronary artery disease

We performed a meta-analysis of 14 genome-wide association studies of coronary artery disease (CAD) comprising 22,233 individuals with CAD (cases) and 64,762 controls of European descent followed by genotyping of top association signals in 56,682 additional individuals. This analysis identified 13 loci newly associated with CAD at P < 5 × 10⁻⁸ and confirmed the association of 10 of 12 previously reported CAD loci. The 13 new loci showed risk allele frequencies ranging from 0.13 to 0.91 and were associated with a 6% to 17% increase in the risk of CAD per allele. Notably, only three of the new loci showed significant association with traditional CAD risk factors and the majority lie in gene regions not previously implicated in the pathogenesis of CAD. Finally, five of the new CAD risk loci appear to have pleiotropic effects, showing strong association with various other human diseases or traits.

Genetic variation in the prostate stem cell antigen gene and upper gastrointestinal cancer in white individuals

BACKGROUND & AIMS

An association between gastric cancer and the rs2294008 (C>T) polymorphism in the prostate stem cell antigen (PSCA) gene has been reported for several Asian populations. We set out to determine whether such an association exists in white individuals.

METHODS

We genotyped 166 relatives of gastric cancer patients, including 43 Helicobacter pylori-infected subjects with hypochlorhydria and gastric atrophy, 65 infected subjects without these abnormalities, 58 H pylori-negative relatives, and 100 population controls. Additionally, a population-based study of chronic atrophic gastritis provided 533 cases and 1054 controls. We then genotyped 2 population-based, case-control studies of upper gastrointestinal cancer: the first included 312 gastric cancer cases and 383 controls; the second included 309 gastric cancer cases, 159 esophageal cancer cases, and 211 controls. Odds ratios were computed from logistic models and adjusted for confounding variables.

RESULTS

Carriage of the risk allele (T) of rs2294008 in PSCA was associated with chronic atrophic gastritis (adjusted odds ratio [OR], 1.5; 95% confidence interval [CI]: 1.1-1.9) and noncardia gastric cancer (OR, 1.9; 95% CI: 1.3-2.8). The association was strongest for the diffuse histologic type (OR, 3.2; 95% CI: 1.2-10.7). An inverse association was observed between carriage of the risk allele and gastric cardia cancer (OR, 0.5; 95% CI: 0.3-0.9), esophageal adenocarcinoma (OR, 0.5; 95% CI: 0.3-0.9), and esophageal squamous cell carcinoma (OR, 0.4; 95% CI: 0.2-0.9).

CONCLUSIONS

The rs2294008 polymorphism in PSCA increases the risk of noncardia gastric cancer and its precursors in white individuals but protects against proximal cancers.

The WST survival assay: an easy and reliable method to screen radiation-sensitive individuals

An easy, fast and reliable method was developed to screen hundreds of Epstein-Barr virus-transformed cell lines (lymphoblastoid cell lines, LCLs) for radiation sensitivity that were generated from lymphocytes isolated from young lung cancer patients. The WST-1 test explores the metabolic activity of the mitochondria as an indicator for the vital status of cells. Cell proliferation as well as indirect cell death can be quantified by this method on a large scale in microtiter plates. Cell survival was measured at 24- and 48-h post-irradiation with 10 Gy ((137)Cs source) by the WST-1 assay and Trypan blue staining. To set up the experimental screening conditions and to establish a positive and a negative control, an ATM-mutated cell line from a radiation-sensitive ATM patient and an ATM proficient cell line from a healthy brother were compared. An optimal differentiation between the two cell lines was demonstrated for 10 Gy and 24- and 48-h cell growth after irradiation. Upon screening 120 LCLs of young lung cancer patients under these conditions, 5 of them were found to be radiation sensitive to a high degree of statistical significance. The results have been confirmed by a second laboratory by means of Trypan blue testing. The WST-1 test represents an efficient and reliable method by means of screening for radiation-sensitive cell lines.

An interferon-induced helicase (IFIH1) gene polymorphism associates with different rates of progression from autoimmunity to type 1 diabetes

OBJECTIVE

Genome-wide association studies have identified gene regions associated with the development of type 1 diabetes. The aim of this study was to determine whether these associations are with the development of autoimmunity and/or progression to diabetes.

RESEARCH DESIGN AND METHODS

Children (n = 1,650) of parents with type 1 diabetes were prospectively followed from birth (median follow-up 10.20 years) for the development of islet autoantibodies, thyroid peroxidase antibodies, tissue transglutaminase antibodies, and diabetes. Genotyping for single-nucleotide polymorphisms of the PTPN22, ERBB3, PTPN2, KIAA0350, CD25, and IFIH1 genes was performed using the MassARRAY system with iPLEX chemistry.

RESULTS

Islet autoantibodies developed in 137 children and diabetes developed in 47 children. Type 1 diabetes risk was associated with the IFIH1 rs2111485 single-nucleotide polymorphism (hazard ratio 2.08; 95% CI 1.16-3.74; P = 0.014). None of the other genes were significantly associated with diabetes development in this cohort. IFIH1 genotypes did not associate with the development of islet autoantibodies (P = 0.80) or autoantibodies against thyroid peroxidase (P = 0.55) and tissue transglutaminase (P = 0.66). Islet autoantibody-positive children with the IFIH1 rs2111485 GG genotype had a faster progression to diabetes (31% within 5 years) than children with the type 1 diabetes protective GA or AA genotypes (11% within 5 years; P = 0.006).

CONCLUSIONS

The findings indicate that IFIH1 genotypes influence progression from autoimmunity to diabetes development, consistent with the notion that protective genotypes downregulate responses to environmental insults after initiation of autoimmunity.

Genetic variants in the FADS gene cluster are associated with arachidonic acid concentrations of human breast milk at 1.5 and 6 mo postpartum and influence the course of milk dodecanoic, tetracosenoic, and trans-9-octadecenoic acid concentrations over the duration of lactation

BACKGROUND

Breastfeeding is considered an optimal nutritional source of n-6 (omega-6) and n-3 (omega-3) fatty acids (FAs) for the proper visual and cognitive development of newborn children. In addition to maternal nutrition as an important regulator of FA concentrations, first results exist on an association of breast-milk FAs with single nucleotide polymorphisms (SNPs) in the FADS gene cluster, which encodes the rate-limiting enzymes in the elongation-desaturation pathway of long-chain polyunsaturated fatty acids (LC-PUFAs).

OBJECTIVE

We analyzed the influence of FADS SNPs on breast-milk FA concentrations and their time course during lactation in the Ulm Birth Cohort study, which comprised 772 nursing mothers at 1.5 mo after giving birth, and in a subset of 463 mothers who were still breastfeeding at 6 mo postpartum.

DESIGN

We conducted linear regression analysis of 8 FADS SNPs with FA concentrations at both time points separately and assessed the genotype effect over time in a longitudinal analysis by using a generalized estimating equation regression model.

RESULTS

We observed significant associations of FADS genotypes with arachidonic acid (AA) concentrations and the 20:4n-6/20:3n-6 ratio at both time points but no association of FADS SNPs with the time course of AA concentrations. A longitudinal analysis of FAs other than LC-PUFAs by genotype over time showed associations for dodecanoic acid, cis-15-tetracosenoic acid, and trans-9-octadecenoic acid.

CONCLUSIONS

Maternal FADS genotypes are associated with breast-milk AA concentrations and might therefore influence the supply of this FA for children. Furthermore, our data indicate an interrelation between the LC-PUFA pathway and saturated and monounsaturated FAs.

Association of a MTNR1B gene variant with fasting glucose and HOMA-B in children and adolescents with high BMI-SDS

CONTEXT

Genome-wide association studies have shown that the melatonin receptor 1B (MTNR1B) gene locus is strongly associated with fasting glucose and β-cell function. However, data are rather limited to the adult population and normal-weight children. So far, little is known whether similar associations are present in overweight and obese children and adolescents.

OBJECTIVE

The aim is to investigate an MTNR1B polymorphism in a sample of 310 overweight and obese children and adolescents (mean body mass index standard deviation score (BMI-SDS)): 2.74 (± 0.55), mean age: 14 (± 2) years), who participated in a short-term weight-loss program based on energy reduction, physical activity, and behavior therapy.

METHODS

We investigated an association between genotype and fasting glucose, fasting insulin, homeostasis model assessment of insulin resistance (HOMA-IR) and of β-cell function (HOMA-B), and anthropometric parameters and their change during intervention.

RESULTS

The minor G allele of polymorphism rs10830963 was significantly associated with increased fasting glucose (0.205  mmol/l, P<0.0001) and decreased HOMA-B (-0.353, P < 0.0001). Categorizing the sample into BMI-SDS groups, these significant associations were abolished in children with BMI-SDS below 2.5 but remained in those with higher BMI-SDS values with stronger β-estimates. The P value for the genotype × BMI-SDS category interaction was 0.012 for fasting glucose and 0.083 for HOMA-B. There was no significant association between genotype and anthropometric parameters and their change during intervention.

CONCLUSIONS

This is the first single study, replicating the association between the MTNR1B locus and diabetes-related traits in overweight and obese children and adolescents. The effect sizes in children and adolescents seem to be stronger than in adults and differed among BMI-SDS categories.

Gaussian graphical modeling reconstructs pathway reactions from high-throughput metabolomics data

Background

With the advent of high-throughput targeted metabolic profiling techniques, the question of how to interpret and analyze the resulting vast amount of data becomes more and more important. In this work we address the reconstruction of metabolic reactions from cross-sectional metabolomics data, that is without the requirement for time-resolved measurements or specific system perturbations. Previous studies in this area mainly focused on Pearson correlation coefficients, which however are generally incapable of distinguishing between direct and indirect metabolic interactions.

Results

In our new approach we propose the application of a Gaussian graphical model (GGM), an undirected probabilistic graphical model estimating the conditional dependence between variables. GGMs are based on partial correlation coefficients, that is pairwise Pearson correlation coefficients conditioned against the correlation with all other metabolites. We first demonstrate the general validity of the method and its advantages over regular correlation networks with computer-simulated reaction systems. Then we estimate a GGM on data from a large human population cohort, covering 1020 fasting blood serum samples with 151 quantified metabolites. The GGM is much sparser than the correlation network, shows a modular structure with respect to metabolite classes, and is stable to the choice of samples in the data set. On the example of human fatty acid metabolism, we demonstrate for the first time that high partial correlation coefficients generally correspond to known metabolic reactions. This feature is evaluated both manually by investigating specific pairs of high-scoring metabolites, and then systematically on a literature-curated model of fatty acid synthesis and degradation. Our method detects many known reactions along with possibly novel pathway interactions, representing candidates for further experimental examination.

Conclusions

In summary, we demonstrate strong signatures of intracellular pathways in blood serum data, and provide a valuable tool for the unbiased reconstruction of metabolic reactions from large-scale metabolomics data sets.

A radiation-induced gene expression signature as a tool to predict acute radiotherapy-induced adverse side effects

The majority of patients tolerate radiotherapy well, but some of them suffer from severe side effects. To find genes possibly predictive for radiosensitivity, mRNA profiles were generated before and 6h after in vitro irradiation with 5Gy. We analyzed lymphocytes from four head and neck and eight breast cancer patients with strong acute radiation toxicity and from 12 matching normal reacting patients in a blind study. Expression was also measured in lymphocyte subpopulations and Epstein-Barr transformed lymphocytes. Radiation response in whole lymphocyte populations was most similar to that of B cells. In peripheral blood lymphocytes of all patients; 153 genes were identified which were statistically significantly altered by a fold change of more than 50% by irradiation. The signatures of radio-responsive genes differed tremendously between primary and transformed cells. Pathway analysis revealed genes involved in p53 signalling, cell cycle control and apoptosis in response to radiation in primary lymphocytes. In these cells, a set of 67 radiation-induced genes was identified capable of differentiating between severe radiosensitive and normal reacting patients. More than one third of such classifying genes belong to the group of apoptosis or cell cycle regulating genes. The classifying potential of the expression signature has now to be validated in further patient cohorts.