Genome-wide analysis of mitochondrial DNA copy number reveals loci implicated in nucleotide metabolism, platelet activation, and megakaryocyte proliferation.

Mitochondrial DNA copy number (mtDNA-CN) measured from blood specimens is a minimally invasive marker of mitochondrial function that exhibits both inter-individual and intercellular variation. To identify genes involved in regulating mitochondrial function, we performed a genome-wide association study (GWAS) in 465,809 White individuals from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium and the UK Biobank (UKB). We identified 133 SNPs with statistically significant, independent effects associated with mtDNA-CN across 100 loci. A combination of fine-mapping, variant annotation, and co-localization analyses was used to prioritize genes within each of the 133 independent sites. Putative causal genes were enriched for known mitochondrial DNA depletion syndromes (p = 3.09 × 10-15) and the gene ontology (GO) terms for mtDNA metabolism (p = 1.43 × 10-8) and mtDNA replication (p = 1.2 × 10-7). A clustering approach leveraged pleiotropy between mtDNA-CN associated SNPs and 41 mtDNA-CN associated phenotypes to identify functional domains, revealing three distinct groups, including platelet activation, megakaryocyte proliferation, and mtDNA metabolism. Finally, using mitochondrial SNPs, we establish causal relationships between mitochondrial function and a variety of blood cell-related traits, kidney function, liver function and overall (p = 0.044) and non-cancer mortality (p = 6.56 × 10-4).

MendelVar: gene prioritization at GWAS loci using phenotypic enrichment of Mendelian disease genes.

MOTIVATION: Gene prioritization at human GWAS loci is challenging due to linkage-disequilibrium and long-range gene regulatory mechanisms. However, identifying the causal gene is crucial to enable identification of potential drug targets and better understanding of molecular mechanisms. Mapping GWAS traits to known phenotypically relevant Mendelian disease genes near a locus is a promising approach to gene prioritization. RESULTS: We present MendelVar, a comprehensive tool that integrates knowledge from four databases on Mendelian disease genes with enrichment testing for a range of associated functional annotations such as Human Phenotype Ontology, Disease Ontology and variants from ClinVar. This open web-based platform enables users to strengthen the case for causal importance of phenotypically matched candidate genes at GWAS loci. We demonstrate the use of MendelVar in post-GWAS gene annotation for type 1 diabetes, type 2 diabetes, blood lipids and atopic dermatitis. AVAILABILITY AND IMPLEMENTATION: MendelVar is freely available at https://mendelvar.mrcieu.ac.uk. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Epigenetic profiling of ADHD symptoms trajectories: a prospective, methylome-wide study.

Attention-deficit/hyperactivity disorder (ADHD) is a prevalent developmental disorder, associated with a range of long-term impairments. Variation in DNA methylation, an epigenetic mechanism, is implicated in both neurobiological functioning and psychiatric health. However, the potential role of DNA methylation in ADHD symptoms is currently unclear. In this study, we examined data from the Avon Longitudinal Study of Parents and Children (ALSPAC)-specifically the subsample forming the Accessible Resource for Integrated Epigenomics Studies (ARIES)-that includes (1) peripheral measures of DNA methylation (Illumina 450k) at birth (n=817, 49% male) and age 7 (n=892, 50% male) and (2) trajectories of ADHD symptoms (7-15 years). We first employed a genome-wide analysis to test whether DNA methylation at birth associates with later ADHD trajectories; and then followed up at age 7 to investigate the stability of associations across early childhood. We found that DNA methylation at birth differentiated ADHD trajectories across multiple genomic locations, including probes annotated to SKI (involved in neural tube development), ZNF544 (previously implicated in ADHD), ST3GAL3 (linked to intellectual disability) and PEX2 (related to perixosomal processes). None of these probes maintained an association with ADHD trajectories at age 7. Findings lend novel insights into the epigenetic landscape of ADHD symptoms, highlighting the potential importance of DNA methylation variation in genes related to neurodevelopmental and peroxisomal processes that play a key role in the maturation and stability of cortical circuits.

DNA methylation and substance-use risk: a prospective, genome-wide study spanning gestation to adolescence.

Epigenetic processes have been implicated in addiction; yet, it remains unclear whether these represent a risk factor and/or a consequence of substance use. Here, we believe we conducted the first genome-wide, longitudinal study to investigate whether DNA methylation patterns in early life prospectively associate with substance use in adolescence. The sample comprised of 244 youth (51% female) from the Avon Longitudinal Study of Parents and Children (ALSPAC), with repeated assessments of DNA methylation (Illumina 450k array; cord blood at birth, whole blood at age 7) and substance use (tobacco, alcohol and cannabis use; age 14-18). We found that, at birth, epigenetic variation across a tightly interconnected genetic network (n=65 loci; q<0.05) associated with greater levels of substance use during adolescence, as well as an earlier age of onset amongst users. Associations were specific to the neonatal period and not observed at age 7. Key annotated genes included PACSIN1, NEUROD4 and NTRK2, implicated in neurodevelopmental processes. Several of the identified loci were associated with known methylation quantitative trait loci, and consequently likely to be under significant genetic control. Collectively, these 65 loci were also found to partially mediate the effect of prenatal maternal tobacco smoking on adolescent substance use. Together, findings lend novel insights into epigenetic correlates of substance use, highlight birth as a potentially sensitive window of biological vulnerability and provide preliminary evidence of an indirect epigenetic pathway linking prenatal tobacco exposure and adolescent substance use.

IGF2BP1, IGF2BP2 and IGF2BP3 genotype, haplotype and genetic model studies in metabolic syndrome traits and diabetes.

OBJECTIVE: Genetic variation at the insulin-like binding protein 2 (IGF2BP2) gene has been associated with type 2 diabetes (T2D) by genome-wide association studies and by replication analyses. Our aim was to explore the underlying genetic model and mechanism of action, factors accounting for non-replications of the associations, and the effect of variation from pathway-related genes IGF2BP1 and IGF2BP3. METHOD: We analysed here the association between T2D (and related traits) and rs4402960 and rs1470579 in IGF2BP2, and rs46522 and rs6949019 (marking IGF2BP1 and IGF2BP3 respectively) from the Age, Gene/Environment Susceptibility (AGES)-Reykjavik Study (N approximately 2500 aged 65-96 years). We undertook a retrospective analysis of the deviations from the multiplicative model in previous studies and the present study. RESULTS: We replicated an association between rs4402960 and T2D status, and reported significant associations with anthropometric traits, fasting insulin, HOMA-IR and HOMA-%B. These associations were also observed for rs1470579, but not for the SNPs marking IGF2BP1 and IGF2BP3. CONCLUSIONS: The lower fasting insulin levels and the impaired beta-cell function associated with IGF2BP2 SNPs are independent of obesity phenotypes. The action of these SNPs on T2D may result from an effect on beta-cell function. This could lead to lower insulin levels, the association with anthropometric traits being secondary. We discuss possible mechanisms of action relating IGF2BP2 with T2D traits. The occurrence of null alleles, the inclusion of T2D patients in analyses of metabolic syndrome risk traits and the genetic model, are possible factors accounting for non-replications of IGF2BP2 associations with T2D.

Refined association mapping for a quantitative trait: weight in the H19-IGF2-INS-TH region.

Previous analyses have provided evidence for one or more loci affecting body weight in the H19-IGF2-INS-TH region on chromosome 11p15. To identify the location of a possible causal locus or loci we applied association analysis by composite likelihood to a large cohort under the Malecot model for body weight. A random sample of 2731 men in the UK were typed for eleven single nucleotide polymorphisms (SNPs) in IGF2, two SNPs in H19, one SNP in INS and one microsatellite marker in the TH genes. Using F tests appropriate to small marker sets, the superiority of regression over correlation was confirmed. All the evidence for association came from IGF2, with P= 0.007 for height-adjusted weight and P= 0.019 for weight additionally adjusted for smoking and alcohol drinking. Although the estimated point location for the suspected causal variant was close to IGF2 ApaI, the 95% confidence and support intervals covered most of IGF2 but none of the other loci. Identification of the causal SNP or SNPs within IGF2 will require typing of more variants in this region.

Genotype of galectin 2 (LGALS2) is associated with insulin-glucose profile in the British Women's Heart and Health Study.

AIMS/HYPOTHESIS: It has been suggested that the gene encoding lymphotoxin-alpha (LTA) is associated with insulin resistance, and genetic association studies in the LTA region offer some support for this. However, LTA is in linkage disequilibrium with both the HLA gene cluster and the gene encoding TNF-alpha, making inferences about causality difficult. In this study, we used the galectin 2 (LGALS2) genotype, which affects LTA secretion but is located on another chromosome than the HLA gene cluster or TNF, to examine the relationship between the LTA pathway and traits of the metabolic syndrome. SUBJECTS: A cross-sectional genetic association study was carried out in 3,272 British women of European origin who were aged 60 to 79 years and were randomly selected from the community. RESULTS: Fasting plasma glucose and serum insulin were statistically significantly associated with LGALS2 rs7291467, with this association being independent of BMI and WHR. The mean difference in fasting insulin per minor allele was -4% (p=0.01 for trend by allele) and the mean per minor allele difference in fasting glucose was -2% (p=0.02 for trend by allele). When women with known diabetes were excluded from the analyses the findings did not differ from those for the whole cohort. CONCLUSIONS/INTERPRETATION: Our findings for the physically unlinked LGALS2, invite further study of LGALS2 specifically and the LTA pathway generally for their influence on glucose-insulin regulation.

Angiotensin II type I receptor gene polymorphism: anthropometric and metabolic syndrome traits.

BACKGROUND: The renin angiotensin system is important in the regulation of vascular tone and fluid and electrolyte balance. The angiotensin converting enzyme gene (ACE) genotype has been shown to affect exercise response and glucose load response dependent on birth weight. Angiotensin II type I receptor (AGTR1) A1166C has previously been associated with the development of hypertension and coronary disease, but its metabolic effects have not been investigated. METHOD: AGTR1 A1166C was genotyped by allele specific PCR in 378 individuals from Hertfordshire, UK, who had been characterised for metabolic syndrome traits. RESULTS: Genotype counts were: AA, 183; AC, 170; CC, 25, consistent with Hardy-Weinberg equilibrium. The CC genotype was associated with significantly lower body mass index (by 1.7 units) in men (p = 0.03), and the same magnitude effect in women with significant lower weight in both genders (p = 0.01), also lower waist circumference and waist-hip ratio (p = 0.01) in men, with a trend for lower waist circumference in women also. Additionally, the CC genotype and/or C allele was associated with lower fasting glucose and insulin, and 30 and 120 min glucose in men (respectively, p = 0.08, 0.04, 0.01, 0.06). Lower means of systolic blood pressure, pulse pressure, cholesterol, and fasting triglyceride were also observed for the CC genotype in both genders though these were not statistically significant. CONCLUSIONS: The AGTR1 1166 CC genotype appears to predispose to favourable anthropometric and metabolic traits, relative to cardiovascular risk.

Positive associations between single nucleotide polymorphisms in the IGF2 gene region and body mass index in adult males.

We previously demonstrated an association between the insulin-like growth factor 2 (IGF2) gene 3'-untranslated region (3'-UTR) ApaI polymorphism and body mass index (BMI) in over 2500 middle-aged Caucasoid males. In the same cohort, we have now tested association with 11 more markers, including seven novel single nucleotide polymorphisms (SNPs), spanning >30 kb across the IGF2 gene. Three SNPs showed significant positive associations with BMI: 6815 A/T in the IGF2 P1 promoter (P = 0.00012, n = 2394) and the newly identified SNPs 1156 C/T in intron 2 (P = 0.017, n = 1567) and 1926 C/G in the 3'-UTR (P = 0.0062, n = 1872). There was strong pairwise linkage disequilibrium (LD) between the ApaI and 1926 C/G sites, whereas LD between ApaI and 6815 A/T, and between ApaI and 1156 T/C, was minimal. Univariately 6815 A/T, 1156 T/C and ApaI explained 1.03, 1.02 and 0.67% of the variation in BMI. Multi-way analysis of variance (ANOVA) models showed that 6815 A/T and 1156 T/C explained a further 0.4 and 0.8% of the variation beyond that accounted for by ApaI and the association of 1926 C/G with BMI disappeared after adjustment. The 6815 A/T, 1156 T/C and ApaI markers in effect constitute independent affirmations of our original hypothesized candidate gene region. In a stepwise multi-way ANOVA model, all three terms were significantly independently associated with BMI. The total proportion of BMI variance explained by this model was 2.25%, strongly suggesting that IGF2 genetic variation is a significant determinant of body weight in middle-aged males.

Linkage analysis of the 5q31-33 candidate region for asthma in 240 UK families.

Atopy and asthma are complex genetic diseases resulting from the interactions of a number of genetic and environmental factors. We had previously reported allelic association between the IL9 marker on chromosome 5q31-33 and atopy. In order to further investigate the role of susceptibility genes on 5q31-33 in the development of atopy and asthma we have studied 240 UK families comprising 131 families selected at random, 60 multiplex families with affected sib pairs, and 49 single proband nuclear families. Polymorphic markers on 5q31-33 were genotyped and both single and multipoint linkage analysis was undertaken using the BETA program. We have used both affection status and quantitative scores for atopy and asthma for phenotypic variables, combining data into scores for asthma and atopy. The strongest suggestion of linkage using multipoint analysis was centred around D5S410 with a maximum Lod of 1.946 at location 171.3 cM and a standard error of 3.3 for the asthma quantitative score. There was no evidence of linkage with atopy, the atopy quantitative score or total serum IgE.

SNP genotyping by combination of 192-well MADGE, ARMS and computerized gel image analysis.

A new modification of the microplate array diagonal gel electrophoresis (MADGE) system accommodates the dual amplification refractory mutation system (ARMS) products of 96 samples on one 192-well gel. Simultaneous electrophoresis of a number of horizontal ARMS-MADGE gels achieves high throughput. Gels are imaged digitally, here using the FluorImager 595 fluorescent scanning system. Customized software by Phoretix enables rapid computerized calling of band patterns in ARMS-MADGE arrays, in which the two wells receiving a pair of allele-specific assays for a single template are juxtaposed to form one virtual track, with genotype data exported directly into Microsoft Excel for statistical analysis. An ARMS assay of the A/T base change at the -23/HphI RFLP in the insulin gene promoter, which initiates from 2.5 ng template DNA, was used here to demonstrate this improved general approach for population SNP analyses.